# Κεφάλαιο 10
# Κατασκευή διαγράμματος διασποράς και υπολογισμός των συντελεστών συσχέτισης Pearson r, Spearman rho και μερικής συσχέτισης στο R

# Φορτώνουμε τα ακόλουθα πακέτα
library(foreign)
library(Hmisc)
library(psych)

# Απόδοση των δεδομένων του αρχείου chapter10_1.sav στο πλαίσιο δεδομένων ch10a
ch10a<-read.spss('chapter10_1.sav', to.data.frame=T)

# Βασικός έλεγχος του πλαισίου δεδομένων
str(ch10a)
summary(ch10a)

# Αυτόματο διάγραμμα διασποράς των δεδομένων
plot(ch10a)
# Διάγραμμα διασποράς με ορισμό των μεταβλητών που συσχετίζονται
plot(esteem~disatisf, ch10a, xlab="Δυσαρέσκεια για το σώμα", ylab="Αυτοπεποίθηση", las=1)

# Υπολογισμός συντελεστή συσχέτισης Pearson r
# Η εντολή υπολογίζει τον συντελεστή Pearson r για υπόθεση δίπλευρου ελέγχου
cor.test(ch10a$disatisf,ch10a$esteem)
# Προσθέτοντας τα ορίσματα alternative και method μπορούμε να παραμετροποιήσουμε τον υπολογισμό ως εξής:
# Το όρισμα alternative παίρνει τις τιμές “two.sided”, η οποία αναφέρεται σε υπόθεση δίπλευρου ελέγχου, “greater”, η οποία αναφέρεται σε θετική συσχέτιση, και “less”, η οποία αναφέρεται σε αρνητική συσχέτιση
# Το όρισμα method δέχεται μία από τις τρεις ακόλουθες επιλογές: “pearson”, “kendall” ή “spearman”
cor.test(ch10a$disatisf,ch10a$esteem, alternative="two.sided", method="pearson")

# Απόδοση των δεδομένων του αρχείου chapter10_2.sav στο πλαίσιο δεδομένων ch10b
ch10b<-read.spss('chapter10_2.sav', to.data.frame=T)

# Υπολογισμός συντελεστή συσχέτισης Spearman rho
cor.test(ch10b$degree,ch10b$unemploy, alternative="less", method="spearman")

# Απόδοση των δεδομένων του αρχείου chapter10_3.sav στο πλαίσιο δεδομένων ch10c
ch10c<-read.spss('chapter10_3.sav', to.data.frame=T)

# Πίνακας συσχετίσεων από τη βιβλιοθήκη Hmisc
rcorr(as.matrix(ch10c))

# Συντελεστής μερικής συσχέτισης από τη βιβλιοθήκη psych
partial.r(ch10c, c(1,2), 3)[1,2]