2017-04-06 48 views
0

我有兩個數據框motivation_on有60個觀察值和motivation_off有146個觀察值,每個數據框包含21個變量值和1個ID列,它位於第一列。R - 如何限制hmisc rcorr的輸出?

現在我想知道瓦爾如何相互關聯,所以我使用:

rcorr(as.matrix(motivation_on[2:ncol(motivation_on)]), type = "spearman") 

rcorr(as.matrix(motivation_off[2:ncol(motivation_off)]), type = "spearman") 

(該子集是爲了擺脫ID列)

現在我想計算在線和離線變量之間的相關性,所以我試過了:

rcorr(as.matrix(motivation_off[2:ncol(motivation_off)]), as.matrix(motivation_on[2:ncol(motivation_on)]) , type = "spearman") 

現在我得到我想要的,但此外它還顯示了motivation_on和之前計算的motivation_off內的變量的所有相關性。這使得輸出非常長。如何完全獲得on_off相關性的rcorr輸出?

編輯澄清: 嘗試以下操作:

x <- as.matrix(mtcars[1:3]) 
y <- as.matrix(mtcars[4:6]) 
rcorr(x,y) 

我想什麼是相關表:MPG,共青團,DISP爲行和HP,DRAT,重量爲列,而不是完整輸出。我目前的工作圍繞:

z <- rcorr(x,y) 
q <- as.data.frame(z[1]) 
q[1:3,4:6] 
+1

快速瀏覽代碼我不認爲這是可能的,但'psych :: corr.test(x,y)'會給你你想要的東西 – user20650

+0

或當然'cor(x,y)'if你不需要pvalues – user20650

+0

不幸的是我也需要p值。 – florian

回答

1

似乎hmisc的rcorr不提供選項提取我想要的東西。我只能手動提取結果,就像我在mtcars例子中演示的那樣。 但是psych::corr.test(x, y)提供了我想要的輸出 - 感謝@ user20650指出了這一點!