Skip to contents

Estimates the parameters of a model created using the slca function.

Usage

estimate(x, ...)

# S3 method for class 'slca'
estimate(x,
   data,
   method = c("em", "hybrid", "nlm"),
   fix2zero = NULL,
   control = slcaControl(), ...)

Arguments

x

an slca object defining the slca model to be estimated.

...

additional arguments passed to the estimation process.

data

a data.frame containing the observed categorical variables included in the model.

method

a character string specifying the estimation method for SLCM parameters. The default is "em", which uses the expectation-maximization (EM) algorithm. The alternative "nlm" employs the Newton-Raphson algorithm via the nlm function, while "hybrid" combines both approaches, starting with EM and finishing with nlm for refined estimates.

fix2zero

a vector specifying parameters to be constrained to zero. See the 'Details' section for further information.

control

a list of control parameters for the estimation procedure. Modify default values using slcaControl().

Value

An object of class slcafit containing the following components:

model

a list describing of the model structure.

method

the estimation method used.

arg

a brief description of the model used during estimation.

mf

the data.frame used for estimation.

par

the log of the estimated paramters.

logit

the log-odds of the estimated parameters.

score

the score function for the estimated parameters.

posterior

a list of posterior probablities for each latent class variable.

convergence

a logical indicator of whether convergence was achieved.

loglikelihood

the loglikelihood value of the estimated model.

control

the control settings used during the estimation process.

The returned object can be further processed using the param function to extract the estimated parameters or their standard errors. The regress function allows for logistic regression analysis using a three-step approach to evaluate the effects of external variables on latent class variables. Additionally, several other methods are available, including predict.slcafit, reorder.slcafit, gof, and others.

Details

The fix2zero argument allows you to constrain specific parameters to zero. Each parameter is associated with a unique index, which can be identified using the param function with the argument index = TRUE. To apply constraints, provide the relevant parameter indices in the fix2zero arguments with vector.

Examples

m <- slca(lc[3] ~ y1 + y2 + y3 + y4)
pi <- rep(1 / 3, 3)
rho <- c(.9, .1, .9, .1, .9, .1, .9, .1, # class 1
         .9, .1, .9, .1, .1, .9, .1, .9, # class 2
         .1, .9, .1, .9, .1, .9, .1, .9) # class 3
dt <- simulate(m, 200, parm = c(pi, rho))
estimate(m, dt$response)
#> Structural Latent Variable Model
#> 
#> Latent variables (Root*):           
#>  Label: lc*
#> nclass: 3  
#> 
#> Measurement model:                            
#>  lc -> { y1, y2, y3, y4 }  a
#> 
#> LogLik: -396.1722

# Several estimation methods
estimate(m, dt$response, method = "em",
         control = slcaControl(verbose = TRUE)) # default
