1

有意性検定を使用cor.ciしてポリコリック相関を取得しようとしていますが、エラー メッセージが表示され続けます。コードは次のとおりです。

install.packages("Hmisc")
library(Hmisc)
mydata <- spss.get("S-IAT for R.sav", use.value.labels=TRUE)

install.packages('psych')
library(psych)
poly.example <- cor.ci(mydata(nvar = 10,n = 100)$items,n.iter = 10,poly = TRUE)
poly.example
print(corr.test(poly.example$rho), short=FALSE)

これが与えるエラーメッセージは次のとおりです。

> library(psych)  
> poly.example <- cor.ci(mydata(nvar = 10,n = 100)$items,n.iter = 10,poly = TRUE)  
Error in cor.ci(mydata(nvar = 10, n = 100)$items, n.iter = 10, poly = TRUE) :  
  could not find function "mydata"  
> poly.example  
Error: object 'poly.example' not found  
> print(corr.test(poly.example$rho), short=FALSE) 
Error in is.data.frame(x) : object 'poly.example' not found

mydata を認識させたり、分析のためにこのデータセットから特定の変数を選択したりするにはどうすればよいですか? ここから上記のコードを取得しました：
Rで有意なポリコリック相関行列

ありがとう！

4

1 に答える 1

0

いくつかの問題があります。1) 前にコメントしたように、mydata を関数として扱っていますが、data.frame として扱う必要があります。したがって、呼び出しは

 poly.example <- cor.ci(mydata,n.iter = 10,poly = TRUE)

最初の 100 個のケースと最初の 10 個の変数だけを取得しようとしている場合は、

poly.example <- cor.ci(mydata[1:10,1:100],n.iter = 10,poly = TRUE)

2) 次に、結果の相関行列に対して corr.test を実行したくない場合。データに対して corr.test を実行する必要があります。

print(corr.test(mydata[1:10,1:100],short=FALSE)

corr.test はピアソン相関をテストしていることに注意してください。

于 2015-05-29T13:37:24.470 に答える