2013-10-21 72 views
2

我一直在急切地尋找一種計算多方相關矩陣的方法,在R中具有顯着性。如果這非常困難,那麼具有顯着性的兩個變量之間的多方相關就足夠了。具有顯着性的多方相關矩陣R

我迄今爲止嘗試:

library(polychor) 
poly <- polychor(var1,var2) 
poly <- polychor(DatM) #where DatM is a DF converted to matrix 

library(polycor) 
hetcor(Dat2) #I am however uncertain hetcor is something I would want if I am looking for polychoric correlation. 

library(psych) 
polychoric(Dat$for2a,smooth=TRUE,global=TRUE,polycor=FALSE, ML = FALSE, std.err=TRUE) 

這些都不放棄的意義。我已經讀過Permutation測試[可以給我意義] [1],用於每個相關對的corr = 0假設。軟件包和lmPerm允許用戶計算排列測試。但是,我不確定如何。

在此先感謝, 海爾吉

回答

4

很抱歉這麼晚纔回復。可以指定poly = TRUE。比爾最近也改變了cor.plot,使其擴展與對比意義的文字大小,但你可能必須從我們的服務器上安裝最新的鬥志包:

install.packages('psych', type = 'source', repos = 'http://personality-project.org/r/') 

下面是一個簡單的例子:

library(psych) 
poly.example <- cor.ci(sim.poly(nvar = 10,n = 100)$items,n.iter = 10,poly = TRUE) 
poly.example 
print(corr.test(poly.example$rho), short=FALSE) 
+0

謝謝您的回答。我不再在這個項目上工作,但下次我需要多邊形相關矩陣時,我一定會嘗試一下。 –

+0

您的意思是corr.p(ply.examples $ rho)嗎?因爲rho已經是一個相關矩陣。 –