# Prepare Data
# -----------------------------------------------------------------------------
require(OpenMx)
contoh<-read.table("D:/SEM/Data/My Data/DATA BARU/data YogyaPabar.csv",header=T,sep=";")
#rataan1<-read.table("D:/SEM/Data/My Data/data rataan.csv",header=T,sep=";")
data1<-cov(contoh)
#rataan<-as.numeric(rataan1)
#names(rataan)<- c("x1", "x2", "x3","x4","x5","x6","x7","x8","x9", "y1", "y2", "y3", "y4", "y5","y6","y7","y8","y9","y10","y11","y12","y13","y14","y15","y16","y17","y18","y19","y20","y21","y22")
#data<-as.numeric(data1)
#names(data)<-(c("x1", "x2", "x3","x4","x5","x6","x7","x8","x9", "y1", "y2", "y3", "y4", "y5","y6","y7","y8","y9","y10","y11","y12","y13","y14","y15","y16","y17","y18","y19","y20","y21","y22"),c("x1", "x2", "x3","x4","x5","x6","x7","x8","x9", "y1", "y2", "y3", "y4", "y5","y6","y7","y8","y9","y10","y11","y12","y13","y14","y15","y16","y17","y18","y19","y20","y21","y22"))
theLat <- c("SKL","Nilai", "Proses","Isi","PTK","Pengelolaan")
FactorModelcontoh<- mxModel("Factor Model contoh Path", type="RAM",
    mxData(
	observed=data1, 
    	type="cov", 
    	numObs=2245
    ),
    manifestVars=c("x1", "x2", "x3","x4","x5","x6","x7","x8","x9", "y1", "y2", "y3", "y4", "y5","y6","y7","y8","y9","y10", "y11","y12","y13","y14","y15","y16","y17","y18","y19","y20","y21","y22"),
    latentVars=theLat,
    mxPath(
    	from=c("x1","x2","x3","x4","x5","x6","x7","x8","x9","y1","y2","y3","y4","y5","y6","y7","y8","y9","y10", "y11","y12","y13","y14","y15","y16","y17","y18","y19","y20","y21","y22"),
        arrows=2,
        free=c(TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE),
        values=c(.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5),
        labels=c("e1","e2","e3","e4","e5","e6","e7","e8","e9","e10","e11","e12","e13","e14","e15","e16","e17","e18","e19","e20","e21","e22","e23","e24","e25","e26","e27","e28","e29","e30","e31")
    ),
	# residual variances
	# -------------------------------------
    mxPath(
    	from=theLat,
        arrows=2,
        free=TRUE,
        values=1,
        labels=paste("var", theLat, sep="")
    ),
	# latent variances and covaraince
	# -------------------------------------
    mxPath(
    	from=c("PTK", "PTK", "Pengelolaan", "Pengelolaan", "Isi", "Isi", "Proses", "Proses", "Nilai"),
        to=c("SKL", "Isi", "SKL", "Isi", "SKL", "Proses", "SKL", "Nilai", "SKL"),
        arrows=1,
        values=.7,
        labels=paste("reg", c("PTK", "PTK", "Pengelolaan", "Pengelolaan", "Isi", "Isi", "Proses", "Proses", "Nilai"), "_to_", c("SKL", "Isi", "SKL", "Isi", "SKL", "Proses", "SKL", "Nilai", "SKL"), sep="")
    ),
	# Latent regressions
	# -------------------------------------
    mxPath(
    	from="SKL",
        to=c("y1","y2","y3","y4","y5"),
        arrows=1,
        free=c(TRUE,TRUE,TRUE,TRUE,TRUE),
        values=c(.5,.5,.5,.5,.5),
        labels=c("l1","l2","l3","l4","l5")
    ),
	#latent variances and covaraince
	# -------------------------------------
    mxPath(
    	from="Nilai",
        to=c("y6","y7","y8","y9","y10"),
        arrows=1,
        free=c(TRUE,TRUE,TRUE,TRUE,TRUE),
        values=c(.4,.4,.4,.4,.4),
        labels=c("l6","l7","l8","l9","l10")
    ),

	#latent variances and covaraince
	# -------------------------------------
    mxPath(
    	from="Proses",
        to=c("y11","y12","y13","y14","y15","y16","y17","y18"),
        arrows=1,
        free=c(TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE),
        values=c(.5,.5,.5,.5,.5,.5,.5,.5),
        labels=c("l11","l12","l13","l14","l15","l16","l17","l18")
    ),
	#latent variances and covaraince
	# -------------------------------------
    mxPath(
    	from="Isi",
        to=c("y19","y20","y21","y22"),
        arrows=1,
        free=c(TRUE,TRUE,TRUE,TRUE),
        values=c(.5,.5,.5,.5),
        labels=c("l19","l20","l21","l22")
    ),

	# factor loadings for x variables
	# -------------------------------------
    mxPath(
    	from="PTK",
        to=c("x1","x2"),
        arrows=1,
        free=c(TRUE,TRUE),
        values=c(.5,.5),
        labels=c("l23","l24")
    ),
	# factor loadings for x variables
	# -------------------------------------
    mxPath(
	from="Pengelolaan",
        to=c("x3","x4","x5","x6","x7","x8","x9"),
        arrows=1,
        free=c(TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE),
        values=c(.5,.5,.5,.5,.5,.5,.5),
        labels=c("l25","l26","l27","l28","l29","l30","l31")
    )
) # close model
# Create an MxModel object
# -----------------------------------------------------------------------------

twoFactorFitcontoh<- mxRun(FactorModelcontoh)

summary(twoFactorFitcontoh)

require(semPlot)
semPaths(twoFactorFitcontoh, layout = "tree2")