×
[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。
library(psych)
library(polycor)
library(foreign)
dat <- read.spss("http://www.indiana.edu/~statmath/stat/all/cfa/values_ord.sav", to.data.frame=T, max.value.labels=99)
pcres <- polychoric(dat)
print(pcres, digits=3)
# データをきちんと変換しておけばhetcor関数のほうが便利かも
x <- dat
dat2 <- data.frame(lapply(x, ordered))
sapply(dat2, class)
res1 <- polychoric(dat2) # エラー
res2 <- hetcor(dat2)
res2
res2$correlations # 相関行列だけとりだす
# 以下参考。感謝
# Rjpwiki: polycorパッケージによる順序相関係数の算出
# 小杉先生の順序尺度水準の相関係数について: http://www.kosugitti.net/labo/polynote.pdf
library(polycor)
library(foreign)
dat <- read.spss("http://www.indiana.edu/~statmath/stat/all/cfa/values_ord.sav", to.data.frame=T, max.value.labels=99)
pcres <- polychoric(dat)
print(pcres, digits=3)
# データをきちんと変換しておけばhetcor関数のほうが便利かも
x <- dat
dat2 <- data.frame(lapply(x, ordered))
sapply(dat2, class)
res1 <- polychoric(dat2) # エラー
res2 <- hetcor(dat2)
res2
res2$correlations # 相関行列だけとりだす
# 以下参考。感謝
# Rjpwiki: polycorパッケージによる順序相関係数の算出
# 小杉先生の順序尺度水準の相関係数について: http://www.kosugitti.net/labo/polynote.pdf
PR
Comment
Trackback
Trackback URL
Comment form
