Rで因子分析、コレスポンデンス分析とかの続き
##########################################
# 因子得点を求める
# scoresを指定しないと出てこない
##########################################
result <-factanal(mydata,6,rotation="promax",scores="Bartlett")
##########################################
# 因子得点をもとにデンドログラムを出す
##########################################
kyori <- dist(result$scores, method = "euclidean")
hcl <- hclust(kyori)
plot(hcl)
##########################################
# 因子得点をもとにクラスタリング
# デンドログラムから、何個のクラスタか決める
##########################################
kmeans(result$scores,5)
【参考文献】
http://rstudio-pubs-static.s3.amazonaws.com/3250_107c2910fec74c68a443a0769e8d22d9.html
※おまけ
ヘッダー&列の名前つきCSV読み込み
read.csv("sample_data.csv",header=T,row.name=1)
【参考文献】
http://cse.naro.affrc.go.jp/takezawa/r-tips/r/40.html
##########################################
# 因子得点を求める
# scoresを指定しないと出てこない
##########################################
result <-factanal(mydata,6,rotation="promax",scores="Bartlett")
##########################################
# 因子得点をもとにデンドログラムを出す
##########################################
kyori <- dist(result$scores, method = "euclidean")
hcl <- hclust(kyori)
plot(hcl)
##########################################
# 因子得点をもとにクラスタリング
# デンドログラムから、何個のクラスタか決める
##########################################
kmeans(result$scores,5)
【参考文献】
http://rstudio-pubs-static.s3.amazonaws.com/3250_107c2910fec74c68a443a0769e8d22d9.html
※おまけ
ヘッダー&列の名前つきCSV読み込み
read.csv("sample_data.csv",header=T,row.name=1)
【参考文献】
http://cse.naro.affrc.go.jp/takezawa/r-tips/r/40.html