Differences

This shows you the differences between two versions of the page.

--- r:correlation [2018/12/26 15:00] – hkimscil
+++ r:correlation [2018/12/26 15:06] – hkimscil
@@ Line 5: / Line 5: @@
 <code>
 data("mtcars")
-my_data <- mtcars[, c(1,3,4,5,6,7)]
+my_data <- mtcars[, c(1,2,3,4,5,6,7)]
 head(my_data)
 </code>
-<code>                   mpg disp  hp drat    wt  qsec
-Mazda RX4         21.0  160 110 3.90 2.620 16.46
+<code>
-Mazda RX4 Wag     21.0  160 110 3.90 2.875 17.02
+                   mpg cyl disp  hp drat    wt  qsec
-Datsun 710        22.8  108  93 3.85 2.320 18.61
+Mazda RX4         21.0   6  160 110 3.90 2.620 16.46
-Hornet 4 Drive    21.4  258 110 3.08 3.215 19.44
+Mazda RX4 Wag     21.0   6  160 110 3.90 2.875 17.02
-Hornet Sportabout 18.7  360 175 3.15 3.440 17.02
+Datsun 710        22.8   4  108  93 3.85 2.320 18.61
-Valiant           18.1  225 105 2.76 3.460 20.22</code>
+Hornet 4 Drive    21.4   6  258 110 3.08 3.215 19.44
+Hornet Sportabout 18.7   8  360 175 3.15 3.440 17.02
+Valiant           18.1   6  225 105 2.76 3.460 20.22
+>
+</code>
 <code>
 install.packages("Hmisc")
@@ Line 22: / Line 27: @@
 <code>
-rcorr(my_data)
+rcorr(as.matrix(my_data))
 </code>
+<code>
+       mpg  disp    hp  drat    wt  qsec
+mpg   1.00 -0.85 -0.78  0.68 -0.87  0.42
+disp -0.85  1.00  0.79 -0.71  0.89 -0.43
+hp   -0.78  0.79  1.00 -0.45  0.66 -0.71
+drat  0.68 -0.71 -0.45  1.00 -0.71  0.09
+wt   -0.87  0.89  0.66 -0.71  1.00 -0.17
+qsec  0.42 -0.43 -0.71  0.09 -0.17  1.00
+n= 32
+P
+     mpg    disp   hp     drat   wt     qsec
+mpg         0.0000 0.0000 0.0000 0.0000 0.0171
+disp 0.0000        0.0000 0.0000 0.0000 0.0131
+hp   0.0000 0.0000        0.0100 0.0000 0.0000
+drat 0.0000 0.0000 0.0100        0.0000 0.6196
+wt   0.0000 0.0000 0.0000 0.0000        0.3389
+qsec 0.0171 0.0131 0.0000 0.6196 0.3389
+> </code>