rm(list=ls()) library(OpenMx) load("C:/Users/u1255320/Dropbox/Phd/Models/OpenMx/thirtyfixationswithdummy.RData") DataRaw <- mxData(observed=test2, type="raw") manifests <- c('F1', 'F2', 'F3', 'F4', 'F5', 'F6', 'F7', 'F8', 'F9', 'F10', 'F11' ,'F12', 'F13', 'F14' , 'F15' , 'F16' , 'F17' ,'F18' ,'F19' ,'F20', 'F21' ,'F22','F23' ,'F24' ,'F25' ,'F26' ,'F27' ,'F28' ,'F29' ,'F30', 'D22','D23' ,'D24' ,'D25' ,'D26' ,'D27' ,'D28' ,'D29' ,'D30') latents <- c('lx1', 'lx2', 'lx3', 'lx4', 'lx5', 'lx6', 'lx7', 'lx8', 'lx9', 'lx10', 'lx11' ,'lx12', 'lx13', 'lx14' , 'lx15' , 'lx16' , 'lx17' ,'lx18' ,'lx19' ,'lx20', 'lx21' ,'lx22','lx23' ,'lx24' ,'lx25' ,'lx26' ,'lx27' ,'lx28' ,'lx29' ,'lx30', 'dx2', 'dx3', 'dx4', 'dx5', 'dx6', 'dx7', 'dx8', 'dx9', 'dx10', 'dx11' ,'dx12', 'dx13', 'dx14' , 'dx15' , 'dx16' , 'dx17' ,'dx18' ,'dx19' ,'dx20', 'dx21' ,'dx22','dx23' ,'dx24' ,'dx25' ,'dx26' ,'dx27' ,'dx28' ,'dx29' ,'dx30', 'b0', 'b1', 'b2') #Residual variances Path1 <- mxPath(from=c('F1', 'F2', 'F3', 'F4', 'F5', 'F6', 'F7', 'F8', 'F9', 'F10', 'F11' ,'F12', 'F13', 'F14' , 'F15' , 'F16' , 'F17' ,'F18' ,'F19' ,'F20', 'F21' ,'F22','F23' ,'F24' ,'F25' ,'F26' ,'F27' ,'F28' ,'F29' ,'F30'), arrows=2, labels=c('ve', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've', 've') ) #Latent True Scores Path2 <- mxPath(from=c('lx1', 'lx2', 'lx3', 'lx4', 'lx5', 'lx6', 'lx7', 'lx8', 'lx9','lx10', 'lx11' ,'lx12', 'lx13', 'lx14' , 'lx15' , 'lx16' , 'lx17' ,'lx18' ,'lx19' ,'lx20', 'lx21' ,'lx22','lx23' ,'lx24' ,'lx25' ,'lx26' ,'lx27' ,'lx28' ,'lx29' ,'lx30'), to=c('F1', 'F2', 'F3', 'F4', 'F5', 'F6', 'F7', 'F8', 'F9', 'F10', 'F11' ,'F12', 'F13', 'F14' , 'F15' , 'F16' , 'F17' ,'F18' ,'F19' ,'F20', 'F21' ,'F22','F23' ,'F24' ,'F25' ,'F26' ,'F27' ,'F28' ,'F29' ,'F30'), free=FALSE, values=1 ) #Autoregressive Paths Path3 <- mxPath(from=c('lx1', 'lx2', 'lx3', 'lx4', 'lx5', 'lx6', 'lx7', 'lx8', 'lx9', 'lx10', 'lx11' ,'lx12', 'lx13', 'lx14' , 'lx15' , 'lx16' , 'lx17' ,'lx18' ,'lx19' ,'lx20', 'lx21' ,'lx22','lx23' ,'lx24' ,'lx25' ,'lx26' ,'lx27' ,'lx28' ,'lx29'), to=c('lx2', 'lx3', 'lx4', 'lx5', 'lx6', 'lx7', 'lx8', 'lx9', 'lx10', 'lx11' ,'lx12', 'lx13', 'lx14' , 'lx15' , 'lx16' , 'lx17' ,'lx18' ,'lx19' ,'lx20', 'lx21' ,'lx22','lx23' ,'lx24' ,'lx25' ,'lx26' ,'lx27' ,'lx28' ,'lx29','lx30'), arrows=1, free=FALSE, values=1 ) #Latent change scores Path4 <- mxPath(from=c('dx2', 'dx3', 'dx4', 'dx5', 'dx6', 'dx7', 'dx8', 'dx9', 'dx10', 'dx11' ,'dx12', 'dx13', 'dx14' , 'dx15' , 'dx16' , 'dx17' ,'dx18' ,'dx19' ,'dx20', 'dx21' ,'dx22','dx23' ,'dx24' ,'dx25' ,'dx26' ,'dx27' ,'dx28' ,'dx29','dx30'), to=c('lx2', 'lx3', 'lx4', 'lx5', 'lx6', 'lx7', 'lx8', 'lx9', 'lx10', 'lx11' ,'lx12', 'lx13', 'lx14' , 'lx15' , 'lx16' , 'lx17' ,'lx18' ,'lx19' ,'lx20', 'lx21' ,'lx22','lx23' ,'lx24' ,'lx25' ,'lx26' ,'lx27' ,'lx28' ,'lx29','lx30'), arrows=1, free=FALSE, values=1 ) #Factor loadings Path5 <- mxPath(from="b0", to="lx1", arrows=1, free=FALSE, values=1) Path6 <- mxPath(from="b1", to=c('dx2', 'dx3', 'dx4', 'dx5', 'dx6', 'dx7', 'dx8', 'dx9', 'dx10', 'dx11' ,'dx12', 'dx13', 'dx14' , 'dx15' , 'dx16' , 'dx17' ,'dx18' ,'dx19' ,'dx20', 'dx21' ,'dx22','dx23' ,'dx24' ,'dx25' ,'dx26' ,'dx27' ,'dx28' ,'dx29' ,'dx30'), arrows=1, free=FALSE, labels=c('Loadb1[1,1]','Loadb1[2,1]','Loadb1[3,1]','Loadb1[4,1]', 'Loadb1[5,1]','Loadb1[6,1]','Loadb1[7,1]','Loadb1[8,1]', 'Loadb1[9,1]','Loadb1[10,1]','Loadb1[11,1]','Loadb1[12,1]', 'Loadb1[13,1]','Loadb1[14,1]','Loadb1[15,1]','Loadb1[16,1]', 'Loadb1[17,1]','Loadb1[18,1]','Loadb1[19,1]','Loadb1[20,1]', 'Loadb1[21,1]','Loadb1[22,1]','Loadb1[23,1]','Loadb1[24,1]', 'Loadb1[25,1]','Loadb1[26,1]','Loadb1[27,1]','Loadb1[28,1]', 'Loadb1[29,1]') ) Path7 <- mxPath(from="b2", to=c('dx2', 'dx3', 'dx4', 'dx5', 'dx6', 'dx7', 'dx8', 'dx9', 'dx10', 'dx11' ,'dx12', 'dx13', 'dx14' , 'dx15' , 'dx16' , 'dx17' ,'dx18' ,'dx19' ,'dx20', 'dx21' ,'dx22','dx23' ,'dx24' ,'dx25' ,'dx26' ,'dx27' ,'dx28' ,'dx29' ,'dx30'), arrows=1, free=FALSE, labels=c('Loadb2[1,1]','Loadb2[2,1]','Loadb2[3,1]','Loadb2[4,1]', 'Loadb2[5,1]','Loadb2[6,1]','Loadb2[7,1]','Loadb2[8,1]', 'Loadb2[9,1]','Loadb2[10,1]','Loadb2[11,1]','Loadb2[12,1]', 'Loadb2[13,1]','Loadb2[14,1]','Loadb2[15,1]','Loadb2[16,1]', 'Loadb2[17,1]','Loadb2[18,1]','Loadb2[19,1]','Loadb2[20,1]', 'Loadb2[21,1]','Loadb2[22,1]','Loadb2[23,1]','Loadb2[24,1]', 'Loadb2[25,1]','Loadb2[26,1]','Loadb2[27,1]','Loadb2[28,1]', 'Loadb2[29,1]') ) #Latent variable Covariances Path8 <- mxPath(from='b0', arrows=2, free=TRUE ) Path9 <- mxPath(from='b0', to='b1', arrows=2, free=TRUE, labels='c_01') Path10 <- mxPath(from='b0', to='b2', arrows=2, free=TRUE, labels='c_02') Path11 <- mxPath(from='b1', arrows=2, free=TRUE, labels='V_1') Path12 <- mxPath(from='b2', arrows=2, free=TRUE, labels='V_2') Path13 <- mxPath(from='b1', to='b2', arrows=2, free=TRUE, labels='c_12') #Latent variables means Path14 <- mxPath(from='one', to=c('b0','b1','b2'), arrows=1, free=c(TRUE,TRUE,FALSE), values=c(0,3,0), labels=c('b_0','mu_1','b_2') ) #mean of measured variables Path15 <- mxPath(from='one', to=c('F1', 'F2', 'F3', 'F4', 'F5', 'F6', 'F7', 'F8', 'F9', 'F10', 'F11' ,'F12', 'F13', 'F14' , 'F15' , 'F16' , 'F17' ,'F18' ,'F19' ,'F20', 'F21' ,'F22','F23' ,'F24' ,'F25' ,'F26' ,'F27' ,'F28' ,'F29' ,'F30'), arrows=1, free=FALSE, values=0 ) #Latent variables means Path16 <- mxPath(from='one', to=c('lx1', 'lx2', 'lx3', 'lx4', 'lx5', 'lx6', 'lx7', 'lx8', 'lx9', 'lx10', 'lx11' ,'lx12', 'lx13', 'lx14' , 'lx15' , 'lx16' , 'lx17' ,'lx18' ,'lx19' ,'lx20', 'lx21' ,'lx22','lx23' ,'lx24' ,'lx25' ,'lx26' ,'lx27' ,'lx28' ,'lx29' ,'lx30'), arrows=1, free=FALSE, values=0 ) Path17 <- mxPath(from='one', to=c('dx2', 'dx3', 'dx4', 'dx5', 'dx6', 'dx7', 'dx8', 'dx9', 'dx10', 'dx11' ,'dx12', 'dx13', 'dx14' , 'dx15' , 'dx16' , 'dx17' ,'dx18' ,'dx19' ,'dx20', 'dx21' ,'dx22','dx23' ,'dx24' ,'dx25' ,'dx26' ,'dx27' ,'dx28' ,'dx29' ,'dx30'), arrows=1, free=FALSE, values=0 ) #Latent variable Covariances Path18 <- mxPath(from=c('lx1', 'lx2', 'lx3', 'lx4', 'lx5', 'lx6', 'lx7', 'lx8', 'lx9', 'lx10', 'lx11' ,'lx12', 'lx13', 'lx14' , 'lx15' , 'lx16' , 'lx17' ,'lx18' ,'lx19' ,'lx20', 'lx21' ,'lx22','lx23' ,'lx24' ,'lx25' ,'lx26' ,'lx27' ,'lx28' ,'lx29' ,'lx30'), arrows=2, free=FALSE, values=0 ) Path19 <- mxPath(from=c('dx2', 'dx3', 'dx4', 'dx5', 'dx6', 'dx7', 'dx8', 'dx9', 'dx10', 'dx11' ,'dx12', 'dx13', 'dx14' , 'dx15' , 'dx16' , 'dx17' ,'dx18' ,'dx19' ,'dx20', 'dx21' ,'dx22','dx23' ,'dx24' ,'dx25' ,'dx26' ,'dx27' ,'dx28' ,'dx29' ,'dx30'), arrows=2, free=FALSE, values=0 ) Path20 <- mxPath(from=c( 'D22','D23' ,'D24' ,'D25' ,'D26' ,'D27' ,'D28' ,'D29' ,'D30'), to=c( 'F22','F23' ,'F24' ,'F25' ,'F26' ,'F27' ,'F28' ,'F29' ,'F30'), arrows=1, free=TRUE) ninefixations <- mxModel("Mplus copy", type="RAM", DataRaw, manifestVars=manifests, latentVars=latents, Path1, Path2, Path3, Path4, Path5, Path6, Path7, Path8, Path9, Path10, Path11, Path12, Path13, Path14, Path15, Path16, Path17, Path18, Path19, Path20, mxMatrix('Full',29,1,free=TRUE,values=3, labels='mu_1', name='mu1'), mxMatrix('Full',29,1,free=TRUE,values=0.6, labels='mu_2', name='mu2'), mxMatrix('Full',29,1,free=FALSE, values=c(1,2,3,4,5,6,7,8,9,10, 11,12,13,14,15,16,17,18,19,20, 21,22,23,24,25,26,27,28,29), name='time' ), mxMatrix('Full',29,1,free=FALSE, values=c(1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1), name='unu' ), mxAlgebra(mu2*exp(-mu2*time),name='Loadb1'), mxAlgebra(mu1*exp(-mu2*time)*(unu-mu2*time),name='Loadb2') ) modelfit <- mxRun(ninefixations) modeltryhard <- mxTryHard(ninefixations) summary(modelfit)