source("http://blue.zero.jp/yokumura/R/haebara/procfactor.txt") dat <- read.csv("http://blue.zero.jp/yokumura/R/haebara/chap10_1.csv",header=F) colnames(dat) <- c("‰·˜a", "—z‹C" , "ŠOŒό“I" ,"eΨ", "ŽΠŒπ“I" ,"‹¦—Ν“I", "Ο‹Ι“I", "‘f’Ό") rownames(dat) <- c("‰·˜a", "—z‹C" , "ŠOŒό“I" ,"eΨ", "ŽΠŒπ“I" ,"‹¦—Ν“I", "Ο‹Ι“I", "‘f’Ό") procfactor(data=dat,n=300,nfact=2) fa1 <- procfactor(data=dat,n=300,nfact=2) fa1$Rota #promax‰ρ“]Œγ‚ΜˆφŽqƒpƒ^ƒ“ fa1$Int #ˆφŽqŠΤ‘ŠŠΦ fa1$Com #‹€’ʐ« fa1$Fit #“K‡“xŽw•W #nfact‚ΕŽw’θ‚΅‚½ˆφŽq”‚πγŒΐ‚ɁC #chisq, df, CAIC, BIC, RMSEA‚πŽZo fa1$Cor #•ͺΝ‚ΕŽg—p‚΅‚½‘ŠŠΦs—ρ