text in bold.
need knitr and markdown
install.packages(“knitr”)
install.packages(“RCurl”)
d = read.table("maps.mixx.epi.bfi.data", header = T)
summary(d)
## epiE epiS epiImp epilie
## Min. : 1.0 Min. : 0.00 Min. :0.00 Min. :0.00
## 1st Qu.:11.0 1st Qu.: 6.00 1st Qu.:3.00 1st Qu.:1.00
## Median :14.0 Median : 8.00 Median :4.00 Median :2.00
## Mean :13.3 Mean : 7.58 Mean :4.37 Mean :2.38
## 3rd Qu.:16.0 3rd Qu.: 9.50 3rd Qu.:6.00 3rd Qu.:3.00
## Max. :22.0 Max. :13.00 Max. :9.00 Max. :7.00
## epiNeur bfagree bfcon bfext bfneur
## Min. : 0.0 Min. : 74 Min. : 53 Min. : 8.0 Min. : 34
## 1st Qu.: 7.0 1st Qu.:112 1st Qu.: 99 1st Qu.: 87.5 1st Qu.: 70
## Median :10.0 Median :126 Median :114 Median :104.0 Median : 90
## Mean :10.4 Mean :125 Mean :113 Mean :102.2 Mean : 88
## 3rd Qu.:14.0 3rd Qu.:136 3rd Qu.:128 3rd Qu.:118.0 3rd Qu.:104
## Max. :23.0 Max. :167 Max. :178 Max. :168.0 Max. :152
## bfopen bdi traitanx stateanx
## Min. : 73 Min. : 0.00 Min. :22 Min. :21.0
## 1st Qu.:110 1st Qu.: 3.00 1st Qu.:32 1st Qu.:32.0
## Median :125 Median : 6.00 Median :38 Median :38.0
## Mean :123 Mean : 6.78 Mean :39 Mean :39.9
## 3rd Qu.:136 3rd Qu.: 9.00 3rd Qu.:44 3rd Qu.:46.5
## Max. :173 Max. :27.00 Max. :71 Max. :79.0
rownames(d)
## [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11"
## [12] "12" "13" "14" "15" "16" "17" "18" "19" "20" "21" "22"
## [23] "23" "24" "25" "26" "27" "28" "29" "30" "31" "32" "33"
## [34] "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44"
## [45] "45" "46" "47" "48" "49" "50" "51" "52" "53" "54" "55"
## [56] "56" "57" "58" "59" "60" "61" "62" "63" "64" "65" "66"
## [67] "67" "68" "69" "70" "71" "72" "73" "74" "75" "76" "77"
## [78] "78" "79" "80" "81" "82" "83" "84" "85" "86" "87" "88"
## [89] "89" "90" "91" "92" "93" "94" "95" "96" "97" "98" "99"
## [100] "100" "101" "102" "103" "104" "105" "106" "107" "108" "109" "110"
## [111] "111" "112" "113" "114" "115" "116" "117" "118" "119" "120" "121"
## [122] "122" "123" "124" "125" "126" "127" "128" "129" "130" "131" "132"
## [133] "133" "134" "135" "136" "137" "138" "139" "140" "141" "142" "143"
## [144] "144" "145" "146" "147" "148" "149" "150" "151" "152" "153" "154"
## [155] "155" "156" "157" "158" "159" "160" "161" "162" "163" "164" "165"
## [166] "166" "167" "168" "169" "170" "171" "172" "173" "174" "175" "176"
## [177] "177" "178" "179" "180" "181" "182" "183" "184" "185" "186" "187"
## [188] "188" "189" "190" "191" "192" "193" "194" "195" "196" "197" "198"
## [199] "199" "200" "201" "202" "203" "204" "205" "206" "207" "208" "209"
## [210] "210" "211" "212" "213" "214" "215" "216" "217" "218" "219" "220"
## [221] "221" "222" "223" "224" "225" "226" "227" "228" "229" "230" "231"
colnames(d)
## [1] "epiE" "epiS" "epiImp" "epilie" "epiNeur" "bfagree"
## [7] "bfcon" "bfext" "bfneur" "bfopen" "bdi" "traitanx"
## [13] "stateanx"
str(d)
## 'data.frame': 231 obs. of 13 variables:
## $ epiE : int 18 16 6 12 14 6 15 18 15 8 ...
## $ epiS : int 10 8 1 6 6 4 9 9 11 5 ...
## $ epiImp : int 7 5 3 4 5 2 4 7 3 2 ...
## $ epilie : int 3 1 2 3 3 5 3 2 3 2 ...
## $ epiNeur : int 9 12 5 15 2 15 12 10 1 10 ...
## $ bfagree : int 138 101 143 104 115 110 109 92 127 74 ...
## $ bfcon : int 96 99 118 106 102 113 58 57 108 100 ...
## $ bfext : int 141 107 38 64 103 61 99 94 108 61 ...
## $ bfneur : int 51 116 68 114 86 54 55 72 35 87 ...
## $ bfopen : int 138 132 90 101 118 149 110 114 86 89 ...
## $ bdi : int 1 7 4 8 8 5 7 0 0 7 ...
## $ traitanx: int 24 41 37 54 39 51 40 32 22 35 ...
## $ stateanx: int 22 40 44 40 67 38 32 41 26 31 ...
attach(d)
plot(bfagree, bfopen, main = "Correlation between Agreeableness and Openness in Big Five",
xlab = "Agreeableness", ylab = "Openness", pch = 20)
abline(lm(bfagree ~ bfopen), col = "blue")
lines(lowess(bfagree, bfopen), col = "purple")
text(165, 120, cor(bfagree, bfopen))