#> EM iteration begin.
#> 1 iterations  logLik: -555.07  diff: inf      
2 iterations  logLik: -511.98  diff: 43.08609      
3 iterations  logLik: -509.90  diff: 2.08325      
4 iterations  logLik: -500.10  diff: 9.80122      
5 iterations  logLik: -470.22  diff: 29.88096      
6 iterations  logLik: -437.56  diff: 32.65298      
7 iterations  logLik: -423.55  diff: 14.01664      
8 iterations  logLik: -411.63  diff: 11.91095      
9 iterations  logLik: -401.67  diff: 9.96603      
10 iterations  logLik: -398.19  diff: 3.48230      
11 iterations  logLik: -397.69  diff: 0.49546      
12 iterations  logLik: -397.64  diff: 0.04873      
13 iterations  logLik: -397.63  diff: 0.01264      
14 iterations  logLik: -397.62  diff: 0.01098      
15 iterations  logLik: -397.61  diff: 0.01176      
16 iterations  logLik: -397.59  diff: 0.01279      
17 iterations  logLik: -397.58  diff: 0.01394      
18 iterations  logLik: -397.56  diff: 0.01520      
19 iterations  logLik: -397.55  diff: 0.01658      
20 iterations  logLik: -397.53  diff: 0.01807      
21 iterations  logLik: -397.51  diff: 0.01969      
22 iterations  logLik: -397.49  diff: 0.02143      
23 iterations  logLik: -397.47  diff: 0.02328      
24 iterations  logLik: -397.44  diff: 0.02523      
25 iterations  logLik: -397.41  diff: 0.02726      
26 iterations  logLik: -397.38  diff: 0.02936      
27 iterations  logLik: -397.35  diff: 0.03149      
28 iterations  logLik: -397.32  diff: 0.03361      
29 iterations  logLik: -397.28  diff: 0.03569      
30 iterations  logLik: -397.25  diff: 0.03767      
31 iterations  logLik: -397.21  diff: 0.03952      
32 iterations  logLik: -397.16  diff: 0.04117      
33 iterations  logLik: -397.12  diff: 0.04259      
34 iterations  logLik: -397.08  diff: 0.04373      
35 iterations  logLik: -397.03  diff: 0.04456      
36 iterations  logLik: -396.99  diff: 0.04507      
37 iterations  logLik: -396.94  diff: 0.04522      
38 iterations  logLik: -396.90  diff: 0.04504      
39 iterations  logLik: -396.85  diff: 0.04453      
40 iterations  logLik: -396.81  diff: 0.04371      
41 iterations  logLik: -396.77  diff: 0.04261      
42 iterations  logLik: -396.73  diff: 0.04127      
43 iterations  logLik: -396.69  diff: 0.03972      
44 iterations  logLik: -396.65  diff: 0.03801      
45 iterations  logLik: -396.61  diff: 0.03617      
46 iterations  logLik: -396.58  diff: 0.03424      
47 iterations  logLik: -396.55  diff: 0.03226      
48 iterations  logLik: -396.52  diff: 0.03027      
49 iterations  logLik: -396.49  diff: 0.02828      
50 iterations  logLik: -396.46  diff: 0.02633      
51 iterations  logLik: -396.44  diff: 0.02443      
52 iterations  logLik: -396.41  diff: 0.02259      
53 iterations  logLik: -396.39  diff: 0.02084      
54 iterations  logLik: -396.37  diff: 0.01918      
55 iterations  logLik: -396.36  diff: 0.01761      
56 iterations  logLik: -396.34  diff: 0.01613      
57 iterations  logLik: -396.33  diff: 0.01476      
58 iterations  logLik: -396.31  diff: 0.01348      
59 iterations  logLik: -396.30  diff: 0.01229      
60 iterations  logLik: -396.29  diff: 0.01120      
61 iterations  logLik: -396.28  diff: 0.01020      
62 iterations  logLik: -396.27  diff: 0.00928      
63 iterations  logLik: -396.26  diff: 0.00843      
64 iterations  logLik: -396.25  diff: 0.00766      
65 iterations  logLik: -396.25  diff: 0.00696      
66 iterations  logLik: -396.24  diff: 0.00632      
67 iterations  logLik: -396.23  diff: 0.00574      
68 iterations  logLik: -396.23  diff: 0.00521      
69 iterations  logLik: -396.22  diff: 0.00473      
70 iterations  logLik: -396.22  diff: 0.00429      
71 iterations  logLik: -396.22  diff: 0.00390      
72 iterations  logLik: -396.21  diff: 0.00354      
73 iterations  logLik: -396.21  diff: 0.00322      
74 iterations  logLik: -396.21  diff: 0.00292      
75 iterations  logLik: -396.20  diff: 0.00266      
76 iterations  logLik: -396.20  diff: 0.00242      
77 iterations  logLik: -396.20  diff: 0.00220      
78 iterations  logLik: -396.20  diff: 0.00200      
79 iterations  logLik: -396.20  diff: 0.00183      
80 iterations  logLik: -396.19  diff: 0.00167      
81 iterations  logLik: -396.19  diff: 0.00152      
82 iterations  logLik: -396.19  diff: 0.00139      
83 iterations  logLik: -396.19  diff: 0.00127      
84 iterations  logLik: -396.19  diff: 0.00116      
85 iterations  logLik: -396.19  diff: 0.00106      
86 iterations  logLik: -396.19  diff: 0.00097      
87 iterations  logLik: -396.19  diff: 0.00089      
88 iterations  logLik: -396.18  diff: 0.00082      
89 iterations  logLik: -396.18  diff: 0.00075      
90 iterations  logLik: -396.18  diff: 0.00069      
91 iterations  logLik: -396.18  diff: 0.00063      
92 iterations  logLik: -396.18  diff: 0.00058      
93 iterations  logLik: -396.18  diff: 0.00054      
94 iterations  logLik: -396.18  diff: 0.00050      
95 iterations  logLik: -396.18  diff: 0.00046      
96 iterations  logLik: -396.18  diff: 0.00042      
97 iterations  logLik: -396.18  diff: 0.00039      
98 iterations  logLik: -396.18  diff: 0.00036      
99 iterations  logLik: -396.18  diff: 0.00034      
100 iterations  logLik: -396.18  diff: 0.00031      
101 iterations  logLik: -396.18  diff: 0.00029      
102 iterations  logLik: -396.18  diff: 0.00027      
103 iterations  logLik: -396.18  diff: 0.00025      
104 iterations  logLik: -396.18  diff: 0.00024      
105 iterations  logLik: -396.18  diff: 0.00022      
106 iterations  logLik: -396.18  diff: 0.00021      
107 iterations  logLik: -396.18  diff: 0.00019      
108 iterations  logLik: -396.18  diff: 0.00018      
109 iterations  logLik: -396.18  diff: 0.00017      
110 iterations  logLik: -396.18  diff: 0.00016      
111 iterations  logLik: -396.18  diff: 0.00015      
112 iterations  logLik: -396.18  diff: 0.00014      
113 iterations  logLik: -396.18  diff: 0.00013      
114 iterations  logLik: -396.18  diff: 0.00013      
115 iterations  logLik: -396.18  diff: 0.00012      
116 iterations  logLik: -396.18  diff: 0.00011      
117 iterations  logLik: -396.18  diff: 0.00011      
118 iterations  logLik: -396.18  diff: 0.00010      
119 iterations  logLik: -396.18  diff: 0.00010      
120 iterations  logLik: -396.18  diff: 0.00009      
121 iterations  logLik: -396.18  diff: 0.00009      
122 iterations  logLik: -396.18  diff: 0.00008      
123 iterations  logLik: -396.17  diff: 0.00008      
124 iterations  logLik: -396.17  diff: 0.00007      
125 iterations  logLik: -396.17  diff: 0.00007      
126 iterations  logLik: -396.17  diff: 0.00007      
127 iterations  logLik: -396.17  diff: 0.00007      
128 iterations  logLik: -396.17  diff: 0.00006      
129 iterations  logLik: -396.17  diff: 0.00006      
130 iterations  logLik: -396.17  diff: 0.00006      
131 iterations  logLik: -396.17  diff: 0.00006      
132 iterations  logLik: -396.17  diff: 0.00005      
133 iterations  logLik: -396.17  diff: 0.00005      
134 iterations  logLik: -396.17  diff: 0.00005      
135 iterations  logLik: -396.17  diff: 0.00005      
136 iterations  logLik: -396.17  diff: 0.00005      
137 iterations  logLik: -396.17  diff: 0.00004      
138 iterations  logLik: -396.17  diff: 0.00004      
139 iterations  logLik: -396.17  diff: 0.00004      
140 iterations  logLik: -396.17  diff: 0.00004      
141 iterations  logLik: -396.17  diff: 0.00004      
142 iterations  logLik: -396.17  diff: 0.00004      
143 iterations  logLik: -396.17  diff: 0.00004      
144 iterations  logLik: -396.17  diff: 0.00003      
145 iterations  logLik: -396.17  diff: 0.00003      
146 iterations  logLik: -396.17  diff: 0.00003      
147 iterations  logLik: -396.17  diff: 0.00003      
148 iterations  logLik: -396.17  diff: 0.00003      
149 iterations  logLik: -396.17  diff: 0.00003      
150 iterations  logLik: -396.17  diff: 0.00003      
151 iterations  logLik: -396.17  diff: 0.00003      
152 iterations  logLik: -396.17  diff: 0.00003      
153 iterations  logLik: -396.17  diff: 0.00003      
154 iterations  logLik: -396.17  diff: 0.00003      
155 iterations  logLik: -396.17  diff: 0.00003      
156 iterations  logLik: -396.17  diff: 0.00002      
157 iterations  logLik: -396.17  diff: 0.00002      
158 iterations  logLik: -396.17  diff: 0.00002      
159 iterations  logLik: -396.17  diff: 0.00002      
160 iterations  logLik: -396.17  diff: 0.00002      
161 iterations  logLik: -396.17  diff: 0.00002      
162 iterations  logLik: -396.17  diff: 0.00002      
163 iterations  logLik: -396.17  diff: 0.00002      
164 iterations  logLik: -396.17  diff: 0.00002      
165 iterations  logLik: -396.17  diff: 0.00002      
166 iterations  logLik: -396.17  diff: 0.00002      
167 iterations  logLik: -396.17  diff: 0.00002      
168 iterations  logLik: -396.17  diff: 0.00002      
169 iterations  logLik: -396.17  diff: 0.00002      
170 iterations  logLik: -396.17  diff: 0.00002      
171 iterations  logLik: -396.17  diff: 0.00002      
172 iterations  logLik: -396.17  diff: 0.00002      
173 iterations  logLik: -396.17  diff: 0.00002      
174 iterations  logLik: -396.17  diff: 0.00002      
175 iterations  logLik: -396.17  diff: 0.00002      
176 iterations  logLik: -396.17  diff: 0.00002      
177 iterations  logLik: -396.17  diff: 0.00002      
178 iterations  logLik: -396.17  diff: 0.00001      
179 iterations  logLik: -396.17  diff: 0.00001      
180 iterations  logLik: -396.17  diff: 0.00001      
181 iterations  logLik: -396.17  diff: 0.00001      
182 iterations  logLik: -396.17  diff: 0.00001      
183 iterations  logLik: -396.17  diff: 0.00001      
184 iterations  logLik: -396.17  diff: 0.00001      
185 iterations  logLik: -396.17  diff: 0.00001      
186 iterations  logLik: -396.17  diff: 0.00001      
187 iterations  logLik: -396.17  diff: 0.00001      
188 iterations  logLik: -396.17  diff: 0.00001      
189 iterations  logLik: -396.17  diff: 0.00001      
190 iterations  logLik: -396.17  diff: 0.00001      
191 iterations  logLik: -396.17  diff: 0.00001      
192 iterations  logLik: -396.17  diff: 0.00001      
193 iterations  logLik: -396.17  diff: 0.00001      
194 iterations  logLik: -396.17  diff: 0.00001      
195 iterations  logLik: -396.17  diff: 0.00001      
196 iterations  logLik: -396.17  diff: 0.00001      
197 iterations  logLik: -396.17  diff: 0.00001      
198 iterations  logLik: -396.17  diff: 0.00001      
199 iterations  logLik: -396.17  diff: 0.00001      
200 iterations  logLik: -396.17  diff: 9.8e-06      
201 iterations  logLik: -396.17  diff: 9.7e-06      
202 iterations  logLik: -396.17  diff: 9.5e-06      
203 iterations  logLik: -396.17  diff: 9.3e-06      
204 iterations  logLik: -396.17  diff: 9.2e-06      
205 iterations  logLik: -396.17  diff: 9.0e-06      
206 iterations  logLik: -396.17  diff: 8.9e-06      
207 iterations  logLik: -396.17  diff: 8.7e-06      
208 iterations  logLik: -396.17  diff: 8.6e-06      
209 iterations  logLik: -396.17  diff: 8.4e-06      
210 iterations  logLik: -396.17  diff: 8.3e-06      
211 iterations  logLik: -396.17  diff: 8.2e-06      
212 iterations  logLik: -396.17  diff: 8.0e-06      
213 iterations  logLik: -396.17  diff: 7.9e-06      
214 iterations  logLik: -396.17  diff: 7.8e-06      
215 iterations  logLik: -396.17  diff: 7.6e-06      
216 iterations  logLik: -396.17  diff: 7.5e-06      
217 iterations  logLik: -396.17  diff: 7.4e-06      
218 iterations  logLik: -396.17  diff: 7.3e-06      
219 iterations  logLik: -396.17  diff: 7.2e-06      
220 iterations  logLik: -396.17  diff: 7.1e-06      
221 iterations  logLik: -396.17  diff: 6.9e-06      
222 iterations  logLik: -396.17  diff: 6.8e-06      
223 iterations  logLik: -396.17  diff: 6.7e-06      
224 iterations  logLik: -396.17  diff: 6.6e-06      
225 iterations  logLik: -396.17  diff: 6.5e-06      
226 iterations  logLik: -396.17  diff: 6.4e-06      
227 iterations  logLik: -396.17  diff: 6.3e-06      
228 iterations  logLik: -396.17  diff: 6.2e-06      
229 iterations  logLik: -396.17  diff: 6.1e-06      
230 iterations  logLik: -396.17  diff: 6.0e-06      
231 iterations  logLik: -396.17  diff: 5.9e-06      
232 iterations  logLik: -396.17  diff: 5.9e-06      
233 iterations  logLik: -396.17  diff: 5.8e-06      
234 iterations  logLik: -396.17  diff: 5.7e-06      
235 iterations  logLik: -396.17  diff: 5.6e-06      
236 iterations  logLik: -396.17  diff: 5.5e-06      
237 iterations  logLik: -396.17  diff: 5.4e-06      
238 iterations  logLik: -396.17  diff: 5.3e-06      
239 iterations  logLik: -396.17  diff: 5.3e-06      
240 iterations  logLik: -396.17  diff: 5.2e-06      
241 iterations  logLik: -396.17  diff: 5.1e-06      
242 iterations  logLik: -396.17  diff: 5.0e-06      
243 iterations  logLik: -396.17  diff: 5.0e-06      
244 iterations  logLik: -396.17  diff: 4.9e-06      
245 iterations  logLik: -396.17  diff: 4.8e-06      
246 iterations  logLik: -396.17  diff: 4.7e-06      
247 iterations  logLik: -396.17  diff: 4.7e-06      
248 iterations  logLik: -396.17  diff: 4.6e-06      
249 iterations  logLik: -396.17  diff: 4.5e-06      
250 iterations  logLik: -396.17  diff: 4.5e-06      
251 iterations  logLik: -396.17  diff: 4.4e-06      
252 iterations  logLik: -396.17  diff: 4.3e-06      
253 iterations  logLik: -396.17  diff: 4.3e-06      
254 iterations  logLik: -396.17  diff: 4.2e-06      
255 iterations  logLik: -396.17  diff: 4.2e-06      
256 iterations  logLik: -396.17  diff: 4.1e-06      
257 iterations  logLik: -396.17  diff: 4.0e-06      
258 iterations  logLik: -396.17  diff: 4.0e-06      
259 iterations  logLik: -396.17  diff: 3.9e-06      
260 iterations  logLik: -396.17  diff: 3.9e-06      
261 iterations  logLik: -396.17  diff: 3.8e-06      
262 iterations  logLik: -396.17  diff: 3.8e-06      
263 iterations  logLik: -396.17  diff: 3.7e-06      
264 iterations  logLik: -396.17  diff: 3.7e-06      
265 iterations  logLik: -396.17  diff: 3.6e-06      
266 iterations  logLik: -396.17  diff: 3.6e-06      
267 iterations  logLik: -396.17  diff: 3.5e-06      
268 iterations  logLik: -396.17  diff: 3.5e-06      
269 iterations  logLik: -396.17  diff: 3.4e-06      
270 iterations  logLik: -396.17  diff: 3.4e-06      
271 iterations  logLik: -396.17  diff: 3.3e-06      
272 iterations  logLik: -396.17  diff: 3.3e-06      
273 iterations  logLik: -396.17  diff: 3.2e-06      
274 iterations  logLik: -396.17  diff: 3.2e-06      
275 iterations  logLik: -396.17  diff: 3.1e-06      
276 iterations  logLik: -396.17  diff: 3.1e-06      
277 iterations  logLik: -396.17  diff: 3.1e-06      
278 iterations  logLik: -396.17  diff: 3.0e-06      
279 iterations  logLik: -396.17  diff: 3.0e-06      
280 iterations  logLik: -396.17  diff: 2.9e-06      
281 iterations  logLik: -396.17  diff: 2.9e-06      
282 iterations  logLik: -396.17  diff: 2.9e-06      
283 iterations  logLik: -396.17  diff: 2.8e-06      
284 iterations  logLik: -396.17  diff: 2.8e-06      
285 iterations  logLik: -396.17  diff: 2.7e-06      
286 iterations  logLik: -396.17  diff: 2.7e-06      
287 iterations  logLik: -396.17  diff: 2.7e-06      
288 iterations  logLik: -396.17  diff: 2.6e-06      
289 iterations  logLik: -396.17  diff: 2.6e-06      
290 iterations  logLik: -396.17  diff: 2.6e-06      
291 iterations  logLik: -396.17  diff: 2.5e-06      
292 iterations  logLik: -396.17  diff: 2.5e-06      
293 iterations  logLik: -396.17  diff: 2.5e-06      
294 iterations  logLik: -396.17  diff: 2.4e-06      
295 iterations  logLik: -396.17  diff: 2.4e-06      
296 iterations  logLik: -396.17  diff: 2.4e-06      
297 iterations  logLik: -396.17  diff: 2.3e-06      
298 iterations  logLik: -396.17  diff: 2.3e-06      
299 iterations  logLik: -396.17  diff: 2.3e-06      
300 iterations  logLik: -396.17  diff: 2.3e-06      
301 iterations  logLik: -396.17  diff: 2.2e-06      
302 iterations  logLik: -396.17  diff: 2.2e-06      
303 iterations  logLik: -396.17  diff: 2.2e-06      
304 iterations  logLik: -396.17  diff: 2.1e-06      
305 iterations  logLik: -396.17  diff: 2.1e-06      
306 iterations  logLik: -396.17  diff: 2.1e-06      
307 iterations  logLik: -396.17  diff: 2.1e-06      
308 iterations  logLik: -396.17  diff: 2.0e-06      
309 iterations  logLik: -396.17  diff: 2.0e-06      
310 iterations  logLik: -396.17  diff: 2.0e-06      
311 iterations  logLik: -396.17  diff: 2.0e-06      
312 iterations  logLik: -396.17  diff: 1.9e-06      
313 iterations  logLik: -396.17  diff: 1.9e-06      
314 iterations  logLik: -396.17  diff: 1.9e-06      
315 iterations  logLik: -396.17  diff: 1.9e-06      
316 iterations  logLik: -396.17  diff: 1.8e-06      
317 iterations  logLik: -396.17  diff: 1.8e-06      
318 iterations  logLik: -396.17  diff: 1.8e-06      
319 iterations  logLik: -396.17  diff: 1.8e-06      
320 iterations  logLik: -396.17  diff: 1.7e-06      
321 iterations  logLik: -396.17  diff: 1.7e-06      
322 iterations  logLik: -396.17  diff: 1.7e-06      
323 iterations  logLik: -396.17  diff: 1.7e-06      
324 iterations  logLik: -396.17  diff: 1.7e-06      
325 iterations  logLik: -396.17  diff: 1.6e-06      
326 iterations  logLik: -396.17  diff: 1.6e-06      
327 iterations  logLik: -396.17  diff: 1.6e-06      
328 iterations  logLik: -396.17  diff: 1.6e-06      
329 iterations  logLik: -396.17  diff: 1.6e-06      
330 iterations  logLik: -396.17  diff: 1.5e-06      
331 iterations  logLik: -396.17  diff: 1.5e-06      
332 iterations  logLik: -396.17  diff: 1.5e-06      
333 iterations  logLik: -396.17  diff: 1.5e-06      
334 iterations  logLik: -396.17  diff: 1.5e-06      
335 iterations  logLik: -396.17  diff: 1.4e-06      
336 iterations  logLik: -396.17  diff: 1.4e-06      
337 iterations  logLik: -396.17  diff: 1.4e-06      
338 iterations  logLik: -396.17  diff: 1.4e-06      
339 iterations  logLik: -396.17  diff: 1.4e-06      
340 iterations  logLik: -396.17  diff: 1.4e-06      
341 iterations  logLik: -396.17  diff: 1.3e-06      
342 iterations  logLik: -396.17  diff: 1.3e-06      
343 iterations  logLik: -396.17  diff: 1.3e-06      
344 iterations  logLik: -396.17  diff: 1.3e-06      
345 iterations  logLik: -396.17  diff: 1.3e-06      
346 iterations  logLik: -396.17  diff: 1.3e-06      
347 iterations  logLik: -396.17  diff: 1.2e-06      
348 iterations  logLik: -396.17  diff: 1.2e-06      
349 iterations  logLik: -396.17  diff: 1.2e-06      
350 iterations  logLik: -396.17  diff: 1.2e-06      
351 iterations  logLik: -396.17  diff: 1.2e-06      
352 iterations  logLik: -396.17  diff: 1.2e-06      
353 iterations  logLik: -396.17  diff: 1.2e-06      
354 iterations  logLik: -396.17  diff: 1.1e-06      
355 iterations  logLik: -396.17  diff: 1.1e-06      
356 iterations  logLik: -396.17  diff: 1.1e-06      
357 iterations  logLik: -396.17  diff: 1.1e-06      
358 iterations  logLik: -396.17  diff: 1.1e-06      
359 iterations  logLik: -396.17  diff: 1.1e-06      
360 iterations  logLik: -396.17  diff: 1.1e-06      
361 iterations  logLik: -396.17  diff: 1.1e-06      
362 iterations  logLik: -396.17  diff: 1.0e-06      
363 iterations  logLik: -396.17  diff: 1.0e-06      
364 iterations  logLik: -396.17  diff: 1.0e-06      
365 iterations  logLik: -396.17  diff: 1.0e-06      
366 iterations  logLik: -396.17  diff: 9.9e-07      
367 iterations  logLik: -396.17  diff: 9.8e-07      
368 iterations  logLik: -396.17  diff: 9.7e-07      
369 iterations  logLik: -396.17  diff: 9.6e-07      
370 iterations  logLik: -396.17  diff: 9.5e-07      
371 iterations  logLik: -396.17  diff: 9.3e-07      
372 iterations  logLik: -396.17  diff: 9.2e-07      
373 iterations  logLik: -396.17  diff: 9.1e-07      
374 iterations  logLik: -396.17  diff: 9.0e-07      
375 iterations  logLik: -396.17  diff: 8.9e-07      
376 iterations  logLik: -396.17  diff: 8.8e-07      
377 iterations  logLik: -396.17  diff: 8.7e-07      
378 iterations  logLik: -396.17  diff: 8.6e-07      
379 iterations  logLik: -396.17  diff: 8.5e-07      
380 iterations  logLik: -396.17  diff: 8.4e-07      
381 iterations  logLik: -396.17  diff: 8.3e-07      
382 iterations  logLik: -396.17  diff: 8.2e-07      
383 iterations  logLik: -396.17  diff: 8.1e-07      
384 iterations  logLik: -396.17  diff: 8.0e-07      
385 iterations  logLik: -396.17  diff: 7.9e-07      
386 iterations  logLik: -396.17  diff: 7.8e-07      
387 iterations  logLik: -396.17  diff: 7.8e-07      
388 iterations  logLik: -396.17  diff: 7.7e-07      
389 iterations  logLik: -396.17  diff: 7.6e-07      
390 iterations  logLik: -396.17  diff: 7.5e-07      
391 iterations  logLik: -396.17  diff: 7.4e-07      
392 iterations  logLik: -396.17  diff: 7.3e-07      
393 iterations  logLik: -396.17  diff: 7.2e-07      
394 iterations  logLik: -396.17  diff: 7.2e-07      
395 iterations  logLik: -396.17  diff: 7.1e-07      
396 iterations  logLik: -396.17  diff: 7.0e-07      
397 iterations  logLik: -396.17  diff: 6.9e-07      
398 iterations  logLik: -396.17  diff: 6.8e-07      
399 iterations  logLik: -396.17  diff: 6.8e-07      
400 iterations  logLik: -396.17  diff: 6.7e-07      
401 iterations  logLik: -396.17  diff: 6.6e-07      
402 iterations  logLik: -396.17  diff: 6.5e-07      
403 iterations  logLik: -396.17  diff: 6.5e-07      
404 iterations  logLik: -396.17  diff: 6.4e-07      
405 iterations  logLik: -396.17  diff: 6.3e-07      
406 iterations  logLik: -396.17  diff: 6.2e-07      
407 iterations  logLik: -396.17  diff: 6.2e-07      
408 iterations  logLik: -396.17  diff: 6.1e-07      
409 iterations  logLik: -396.17  diff: 6.0e-07      
410 iterations  logLik: -396.17  diff: 6.0e-07      
411 iterations  logLik: -396.17  diff: 5.9e-07      
412 iterations  logLik: -396.17  diff: 5.8e-07      
413 iterations  logLik: -396.17  diff: 5.8e-07      
414 iterations  logLik: -396.17  diff: 5.7e-07      
415 iterations  logLik: -396.17  diff: 5.6e-07      
416 iterations  logLik: -396.17  diff: 5.6e-07      
417 iterations  logLik: -396.17  diff: 5.5e-07      
418 iterations  logLik: -396.17  diff: 5.4e-07      
419 iterations  logLik: -396.17  diff: 5.4e-07      
420 iterations  logLik: -396.17  diff: 5.3e-07      
421 iterations  logLik: -396.17  diff: 5.3e-07      
422 iterations  logLik: -396.17  diff: 5.2e-07      
423 iterations  logLik: -396.17  diff: 5.2e-07      
424 iterations  logLik: -396.17  diff: 5.1e-07      
425 iterations  logLik: -396.17  diff: 5.0e-07      
426 iterations  logLik: -396.17  diff: 5.0e-07      
427 iterations  logLik: -396.17  diff: 4.9e-07      
428 iterations  logLik: -396.17  diff: 4.9e-07      
429 iterations  logLik: -396.17  diff: 4.8e-07      
430 iterations  logLik: -396.17  diff: 4.8e-07      
431 iterations  logLik: -396.17  diff: 4.7e-07      
432 iterations  logLik: -396.17  diff: 4.7e-07      
433 iterations  logLik: -396.17  diff: 4.6e-07      
434 iterations  logLik: -396.17  diff: 4.6e-07      
435 iterations  logLik: -396.17  diff: 4.5e-07      
436 iterations  logLik: -396.17  diff: 4.5e-07      
437 iterations  logLik: -396.17  diff: 4.4e-07      
438 iterations  logLik: -396.17  diff: 4.4e-07      
439 iterations  logLik: -396.17  diff: 4.3e-07      
440 iterations  logLik: -396.17  diff: 4.3e-07      
441 iterations  logLik: -396.17  diff: 4.2e-07      
442 iterations  logLik: -396.17  diff: 4.2e-07      
443 iterations  logLik: -396.17  diff: 4.1e-07      
444 iterations  logLik: -396.17  diff: 4.1e-07      
445 iterations  logLik: -396.17  diff: 4.0e-07      
446 iterations  logLik: -396.17  diff: 4.0e-07      
447 iterations  logLik: -396.17  diff: 3.9e-07      
448 iterations  logLik: -396.17  diff: 3.9e-07      
449 iterations  logLik: -396.17  diff: 3.9e-07      
450 iterations  logLik: -396.17  diff: 3.8e-07      
451 iterations  logLik: -396.17  diff: 3.8e-07      
452 iterations  logLik: -396.17  diff: 3.7e-07      
453 iterations  logLik: -396.17  diff: 3.7e-07      
454 iterations  logLik: -396.17  diff: 3.7e-07      
455 iterations  logLik: -396.17  diff: 3.6e-07      
456 iterations  logLik: -396.17  diff: 3.6e-07      
457 iterations  logLik: -396.17  diff: 3.5e-07      
458 iterations  logLik: -396.17  diff: 3.5e-07      
459 iterations  logLik: -396.17  diff: 3.5e-07      
460 iterations  logLik: -396.17  diff: 3.4e-07      
461 iterations  logLik: -396.17  diff: 3.4e-07      
462 iterations  logLik: -396.17  diff: 3.4e-07      
463 iterations  logLik: -396.17  diff: 3.3e-07      
464 iterations  logLik: -396.17  diff: 3.3e-07      
465 iterations  logLik: -396.17  diff: 3.2e-07      
466 iterations  logLik: -396.17  diff: 3.2e-07      
467 iterations  logLik: -396.17  diff: 3.2e-07      
468 iterations  logLik: -396.17  diff: 3.1e-07      
469 iterations  logLik: -396.17  diff: 3.1e-07      
470 iterations  logLik: -396.17  diff: 3.1e-07      
471 iterations  logLik: -396.17  diff: 3.0e-07      
472 iterations  logLik: -396.17  diff: 3.0e-07      
473 iterations  logLik: -396.17  diff: 3.0e-07      
474 iterations  logLik: -396.17  diff: 2.9e-07      
475 iterations  logLik: -396.17  diff: 2.9e-07      
476 iterations  logLik: -396.17  diff: 2.9e-07      
477 iterations  logLik: -396.17  diff: 2.9e-07      
478 iterations  logLik: -396.17  diff: 2.8e-07      
479 iterations  logLik: -396.17  diff: 2.8e-07      
480 iterations  logLik: -396.17  diff: 2.8e-07      
481 iterations  logLik: -396.17  diff: 2.7e-07      
482 iterations  logLik: -396.17  diff: 2.7e-07      
483 iterations  logLik: -396.17  diff: 2.7e-07      
484 iterations  logLik: -396.17  diff: 2.6e-07      
485 iterations  logLik: -396.17  diff: 2.6e-07      
486 iterations  logLik: -396.17  diff: 2.6e-07      
487 iterations  logLik: -396.17  diff: 2.6e-07      
488 iterations  logLik: -396.17  diff: 2.5e-07      
489 iterations  logLik: -396.17  diff: 2.5e-07      
490 iterations  logLik: -396.17  diff: 2.5e-07      
491 iterations  logLik: -396.17  diff: 2.5e-07      
492 iterations  logLik: -396.17  diff: 2.4e-07      
493 iterations  logLik: -396.17  diff: 2.4e-07      
494 iterations  logLik: -396.17  diff: 2.4e-07      
495 iterations  logLik: -396.17  diff: 2.4e-07      
496 iterations  logLik: -396.17  diff: 2.3e-07      
497 iterations  logLik: -396.17  diff: 2.3e-07      
498 iterations  logLik: -396.17  diff: 2.3e-07      
499 iterations  logLik: -396.17  diff: 2.3e-07      
500 iterations  logLik: -396.17  diff: 2.2e-07      
501 iterations  logLik: -396.17  diff: 2.2e-07      
502 iterations  logLik: -396.17  diff: 2.2e-07      
503 iterations  logLik: -396.17  diff: 2.2e-07      
504 iterations  logLik: -396.17  diff: 2.1e-07      
505 iterations  logLik: -396.17  diff: 2.1e-07      
506 iterations  logLik: -396.17  diff: 2.1e-07      
507 iterations  logLik: -396.17  diff: 2.1e-07      
508 iterations  logLik: -396.17  diff: 2.1e-07      
509 iterations  logLik: -396.17  diff: 2.0e-07      
510 iterations  logLik: -396.17  diff: 2.0e-07      
511 iterations  logLik: -396.17  diff: 2.0e-07      
512 iterations  logLik: -396.17  diff: 2.0e-07      
513 iterations  logLik: -396.17  diff: 1.9e-07      
514 iterations  logLik: -396.17  diff: 1.9e-07      
515 iterations  logLik: -396.17  diff: 1.9e-07      
516 iterations  logLik: -396.17  diff: 1.9e-07      
517 iterations  logLik: -396.17  diff: 1.9e-07      
518 iterations  logLik: -396.17  diff: 1.9e-07      
519 iterations  logLik: -396.17  diff: 1.8e-07      
520 iterations  logLik: -396.17  diff: 1.8e-07      
521 iterations  logLik: -396.17  diff: 1.8e-07      
522 iterations  logLik: -396.17  diff: 1.8e-07      
523 iterations  logLik: -396.17  diff: 1.8e-07      
524 iterations  logLik: -396.17  diff: 1.7e-07      
525 iterations  logLik: -396.17  diff: 1.7e-07      
526 iterations  logLik: -396.17  diff: 1.7e-07      
527 iterations  logLik: -396.17  diff: 1.7e-07      
528 iterations  logLik: -396.17  diff: 1.7e-07      
529 iterations  logLik: -396.17  diff: 1.7e-07      
530 iterations  logLik: -396.17  diff: 1.6e-07      
531 iterations  logLik: -396.17  diff: 1.6e-07      
532 iterations  logLik: -396.17  diff: 1.6e-07      
533 iterations  logLik: -396.17  diff: 1.6e-07      
534 iterations  logLik: -396.17  diff: 1.6e-07      
535 iterations  logLik: -396.17  diff: 1.6e-07      
536 iterations  logLik: -396.17  diff: 1.5e-07      
537 iterations  logLik: -396.17  diff: 1.5e-07      
538 iterations  logLik: -396.17  diff: 1.5e-07      
539 iterations  logLik: -396.17  diff: 1.5e-07      
540 iterations  logLik: -396.17  diff: 1.5e-07      
541 iterations  logLik: -396.17  diff: 1.5e-07      
542 iterations  logLik: -396.17  diff: 1.4e-07      
543 iterations  logLik: -396.17  diff: 1.4e-07      
544 iterations  logLik: -396.17  diff: 1.4e-07      
545 iterations  logLik: -396.17  diff: 1.4e-07      
546 iterations  logLik: -396.17  diff: 1.4e-07      
547 iterations  logLik: -396.17  diff: 1.4e-07      
548 iterations  logLik: -396.17  diff: 1.4e-07      
549 iterations  logLik: -396.17  diff: 1.3e-07      
550 iterations  logLik: -396.17  diff: 1.3e-07      
551 iterations  logLik: -396.17  diff: 1.3e-07      
552 iterations  logLik: -396.17  diff: 1.3e-07      
553 iterations  logLik: -396.17  diff: 1.3e-07      
554 iterations  logLik: -396.17  diff: 1.3e-07      
555 iterations  logLik: -396.17  diff: 1.3e-07      
556 iterations  logLik: -396.17  diff: 1.2e-07      
557 iterations  logLik: -396.17  diff: 1.2e-07      
558 iterations  logLik: -396.17  diff: 1.2e-07      
559 iterations  logLik: -396.17  diff: 1.2e-07      
560 iterations  logLik: -396.17  diff: 1.2e-07      
561 iterations  logLik: -396.17  diff: 1.2e-07      
562 iterations  logLik: -396.17  diff: 1.2e-07      
563 iterations  logLik: -396.17  diff: 1.2e-07      
564 iterations  logLik: -396.17  diff: 1.2e-07      
565 iterations  logLik: -396.17  diff: 1.1e-07      
566 iterations  logLik: -396.17  diff: 1.1e-07      
567 iterations  logLik: -396.17  diff: 1.1e-07      
568 iterations  logLik: -396.17  diff: 1.1e-07      
569 iterations  logLik: -396.17  diff: 1.1e-07      
570 iterations  logLik: -396.17  diff: 1.1e-07      
571 iterations  logLik: -396.17  diff: 1.1e-07      
572 iterations  logLik: -396.17  diff: 1.1e-07      
573 iterations  logLik: -396.17  diff: 1.1e-07      
574 iterations  logLik: -396.17  diff: 1.0e-07      
575 iterations  logLik: -396.17  diff: 1.0e-07      
576 iterations  logLik: -396.17  diff: 1.0e-07      
577 iterations  logLik: -396.17  diff: 1.0e-07      
578 iterations  logLik: -396.17  diff: 1.0e-07      
579 iterations  logLik: -396.17  diff: 9.9e-08      
580 iterations  logLik: -396.17  diff: 9.8e-08      
581 iterations  logLik: -396.17  diff: 9.7e-08      
582 iterations  logLik: -396.17  diff: 9.6e-08      
583 iterations  logLik: -396.17  diff: 9.5e-08      
584 iterations  logLik: -396.17  diff: 9.4e-08      
585 iterations  logLik: -396.17  diff: 9.3e-08      
586 iterations  logLik: -396.17  diff: 9.2e-08      
587 iterations  logLik: -396.17  diff: 9.1e-08      
588 iterations  logLik: -396.17  diff: 9.0e-08      
589 iterations  logLik: -396.17  diff: 8.9e-08      
590 iterations  logLik: -396.17  diff: 8.8e-08      
591 iterations  logLik: -396.17  diff: 8.8e-08      
592 iterations  logLik: -396.17  diff: 8.7e-08      
593 iterations  logLik: -396.17  diff: 8.6e-08      
594 iterations  logLik: -396.17  diff: 8.5e-08      
595 iterations  logLik: -396.17  diff: 8.4e-08      
596 iterations  logLik: -396.17  diff: 8.3e-08      
597 iterations  logLik: -396.17  diff: 8.2e-08      
598 iterations  logLik: -396.17  diff: 8.2e-08      
599 iterations  logLik: -396.17  diff: 8.1e-08      
600 iterations  logLik: -396.17  diff: 8.0e-08      
601 iterations  logLik: -396.17  diff: 7.9e-08      
602 iterations  logLik: -396.17  diff: 7.8e-08      
603 iterations  logLik: -396.17  diff: 7.8e-08      
604 iterations  logLik: -396.17  diff: 7.7e-08      
605 iterations  logLik: -396.17  diff: 7.6e-08      
606 iterations  logLik: -396.17  diff: 7.5e-08      
607 iterations  logLik: -396.17  diff: 7.5e-08      
608 iterations  logLik: -396.17  diff: 7.4e-08      
609 iterations  logLik: -396.17  diff: 7.3e-08      
610 iterations  logLik: -396.17  diff: 7.2e-08      
611 iterations  logLik: -396.17  diff: 7.2e-08      
612 iterations  logLik: -396.17  diff: 7.1e-08      
613 iterations  logLik: -396.17  diff: 7.0e-08      
614 iterations  logLik: -396.17  diff: 6.9e-08      
615 iterations  logLik: -396.17  diff: 6.9e-08      
616 iterations  logLik: -396.17  diff: 6.8e-08      
617 iterations  logLik: -396.17  diff: 6.7e-08      
618 iterations  logLik: -396.17  diff: 6.7e-08      
619 iterations  logLik: -396.17  diff: 6.6e-08      
620 iterations  logLik: -396.17  diff: 6.5e-08      
621 iterations  logLik: -396.17  diff: 6.5e-08      
622 iterations  logLik: -396.17  diff: 6.4e-08      
623 iterations  logLik: -396.17  diff: 6.4e-08      
624 iterations  logLik: -396.17  diff: 6.3e-08      
625 iterations  logLik: -396.17  diff: 6.2e-08      
626 iterations  logLik: -396.17  diff: 6.2e-08      
627 iterations  logLik: -396.17  diff: 6.1e-08      
628 iterations  logLik: -396.17  diff: 6.0e-08      
629 iterations  logLik: -396.17  diff: 6.0e-08      
630 iterations  logLik: -396.17  diff: 5.9e-08      
631 iterations  logLik: -396.17  diff: 5.9e-08      
632 iterations  logLik: -396.17  diff: 5.8e-08      
633 iterations  logLik: -396.17  diff: 5.8e-08      
634 iterations  logLik: -396.17  diff: 5.7e-08      
635 iterations  logLik: -396.17  diff: 5.6e-08      
636 iterations  logLik: -396.17  diff: 5.6e-08      
637 iterations  logLik: -396.17  diff: 5.5e-08      
638 iterations  logLik: -396.17  diff: 5.5e-08      
639 iterations  logLik: -396.17  diff: 5.4e-08      
640 iterations  logLik: -396.17  diff: 5.4e-08      
641 iterations  logLik: -396.17  diff: 5.3e-08      
642 iterations  logLik: -396.17  diff: 5.3e-08      
643 iterations  logLik: -396.17  diff: 5.2e-08      
644 iterations  logLik: -396.17  diff: 5.2e-08      
645 iterations  logLik: -396.17  diff: 5.1e-08      
646 iterations  logLik: -396.17  diff: 5.1e-08      
647 iterations  logLik: -396.17  diff: 5.0e-08      
648 iterations  logLik: -396.17  diff: 5.0e-08      
649 iterations  logLik: -396.17  diff: 4.9e-08      
650 iterations  logLik: -396.17  diff: 4.9e-08      
651 iterations  logLik: -396.17  diff: 4.8e-08      
652 iterations  logLik: -396.17  diff: 4.8e-08      
653 iterations  logLik: -396.17  diff: 4.7e-08      
654 iterations  logLik: -396.17  diff: 4.7e-08      
655 iterations  logLik: -396.17  diff: 4.6e-08      
656 iterations  logLik: -396.17  diff: 4.6e-08      
657 iterations  logLik: -396.17  diff: 4.5e-08      
658 iterations  logLik: -396.17  diff: 4.5e-08      
659 iterations  logLik: -396.17  diff: 4.4e-08      
660 iterations  logLik: -396.17  diff: 4.4e-08      
661 iterations  logLik: -396.17  diff: 4.4e-08      
662 iterations  logLik: -396.17  diff: 4.3e-08      
663 iterations  logLik: -396.17  diff: 4.3e-08      
664 iterations  logLik: -396.17  diff: 4.2e-08      
665 iterations  logLik: -396.17  diff: 4.2e-08      
666 iterations  logLik: -396.17  diff: 4.2e-08      
667 iterations  logLik: -396.17  diff: 4.1e-08      
668 iterations  logLik: -396.17  diff: 4.1e-08      
669 iterations  logLik: -396.17  diff: 4.0e-08      
670 iterations  logLik: -396.17  diff: 4.0e-08      
671 iterations  logLik: -396.17  diff: 4.0e-08      
672 iterations  logLik: -396.17  diff: 3.9e-08      
673 iterations  logLik: -396.17  diff: 3.9e-08      
674 iterations  logLik: -396.17  diff: 3.8e-08      
675 iterations  logLik: -396.17  diff: 3.8e-08      
676 iterations  logLik: -396.17  diff: 3.8e-08      
677 iterations  logLik: -396.17  diff: 3.7e-08      
678 iterations  logLik: -396.17  diff: 3.7e-08      
679 iterations  logLik: -396.17  diff: 3.7e-08      
680 iterations  logLik: -396.17  diff: 3.6e-08      
681 iterations  logLik: -396.17  diff: 3.6e-08      
682 iterations  logLik: -396.17  diff: 3.5e-08      
683 iterations  logLik: -396.17  diff: 3.5e-08      
684 iterations  logLik: -396.17  diff: 3.5e-08      
685 iterations  logLik: -396.17  diff: 3.4e-08      
686 iterations  logLik: -396.17  diff: 3.4e-08      
687 iterations  logLik: -396.17  diff: 3.4e-08      
688 iterations  logLik: -396.17  diff: 3.3e-08      
689 iterations  logLik: -396.17  diff: 3.3e-08      
690 iterations  logLik: -396.17  diff: 3.3e-08      
691 iterations  logLik: -396.17  diff: 3.2e-08      
692 iterations  logLik: -396.17  diff: 3.2e-08      
693 iterations  logLik: -396.17  diff: 3.2e-08      
694 iterations  logLik: -396.17  diff: 3.2e-08      
695 iterations  logLik: -396.17  diff: 3.1e-08      
696 iterations  logLik: -396.17  diff: 3.1e-08      
697 iterations  logLik: -396.17  diff: 3.1e-08      
698 iterations  logLik: -396.17  diff: 3.0e-08      
699 iterations  logLik: -396.17  diff: 3.0e-08      
700 iterations  logLik: -396.17  diff: 3.0e-08      
701 iterations  logLik: -396.17  diff: 2.9e-08      
702 iterations  logLik: -396.17  diff: 2.9e-08      
703 iterations  logLik: -396.17  diff: 2.9e-08      
704 iterations  logLik: -396.17  diff: 2.9e-08      
705 iterations  logLik: -396.17  diff: 2.8e-08      
706 iterations  logLik: -396.17  diff: 2.8e-08      
707 iterations  logLik: -396.17  diff: 2.8e-08      
708 iterations  logLik: -396.17  diff: 2.8e-08      
709 iterations  logLik: -396.17  diff: 2.7e-08      
710 iterations  logLik: -396.17  diff: 2.7e-08      
711 iterations  logLik: -396.17  diff: 2.7e-08      
712 iterations  logLik: -396.17  diff: 2.6e-08      
713 iterations  logLik: -396.17  diff: 2.6e-08      
714 iterations  logLik: -396.17  diff: 2.6e-08      
715 iterations  logLik: -396.17  diff: 2.6e-08      
716 iterations  logLik: -396.17  diff: 2.5e-08      
717 iterations  logLik: -396.17  diff: 2.5e-08      
718 iterations  logLik: -396.17  diff: 2.5e-08      
719 iterations  logLik: -396.17  diff: 2.5e-08      
720 iterations  logLik: -396.17  diff: 2.4e-08      
721 iterations  logLik: -396.17  diff: 2.4e-08      
722 iterations  logLik: -396.17  diff: 2.4e-08      
723 iterations  logLik: -396.17  diff: 2.4e-08      
724 iterations  logLik: -396.17  diff: 2.4e-08      
725 iterations  logLik: -396.17  diff: 2.3e-08      
726 iterations  logLik: -396.17  diff: 2.3e-08      
727 iterations  logLik: -396.17  diff: 2.3e-08      
728 iterations  logLik: -396.17  diff: 2.3e-08      
729 iterations  logLik: -396.17  diff: 2.2e-08      
730 iterations  logLik: -396.17  diff: 2.2e-08      
731 iterations  logLik: -396.17  diff: 2.2e-08      
732 iterations  logLik: -396.17  diff: 2.2e-08      
733 iterations  logLik: -396.17  diff: 2.2e-08      
734 iterations  logLik: -396.17  diff: 2.1e-08      
735 iterations  logLik: -396.17  diff: 2.1e-08      
736 iterations  logLik: -396.17  diff: 2.1e-08      
737 iterations  logLik: -396.17  diff: 2.1e-08      
738 iterations  logLik: -396.17  diff: 2.1e-08      
739 iterations  logLik: -396.17  diff: 2.0e-08      
740 iterations  logLik: -396.17  diff: 2.0e-08      
741 iterations  logLik: -396.17  diff: 2.0e-08      
742 iterations  logLik: -396.17  diff: 2.0e-08      
743 iterations  logLik: -396.17  diff: 2.0e-08      
744 iterations  logLik: -396.17  diff: 1.9e-08      
745 iterations  logLik: -396.17  diff: 1.9e-08      
746 iterations  logLik: -396.17  diff: 1.9e-08      
747 iterations  logLik: -396.17  diff: 1.9e-08      
748 iterations  logLik: -396.17  diff: 1.9e-08      
749 iterations  logLik: -396.17  diff: 1.8e-08      
750 iterations  logLik: -396.17  diff: 1.8e-08      
751 iterations  logLik: -396.17  diff: 1.8e-08      
752 iterations  logLik: -396.17  diff: 1.8e-08      
753 iterations  logLik: -396.17  diff: 1.8e-08      
754 iterations  logLik: -396.17  diff: 1.8e-08      
755 iterations  logLik: -396.17  diff: 1.7e-08      
756 iterations  logLik: -396.17  diff: 1.7e-08      
757 iterations  logLik: -396.17  diff: 1.7e-08      
758 iterations  logLik: -396.17  diff: 1.7e-08      
759 iterations  logLik: -396.17  diff: 1.7e-08      
760 iterations  logLik: -396.17  diff: 1.7e-08      
761 iterations  logLik: -396.17  diff: 1.6e-08      
762 iterations  logLik: -396.17  diff: 1.6e-08      
763 iterations  logLik: -396.17  diff: 1.6e-08      
764 iterations  logLik: -396.17  diff: 1.6e-08      
765 iterations  logLik: -396.17  diff: 1.6e-08      
766 iterations  logLik: -396.17  diff: 1.6e-08      
767 iterations  logLik: -396.17  diff: 1.6e-08      
768 iterations  logLik: -396.17  diff: 1.5e-08      
769 iterations  logLik: -396.17  diff: 1.5e-08      
770 iterations  logLik: -396.17  diff: 1.5e-08      
771 iterations  logLik: -396.17  diff: 1.5e-08      
772 iterations  logLik: -396.17  diff: 1.5e-08      
773 iterations  logLik: -396.17  diff: 1.5e-08      
774 iterations  logLik: -396.17  diff: 1.5e-08      
775 iterations  logLik: -396.17  diff: 1.4e-08      
776 iterations  logLik: -396.17  diff: 1.4e-08      
777 iterations  logLik: -396.17  diff: 1.4e-08      
778 iterations  logLik: -396.17  diff: 1.4e-08      
779 iterations  logLik: -396.17  diff: 1.4e-08      
780 iterations  logLik: -396.17  diff: 1.4e-08      
781 iterations  logLik: -396.17  diff: 1.4e-08      
782 iterations  logLik: -396.17  diff: 1.3e-08      
783 iterations  logLik: -396.17  diff: 1.3e-08      
784 iterations  logLik: -396.17  diff: 1.3e-08      
785 iterations  logLik: -396.17  diff: 1.3e-08      
786 iterations  logLik: -396.17  diff: 1.3e-08      
787 iterations  logLik: -396.17  diff: 1.3e-08      
788 iterations  logLik: -396.17  diff: 1.3e-08      
789 iterations  logLik: -396.17  diff: 1.3e-08      
790 iterations  logLik: -396.17  diff: 1.2e-08      
791 iterations  logLik: -396.17  diff: 1.2e-08      
792 iterations  logLik: -396.17  diff: 1.2e-08      
793 iterations  logLik: -396.17  diff: 1.2e-08      
794 iterations  logLik: -396.17  diff: 1.2e-08      
795 iterations  logLik: -396.17  diff: 1.2e-08      
796 iterations  logLik: -396.17  diff: 1.2e-08      
797 iterations  logLik: -396.17  diff: 1.2e-08      
798 iterations  logLik: -396.17  diff: 1.2e-08      
799 iterations  logLik: -396.17  diff: 1.1e-08      
800 iterations  logLik: -396.17  diff: 1.1e-08      
801 iterations  logLik: -396.17  diff: 1.1e-08      
802 iterations  logLik: -396.17  diff: 1.1e-08      
803 iterations  logLik: -396.17  diff: 1.1e-08      
804 iterations  logLik: -396.17  diff: 1.1e-08      
805 iterations  logLik: -396.17  diff: 1.1e-08      
806 iterations  logLik: -396.17  diff: 1.1e-08      
807 iterations  logLik: -396.17  diff: 1.1e-08      
808 iterations  logLik: -396.17  diff: 1.0e-08      
809 iterations  logLik: -396.17  diff: 1.0e-08      
810 iterations  logLik: -396.17  diff: 1.0e-08      
811 iterations  logLik: -396.17  diff: 1.0e-08      
812 iterations  logLik: -396.17  diff: 1.0e-08      
813 iterations  logLik: -396.17  diff: 1.0e-08      

