library(psych)
install.packages("psych")
library(sem)
ls()
data(bfi)
head(bfi)

bfi[1:5,1:2]

bfi[,1:25]		# výběr dat

BF.upr<-na.omit(bfi[,1:25])	# vynechání chybějících hodnot

BF.c<-cov(BF.upr)

model.bfi<-specifyModel()
agree->A1, fa1, NA
agree->A2, fa2, NA
agree->A3, fa3, NA
agree->A4, fa4, NA
agree->A5, fa5, NA
consc->C1, fc1, NA
consc->C2, fc2, NA
consc->C3, fc3, NA
consc->C4, fc4, NA
consc->C5, fc5, NA
extra->E1, fe1, NA
extra->E2, fe2, NA
extra->E3, fe3, NA
extra->E4, fe4, NA
extra->E5, fe5, NA
neuro->N1, fn1, NA
neuro->N2, fn2, NA
neuro->N3, fn3, NA
neuro->N4, fn4, NA
neuro->N5, fn5, NA
opene->O1, fo1, NA
opene->O2, fo2, NA
opene->O3, fo3, NA
opene->O4, fo4, NA
opene->O5, fo5, NA
agree<->agree, NA, 1
consc<->consc, NA, 1
extra<->extra, NA, 1
neuro<->neuro, NA, 1
opene<->opene, NA, 1

BF.fit<-sem(model.bfi,BF.c,nrow(BF.upr))

summary(BF.fit)
stdCoef(BF.fit)
summary(BF.fit,fit.indices=c("GFI", "AGFI", "RMSEA", "NFI", "NNFI", "CFI", "RNI", "IFI", "SRMR", "AIC", "AICc", "BIC", "CAIC"))
modIndices(BF.fit)		# provést modifikaci modelu

######################################

hp<-read.table("HPSQ.txt",sep="\t",head=T)
head(hs)
hp$Grade.1
names(hp)
pc<-princomp(hp[,2:11],cor=T)
plot(pc,type="lines")
pc<-principal(hp[,2:11],nfactors=3,rotate="varimax")

pt<-princomp(hp[,19:28],cor=T)
plot(pt,type="lines")
pt<-principal(hp[,19:28],nfactors=3,rotate="varimax")

library(sem)
hpc.c<-cov(hp[,2:11])

mod.hpc<-specifyModel()
leg->p1HPSQC, fa1, NA
leg->p2HPSQC, fa2, NA
pet->p3HPSQC, fa3, NA
pet->p4HPSQC, fa4, NA
pew->p5HPSQC, fa5, NA
pew->p6HPSQC, fa6, NA
pew->p7HPSQC, fa7, NA
pew->p8HPSQC, fa8, NA
pet->p9HPSQC, fa9, NA
leg->p10HPSQC, fa10, NA
leg<->leg, NA, 1
pet<->pet, NA, 1
pew<->pew, NA, 1

hpc.fit<-sem(mod.hpc,hpc.c,nrow(hp))

summary(hpc.fit)
stdCoef(hpc.fit)
summary(hpc.fit,fit.indices=c("GFI", "AGFI", "RMSEA", "NFI", "NNFI", "CFI", "RNI", "IFI", "SRMR", "AIC", "AICc", "BIC", "CAIC"))
modIndices(hpc.fit)		# provést modifikaci modelu