#########################################################
##excel vba program for faccon.exe コバンザメアプリ v1.6の例題 
#http://www.ec.kagawa-u.ac.jp/~hori/delphistat/hattori.html
#########################################################
source("http://blue.zero.jp/yokumura/R/testtheory/ctt.txt")


#------------------------------------------------------#
#est.faccon(
# data = 生データ/相関行列
# parallel = 100 (平行分析の時の試行数)
# ns       = 相関行列を指定した場合の標本サイズ 
# nfactors = 1:3 (1因子-3因子までを分析。ただし，ledermanの限界は超えないよう分析される)
# method = ML/ULS (最尤推定法，最小二乗法)
# data.type = raw/cor (生データ/相関行列)
# path = faccon.exeがインストールされているフォルダ
# sort = F (因子パターンをソートするか否か)
#
# $fitstat...適合度指標
# $nf1-nf因子数...当該因子数の因子パタン (プロマックス回転後)，共通性，独自性
# $nf1corr-nf因子数corr...当該因子数の因子間相関
#------------------------------------------------------#

#------------------------------------#
#忍者ハットリ君 thurstone21 
#------------------------------------#
dat1 <- as.matrix(read.table("http://blue.zero.jp/yokumura/R/factor/thurstone21.csv", header=T, sep=",", row.names="X"))
colnames(dat1) <- rownames(dat1)
dat1 <- lower2full(dat1)
dat1
est.faccon(data=dat1, nfactors=1:10, data.type="cor", method="ML", ns=437, path="Z:/data/factor/faccon/")
est.faccon(data=dat1, nfactors=1:10, data.type="cor", method="ML", ns=437, sort=T, path="Z:/data/factor/faccon/")
est.faccon(data=dat1, nfactors=1:10, data.type="cor", method="ULS", ns=437, sort=T, path="Z:/data/factor/faccon/")



#------------------------------------#
#忍者ハットリ君 Thurstone因子間 
#------------------------------------#
dat2 <- as.matrix(read.table("http://blue.zero.jp/yokumura/R/factor/thurstone.csv", header=T, sep=","))
rownames(dat2) <- colnames(dat2)
dat2
est.faccon(data=dat2, nfactors=1:4, data.type="cor", method="ML", ns=437, path="Z:/data/factor/faccon/")
est.faccon(data=dat2, nfactors=1:4, data.type="cor", method="ML", ns=437, sort=T, path="Z:/data/factor/faccon/")
est.faccon(data=dat2, nfactors=1:4, data.type="cor", method="ULS", ns=437, sort=T, path="Z:/data/factor/faccon/")


#------------------------------------#
#忍者ハットリ君 holizinger P 
#------------------------------------#
dat3 <- read.table("http://blue.zero.jp/yokumura/R/factor/holizinger.csv", header=T, sep=",")
dat3
est.faccon(data=dat3[,-1], nfactors=1:11, data.type="raw", method="ML", path="Z:/data/factor/faccon/")
est.faccon(data=dat3[,-1], nfactors=1:11, data.type="raw", method="ML", sort=T, path="Z:/data/factor/faccon/")
est.faccon(data=dat3[,-1], nfactors=1:11, data.type="raw", method="ULS", sort=T, path="Z:/data/factor/faccon/")


#------------------------------------#
#忍者ハットリ君 服部a 
#------------------------------------#
dat4 <- read.table("http://blue.zero.jp/yokumura/R/factor/hattori.csv", header=T, sep=",")
dat4
est.faccon(data=dat4, nfactors=1:16, data.type="raw", method="ML", path="Z:/data/factor/faccon/")
est.faccon(data=dat4, nfactors=1:16, data.type="raw", method="ML", sort=T, path="Z:/data/factor/faccon/")
est.faccon(data=dat4, nfactors=1:16, data.type="raw", method="ULS", sort=T, path="Z:/data/factor/faccon/")