2017-04-06 15 views
0

私は60の観測とmotivation_offと146の観測との2つのデータフレームmotif_onを持っています。それぞれ21のvarsと1つのID列から成ります。これは最初の列にあります。R - hmisc rcorrの出力を制限する方法は?

今私はVARSは、私が使用して互いに相関方法を知りたい:

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内のすべてのvarsの相関を表示します。これにより、出力が極端に長くなります。どのようにon_offの相関関係のためのrcorr出力を排他的に得るか?明確化のため

編集: 次のことを試してみてください。

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

私が望んでいたこととの相関表である:カラムといっぱいになっていないなどの重量のmpg、CYL、行やHPなどDISP、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のおかげで!

関連する問題