2014-10-01 99 views
-4

我有一個包含13個屬性的數據集。我使用R.散點圖和R,難度很大

  1. 我想以產生用於每對屬性的散點圖(帶有適當標記x和y軸),這將會給我共有156個散點圖。

  2. 我希望所有這些156散點圖以1個JPEG/PDF文件顯示。

  3. 我也想找到相關係數最高的一對。

我該如何去做這件事?我很困擾。

+0

您是否試圖爲變量之間的關係創建散點圖?這個情節看起來很亂,很難解釋。 – 2014-10-01 18:19:53

+2

您能否向我們提供一些樣本數據?你試過什麼了? – 2014-10-01 18:20:01

+0

你是否使用過Google? ['雙'功能這樣做。](http://astrostatistics.psu.edu/su07/R/html/graphics/html/pairs.html) – 2014-10-01 18:20:21

回答

2

答案其實很簡單。

plot(mtcars) 

enter image description here

,並保存它:

png() 
plot(mtcars) 
dev.off() 

而獲得的相關性:

cor(mtcars) 
      mpg  cyl  disp   hp  drat   wt 
mpg 1.0000000 -0.8521620 -0.8475514 -0.7761684 0.68117191 -0.8676594 
cyl -0.8521620 1.0000000 0.9020329 0.8324475 -0.69993811 0.7824958 
disp -0.8475514 0.9020329 1.0000000 0.7909486 -0.71021393 0.8879799 
hp -0.7761684 0.8324475 0.7909486 1.0000000 -0.44875912 0.6587479 
drat 0.6811719 -0.6999381 -0.7102139 -0.4487591 1.00000000 -0.7124406 
wt -0.8676594 0.7824958 0.8879799 0.6587479 -0.71244065 1.0000000 
qsec 0.4186840 -0.5912421 -0.4336979 -0.7082234 0.09120476 -0.1747159 
vs 0.6640389 -0.8108118 -0.7104159 -0.7230967 0.44027846 -0.5549157 
am 0.5998324 -0.5226070 -0.5912270 -0.2432043 0.71271113 -0.6924953 
gear 0.4802848 -0.4926866 -0.5555692 -0.1257043 0.69961013 -0.5832870 
carb -0.5509251 0.5269883 0.3949769 0.7498125 -0.09078980 0.4276059 
      qsec   vs   am  gear  carb 
mpg 0.41868403 0.6640389 0.59983243 0.4802848 -0.55092507 
cyl -0.59124207 -0.8108118 -0.52260705 -0.4926866 0.52698829 
disp -0.43369788 -0.7104159 -0.59122704 -0.5555692 0.39497686 
hp -0.70822339 -0.7230967 -0.24320426 -0.1257043 0.74981247 
drat 0.09120476 0.4402785 0.71271113 0.6996101 -0.09078980 
wt -0.17471588 -0.5549157 -0.69249526 -0.5832870 0.42760594 
qsec 1.00000000 0.7445354 -0.22986086 -0.2126822 -0.65624923 
vs 0.74453544 1.0000000 0.16834512 0.2060233 -0.56960714 
am -0.22986086 0.1683451 1.00000000 0.7940588 0.05753435 
gear -0.21268223 0.2060233 0.79405876 1.0000000 0.27407284 
carb -0.65624923 -0.5696071 0.05753435 0.2740728 1.00000000 

而且,以獲得最高的相關性是一個比較棘手:

cor_values = cor(mtcars) 
diag(cor_values) = NA # Set diagonal to NA or else that will be the highest correlation 
which(cor_values == max(cor_values, na.rm = TRUE), arr.ind = TRUE) 
    row col 
disp 3 2 
cyl 2 3 

在這種情況下,涉及dispcyl變量。

+1

我喜歡用'plot(mtcars,gap = 0)' – 2014-10-01 21:50:26

+0

@BenBolker,這確實是一個非常不錯的添加。 – 2014-10-02 06:36:49