myscratchpad

spssでは標準で出力される (らしい) Kaiser-Meyer-Olkinの標本妥当性と、Bartlettの球面性検定の関数を見つけたのでメモしておく

# Kaiser-Meyer-Olkinの標本妥当性
青木先生のページから。感謝。
kmo(dat) # datは数値のみのデータフレームか行列

全く同じものがこのページにもあった
kmo.test <- function(df)
{
        cor.sq = cor(df)^2
    cor.sumsq = (sum(cor.sq)-dim(cor.sq)[1])/2
    library(corpcor)
    pcor.sq = cor2pcor(cor(df))^2
    pcor.sumsq = (sum(pcor.sq)-dim(pcor.sq)[1])/2
    kmo = cor.sumsq/(cor.sumsq+pcor.sumsq)
    return(kmo)
}


# Bartlettの球面性
# psychパッケージは偉大すぎる
library(psych)
cortest.bartlettdat)
#cortest.bartlett(cor(dat)) # 相関行列を与えること、とマニュアルにはあるが、spssの結果と一致するのはローデータを与えたときのようだ