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 theslca
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 thenlm
function, while"hybrid"
combines both approaches, starting with EM and finishing withnlm
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 usingslcaControl()
.
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