# Prepare Data # ----------------------------------------------------------------------------- require(OpenMx) dataskripsi<-read.table("D:/SEM/Data/My Data/DATA FIX/DATA SKRIPSI.csv",header=T,sep=";") data1<-cov(dataskripsi) theLatent <- c("SKL","Nilai", "Proses","Isi","PTK","Pengelolaan") FactorModelskripsi<- mxModel("Factor Model Skripsi Path", type="RAM", mxData( observed=data1, type="cov", numObs=3953 ), manifestVars=c ("y1","y2","y3","y4","y5","y6","y7","y8","y9","y10","y11","y12","y13","y14","y15","y16","y17","y18","y19","y20","y21","y22","y23","y24","x1","x2","x3","x4","x5","x6","x7"), latentVars=theLatent, mxPath( from=c ("y1","y2","y3","y4","y5","y6","y7","y8","y9","y10","y11","y12","y13","y14","y15","y16","y17","y18","y19","y20","y21","y22","y23","y24","x1","x2","x3","x4","x5","x6","x7"), 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), 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 laten # ------------------------------------- mxPath( from=theLatent, arrows=2, free=TRUE, values=1, labels=paste("var", theLatent, sep="") ), # gamma dan beta # ------------------------------------- mxPath( from=c("Pengelolaan","PTK", "PTK","PTK","Isi","Proses","Isi","Nilai","Proses"), to=c("PTK","Isi","Proses","Nilai","Proses","Nilai","SKL","SKL","SKL"), arrows=1, values=.9, labels=paste("reg", c("Pengelolaan","PTK", "PTK","PTK","Isi","Proses","Isi","Nilai","Proses"), "_to_", c ("PTK","Isi","Proses","Nilai","Proses","Nilai","SKL","SKL","SKL"), sep="") ), # lambda # ------------------------------------- mxPath( from="SKL", to=c("y1","y2","y3","y4","y5"), arrows=1, free=c(TRUE,TRUE,TRUE,TRUE,TRUE), values=c(1,0,0,0,0), labels=c("l1","l2","l3","l4","l5") ), #latent variances and covaraince # ------------------------------------- mxPath( from="Isi", to=c("y6","y7","y8","y9"), arrows=1, free=c(TRUE,TRUE,TRUE,TRUE), values=c(1,0,0,0), labels=c("l6","l7","l8","l9") ), #latent variances and covaraince # ------------------------------------- mxPath( from="Proses", to=c("y10","y11","y12","y13","y14","y15","y16","y17"), arrows=1, free=c(TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE), values=c(1,0,0,0,0,0,0,0), labels=c("l10","l11","l12","l13","l14","l15","l16","l17") ), #latent variances and covaraince # ------------------------------------- mxPath( from="Nilai", to=c("y18","y19","y20","y21","y22"), arrows=1, free=c(TRUE,TRUE,TRUE,TRUE,TRUE), values=c(1,0,0,0,0), labels=c("l18","l19","l20","l21","l22") ), # factor loadings for y variables # ------------------------------------- mxPath( from="PTK", to=c("y23","y24"), arrows=1, free=c(TRUE,TRUE), values=c(1,0), labels=c("l23","l24") ), # factor loadings for x variables # ------------------------------------- mxPath( from="Pengelolaan", to=c("x1","x2","x3","x4","x5","x6","x7"), arrows=1, free=c(TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,TRUE), values=c(1,0,0,0,0,0,0), labels=c("l25","l26","l27","l28","l29","l30","l31") ) ) # close model # Create an MxModel object # ----------------------------------------------------------------------------- twoFactorFitskripsi<- mxRun(FactorModelskripsi) summary(twoFactorFitskripsi) require(semPlot) semPaths(twoFactorFitskripsi, layout = "tree2") #require(semTools) #fitMeasuresMx(twoFactorFitskripsi) output<-summary(twoFactorFitskripsi) name<-output[1]$parameters[,1] col<-output[1]$parameters[,4] est<-output[1]$parameters[,5] ste<-output[1]$parameters[,6] tvalue<-est/ste N<-2916 pvalue<-2*(1-pt(abs(tvalue),N-1)) newtable<-cbind(est,ste,tvalue,pvalue) rownames(newtable)<-name newtable