自分へのメモ、やり方のみ
##############################################
#
#因子分析は、psychパッケージのfactanal
#
##############################################
install.packages("psych",dependencies=TRUE)
library(psych)
#今回データはbfiを使う
summary(bfi)
describe(bfi)
mydata<-na.omit(bfi)
#スクリープロットを作図
VSS.scree(mydata)
#因子分析
factanal(mydata,6,rotation="promax")
#相関を使って求める
mycor<-cor(mydata)
factanal(covmat=mycor,factors=6,rotation="promax")
##############################################
#
#コレポンは、MASSパッケージのcorresp
#
##############################################
library(MASS)
myc<-corresp(mydata,nf=3)
#1軸、2軸の表示
biploat(myc)
##############################################
#
#おまけ。回帰。式だけ。出てくる値に意味なさそう
#
##############################################
lm(education~.,data=mydata)
<HR>
■参考サイト
http://mizumot.com/handbook/?page_id=545
http://www.slideshare.net/mitsuoshimohata/ss-24419059
http://www1.doshisha.ac.jp/~mjin/R/26/26.html
##############################################
#
#因子分析は、psychパッケージのfactanal
#
##############################################
install.packages("psych",dependencies=TRUE)
library(psych)
#今回データはbfiを使う
summary(bfi)
describe(bfi)
mydata<-na.omit(bfi)
#スクリープロットを作図
VSS.scree(mydata)
#因子分析
factanal(mydata,6,rotation="promax")
#相関を使って求める
mycor<-cor(mydata)
factanal(covmat=mycor,factors=6,rotation="promax")
##############################################
#
#コレポンは、MASSパッケージのcorresp
#
##############################################
library(MASS)
myc<-corresp(mydata,nf=3)
#1軸、2軸の表示
biploat(myc)
##############################################
#
#おまけ。回帰。式だけ。出てくる値に意味なさそう
#
##############################################
lm(education~.,data=mydata)
<HR>
■参考サイト
http://mizumot.com/handbook/?page_id=545
http://www.slideshare.net/mitsuoshimohata/ss-24419059
http://www1.doshisha.ac.jp/~mjin/R/26/26.html