#> .. done.
#> Structural Latent Variable Model
#> 
#> Latent variables (Root*):           
#>  Label: lc*
#> nclass: 3  
#> 
#> Measurement model:                            
#>  lc -> { y1, y2, y3, y4 }  a
#> 
#> LogLik: -396.1722
estimate(m, dt$response, method = "nlm",
         control = slcaControl(verbose = TRUE))
#> nlm iteration begin.
#> .. done.
#> Structural Latent Variable Model
#> 
#> Latent variables (Root*):           
#>  Label: lc*
#> nclass: 3  
#> 
#> Measurement model:                            
#>  lc -> { y1, y2, y3, y4 }  a
#> 
#> LogLik: -396.1722
estimate(m, dt$response, method = "hybrid",
         control = slcaControl(verbose = TRUE))
#> EM iteration begin.
#> 1 iterations  logLik: -555.41  diff: inf      
2 iterations  logLik: -512.40  diff: 43.00670      
3 iterations  logLik: -512.18  diff: 0.22563      
4 iterations  logLik: -510.99  diff: 1.18797      
5 iterations  logLik: -505.07  diff: 5.92316      
6 iterations  logLik: -482.78  diff: 22.28088      
7 iterations  logLik: -446.04  diff: 36.74039      
8 iterations  logLik: -428.90  diff: 17.14771      
9 iterations  logLik: -419.91  diff: 8.98238      
10 iterations  logLik: -408.31  diff: 11.60348      
11 iterations  logLik: -400.08  diff: 8.23315      
12 iterations  logLik: -397.76  diff: 2.32009      
13 iterations  logLik: -397.45  diff: 0.30874      
14 iterations  logLik: -397.40  diff: 0.04823      
15 iterations  logLik: -397.37  diff: 0.02836      
16 iterations  logLik: -397.34  diff: 0.02808      
17 iterations  logLik: -397.32  diff: 0.02913      
18 iterations  logLik: -397.28  diff: 0.03025      
19 iterations  logLik: -397.25  diff: 0.03136      
20 iterations  logLik: -397.22  diff: 0.03243      
21 iterations  logLik: -397.19  diff: 0.03346      
22 iterations  logLik: -397.15  diff: 0.03442      
23 iterations  logLik: -397.12  diff: 0.03530      
24 iterations  logLik: -397.08  diff: 0.03610      
25 iterations  logLik: -397.05  diff: 0.03678      
26 iterations  logLik: -397.01  diff: 0.03733      
27 iterations  logLik: -396.97  diff: 0.03773      
28 iterations  logLik: -396.93  diff: 0.03794      
29 iterations  logLik: -396.89  diff: 0.03796      
30 iterations  logLik: -396.86  diff: 0.03776      
31 iterations  logLik: -396.82  diff: 0.03732      
32 iterations  logLik: -396.78  diff: 0.03665      
33 iterations  logLik: -396.75  diff: 0.03573      
34 iterations  logLik: -396.71  diff: 0.03460      
35 iterations  logLik: -396.68  diff: 0.03326      
36 iterations  logLik: -396.65  diff: 0.03176      
37 iterations  logLik: -396.62  diff: 0.03013      
38 iterations  logLik: -396.59  diff: 0.02841      
39 iterations  logLik: -396.56  diff: 0.02664      
40 iterations  logLik: -396.54  diff: 0.02487      
41 iterations  logLik: -396.51  diff: 0.02311      
42 iterations  logLik: -396.49  diff: 0.02141      
43 iterations  logLik: -396.47  diff: 0.01977      
44 iterations  logLik: -396.45  diff: 0.01822      
45 iterations  logLik: -396.44  diff: 0.01677      
46 iterations  logLik: -396.42  diff: 0.01542      
47 iterations  logLik: -396.41  diff: 0.01417      
48 iterations  logLik: -396.40  diff: 0.01302      
49 iterations  logLik: -396.38  diff: 0.01198      
50 iterations  logLik: -396.37  diff: 0.01102      
51 iterations  logLik: -396.36  diff: 0.01016      
52 iterations  logLik: -396.35  diff: 0.00937      
53 iterations  logLik: -396.34  diff: 0.00866      
54 iterations  logLik: -396.34  diff: 0.00803      
55 iterations  logLik: -396.33  diff: 0.00745      
56 iterations  logLik: -396.32  diff: 0.00693      
57 iterations  logLik: -396.31  diff: 0.00646      
58 iterations  logLik: -396.31  diff: 0.00604      
59 iterations  logLik: -396.30  diff: 0.00565      
60 iterations  logLik: -396.30  diff: 0.00531      
61 iterations  logLik: -396.29  diff: 0.00499      
62 iterations  logLik: -396.29  diff: 0.00471      
63 iterations  logLik: -396.28  diff: 0.00445      
64 iterations  logLik: -396.28  diff: 0.00421      
65 iterations  logLik: -396.28  diff: 0.00400      
66 iterations  logLik: -396.27  diff: 0.00380      
67 iterations  logLik: -396.27  diff: 0.00362      
68 iterations  logLik: -396.26  diff: 0.00345      
69 iterations  logLik: -396.26  diff: 0.00330      
70 iterations  logLik: -396.26  diff: 0.00315      
71 iterations  logLik: -396.26  diff: 0.00302      
72 iterations  logLik: -396.25  diff: 0.00290      
73 iterations  logLik: -396.25  diff: 0.00278      
74 iterations  logLik: -396.25  diff: 0.00267      
75 iterations  logLik: -396.24  diff: 0.00257      
76 iterations  logLik: -396.24  diff: 0.00247      
77 iterations  logLik: -396.24  diff: 0.00238      
78 iterations  logLik: -396.24  diff: 0.00230      
79 iterations  logLik: -396.24  diff: 0.00222      
80 iterations  logLik: -396.23  diff: 0.00214      
81 iterations  logLik: -396.23  diff: 0.00206      
82 iterations  logLik: -396.23  diff: 0.00199      
83 iterations  logLik: -396.23  diff: 0.00192      
84 iterations  logLik: -396.23  diff: 0.00186      
85 iterations  logLik: -396.22  diff: 0.00180      
86 iterations  logLik: -396.22  diff: 0.00174      
87 iterations  logLik: -396.22  diff: 0.00168      
88 iterations  logLik: -396.22  diff: 0.00162      
89 iterations  logLik: -396.22  diff: 0.00157      
90 iterations  logLik: -396.22  diff: 0.00152      
91 iterations  logLik: -396.21  diff: 0.00147      
92 iterations  logLik: -396.21  diff: 0.00142      
93 iterations  logLik: -396.21  diff: 0.00137      
94 iterations  logLik: -396.21  diff: 0.00132      
95 iterations  logLik: -396.21  diff: 0.00128      
96 iterations  logLik: -396.21  diff: 0.00124      
97 iterations  logLik: -396.21  diff: 0.00119      
98 iterations  logLik: -396.20  diff: 0.00115      
99 iterations  logLik: -396.20  diff: 0.00111      
100 iterations  logLik: -396.20  diff: 0.00108      
101 iterations  logLik: -396.20  diff: 0.00104      
102 iterations  logLik: -396.20  diff: 0.00100      
103 iterations  logLik: -396.20  diff: 0.00097      
104 iterations  logLik: -396.20  diff: 0.00094      
105 iterations  logLik: -396.20  diff: 0.00090      
106 iterations  logLik: -396.20  diff: 0.00087      
107 iterations  logLik: -396.20  diff: 0.00084      
108 iterations  logLik: -396.20  diff: 0.00081      
109 iterations  logLik: -396.19  diff: 0.00078      
110 iterations  logLik: -396.19  diff: 0.00075      
111 iterations  logLik: -396.19  diff: 0.00073      
112 iterations  logLik: -396.19  diff: 0.00070      
113 iterations  logLik: -396.19  diff: 0.00068      
114 iterations  logLik: -396.19  diff: 0.00065      
115 iterations  logLik: -396.19  diff: 0.00063      
116 iterations  logLik: -396.19  diff: 0.00061      
117 iterations  logLik: -396.19  diff: 0.00058      
118 iterations  logLik: -396.19  diff: 0.00056      
119 iterations  logLik: -396.19  diff: 0.00054      
120 iterations  logLik: -396.19  diff: 0.00052      
121 iterations  logLik: -396.19  diff: 0.00050      
122 iterations  logLik: -396.19  diff: 0.00049      
123 iterations  logLik: -396.19  diff: 0.00047      
124 iterations  logLik: -396.19  diff: 0.00045      
125 iterations  logLik: -396.19  diff: 0.00043      
126 iterations  logLik: -396.18  diff: 0.00042      
127 iterations  logLik: -396.18  diff: 0.00040      
128 iterations  logLik: -396.18  diff: 0.00039      
129 iterations  logLik: -396.18  diff: 0.00037      
130 iterations  logLik: -396.18  diff: 0.00036      
131 iterations  logLik: -396.18  diff: 0.00035      
132 iterations  logLik: -396.18  diff: 0.00033      
133 iterations  logLik: -396.18  diff: 0.00032      
134 iterations  logLik: -396.18  diff: 0.00031      
135 iterations  logLik: -396.18  diff: 0.00030      
136 iterations  logLik: -396.18  diff: 0.00029      
137 iterations  logLik: -396.18  diff: 0.00028      
138 iterations  logLik: -396.18  diff: 0.00027      
139 iterations  logLik: -396.18  diff: 0.00026      
140 iterations  logLik: -396.18  diff: 0.00025      
141 iterations  logLik: -396.18  diff: 0.00024      
142 iterations  logLik: -396.18  diff: 0.00023      
143 iterations  logLik: -396.18  diff: 0.00022      
144 iterations  logLik: -396.18  diff: 0.00021      
145 iterations  logLik: -396.18  diff: 0.00020      
146 iterations  logLik: -396.18  diff: 0.00020      
147 iterations  logLik: -396.18  diff: 0.00019      
148 iterations  logLik: -396.18  diff: 0.00018      
149 iterations  logLik: -396.18  diff: 0.00018      
150 iterations  logLik: -396.18  diff: 0.00017      
151 iterations  logLik: -396.18  diff: 0.00016      
152 iterations  logLik: -396.18  diff: 0.00016      
153 iterations  logLik: -396.18  diff: 0.00015      
154 iterations  logLik: -396.18  diff: 0.00015      
155 iterations  logLik: -396.18  diff: 0.00014      
156 iterations  logLik: -396.18  diff: 0.00014      
157 iterations  logLik: -396.18  diff: 0.00013      
158 iterations  logLik: -396.18  diff: 0.00013      
159 iterations  logLik: -396.18  diff: 0.00012      
160 iterations  logLik: -396.18  diff: 0.00012      
161 iterations  logLik: -396.18  diff: 0.00012      
162 iterations  logLik: -396.18  diff: 0.00011      
163 iterations  logLik: -396.18  diff: 0.00011      
164 iterations  logLik: -396.18  diff: 0.00010      
165 iterations  logLik: -396.18  diff: 0.00010      
166 iterations  logLik: -396.18  diff: 0.00010      
167 iterations  logLik: -396.18  diff: 0.00009      
168 iterations  logLik: -396.18  diff: 0.00009      
169 iterations  logLik: -396.18  diff: 0.00009      
170 iterations  logLik: -396.18  diff: 0.00009      
171 iterations  logLik: -396.18  diff: 0.00008      
172 iterations  logLik: -396.18  diff: 0.00008      
173 iterations  logLik: -396.18  diff: 0.00008      
174 iterations  logLik: -396.18  diff: 0.00007      
175 iterations  logLik: -396.18  diff: 0.00007      
176 iterations  logLik: -396.18  diff: 0.00007      
177 iterations  logLik: -396.18  diff: 0.00007      
178 iterations  logLik: -396.18  diff: 0.00007      
179 iterations  logLik: -396.18  diff: 0.00006      
180 iterations  logLik: -396.18  diff: 0.00006      
181 iterations  logLik: -396.17  diff: 0.00006      
182 iterations  logLik: -396.17  diff: 0.00006      
183 iterations  logLik: -396.17  diff: 0.00006      
184 iterations  logLik: -396.17  diff: 0.00006      
185 iterations  logLik: -396.17  diff: 0.00005      
186 iterations  logLik: -396.17  diff: 0.00005      
187 iterations  logLik: -396.17  diff: 0.00005      
188 iterations  logLik: -396.17  diff: 0.00005      
189 iterations  logLik: -396.17  diff: 0.00005      
190 iterations  logLik: -396.17  diff: 0.00005      
191 iterations  logLik: -396.17  diff: 0.00005      
192 iterations  logLik: -396.17  diff: 0.00004      
193 iterations  logLik: -396.17  diff: 0.00004      
194 iterations  logLik: -396.17  diff: 0.00004      
195 iterations  logLik: -396.17  diff: 0.00004      
196 iterations  logLik: -396.17  diff: 0.00004      
197 iterations  logLik: -396.17  diff: 0.00004      
198 iterations  logLik: -396.17  diff: 0.00004      
199 iterations  logLik: -396.17  diff: 0.00004      
200 iterations  logLik: -396.17  diff: 0.00004      
201 iterations  logLik: -396.17  diff: 0.00003      
202 iterations  logLik: -396.17  diff: 0.00003      
203 iterations  logLik: -396.17  diff: 0.00003      
204 iterations  logLik: -396.17  diff: 0.00003      
205 iterations  logLik: -396.17  diff: 0.00003      
206 iterations  logLik: -396.17  diff: 0.00003      
207 iterations  logLik: -396.17  diff: 0.00003      
208 iterations  logLik: -396.17  diff: 0.00003      
209 iterations  logLik: -396.17  diff: 0.00003      
210 iterations  logLik: -396.17  diff: 0.00003      
211 iterations  logLik: -396.17  diff: 0.00003      
212 iterations  logLik: -396.17  diff: 0.00003      
213 iterations  logLik: -396.17  diff: 0.00003      
214 iterations  logLik: -396.17  diff: 0.00003      
215 iterations  logLik: -396.17  diff: 0.00002      
216 iterations  logLik: -396.17  diff: 0.00002      
217 iterations  logLik: -396.17  diff: 0.00002      
218 iterations  logLik: -396.17  diff: 0.00002      
219 iterations  logLik: -396.17  diff: 0.00002      
220 iterations  logLik: -396.17  diff: 0.00002      
221 iterations  logLik: -396.17  diff: 0.00002      
222 iterations  logLik: -396.17  diff: 0.00002      
223 iterations  logLik: -396.17  diff: 0.00002      
224 iterations  logLik: -396.17  diff: 0.00002      
225 iterations  logLik: -396.17  diff: 0.00002      
226 iterations  logLik: -396.17  diff: 0.00002      
227 iterations  logLik: -396.17  diff: 0.00002      
228 iterations  logLik: -396.17  diff: 0.00002      
229 iterations  logLik: -396.17  diff: 0.00002      
230 iterations  logLik: -396.17  diff: 0.00002      
231 iterations  logLik: -396.17  diff: 0.00002      
232 iterations  logLik: -396.17  diff: 0.00002      
233 iterations  logLik: -396.17  diff: 0.00002      
234 iterations  logLik: -396.17  diff: 0.00002      
235 iterations  logLik: -396.17  diff: 0.00002      
236 iterations  logLik: -396.17  diff: 0.00002      
237 iterations  logLik: -396.17  diff: 0.00002      
238 iterations  logLik: -396.17  diff: 0.00002      
239 iterations  logLik: -396.17  diff: 0.00002      
240 iterations  logLik: -396.17  diff: 0.00001      
241 iterations  logLik: -396.17  diff: 0.00001      
242 iterations  logLik: -396.17  diff: 0.00001      
243 iterations  logLik: -396.17  diff: 0.00001      
244 iterations  logLik: -396.17  diff: 0.00001      
245 iterations  logLik: -396.17  diff: 0.00001      
246 iterations  logLik: -396.17  diff: 0.00001      
247 iterations  logLik: -396.17  diff: 0.00001      
248 iterations  logLik: -396.17  diff: 0.00001      
249 iterations  logLik: -396.17  diff: 0.00001      
250 iterations  logLik: -396.17  diff: 0.00001      
251 iterations  logLik: -396.17  diff: 0.00001      
252 iterations  logLik: -396.17  diff: 0.00001      
253 iterations  logLik: -396.17  diff: 0.00001      
254 iterations  logLik: -396.17  diff: 0.00001      
255 iterations  logLik: -396.17  diff: 0.00001      
256 iterations  logLik: -396.17  diff: 0.00001      
257 iterations  logLik: -396.17  diff: 0.00001      
258 iterations  logLik: -396.17  diff: 0.00001      
259 iterations  logLik: -396.17  diff: 0.00001      
260 iterations  logLik: -396.17  diff: 0.00001      
261 iterations  logLik: -396.17  diff: 0.00001      
262 iterations  logLik: -396.17  diff: 0.00001      
263 iterations  logLik: -396.17  diff: 9.9e-06      
264 iterations  logLik: -396.17  diff: 9.7e-06      
265 iterations  logLik: -396.17  diff: 9.5e-06      
266 iterations  logLik: -396.17  diff: 9.4e-06      
267 iterations  logLik: -396.17  diff: 9.2e-06      
268 iterations  logLik: -396.17  diff: 9.1e-06      
269 iterations  logLik: -396.17  diff: 8.9e-06      
270 iterations  logLik: -396.17  diff: 8.8e-06      
271 iterations  logLik: -396.17  diff: 8.6e-06      
272 iterations  logLik: -396.17  diff: 8.5e-06      
273 iterations  logLik: -396.17  diff: 8.4e-06      
274 iterations  logLik: -396.17  diff: 8.2e-06      
275 iterations  logLik: -396.17  diff: 8.1e-06      
276 iterations  logLik: -396.17  diff: 8.0e-06      
277 iterations  logLik: -396.17  diff: 7.8e-06      
278 iterations  logLik: -396.17  diff: 7.7e-06      
279 iterations  logLik: -396.17  diff: 7.6e-06      
280 iterations  logLik: -396.17  diff: 7.5e-06      
281 iterations  logLik: -396.17  diff: 7.3e-06      
282 iterations  logLik: -396.17  diff: 7.2e-06      
283 iterations  logLik: -396.17  diff: 7.1e-06      
284 iterations  logLik: -396.17  diff: 7.0e-06      
285 iterations  logLik: -396.17  diff: 6.9e-06      
286 iterations  logLik: -396.17  diff: 6.8e-06      
287 iterations  logLik: -396.17  diff: 6.7e-06      
288 iterations  logLik: -396.17  diff: 6.6e-06      
289 iterations  logLik: -396.17  diff: 6.5e-06      
290 iterations  logLik: -396.17  diff: 6.4e-06      
291 iterations  logLik: -396.17  diff: 6.3e-06      
292 iterations  logLik: -396.17  diff: 6.2e-06      
293 iterations  logLik: -396.17  diff: 6.1e-06      
294 iterations  logLik: -396.17  diff: 6.0e-06      
295 iterations  logLik: -396.17  diff: 5.9e-06      
296 iterations  logLik: -396.17  diff: 5.8e-06      
297 iterations  logLik: -396.17  diff: 5.7e-06      
298 iterations  logLik: -396.17  diff: 5.6e-06      
299 iterations  logLik: -396.17  diff: 5.6e-06      
300 iterations  logLik: -396.17  diff: 5.5e-06      
301 iterations  logLik: -396.17  diff: 5.4e-06      
302 iterations  logLik: -396.17  diff: 5.3e-06      
303 iterations  logLik: -396.17  diff: 5.2e-06      
304 iterations  logLik: -396.17  diff: 5.2e-06      
305 iterations  logLik: -396.17  diff: 5.1e-06      
306 iterations  logLik: -396.17  diff: 5.0e-06      
307 iterations  logLik: -396.17  diff: 4.9e-06      
308 iterations  logLik: -396.17  diff: 4.9e-06      
309 iterations  logLik: -396.17  diff: 4.8e-06      
310 iterations  logLik: -396.17  diff: 4.7e-06      
311 iterations  logLik: -396.17  diff: 4.7e-06      
312 iterations  logLik: -396.17  diff: 4.6e-06      
313 iterations  logLik: -396.17  diff: 4.5e-06      
314 iterations  logLik: -396.17  diff: 4.5e-06      
315 iterations  logLik: -396.17  diff: 4.4e-06      
316 iterations  logLik: -396.17  diff: 4.3e-06      
317 iterations  logLik: -396.17  diff: 4.3e-06      
318 iterations  logLik: -396.17  diff: 4.2e-06      
319 iterations  logLik: -396.17  diff: 4.1e-06      
320 iterations  logLik: -396.17  diff: 4.1e-06      
321 iterations  logLik: -396.17  diff: 4.0e-06      
322 iterations  logLik: -396.17  diff: 4.0e-06      
323 iterations  logLik: -396.17  diff: 3.9e-06      
324 iterations  logLik: -396.17  diff: 3.9e-06      
325 iterations  logLik: -396.17  diff: 3.8e-06      
326 iterations  logLik: -396.17  diff: 3.8e-06      
327 iterations  logLik: -396.17  diff: 3.7e-06      
328 iterations  logLik: -396.17  diff: 3.7e-06      
329 iterations  logLik: -396.17  diff: 3.6e-06      
330 iterations  logLik: -396.17  diff: 3.6e-06      
331 iterations  logLik: -396.17  diff: 3.5e-06      
332 iterations  logLik: -396.17  diff: 3.5e-06      
333 iterations  logLik: -396.17  diff: 3.4e-06      
334 iterations  logLik: -396.17  diff: 3.4e-06      
335 iterations  logLik: -396.17  diff: 3.3e-06      
336 iterations  logLik: -396.17  diff: 3.3e-06      
337 iterations  logLik: -396.17  diff: 3.2e-06      
338 iterations  logLik: -396.17  diff: 3.2e-06      
339 iterations  logLik: -396.17  diff: 3.1e-06      
340 iterations  logLik: -396.17  diff: 3.1e-06      
341 iterations  logLik: -396.17  diff: 3.1e-06      
342 iterations  logLik: -396.17  diff: 3.0e-06      
343 iterations  logLik: -396.17  diff: 3.0e-06      
344 iterations  logLik: -396.17  diff: 2.9e-06      
345 iterations  logLik: -396.17  diff: 2.9e-06      
346 iterations  logLik: -396.17  diff: 2.9e-06      
347 iterations  logLik: -396.17  diff: 2.8e-06      
348 iterations  logLik: -396.17  diff: 2.8e-06      
349 iterations  logLik: -396.17  diff: 2.7e-06      
350 iterations  logLik: -396.17  diff: 2.7e-06      
351 iterations  logLik: -396.17  diff: 2.7e-06      
352 iterations  logLik: -396.17  diff: 2.6e-06      
353 iterations  logLik: -396.17  diff: 2.6e-06      
354 iterations  logLik: -396.17  diff: 2.6e-06      
355 iterations  logLik: -396.17  diff: 2.5e-06      
356 iterations  logLik: -396.17  diff: 2.5e-06      
357 iterations  logLik: -396.17  diff: 2.5e-06      
358 iterations  logLik: -396.17  diff: 2.4e-06      
359 iterations  logLik: -396.17  diff: 2.4e-06      
360 iterations  logLik: -396.17  diff: 2.4e-06      
361 iterations  logLik: -396.17  diff: 2.3e-06      
362 iterations  logLik: -396.17  diff: 2.3e-06      
363 iterations  logLik: -396.17  diff: 2.3e-06      
364 iterations  logLik: -396.17  diff: 2.2e-06      
365 iterations  logLik: -396.17  diff: 2.2e-06      
366 iterations  logLik: -396.17  diff: 2.2e-06      
367 iterations  logLik: -396.17  diff: 2.2e-06      
368 iterations  logLik: -396.17  diff: 2.1e-06      
369 iterations  logLik: -396.17  diff: 2.1e-06      
370 iterations  logLik: -396.17  diff: 2.1e-06      
371 iterations  logLik: -396.17  diff: 2.1e-06      
372 iterations  logLik: -396.17  diff: 2.0e-06      
373 iterations  logLik: -396.17  diff: 2.0e-06      
374 iterations  logLik: -396.17  diff: 2.0e-06      
375 iterations  logLik: -396.17  diff: 1.9e-06      
376 iterations  logLik: -396.17  diff: 1.9e-06      
377 iterations  logLik: -396.17  diff: 1.9e-06      
378 iterations  logLik: -396.17  diff: 1.9e-06      
379 iterations  logLik: -396.17  diff: 1.9e-06      
380 iterations  logLik: -396.17  diff: 1.8e-06      
381 iterations  logLik: -396.17  diff: 1.8e-06      
382 iterations  logLik: -396.17  diff: 1.8e-06      
383 iterations  logLik: -396.17  diff: 1.8e-06      
384 iterations  logLik: -396.17  diff: 1.7e-06      
385 iterations  logLik: -396.17  diff: 1.7e-06      
386 iterations  logLik: -396.17  diff: 1.7e-06      
387 iterations  logLik: -396.17  diff: 1.7e-06      
388 iterations  logLik: -396.17  diff: 1.7e-06      
389 iterations  logLik: -396.17  diff: 1.6e-06      
390 iterations  logLik: -396.17  diff: 1.6e-06      
391 iterations  logLik: -396.17  diff: 1.6e-06      
392 iterations  logLik: -396.17  diff: 1.6e-06      
393 iterations  logLik: -396.17  diff: 1.6e-06      
394 iterations  logLik: -396.17  diff: 1.5e-06      
395 iterations  logLik: -396.17  diff: 1.5e-06      
396 iterations  logLik: -396.17  diff: 1.5e-06      
397 iterations  logLik: -396.17  diff: 1.5e-06      
398 iterations  logLik: -396.17  diff: 1.5e-06      
399 iterations  logLik: -396.17  diff: 1.4e-06      
400 iterations  logLik: -396.17  diff: 1.4e-06      
401 iterations  logLik: -396.17  diff: 1.4e-06      
402 iterations  logLik: -396.17  diff: 1.4e-06      
403 iterations  logLik: -396.17  diff: 1.4e-06      
404 iterations  logLik: -396.17  diff: 1.4e-06      
405 iterations  logLik: -396.17  diff: 1.3e-06      
406 iterations  logLik: -396.17  diff: 1.3e-06      
407 iterations  logLik: -396.17  diff: 1.3e-06      
408 iterations  logLik: -396.17  diff: 1.3e-06      
409 iterations  logLik: -396.17  diff: 1.3e-06      
410 iterations  logLik: -396.17  diff: 1.3e-06      
411 iterations  logLik: -396.17  diff: 1.2e-06      
412 iterations  logLik: -396.17  diff: 1.2e-06      
413 iterations  logLik: -396.17  diff: 1.2e-06      
414 iterations  logLik: -396.17  diff: 1.2e-06      
415 iterations  logLik: -396.17  diff: 1.2e-06      
416 iterations  logLik: -396.17  diff: 1.2e-06      
417 iterations  logLik: -396.17  diff: 1.2e-06      
418 iterations  logLik: -396.17  diff: 1.1e-06      
419 iterations  logLik: -396.17  diff: 1.1e-06      
420 iterations  logLik: -396.17  diff: 1.1e-06      
421 iterations  logLik: -396.17  diff: 1.1e-06      
422 iterations  logLik: -396.17  diff: 1.1e-06      
423 iterations  logLik: -396.17  diff: 1.1e-06      
424 iterations  logLik: -396.17  diff: 1.1e-06      
425 iterations  logLik: -396.17  diff: 1.1e-06      
426 iterations  logLik: -396.17  diff: 1.0e-06      
427 iterations  logLik: -396.17  diff: 1.0e-06      
428 iterations  logLik: -396.17  diff: 1.0e-06      
429 iterations  logLik: -396.17  diff: 1.0e-06      
430 iterations  logLik: -396.17  diff: 9.9e-07      
431 iterations  logLik: -396.17  diff: 9.8e-07      
432 iterations  logLik: -396.17  diff: 9.7e-07      
433 iterations  logLik: -396.17  diff: 9.6e-07      
434 iterations  logLik: -396.17  diff: 9.4e-07      
435 iterations  logLik: -396.17  diff: 9.3e-07      
436 iterations  logLik: -396.17  diff: 9.2e-07      
437 iterations  logLik: -396.17  diff: 9.1e-07      
438 iterations  logLik: -396.17  diff: 9.0e-07      
439 iterations  logLik: -396.17  diff: 8.9e-07      
440 iterations  logLik: -396.17  diff: 8.8e-07      
441 iterations  logLik: -396.17  diff: 8.7e-07      
442 iterations  logLik: -396.17  diff: 8.6e-07      
443 iterations  logLik: -396.17  diff: 8.5e-07      
444 iterations  logLik: -396.17  diff: 8.4e-07      
445 iterations  logLik: -396.17  diff: 8.3e-07      
446 iterations  logLik: -396.17  diff: 8.2e-07      
447 iterations  logLik: -396.17  diff: 8.1e-07      
448 iterations  logLik: -396.17  diff: 8.0e-07      
449 iterations  logLik: -396.17  diff: 7.9e-07      
450 iterations  logLik: -396.17  diff: 7.8e-07      
451 iterations  logLik: -396.17  diff: 7.7e-07      
452 iterations  logLik: -396.17  diff: 7.6e-07      
453 iterations  logLik: -396.17  diff: 7.6e-07      
454 iterations  logLik: -396.17  diff: 7.5e-07      
455 iterations  logLik: -396.17  diff: 7.4e-07      
456 iterations  logLik: -396.17  diff: 7.3e-07      
457 iterations  logLik: -396.17  diff: 7.2e-07      
458 iterations  logLik: -396.17  diff: 7.1e-07      
459 iterations  logLik: -396.17  diff: 7.1e-07      
460 iterations  logLik: -396.17  diff: 7.0e-07      
461 iterations  logLik: -396.17  diff: 6.9e-07      
462 iterations  logLik: -396.17  diff: 6.8e-07      
463 iterations  logLik: -396.17  diff: 6.7e-07      
464 iterations  logLik: -396.17  diff: 6.7e-07      
465 iterations  logLik: -396.17  diff: 6.6e-07      
466 iterations  logLik: -396.17  diff: 6.5e-07      
467 iterations  logLik: -396.17  diff: 6.4e-07      
468 iterations  logLik: -396.17  diff: 6.4e-07      
469 iterations  logLik: -396.17  diff: 6.3e-07      
470 iterations  logLik: -396.17  diff: 6.2e-07      
471 iterations  logLik: -396.17  diff: 6.2e-07      
472 iterations  logLik: -396.17  diff: 6.1e-07      
473 iterations  logLik: -396.17  diff: 6.0e-07      
474 iterations  logLik: -396.17  diff: 5.9e-07      
475 iterations  logLik: -396.17  diff: 5.9e-07      
476 iterations  logLik: -396.17  diff: 5.8e-07      
477 iterations  logLik: -396.17  diff: 5.7e-07      
478 iterations  logLik: -396.17  diff: 5.7e-07      
479 iterations  logLik: -396.17  diff: 5.6e-07      
480 iterations  logLik: -396.17  diff: 5.6e-07      
481 iterations  logLik: -396.17  diff: 5.5e-07      
482 iterations  logLik: -396.17  diff: 5.4e-07      
483 iterations  logLik: -396.17  diff: 5.4e-07      
484 iterations  logLik: -396.17  diff: 5.3e-07      
485 iterations  logLik: -396.17  diff: 5.3e-07      
486 iterations  logLik: -396.17  diff: 5.2e-07      
487 iterations  logLik: -396.17  diff: 5.1e-07      
488 iterations  logLik: -396.17  diff: 5.1e-07      
489 iterations  logLik: -396.17  diff: 5.0e-07      
490 iterations  logLik: -396.17  diff: 5.0e-07      
491 iterations  logLik: -396.17  diff: 4.9e-07      
492 iterations  logLik: -396.17  diff: 4.9e-07      
493 iterations  logLik: -396.17  diff: 4.8e-07      
494 iterations  logLik: -396.17  diff: 4.8e-07      
495 iterations  logLik: -396.17  diff: 4.7e-07      
496 iterations  logLik: -396.17  diff: 4.6e-07      
497 iterations  logLik: -396.17  diff: 4.6e-07      
498 iterations  logLik: -396.17  diff: 4.5e-07      
499 iterations  logLik: -396.17  diff: 4.5e-07      
500 iterations  logLik: -396.17  diff: 4.4e-07      
501 iterations  logLik: -396.17  diff: 4.4e-07      
502 iterations  logLik: -396.17  diff: 4.4e-07      
503 iterations  logLik: -396.17  diff: 4.3e-07      
504 iterations  logLik: -396.17  diff: 4.3e-07      
505 iterations  logLik: -396.17  diff: 4.2e-07      
506 iterations  logLik: -396.17  diff: 4.2e-07      
507 iterations  logLik: -396.17  diff: 4.1e-07      
508 iterations  logLik: -396.17  diff: 4.1e-07      
509 iterations  logLik: -396.17  diff: 4.0e-07      
510 iterations  logLik: -396.17  diff: 4.0e-07      
511 iterations  logLik: -396.17  diff: 3.9e-07      
512 iterations  logLik: -396.17  diff: 3.9e-07      
513 iterations  logLik: -396.17  diff: 3.9e-07      
514 iterations  logLik: -396.17  diff: 3.8e-07      
515 iterations  logLik: -396.17  diff: 3.8e-07      
516 iterations  logLik: -396.17  diff: 3.7e-07      
517 iterations  logLik: -396.17  diff: 3.7e-07      
518 iterations  logLik: -396.17  diff: 3.7e-07      
519 iterations  logLik: -396.17  diff: 3.6e-07      
520 iterations  logLik: -396.17  diff: 3.6e-07      
521 iterations  logLik: -396.17  diff: 3.5e-07      
522 iterations  logLik: -396.17  diff: 3.5e-07      
523 iterations  logLik: -396.17  diff: 3.5e-07      
524 iterations  logLik: -396.17  diff: 3.4e-07      
525 iterations  logLik: -396.17  diff: 3.4e-07      
526 iterations  logLik: -396.17  diff: 3.3e-07      
527 iterations  logLik: -396.17  diff: 3.3e-07      
528 iterations  logLik: -396.17  diff: 3.3e-07      
529 iterations  logLik: -396.17  diff: 3.2e-07      
530 iterations  logLik: -396.17  diff: 3.2e-07      
531 iterations  logLik: -396.17  diff: 3.2e-07      
532 iterations  logLik: -396.17  diff: 3.1e-07      
533 iterations  logLik: -396.17  diff: 3.1e-07      
534 iterations  logLik: -396.17  diff: 3.1e-07      
535 iterations  logLik: -396.17  diff: 3.0e-07      
536 iterations  logLik: -396.17  diff: 3.0e-07      
537 iterations  logLik: -396.17  diff: 3.0e-07      
538 iterations  logLik: -396.17  diff: 2.9e-07      
539 iterations  logLik: -396.17  diff: 2.9e-07      
540 iterations  logLik: -396.17  diff: 2.9e-07      
541 iterations  logLik: -396.17  diff: 2.8e-07      
542 iterations  logLik: -396.17  diff: 2.8e-07      
543 iterations  logLik: -396.17  diff: 2.8e-07      
544 iterations  logLik: -396.17  diff: 2.8e-07      
545 iterations  logLik: -396.17  diff: 2.7e-07      
546 iterations  logLik: -396.17  diff: 2.7e-07      
547 iterations  logLik: -396.17  diff: 2.7e-07      
548 iterations  logLik: -396.17  diff: 2.6e-07      
549 iterations  logLik: -396.17  diff: 2.6e-07      
550 iterations  logLik: -396.17  diff: 2.6e-07      
551 iterations  logLik: -396.17  diff: 2.6e-07      
552 iterations  logLik: -396.17  diff: 2.5e-07      
553 iterations  logLik: -396.17  diff: 2.5e-07      
554 iterations  logLik: -396.17  diff: 2.5e-07      
555 iterations  logLik: -396.17  diff: 2.5e-07      
556 iterations  logLik: -396.17  diff: 2.4e-07      
557 iterations  logLik: -396.17  diff: 2.4e-07      
558 iterations  logLik: -396.17  diff: 2.4e-07      
559 iterations  logLik: -396.17  diff: 2.4e-07      
560 iterations  logLik: -396.17  diff: 2.3e-07      
561 iterations  logLik: -396.17  diff: 2.3e-07      
562 iterations  logLik: -396.17  diff: 2.3e-07      
563 iterations  logLik: -396.17  diff: 2.3e-07      
564 iterations  logLik: -396.17  diff: 2.2e-07      
565 iterations  logLik: -396.17  diff: 2.2e-07      
566 iterations  logLik: -396.17  diff: 2.2e-07      
567 iterations  logLik: -396.17  diff: 2.2e-07      
568 iterations  logLik: -396.17  diff: 2.1e-07      
569 iterations  logLik: -396.17  diff: 2.1e-07      
570 iterations  logLik: -396.17  diff: 2.1e-07      
571 iterations  logLik: -396.17  diff: 2.1e-07      
572 iterations  logLik: -396.17  diff: 2.1e-07      
573 iterations  logLik: -396.17  diff: 2.0e-07      
574 iterations  logLik: -396.17  diff: 2.0e-07      
575 iterations  logLik: -396.17  diff: 2.0e-07      
576 iterations  logLik: -396.17  diff: 2.0e-07      
577 iterations  logLik: -396.17  diff: 1.9e-07      
578 iterations  logLik: -396.17  diff: 1.9e-07      
579 iterations  logLik: -396.17  diff: 1.9e-07      
580 iterations  logLik: -396.17  diff: 1.9e-07      
581 iterations  logLik: -396.17  diff: 1.9e-07      
582 iterations  logLik: -396.17  diff: 1.8e-07      
583 iterations  logLik: -396.17  diff: 1.8e-07      
584 iterations  logLik: -396.17  diff: 1.8e-07      
585 iterations  logLik: -396.17  diff: 1.8e-07      
586 iterations  logLik: -396.17  diff: 1.8e-07      
587 iterations  logLik: -396.17  diff: 1.8e-07      
588 iterations  logLik: -396.17  diff: 1.7e-07      
589 iterations  logLik: -396.17  diff: 1.7e-07      
590 iterations  logLik: -396.17  diff: 1.7e-07      
591 iterations  logLik: -396.17  diff: 1.7e-07      
592 iterations  logLik: -396.17  diff: 1.7e-07      
593 iterations  logLik: -396.17  diff: 1.6e-07      
594 iterations  logLik: -396.17  diff: 1.6e-07      
595 iterations  logLik: -396.17  diff: 1.6e-07      
596 iterations  logLik: -396.17  diff: 1.6e-07      
597 iterations  logLik: -396.17  diff: 1.6e-07      
598 iterations  logLik: -396.17  diff: 1.6e-07      
599 iterations  logLik: -396.17  diff: 1.5e-07      
600 iterations  logLik: -396.17  diff: 1.5e-07      
601 iterations  logLik: -396.17  diff: 1.5e-07      
602 iterations  logLik: -396.17  diff: 1.5e-07      
603 iterations  logLik: -396.17  diff: 1.5e-07      
604 iterations  logLik: -396.17  diff: 1.5e-07      
605 iterations  logLik: -396.17  diff: 1.5e-07      
606 iterations  logLik: -396.17  diff: 1.4e-07      
607 iterations  logLik: -396.17  diff: 1.4e-07      
608 iterations  logLik: -396.17  diff: 1.4e-07      
609 iterations  logLik: -396.17  diff: 1.4e-07      
610 iterations  logLik: -396.17  diff: 1.4e-07      
611 iterations  logLik: -396.17  diff: 1.4e-07      
612 iterations  logLik: -396.17  diff: 1.4e-07      
613 iterations  logLik: -396.17  diff: 1.3e-07      
614 iterations  logLik: -396.17  diff: 1.3e-07      
615 iterations  logLik: -396.17  diff: 1.3e-07      
616 iterations  logLik: -396.17  diff: 1.3e-07      
617 iterations  logLik: -396.17  diff: 1.3e-07      
618 iterations  logLik: -396.17  diff: 1.3e-07      
619 iterations  logLik: -396.17  diff: 1.3e-07      
620 iterations  logLik: -396.17  diff: 1.2e-07      
621 iterations  logLik: -396.17  diff: 1.2e-07      
622 iterations  logLik: -396.17  diff: 1.2e-07      
623 iterations  logLik: -396.17  diff: 1.2e-07      
624 iterations  logLik: -396.17  diff: 1.2e-07      
625 iterations  logLik: -396.17  diff: 1.2e-07      
626 iterations  logLik: -396.17  diff: 1.2e-07      
627 iterations  logLik: -396.17  diff: 1.2e-07      
628 iterations  logLik: -396.17  diff: 1.1e-07      
629 iterations  logLik: -396.17  diff: 1.1e-07      
630 iterations  logLik: -396.17  diff: 1.1e-07      
631 iterations  logLik: -396.17  diff: 1.1e-07      
632 iterations  logLik: -396.17  diff: 1.1e-07      
633 iterations  logLik: -396.17  diff: 1.1e-07      
634 iterations  logLik: -396.17  diff: 1.1e-07      
635 iterations  logLik: -396.17  diff: 1.1e-07      
636 iterations  logLik: -396.17  diff: 1.1e-07      
637 iterations  logLik: -396.17  diff: 1.0e-07      
638 iterations  logLik: -396.17  diff: 1.0e-07      
639 iterations  logLik: -396.17  diff: 1.0e-07      
640 iterations  logLik: -396.17  diff: 1.0e-07      
641 iterations  logLik: -396.17  diff: 1.0e-07      
642 iterations  logLik: -396.17  diff: 1.0e-07      
643 iterations  logLik: -396.17  diff: 9.9e-08      
644 iterations  logLik: -396.17  diff: 9.8e-08      
645 iterations  logLik: -396.17  diff: 9.7e-08      
646 iterations  logLik: -396.17  diff: 9.6e-08      
647 iterations  logLik: -396.17  diff: 9.5e-08      
648 iterations  logLik: -396.17  diff: 9.4e-08      
649 iterations  logLik: -396.17  diff: 9.3e-08      
650 iterations  logLik: -396.17  diff: 9.2e-08      
651 iterations  logLik: -396.17  diff: 9.1e-08      
652 iterations  logLik: -396.17  diff: 9.0e-08      
653 iterations  logLik: -396.17  diff: 8.9e-08      
654 iterations  logLik: -396.17  diff: 8.8e-08      
655 iterations  logLik: -396.17  diff: 8.7e-08      
656 iterations  logLik: -396.17  diff: 8.6e-08      
657 iterations  logLik: -396.17  diff: 8.6e-08      
658 iterations  logLik: -396.17  diff: 8.5e-08      
659 iterations  logLik: -396.17  diff: 8.4e-08      
660 iterations  logLik: -396.17  diff: 8.3e-08      
661 iterations  logLik: -396.17  diff: 8.2e-08      
662 iterations  logLik: -396.17  diff: 8.1e-08      
663 iterations  logLik: -396.17  diff: 8.1e-08      
664 iterations  logLik: -396.17  diff: 8.0e-08      
665 iterations  logLik: -396.17  diff: 7.9e-08      
666 iterations  logLik: -396.17  diff: 7.8e-08      
667 iterations  logLik: -396.17  diff: 7.7e-08      
668 iterations  logLik: -396.17  diff: 7.7e-08      
669 iterations  logLik: -396.17  diff: 7.6e-08      
670 iterations  logLik: -396.17  diff: 7.5e-08      
671 iterations  logLik: -396.17  diff: 7.4e-08      
672 iterations  logLik: -396.17  diff: 7.4e-08      
673 iterations  logLik: -396.17  diff: 7.3e-08      
674 iterations  logLik: -396.17  diff: 7.2e-08      
675 iterations  logLik: -396.17  diff: 7.1e-08      
676 iterations  logLik: -396.17  diff: 7.1e-08      
677 iterations  logLik: -396.17  diff: 7.0e-08      
678 iterations  logLik: -396.17  diff: 6.9e-08      
679 iterations  logLik: -396.17  diff: 6.9e-08      
680 iterations  logLik: -396.17  diff: 6.8e-08      
681 iterations  logLik: -396.17  diff: 6.7e-08      
682 iterations  logLik: -396.17  diff: 6.7e-08      
683 iterations  logLik: -396.17  diff: 6.6e-08      
684 iterations  logLik: -396.17  diff: 6.5e-08      
685 iterations  logLik: -396.17  diff: 6.5e-08      
686 iterations  logLik: -396.17  diff: 6.4e-08      
687 iterations  logLik: -396.17  diff: 6.3e-08      
688 iterations  logLik: -396.17  diff: 6.3e-08      
689 iterations  logLik: -396.17  diff: 6.2e-08      
690 iterations  logLik: -396.17  diff: 6.2e-08      
691 iterations  logLik: -396.17  diff: 6.1e-08      
692 iterations  logLik: -396.17  diff: 6.0e-08      
693 iterations  logLik: -396.17  diff: 6.0e-08      
694 iterations  logLik: -396.17  diff: 5.9e-08      
695 iterations  logLik: -396.17  diff: 5.9e-08      
696 iterations  logLik: -396.17  diff: 5.8e-08      
697 iterations  logLik: -396.17  diff: 5.7e-08      
698 iterations  logLik: -396.17  diff: 5.7e-08      
699 iterations  logLik: -396.17  diff: 5.6e-08      
700 iterations  logLik: -396.17  diff: 5.6e-08      
701 iterations  logLik: -396.17  diff: 5.5e-08      
702 iterations  logLik: -396.17  diff: 5.5e-08      
703 iterations  logLik: -396.17  diff: 5.4e-08      
704 iterations  logLik: -396.17  diff: 5.4e-08      
705 iterations  logLik: -396.17  diff: 5.3e-08      
706 iterations  logLik: -396.17  diff: 5.2e-08      
707 iterations  logLik: -396.17  diff: 5.2e-08      
708 iterations  logLik: -396.17  diff: 5.1e-08      
709 iterations  logLik: -396.17  diff: 5.1e-08      
710 iterations  logLik: -396.17  diff: 5.0e-08      
711 iterations  logLik: -396.17  diff: 5.0e-08      
712 iterations  logLik: -396.17  diff: 4.9e-08      
713 iterations  logLik: -396.17  diff: 4.9e-08      
714 iterations  logLik: -396.17  diff: 4.8e-08      
715 iterations  logLik: -396.17  diff: 4.8e-08      
716 iterations  logLik: -396.17  diff: 4.8e-08      
717 iterations  logLik: -396.17  diff: 4.7e-08      
718 iterations  logLik: -396.17  diff: 4.7e-08      
719 iterations  logLik: -396.17  diff: 4.6e-08      
720 iterations  logLik: -396.17  diff: 4.6e-08      
721 iterations  logLik: -396.17  diff: 4.5e-08      
722 iterations  logLik: -396.17  diff: 4.5e-08      
723 iterations  logLik: -396.17  diff: 4.4e-08      
724 iterations  logLik: -396.17  diff: 4.4e-08      
725 iterations  logLik: -396.17  diff: 4.4e-08      
726 iterations  logLik: -396.17  diff: 4.3e-08      
727 iterations  logLik: -396.17  diff: 4.3e-08      
728 iterations  logLik: -396.17  diff: 4.2e-08      
729 iterations  logLik: -396.17  diff: 4.2e-08      
730 iterations  logLik: -396.17  diff: 4.1e-08      
731 iterations  logLik: -396.17  diff: 4.1e-08      
732 iterations  logLik: -396.17  diff: 4.1e-08      
733 iterations  logLik: -396.17  diff: 4.0e-08      
734 iterations  logLik: -396.17  diff: 4.0e-08      
735 iterations  logLik: -396.17  diff: 3.9e-08      
736 iterations  logLik: -396.17  diff: 3.9e-08      
737 iterations  logLik: -396.17  diff: 3.9e-08      
738 iterations  logLik: -396.17  diff: 3.8e-08      
739 iterations  logLik: -396.17  diff: 3.8e-08      
740 iterations  logLik: -396.17  diff: 3.8e-08      
741 iterations  logLik: -396.17  diff: 3.7e-08      
742 iterations  logLik: -396.17  diff: 3.7e-08      
743 iterations  logLik: -396.17  diff: 3.6e-08      
744 iterations  logLik: -396.17  diff: 3.6e-08      
745 iterations  logLik: -396.17  diff: 3.6e-08      
746 iterations  logLik: -396.17  diff: 3.5e-08      
747 iterations  logLik: -396.17  diff: 3.5e-08      
748 iterations  logLik: -396.17  diff: 3.5e-08      
749 iterations  logLik: -396.17  diff: 3.4e-08      
750 iterations  logLik: -396.17  diff: 3.4e-08      
751 iterations  logLik: -396.17  diff: 3.4e-08      
752 iterations  logLik: -396.17  diff: 3.3e-08      
753 iterations  logLik: -396.17  diff: 3.3e-08      
754 iterations  logLik: -396.17  diff: 3.3e-08      
755 iterations  logLik: -396.17  diff: 3.2e-08      
756 iterations  logLik: -396.17  diff: 3.2e-08      
757 iterations  logLik: -396.17  diff: 3.2e-08      
758 iterations  logLik: -396.17  diff: 3.1e-08      
759 iterations  logLik: -396.17  diff: 3.1e-08      
760 iterations  logLik: -396.17  diff: 3.1e-08      
761 iterations  logLik: -396.17  diff: 3.1e-08      
762 iterations  logLik: -396.17  diff: 3.0e-08      
763 iterations  logLik: -396.17  diff: 3.0e-08      
764 iterations  logLik: -396.17  diff: 3.0e-08      
765 iterations  logLik: -396.17  diff: 2.9e-08      
766 iterations  logLik: -396.17  diff: 2.9e-08      
767 iterations  logLik: -396.17  diff: 2.9e-08      
768 iterations  logLik: -396.17  diff: 2.9e-08      
769 iterations  logLik: -396.17  diff: 2.8e-08      
770 iterations  logLik: -396.17  diff: 2.8e-08      
771 iterations  logLik: -396.17  diff: 2.8e-08      
772 iterations  logLik: -396.17  diff: 2.7e-08      
773 iterations  logLik: -396.17  diff: 2.7e-08      
774 iterations  logLik: -396.17  diff: 2.7e-08      
775 iterations  logLik: -396.17  diff: 2.7e-08      
776 iterations  logLik: -396.17  diff: 2.6e-08      
777 iterations  logLik: -396.17  diff: 2.6e-08      
778 iterations  logLik: -396.17  diff: 2.6e-08      
779 iterations  logLik: -396.17  diff: 2.6e-08      
780 iterations  logLik: -396.17  diff: 2.5e-08      
781 iterations  logLik: -396.17  diff: 2.5e-08      
782 iterations  logLik: -396.17  diff: 2.5e-08      
783 iterations  logLik: -396.17  diff: 2.5e-08      
784 iterations  logLik: -396.17  diff: 2.4e-08      
785 iterations  logLik: -396.17  diff: 2.4e-08      
786 iterations  logLik: -396.17  diff: 2.4e-08      
787 iterations  logLik: -396.17  diff: 2.4e-08      
788 iterations  logLik: -396.17  diff: 2.4e-08      
789 iterations  logLik: -396.17  diff: 2.3e-08      
790 iterations  logLik: -396.17  diff: 2.3e-08      
791 iterations  logLik: -396.17  diff: 2.3e-08      
792 iterations  logLik: -396.17  diff: 2.3e-08      
793 iterations  logLik: -396.17  diff: 2.2e-08      
794 iterations  logLik: -396.17  diff: 2.2e-08      
795 iterations  logLik: -396.17  diff: 2.2e-08      
796 iterations  logLik: -396.17  diff: 2.2e-08      
797 iterations  logLik: -396.17  diff: 2.2e-08      
798 iterations  logLik: -396.17  diff: 2.1e-08      
799 iterations  logLik: -396.17  diff: 2.1e-08      
800 iterations  logLik: -396.17  diff: 2.1e-08      
801 iterations  logLik: -396.17  diff: 2.1e-08      
802 iterations  logLik: -396.17  diff: 2.1e-08      
803 iterations  logLik: -396.17  diff: 2.0e-08      
804 iterations  logLik: -396.17  diff: 2.0e-08      
805 iterations  logLik: -396.17  diff: 2.0e-08      
806 iterations  logLik: -396.17  diff: 2.0e-08      
807 iterations  logLik: -396.17  diff: 2.0e-08      
808 iterations  logLik: -396.17  diff: 1.9e-08      
809 iterations  logLik: -396.17  diff: 1.9e-08      
810 iterations  logLik: -396.17  diff: 1.9e-08      
811 iterations  logLik: -396.17  diff: 1.9e-08      
812 iterations  logLik: -396.17  diff: 1.9e-08      
813 iterations  logLik: -396.17  diff: 1.8e-08      
814 iterations  logLik: -396.17  diff: 1.8e-08      
815 iterations  logLik: -396.17  diff: 1.8e-08      
816 iterations  logLik: -396.17  diff: 1.8e-08      
817 iterations  logLik: -396.17  diff: 1.8e-08      
818 iterations  logLik: -396.17  diff: 1.8e-08      
819 iterations  logLik: -396.17  diff: 1.7e-08      
820 iterations  logLik: -396.17  diff: 1.7e-08      
821 iterations  logLik: -396.17  diff: 1.7e-08      
822 iterations  logLik: -396.17  diff: 1.7e-08      
823 iterations  logLik: -396.17  diff: 1.7e-08      
824 iterations  logLik: -396.17  diff: 1.7e-08      
825 iterations  logLik: -396.17  diff: 1.6e-08      
826 iterations  logLik: -396.17  diff: 1.6e-08      
827 iterations  logLik: -396.17  diff: 1.6e-08      
828 iterations  logLik: -396.17  diff: 1.6e-08      
829 iterations  logLik: -396.17  diff: 1.6e-08      
830 iterations  logLik: -396.17  diff: 1.6e-08      
831 iterations  logLik: -396.17  diff: 1.6e-08      
832 iterations  logLik: -396.17  diff: 1.5e-08      
833 iterations  logLik: -396.17  diff: 1.5e-08      
834 iterations  logLik: -396.17  diff: 1.5e-08      
835 iterations  logLik: -396.17  diff: 1.5e-08      
836 iterations  logLik: -396.17  diff: 1.5e-08      
837 iterations  logLik: -396.17  diff: 1.5e-08      
838 iterations  logLik: -396.17  diff: 1.5e-08      
839 iterations  logLik: -396.17  diff: 1.4e-08      
840 iterations  logLik: -396.17  diff: 1.4e-08      
841 iterations  logLik: -396.17  diff: 1.4e-08      
842 iterations  logLik: -396.17  diff: 1.4e-08      
843 iterations  logLik: -396.17  diff: 1.4e-08      
844 iterations  logLik: -396.17  diff: 1.4e-08      
845 iterations  logLik: -396.17  diff: 1.4e-08      
846 iterations  logLik: -396.17  diff: 1.3e-08      
847 iterations  logLik: -396.17  diff: 1.3e-08      
848 iterations  logLik: -396.17  diff: 1.3e-08      
849 iterations  logLik: -396.17  diff: 1.3e-08      
850 iterations  logLik: -396.17  diff: 1.3e-08      
851 iterations  logLik: -396.17  diff: 1.3e-08      
852 iterations  logLik: -396.17  diff: 1.3e-08      
853 iterations  logLik: -396.17  diff: 1.3e-08      
854 iterations  logLik: -396.17  diff: 1.2e-08      
855 iterations  logLik: -396.17  diff: 1.2e-08      
856 iterations  logLik: -396.17  diff: 1.2e-08      
857 iterations  logLik: -396.17  diff: 1.2e-08      
858 iterations  logLik: -396.17  diff: 1.2e-08      
859 iterations  logLik: -396.17  diff: 1.2e-08      
860 iterations  logLik: -396.17  diff: 1.2e-08      
861 iterations  logLik: -396.17  diff: 1.2e-08      
862 iterations  logLik: -396.17  diff: 1.2e-08      
863 iterations  logLik: -396.17  diff: 1.1e-08      
864 iterations  logLik: -396.17  diff: 1.1e-08      
865 iterations  logLik: -396.17  diff: 1.1e-08      
866 iterations  logLik: -396.17  diff: 1.1e-08      
867 iterations  logLik: -396.17  diff: 1.1e-08      
868 iterations  logLik: -396.17  diff: 1.1e-08      
869 iterations  logLik: -396.17  diff: 1.1e-08      
870 iterations  logLik: -396.17  diff: 1.1e-08      
871 iterations  logLik: -396.17  diff: 1.1e-08      
872 iterations  logLik: -396.17  diff: 1.0e-08      
873 iterations  logLik: -396.17  diff: 1.0e-08      
874 iterations  logLik: -396.17  diff: 1.0e-08      
875 iterations  logLik: -396.17  diff: 1.0e-08      
876 iterations  logLik: -396.17  diff: 1.0e-08      
877 iterations  logLik: -396.17  diff: 1.0e-08      

