Rで因子分析: KMOとBartlett | myscratchpadspssでは標準で出力される (らしい) 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
