//LG5.0//
version = 5.0
infile 'rep1_N150_D7_t8.sim'

model
options
   algorithm 
      tolerance=1e-008 emtolerance=0.01 emiterations=300 nriterations=200;
   startvalues
      seed=0 sets=100 tolerance=1e-005 iterations=50;
   bayes
      categorical=1 variances=1 latent=0 poisson=1;
   montecarlo
      seed=0 replicates=500 tolerance=1e-008;
   quadrature  nodes=10;
   missing  includeall;
   output      
      parameters=first standarderrors write  'rep1_N150_D7_t8_CT.csv';
variables
  caseid ID;  
  timeinterval Time100;  
  dependent well ordinal 3, good ordinal 3;  
independent Newday nominal, 
Notnewday nominal, 
Firsthour nominal, 
Laterhour nominal,
Time100 inactive;
latent DClass nominal dynamic 2, 
State nominal dynamic 3;
equations
DClass[=0] <- (a) 1;  
State[=0] <- (b) 1; 
well <- (c~err) 1 | State + DClass;	
good <- (d~err) 1 | State + DClass;	
DClass <- (e~tra) newday | DClass[-1] + (f~tra) notnewday | DClass[-1];	
State <- (g~tra) newday | State[-1] + (h~tra) firsthour | State[-1] DClass + (i~tra~int) laterhour | State[-1] DClass; 	

b = -15;
c[1,2] = -15;
c[3,1] = -15;
d[1,2] = -15;
d[3,1] = -15;
f = -15;	
g = -15;   	
g[2,1] = 15;	
g[3,1] = 15;	
h[3] = -15;	
h[4] = -15;	
h[5] = -15;	
h[6] = -15;	
i[1,2] = -15;
i[2,2] = -15;
i[5,1] = -15;
i[6,1] = -15;	 

{
-0.72
-2.37
-2.93
-3.6
-0.56
1.15
-2.72
-3.56
-3.59
-0.31
1.91
-2.19
-2.28
0.48
-2.1
0.2
-0.12
-0.35
0.47
-2
-1.97
-0.89
-0.68
-0.27
-0.88
}
end model