#> .. done. 
#> nlm iteration begin.
#> .. done.
#> Structural Latent Variable Model
#> 
#> Latent variables (Root*):           
#>  Label: lc*
#> nclass: 3  
#> 
#> Measurement model:                            
#>  lc -> { y1, y2, y3, y4 }  a
#> 
#> LogLik: -396.1722

# Parameter restriction
mf <- estimate(m, dt$response)
param(mf, index = TRUE)
#> PI :
#> (lc)
#>   class
#>              1           2           3
#>     0.3123 (1)  0.3295 (2)  0.3582 (3)
#> 
#> RHO :
#> (a)
#>         class
#> response            1            2            3
#>    1(V1)   0.1678 (4)  0.9689 (12)  0.8329 (20)
#>    2       0.8322 (5)  0.0311 (13)  0.1671 (21)
#>    1(V2)   0.0000 (6)  0.9463 (14)  1.0000 (22)
#>    2       1.0000 (7)  0.0537 (15)  0.0000 (23)
#>    1(V3)   0.0492 (8)  0.8358 (16)  0.1235 (24)
#>    2       0.9508 (9)  0.1642 (17)  0.8765 (25)
#>    1(V4)  0.1357 (10)  0.9791 (18)  0.0000 (26)
#>    2      0.8643 (11)  0.0209 (19)  1.0000 (27)
#> 
#>    V1 V2 V3 V4
#> lc y1 y2 y3 y4
mf0 <- estimate(mf, fix2zero = c(4, 6, 8, 10))
param(mf0)
#> PI :
#> (lc)
#>   class
#>          1       2       3
#>     0.2005  0.4529  0.3466
#> 
#> RHO :
#> (a)
#>         class
#> response       1       2       3
#>    1(V1)  0.0000  0.8488  0.8239
#>    2      1.0000  0.1512  0.1761
#>    1(V2)  0.0000  0.8256  0.8543
#>    2      1.0000  0.1744  0.1457
#>    1(V3)  0.0000  0.7397  0.0000
#>    2      1.0000  0.2603  1.0000
#>    1(V4)  0.0000  0.8060  0.0000
#>    2      1.0000  0.1940  1.0000
#> 
#>    V1 V2 V3 V4
#> lc y1 y2 y3 y4

# Estimation control
estimate(m, dt$response, control = slcaControl(nrep = 3, verbose = TRUE))
#> Inital parameter test: 
#> 1 / 3  logLik: -396.1722 
#> 2 / 3  logLik: -396.1722 
#> 3 / 3  logLik: -396.1722 
#> 
#> 1th parameter set has been selected.
#> 
#> EM iteration begin.
#> 1 iterations  logLik: -396.17  diff: inf      
2 iterations  logLik: -396.17  diff: 8.3e-09      

#> .. done.
#> Structural Latent Variable Model
#> 
#> Latent variables (Root*):           
#>  Label: lc*
#> nclass: 3  
#> 
#> Measurement model:                            
#>  lc -> { y1, y2, y3, y4 }  a
#> 
#> LogLik: -396.1722
estimate(m, dt$response, control = slcaControl(init.param = c(pi, rho)))
#> Structural Latent Variable Model
#> 
#> Latent variables (Root*):           
#>  Label: lc*
#> nclass: 3  
#> 
#> Measurement model:                            
#>  lc -> { y1, y2, y3, y4 }  a
#> 
#> LogLik: -396.1722