R version 3.6.3 (2020-02-29) -- "Holding the Windsock"
Copyright (C) 2020 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> ############################################################################
> #     MLwiN MCMC Manual
> #
> # 19  Mixed Response Models and Correlated Residuals . . . . . . . . . . 287
> #
> #     Browne, W.J. (2009) MCMC Estimation in MLwiN, v2.13. Centre for
> #     Multilevel Modelling, University of Bristol.
> ############################################################################
> #     R script to replicate all analyses using R2MLwiN
> #
> #     Zhang, Z., Charlton, C., Parker, R, Leckie, G., and Browne, W.J.
> #     Centre for Multilevel Modelling, 2012
> #     http://www.bristol.ac.uk/cmm/software/R2MLwiN/
> ############################################################################
> 
> # 19.1 Mixed response models . . . . . . . . . . . . . . . . . . . . . . 287
> 
> # 19.2 The JSP mixed response example . . . . . . . . . . . . . . . . . .289
> 
> library(R2MLwiN)

R2MLwiN: A package to run models implemented in MLwiN from R
Copyright 2013-2017 Zhengzheng Zhang, Christopher M. J. Charlton, Richard M. A. Parker, William J. Browne and George Leckie

Support provided by the Economic and Social Research Council (ESRC)
(Grants RES-149-25-1084, RES-576-25-0032 and ES/K007246/1)

To cite R2MLwiN in publications use:

  Zhengzheng Zhang, Richard M. A. Parker, Christopher M. J. Charlton,
  George Leckie, William J. Browne (2016). R2MLwiN: A Package to Run
  MLwiN from within R. Journal of Statistical Software, 72(10), 1-43.
  doi:10.18637/jss.v072.i10

A BibTeX entry for LaTeX users is

  @Article{,
    title = {{R2MLwiN}: A Package to Run {MLwiN} from within {R}},
    author = {Zhengzheng Zhang and Richard M. A. Parker and Christopher M. J. Charlton and George Leckie and William J. Browne},
    journal = {Journal of Statistical Software},
    year = {2016},
    volume = {72},
    number = {10},
    pages = {1--43},
    doi = {10.18637/jss.v072.i10},
  }

The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.05/
To change this use: options(MLwiN_path="<path to MLwiN>")
> # MLwiN folder
> mlwin <- getOption("MLwiN_path")
> while (!file.access(mlwin, mode = 1) == 0) {
+   cat("Please specify the root MLwiN folder or the full path to the MLwiN executable:\n")
+   mlwin <- scan(what = character(0), sep = "\n")
+   mlwin <- gsub("\\", "/", mlwin, fixed = TRUE)
+ }
> options(MLwiN_path = mlwin)
> 
> # User's input if necessary
> 
> ## Read jspmix1 data
> data(jspmix1, package = "R2MLwiN")
> 
> tab1 <- matrix(, 3, 3)
> colnames(tab1) <- c("0", "1", "TOTALS")
> rownames(tab1) <- c("N", "MEANS", "SDs")
> tab1[1, 1:2] <- colSums(table(jspmix1$english, jspmix1$behaviour))
> tab1[1, 3] <- sum(tab1[1, 1:2])
> tab1[2, 1] <- mean(jspmix1$english[jspmix1$behaviour == 0])
> tab1[2, 2] <- mean(jspmix1$english[jspmix1$behaviour == 1])
> tab1[2, 3] <- mean(jspmix1$english)
> tab1[3, 1] <- sd(jspmix1$english[jspmix1$behaviour == 0])
> tab1[3, 2] <- sd(jspmix1$english[jspmix1$behaviour == 1])
> tab1[3, 3] <- sd(jspmix1$english)
> formatC(tab1)
      0     1     TOTALS 
N     "248" "871" "1119" 
MEANS "NaN" "NaN" "41.47"
SDs   " NA" " NA" "21.56"
> 
> vars <- cbind(as.numeric(jspmix1$sex) - 1, jspmix1$fluent, jspmix1$ravens, jspmix1$english, as.numeric(jspmix1$behaviour) - 
+   1)
> colnames(vars) <- c("sex", "fluent", "ravens", "english", "behaviour")
> round(cor(vars), 4)
              sex  fluent ravens english behaviour
sex        1.0000 -0.0229 0.0341 -0.1479   -0.1272
fluent    -0.0229  1.0000 0.1705  0.2054   -0.0038
ravens     0.0341  0.1705 1.0000  0.5042    0.2181
english   -0.1479  0.2054 0.5042  1.0000    0.3122
behaviour -0.1272 -0.0038 0.2181  0.3122    1.0000
> 
> # 19.3 Setting up a single level mixed response model . . . . . . . . . .291
> 
> (mymodel <- runMLwiN(c(english, probit(behaviour)) ~ 1 + sex + ravens + fluent[1] + (1[1] | id), D = c("Mixed", 
+   "Normal", "Binomial"), estoptions = list(EstM = 1, mcmcMeth = list(fixM = 1, residM = 1, Lev1VarM = 1), sort.ignore = TRUE), 
+   data = jspmix1))
MLwiN is running, please wait......
/nogui option ignored
ECHO    0


Echoing is ON
BATC 1

Batch mode is ON
MAXI 2
STAR
iteration 0
iteration 1

Convergence not achieved
TOLE 2
MAXI 20
BATC 1

Batch mode is ON
NEXT
iteration 2
iteration 3
iteration 4

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 5
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 12805.994829
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 12701.732578
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 12724.390456
Adapting finished and took 400 iterations
Adapting took 400 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 12644.886667
Burning in for 100 iterations out of 500
-2 * Loglike = 12677.914744
Burning in for 150 iterations out of 500
-2 * Loglike = 12519.580841
Burning in for 200 iterations out of 500
-2 * Loglike = 12667.799806
Burning in for 250 iterations out of 500
-2 * Loglike = 12634.161580
Burning in for 300 iterations out of 500
-2 * Loglike = 12694.642713
Burning in for 350 iterations out of 500
-2 * Loglike = 12719.620791
Burning in for 400 iterations out of 500
-2 * Loglike = 12658.786407
Burning in for 450 iterations out of 500
-2 * Loglike = 12592.747762
Burning in for 500 iterations out of 500
-2 * Loglike = 12746.360304

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 5
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 12796.522205
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 12656.064815
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 12694.258119
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 12606.728281
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 12577.775711
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 12698.818511
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 12569.306425
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 12789.942996
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 12610.843198
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 12682.706838
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 12597.893347
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 12559.048824
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 12568.441026
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 12765.174713
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 12658.175428
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 12683.844532
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 12610.634044
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 12752.638928
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 12649.882380
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 12646.046502
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 12647.973320
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 12747.748707
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 12657.026643
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 12705.140004
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 12675.142654
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 12717.794227
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 12667.445077
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 12630.563382
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 12705.050444
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 12676.193995
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 12742.408087
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 12753.110196
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 12704.490821
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 12667.289229
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 12773.461726
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 12801.352394
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 12695.444529
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 12796.098093
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 12682.202371
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 12681.172716
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 12698.496149
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 12671.863176
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 12663.311681
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 12727.171161
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 12666.202639
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 12652.323066
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 12684.371903
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 12771.840974
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 12703.329278
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 12688.775234
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 12664.284521
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 12640.038178
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 12699.283669
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 12649.062300
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 12658.294228
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 12652.610130
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 12613.582467
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 12637.030677
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 12671.673168
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 12660.789255
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 12696.443953
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 12745.679444
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 12667.638189
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 12724.881927
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 12707.733312
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 12681.875177
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 12678.754785
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 12732.905554
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 12777.514863
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 12685.130835
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 12737.719115
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 12718.772062
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 12696.546029
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 12777.511165
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 12747.208115
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 12651.180201
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 12768.306785
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 12760.771298
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 12784.285065
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 12726.061394
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 12795.379811
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 12642.817884
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 12787.906320
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 12726.661034
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 12575.547200
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 12668.978327
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 12666.411092
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 12622.720101
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 12621.221227
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 12753.835390
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 12773.303382
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 12662.229789
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 12712.203349
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 12726.830417
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 12741.282427
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 12737.291126
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 12730.633125
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 12815.318382
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 12666.945986
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 12778.902232
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          12698.
S.D.    =          54.404
S.E.M.  =         0.76939

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  multilevel model (Mixed) 
Estimation algorithm:  MCMC      Elapsed time : 3.03s 
Number of obs:  1119 (from total 1119)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Deviance statistic:  NA 
--------------------------------------------------------------------------------------------------- 
The model formula:
c(english, probit(behaviour)) ~ 1 + sex + ravens + fluent[1] + 
    (1[1] | id)
Level 1: id      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
                        Coef.   Std. Err.      z    Pr(>|z|)       [95% Cred.   Interval]    ESS 
sexmale_english      -6.96715     1.09320  -6.37   1.851e-10  ***    -9.09428    -4.82603   5000 
sexmale_behaviour    -0.39743     0.08626  -4.61   4.079e-06  ***    -0.56673    -0.22984   1497 
ravens_english        1.74997     0.09191  19.04   8.063e-81  ***     1.56448     1.92459   5000 
ravens_behaviour      0.05622     0.00727   7.73   1.066e-14  ***     0.04208     0.07034   1382 
fluent_1              6.77051     1.31461   5.15   2.602e-07  ***     4.27870     9.34429   4315 
Intercept_1         -11.69678     3.19541  -3.66   0.0002517  ***   -17.87894    -5.36995   4601 
Intercept_2          -0.38432     0.17993  -2.14     0.03269  *      -0.72889    -0.03304   1715 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the id level: 
                              Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_bcons_2                 1.00000     0.00001      1.00000     1.00000   5000 
cov_bcons_2_Intercept_1     6.22678     0.71826      4.80217     7.68186    371 
var_Intercept_1           329.54290    13.85086    303.56265   357.43881    886 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 19.4 Multilevel mixed response model . . . . . . . . . . . . . . . . . 294
> 
> (mymodel <- runMLwiN(c(english, probit(behaviour)) ~ 1 + sex + ravens + fluent[1] + (1 | school) + (1[1] | id), 
+   D = c("Mixed", "Normal", "Binomial"), estoptions = list(EstM = 1, mcmcMeth = list(fixM = 1, residM = 1, Lev1VarM = 1)), 
+   data = jspmix1))
MLwiN is running, please wait......
/nogui option ignored
ECHO    0


Echoing is ON
BATC 1

Batch mode is ON
MAXI 2
STAR
iteration 0
iteration 1

Convergence not achieved
TOLE 2
MAXI 20
BATC 1

Batch mode is ON
NEXT
iteration 2
iteration 3
iteration 4
iteration 5
iteration 6
iteration 7
iteration 8

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 5
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 12525.658322
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 12593.457854
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 12542.630286
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 12467.925995
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 12569.865026
Adapting for 600 iterations (Maximum 5000)
-2 * Loglike = 12561.878881
Adapting for 700 iterations (Maximum 5000)
-2 * Loglike = 12485.786776
Adapting finished and took 800 iterations
Adapting took 800 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 12561.833984
Burning in for 100 iterations out of 500
-2 * Loglike = 12590.811909
Burning in for 150 iterations out of 500
-2 * Loglike = 12484.819059
Burning in for 200 iterations out of 500
-2 * Loglike = 12582.254061
Burning in for 250 iterations out of 500
-2 * Loglike = 12572.319424
Burning in for 300 iterations out of 500
-2 * Loglike = 12568.440828
Burning in for 350 iterations out of 500
-2 * Loglike = 12566.051808
Burning in for 400 iterations out of 500
-2 * Loglike = 12551.209439
Burning in for 450 iterations out of 500
-2 * Loglike = 12523.790179
Burning in for 500 iterations out of 500
-2 * Loglike = 12562.427806

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 5
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 12588.808395
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 12510.453167
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 12664.435852
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 12542.722437
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 12522.155920
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 12578.344955
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 12566.853411
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 12465.874142
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 12584.190993
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 12602.940999
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 12445.486613
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 12537.055338
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 12617.291428
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 12462.995124
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 12545.681692
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 12506.899370
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 12474.618203
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 12554.530305
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 12566.010778
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 12519.719816
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 12450.965390
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 12564.980219
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 12499.566213
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 12562.123366
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 12560.872773
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 12468.739685
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 12508.356907
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 12477.288177
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 12422.673899
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 12602.517576
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 12581.578312
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 12585.553877
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 12549.691145
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 12626.643663
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 12586.903014
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 12484.665680
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 12486.328551
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 12497.814341
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 12487.421476
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 12488.144779
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 12465.074512
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 12546.661222
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 12507.636723
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 12608.623492
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 12500.235270
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 12605.135751
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 12643.360000
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 12541.408777
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 12646.548957
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 12598.267247
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 12456.516485
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 12486.544571
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 12523.251366
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 12630.881340
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 12576.997843
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 12579.945712
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 12506.448137
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 12510.190371
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 12472.260589
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 12464.480791
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 12393.313306
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 12533.290687
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 12562.288463
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 12499.309589
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 12526.616523
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 12581.839372
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 12487.938207
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 12479.313460
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 12545.773577
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 12505.310953
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 12452.990205
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 12570.980498
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 12647.585685
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 12569.762979
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 12522.044375
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 12568.953657
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 12636.900420
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 12502.428452
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 12440.901754
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 12656.500634
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 12610.232973
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 12619.856266
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 12544.921528
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 12562.320164
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 12549.874046
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 12546.631284
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 12468.626822
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 12653.896506
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 12531.798163
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 12566.774518
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 12398.839179
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 12518.103460
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 12475.833195
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 12485.950731
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 12503.389621
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 12641.030255
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 12471.456858
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 12560.844191
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 12515.566571
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 12559.674437
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          12540.
S.D.    =          55.432
S.E.M.  =         0.78393

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  multilevel model (Mixed) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 47   7 23.80851  76         47            7      23.80851           76
Estimation algorithm:  MCMC      Elapsed time : 4.44s 
Number of obs:  1119 (from total 1119)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Deviance statistic:  NA 
--------------------------------------------------------------------------------------------------- 
The model formula:
c(english, probit(behaviour)) ~ 1 + sex + ravens + fluent[1] + 
    (1 | school) + (1[1] | id)
Level 2: school     Level 1: id      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
                       Coef.   Std. Err.      z    Pr(>|z|)       [95% Cred.   Interval]    ESS 
sexmale_english     -6.25128     1.04875  -5.96   2.512e-09  ***    -8.24857    -4.14073   5000 
sexmale_behaviour   -0.41709     0.08993  -4.64   3.515e-06  ***    -0.59170    -0.23841   1362 
ravens_english       1.66182     0.09105  18.25   2.026e-74  ***     1.48512     1.83777   4083 
ravens_behaviour     0.05683     0.00772   7.36   1.791e-13  ***     0.04202     0.07219   1155 
fluent_1             6.37932     1.28503   4.96   6.894e-07  ***     3.90506     8.87736   4348 
Intercept_1         -9.27886     3.29575  -2.82    0.004872  **    -15.81773    -2.67064   2931 
Intercept_2         -0.35566     0.19942  -1.78     0.07451  .      -0.74788     0.03255   1185 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the school level: 
                                 Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept_1               40.87154    11.81184     22.55033    68.01690   1700 
cov_Intercept_1_Intercept_2    0.03318     0.42400     -0.79115     0.91364    694 
var_Intercept_2                0.06268     0.03102      0.01903     0.13664    363 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the id level: 
                              Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_bcons_2                 1.00000     0.00001      1.00000     1.00000   5000 
cov_bcons_2_Intercept_1     6.22421     0.70235      4.85780     7.59930    345 
var_Intercept_1           291.33014    12.34688    268.03250   316.60036    782 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 19.5 Rats dataset . . . . . . . . . . . . . . . . . . . . . . . . . . .295
> 
> ## Read rats data
> data(rats, package = "R2MLwiN")
> 
> (mymodel <- runMLwiN(c(y8, y15, y22, y29, y36) ~ 1 + (1 | rat), D = "Multivariate Normal", estoptions = list(EstM = 1), 
+   data = rats))
MLwiN is running, please wait......
/nogui option ignored
ECHO    0


Echoing is ON
BATC 1

Batch mode is ON
MAXI 2
STAR
iteration 0
iteration 1

Convergence achieved
TOLE 2
MAXI 20
BATC 1

Batch mode is ON
NEXT
iteration 2

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 4
Burning in for 50 iterations out of 500
-2 * Loglike = 1009.437068
Burning in for 100 iterations out of 500
-2 * Loglike = 1006.298254
Burning in for 150 iterations out of 500
-2 * Loglike = 1008.607885
Burning in for 200 iterations out of 500
-2 * Loglike = 1005.985450
Burning in for 250 iterations out of 500
-2 * Loglike = 1014.145045
Burning in for 300 iterations out of 500
-2 * Loglike = 1004.938369
Burning in for 350 iterations out of 500
-2 * Loglike = 1025.544161
Burning in for 400 iterations out of 500
-2 * Loglike = 1011.374624
Burning in for 450 iterations out of 500
-2 * Loglike = 1005.990961
Burning in for 500 iterations out of 500
-2 * Loglike = 1013.121401

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 4
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 1005.847619
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 1000.383434
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 1004.012006
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 1007.236221
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 999.010435
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 1017.094345
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 995.171524
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 1003.569097
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 1004.595807
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 1006.079451
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 1020.354743
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 1001.352172
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 1005.268099
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 1003.690802
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 1016.582578
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 1014.114350
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 1006.188743
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 1016.805220
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 1007.322041
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 1014.086949
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 1011.356589
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 1011.278404
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 1005.294994
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 1012.864830
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 1018.141424
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 1002.441202
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 1007.791599
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 1006.282787
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 1002.465128
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 1008.373744
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 998.012181
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 1011.666889
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 1012.384353
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 1007.227055
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 1005.776334
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 998.252170
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 1013.941667
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 1000.592760
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 1043.017936
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 1008.574665
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 1008.853444
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 1007.879983
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 998.853832
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 997.792486
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 1009.573882
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 1005.072069
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 1014.309930
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 1009.381119
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 998.066328
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 1003.244556
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 1012.617477
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 1011.178643
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 999.873629
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 1000.019906
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 1002.187112
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 1014.454894
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 1003.814192
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 1011.189844
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 1009.392813
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 1012.664242
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 1003.039396
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 1000.360655
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 1005.715447
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 1002.291413
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 1017.965747
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 1001.340079
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 1012.673352
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 1006.586695
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 1004.155966
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 1012.705577
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 1016.562440
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 1002.849180
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 1007.336922
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 1012.578390
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 1016.766568
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 1015.016906
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 1014.993475
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 1006.696834
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 1013.314038
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 1012.698703
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 1019.858500
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 997.735015
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 1009.535942
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 1009.256408
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 1011.503094
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 1017.411966
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 1011.252261
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 1005.894385
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 1003.128223
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 1005.864592
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 995.244438
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 1001.822149
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 1015.106792
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 1017.796146
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 1009.323841
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 1016.696885
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 1018.440639
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 1007.334206
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 1005.905818
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 1004.216268
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          1008.8
S.D.    =          6.0605
S.E.M.  =        0.085709

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  multilevel model (Multivariate Normal) 
Estimation algorithm:  MCMC      Elapsed time : 0.37s 
Number of obs:  30 (from total 30)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
1008.798   993.349    15.449     1024.246   
--------------------------------------------------------------------------------------------------- 
The model formula:
c(y8, y15, y22, y29, y36) ~ 1 + (1 | rat)
Level 1: rat      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
                    Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept_y8    152.15244     2.20921  68.87          0  ***   147.73541   156.45141   5418 
Intercept_y15   201.74869     2.51077  80.35          0  ***   196.71465   206.70704   5095 
Intercept_y22   244.98005     3.02945  80.87          0  ***   238.86493   250.89769   5057 
Intercept_y29   289.42062     3.79411  76.28          0  ***   281.92528   296.81249   5224 
Intercept_y36   324.70046     3.86453  84.02          0  ***   317.22035   332.38677   5222 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the rat level: 
                                      Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept_y8                  150.54482    41.72600     89.02023   252.82079   4670 
cov_Intercept_y8_Intercept_y15    156.76883    45.79837     90.00509   268.72897   4544 
var_Intercept_y15                 194.00095    54.64486    115.22264   328.33698   3256 
cov_Intercept_y8_Intercept_y22    170.61354    53.13663     93.99026   301.82521   4549 
cov_Intercept_y15_Intercept_y22   216.67225    63.84111    124.67072   375.76304   4513 
var_Intercept_y22                 284.80724    80.51540    169.55108   482.66966   4547 
cov_Intercept_y8_Intercept_y29    187.41159    63.35202     94.61906   342.49951   4554 
cov_Intercept_y15_Intercept_y29   247.17615    76.65214    136.53310   438.54704   4586 
cov_Intercept_y22_Intercept_y29   336.73233    97.73940    195.26309   576.76341   4616 
var_Intercept_y29                 443.26810   125.17855    262.92079   743.86706   4695 
cov_Intercept_y8_Intercept_y36    161.00551    60.74671     67.94944   309.01290   4597 
cov_Intercept_y15_Intercept_y36   208.38045    72.35116    102.60791   385.91683   4613 
cov_Intercept_y22_Intercept_y36   309.57136    94.68515    171.94331   537.56497   4645 
cov_Intercept_y29_Intercept_y36   413.91606   121.63624    238.59326   698.03954   4706 
var_Intercept_y36                 457.05567   128.41577    275.64905   756.37108   4770 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> sixway(mymodel@chains[, "RP1_var_Intercept_y8", drop = FALSE], "sigma2u0")
> 
> covM1 <- matrix(, 5, 5)
> colnames(covM1) <- rownames(covM1) <- c("cons.y8", "cons.y15", "cons.y22", "cons.y29", "cons.y36")
> covM1[upper.tri(covM1, diag = TRUE)] <- mymodel@RP
> # covM1[lower.tri(covM1)] <- t(covM1)[lower.tri(covM1)]
> round(t(covM1), 3)
         cons.y8 cons.y15 cons.y22 cons.y29 cons.y36
cons.y8  150.545       NA       NA       NA       NA
cons.y15 156.769  194.001       NA       NA       NA
cons.y22 170.614  216.672  284.807       NA       NA
cons.y29 187.412  247.176  336.732  443.268       NA
cons.y36 161.006  208.380  309.571  413.916  457.056
> round(cov2cor(t(covM1)), 3)
         cons.y8 cons.y15 cons.y22 cons.y29 cons.y36
cons.y8    1.000       NA       NA       NA       NA
cons.y15   0.917    1.000       NA       NA       NA
cons.y22   0.824    0.922    1.000       NA       NA
cons.y29   0.725    0.843    0.948     1.00       NA
cons.y36   0.614    0.700    0.858     0.92        1
> 
> # 19.6 Fitting an autoregressive structure to the variance matrix . . . .298
> 
> (mymodel <- runMLwiN(c(y8, y15, y22, y29, y36) ~ 1 + (1 | rat), D = "Multivariate Normal", estoptions = list(EstM = 1, 
+   mcmcMeth = list(iterations = 50000), mcmcOptions = list(mcco = 4)), data = rats))
MLwiN is running, please wait......
/nogui option ignored
ECHO    0


Echoing is ON
BATC 1

Batch mode is ON
MAXI 2
STAR
iteration 0
iteration 1

Convergence achieved
TOLE 2
MAXI 20
BATC 1

Batch mode is ON
NEXT
iteration 2

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 4
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 1016.624698
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 1008.605908
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 1009.070767
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 1018.927028
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 1007.533016
Adapting for 600 iterations (Maximum 5000)
-2 * Loglike = 1020.805994
Adapting for 700 iterations (Maximum 5000)
-2 * Loglike = 1011.440789
Adapting for 800 iterations (Maximum 5000)
-2 * Loglike = 1011.484499
Adapting for 900 iterations (Maximum 5000)
-2 * Loglike = 1011.513250
Adapting for 1000 iterations (Maximum 5000)
-2 * Loglike = 1015.574310
Adapting for 1100 iterations (Maximum 5000)
-2 * Loglike = 1017.289232
Adapting for 1200 iterations (Maximum 5000)
-2 * Loglike = 1020.070828
Adapting for 1300 iterations (Maximum 5000)
-2 * Loglike = 1009.814970
Adapting for 1400 iterations (Maximum 5000)
-2 * Loglike = 1024.969594
Adapting finished and took 1500 iterations
Adapting took 1500 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 1010.566090
Burning in for 100 iterations out of 500
-2 * Loglike = 1008.441573
Burning in for 150 iterations out of 500
-2 * Loglike = 1008.100849
Burning in for 200 iterations out of 500
-2 * Loglike = 1020.348375
Burning in for 250 iterations out of 500
-2 * Loglike = 1020.824615
Burning in for 300 iterations out of 500
-2 * Loglike = 1014.750780
Burning in for 350 iterations out of 500
-2 * Loglike = 1029.249535
Burning in for 400 iterations out of 500
-2 * Loglike = 1017.027205
Burning in for 450 iterations out of 500
-2 * Loglike = 1017.263181
Burning in for 500 iterations out of 500
-2 * Loglike = 1020.356006

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 50000 1 c1090 c1091 c1003 c1004 1 4
Actual update 50 of 50000, Stored update 50 of 50000
-2 * Loglike = 1016.441678
Actual update 100 of 50000, Stored update 100 of 50000
-2 * Loglike = 1010.202563
Actual update 150 of 50000, Stored update 150 of 50000
-2 * Loglike = 1011.831392
Actual update 200 of 50000, Stored update 200 of 50000
-2 * Loglike = 1008.250494
Actual update 250 of 50000, Stored update 250 of 50000
-2 * Loglike = 1006.935569
Actual update 300 of 50000, Stored update 300 of 50000
-2 * Loglike = 1017.039697
Actual update 350 of 50000, Stored update 350 of 50000
-2 * Loglike = 1021.248399
Actual update 400 of 50000, Stored update 400 of 50000
-2 * Loglike = 1008.783608
Actual update 450 of 50000, Stored update 450 of 50000
-2 * Loglike = 1021.509422
Actual update 500 of 50000, Stored update 500 of 50000
-2 * Loglike = 1017.593841
Actual update 550 of 50000, Stored update 550 of 50000
-2 * Loglike = 1006.595263
Actual update 600 of 50000, Stored update 600 of 50000
-2 * Loglike = 1012.804379
Actual update 650 of 50000, Stored update 650 of 50000
-2 * Loglike = 1021.202253
Actual update 700 of 50000, Stored update 700 of 50000
-2 * Loglike = 1012.028717
Actual update 750 of 50000, Stored update 750 of 50000
-2 * Loglike = 1022.903985
Actual update 800 of 50000, Stored update 800 of 50000
-2 * Loglike = 1022.289997
Actual update 850 of 50000, Stored update 850 of 50000
-2 * Loglike = 1014.408990
Actual update 900 of 50000, Stored update 900 of 50000
-2 * Loglike = 1007.695558
Actual update 950 of 50000, Stored update 950 of 50000
-2 * Loglike = 1015.918717
Actual update 1000 of 50000, Stored update 1000 of 50000
-2 * Loglike = 1012.241524
Actual update 1050 of 50000, Stored update 1050 of 50000
-2 * Loglike = 1005.727757
Actual update 1100 of 50000, Stored update 1100 of 50000
-2 * Loglike = 1006.808714
Actual update 1150 of 50000, Stored update 1150 of 50000
-2 * Loglike = 1012.472998
Actual update 1200 of 50000, Stored update 1200 of 50000
-2 * Loglike = 1007.903318
Actual update 1250 of 50000, Stored update 1250 of 50000
-2 * Loglike = 1008.337983
Actual update 1300 of 50000, Stored update 1300 of 50000
-2 * Loglike = 1009.058194
Actual update 1350 of 50000, Stored update 1350 of 50000
-2 * Loglike = 1018.889231
Actual update 1400 of 50000, Stored update 1400 of 50000
-2 * Loglike = 1006.795302
Actual update 1450 of 50000, Stored update 1450 of 50000
-2 * Loglike = 1008.457421
Actual update 1500 of 50000, Stored update 1500 of 50000
-2 * Loglike = 1007.494342
Actual update 1550 of 50000, Stored update 1550 of 50000
-2 * Loglike = 1020.627549
Actual update 1600 of 50000, Stored update 1600 of 50000
-2 * Loglike = 1014.934025
Actual update 1650 of 50000, Stored update 1650 of 50000
-2 * Loglike = 1016.502796
Actual update 1700 of 50000, Stored update 1700 of 50000
-2 * Loglike = 1020.007771
Actual update 1750 of 50000, Stored update 1750 of 50000
-2 * Loglike = 1017.572431
Actual update 1800 of 50000, Stored update 1800 of 50000
-2 * Loglike = 1015.713737
Actual update 1850 of 50000, Stored update 1850 of 50000
-2 * Loglike = 1010.557916
Actual update 1900 of 50000, Stored update 1900 of 50000
-2 * Loglike = 1014.360348
Actual update 1950 of 50000, Stored update 1950 of 50000
-2 * Loglike = 1005.112777
Actual update 2000 of 50000, Stored update 2000 of 50000
-2 * Loglike = 1020.030861
Actual update 2050 of 50000, Stored update 2050 of 50000
-2 * Loglike = 1014.048123
Actual update 2100 of 50000, Stored update 2100 of 50000
-2 * Loglike = 1020.780800
Actual update 2150 of 50000, Stored update 2150 of 50000
-2 * Loglike = 1016.623249
Actual update 2200 of 50000, Stored update 2200 of 50000
-2 * Loglike = 1015.047747
Actual update 2250 of 50000, Stored update 2250 of 50000
-2 * Loglike = 1014.432035
Actual update 2300 of 50000, Stored update 2300 of 50000
-2 * Loglike = 1020.842533
Actual update 2350 of 50000, Stored update 2350 of 50000
-2 * Loglike = 1008.342615
Actual update 2400 of 50000, Stored update 2400 of 50000
-2 * Loglike = 1019.480506
Actual update 2450 of 50000, Stored update 2450 of 50000
-2 * Loglike = 1023.712967
Actual update 2500 of 50000, Stored update 2500 of 50000
-2 * Loglike = 1029.615857
Actual update 2550 of 50000, Stored update 2550 of 50000
-2 * Loglike = 1015.573914
Actual update 2600 of 50000, Stored update 2600 of 50000
-2 * Loglike = 1010.574522
Actual update 2650 of 50000, Stored update 2650 of 50000
-2 * Loglike = 1024.022687
Actual update 2700 of 50000, Stored update 2700 of 50000
-2 * Loglike = 1016.806892
Actual update 2750 of 50000, Stored update 2750 of 50000
-2 * Loglike = 1017.268107
Actual update 2800 of 50000, Stored update 2800 of 50000
-2 * Loglike = 1017.498462
Actual update 2850 of 50000, Stored update 2850 of 50000
-2 * Loglike = 1020.170940
Actual update 2900 of 50000, Stored update 2900 of 50000
-2 * Loglike = 1017.742886
Actual update 2950 of 50000, Stored update 2950 of 50000
-2 * Loglike = 1022.324671
Actual update 3000 of 50000, Stored update 3000 of 50000
-2 * Loglike = 1016.436711
Actual update 3050 of 50000, Stored update 3050 of 50000
-2 * Loglike = 1015.752862
Actual update 3100 of 50000, Stored update 3100 of 50000
-2 * Loglike = 1022.197505
Actual update 3150 of 50000, Stored update 3150 of 50000
-2 * Loglike = 1005.969966
Actual update 3200 of 50000, Stored update 3200 of 50000
-2 * Loglike = 1020.390297
Actual update 3250 of 50000, Stored update 3250 of 50000
-2 * Loglike = 1013.271658
Actual update 3300 of 50000, Stored update 3300 of 50000
-2 * Loglike = 1030.192388
Actual update 3350 of 50000, Stored update 3350 of 50000
-2 * Loglike = 1007.521462
Actual update 3400 of 50000, Stored update 3400 of 50000
-2 * Loglike = 1008.637029
Actual update 3450 of 50000, Stored update 3450 of 50000
-2 * Loglike = 1016.729922
Actual update 3500 of 50000, Stored update 3500 of 50000
-2 * Loglike = 1015.986462
Actual update 3550 of 50000, Stored update 3550 of 50000
-2 * Loglike = 1013.293208
Actual update 3600 of 50000, Stored update 3600 of 50000
-2 * Loglike = 1009.965805
Actual update 3650 of 50000, Stored update 3650 of 50000
-2 * Loglike = 1015.833158
Actual update 3700 of 50000, Stored update 3700 of 50000
-2 * Loglike = 1023.063515
Actual update 3750 of 50000, Stored update 3750 of 50000
-2 * Loglike = 1013.586550
Actual update 3800 of 50000, Stored update 3800 of 50000
-2 * Loglike = 1005.514443
Actual update 3850 of 50000, Stored update 3850 of 50000
-2 * Loglike = 1010.019698
Actual update 3900 of 50000, Stored update 3900 of 50000
-2 * Loglike = 1013.028663
Actual update 3950 of 50000, Stored update 3950 of 50000
-2 * Loglike = 1015.275639
Actual update 4000 of 50000, Stored update 4000 of 50000
-2 * Loglike = 1016.561760
Actual update 4050 of 50000, Stored update 4050 of 50000
-2 * Loglike = 1018.456227
Actual update 4100 of 50000, Stored update 4100 of 50000
-2 * Loglike = 1014.824768
Actual update 4150 of 50000, Stored update 4150 of 50000
-2 * Loglike = 1009.189601
Actual update 4200 of 50000, Stored update 4200 of 50000
-2 * Loglike = 1022.692343
Actual update 4250 of 50000, Stored update 4250 of 50000
-2 * Loglike = 1009.042323
Actual update 4300 of 50000, Stored update 4300 of 50000
-2 * Loglike = 1011.655352
Actual update 4350 of 50000, Stored update 4350 of 50000
-2 * Loglike = 1008.933639
Actual update 4400 of 50000, Stored update 4400 of 50000
-2 * Loglike = 1015.695589
Actual update 4450 of 50000, Stored update 4450 of 50000
-2 * Loglike = 1013.718937
Actual update 4500 of 50000, Stored update 4500 of 50000
-2 * Loglike = 1018.093306
Actual update 4550 of 50000, Stored update 4550 of 50000
-2 * Loglike = 1020.682127
Actual update 4600 of 50000, Stored update 4600 of 50000
-2 * Loglike = 1024.431097
Actual update 4650 of 50000, Stored update 4650 of 50000
-2 * Loglike = 1011.192936
Actual update 4700 of 50000, Stored update 4700 of 50000
-2 * Loglike = 1019.419434
Actual update 4750 of 50000, Stored update 4750 of 50000
-2 * Loglike = 1013.498805
Actual update 4800 of 50000, Stored update 4800 of 50000
-2 * Loglike = 1016.984248
Actual update 4850 of 50000, Stored update 4850 of 50000
-2 * Loglike = 1021.163006
Actual update 4900 of 50000, Stored update 4900 of 50000
-2 * Loglike = 1025.831540
Actual update 4950 of 50000, Stored update 4950 of 50000
-2 * Loglike = 1028.636337
Actual update 5000 of 50000, Stored update 5000 of 50000
-2 * Loglike = 1047.021764
Actual update 5050 of 50000, Stored update 5050 of 50000
-2 * Loglike = 1020.401698
Actual update 5100 of 50000, Stored update 5100 of 50000
-2 * Loglike = 1017.551938
Actual update 5150 of 50000, Stored update 5150 of 50000
-2 * Loglike = 1023.805317
Actual update 5200 of 50000, Stored update 5200 of 50000
-2 * Loglike = 1020.707470
Actual update 5250 of 50000, Stored update 5250 of 50000
-2 * Loglike = 1020.512011
Actual update 5300 of 50000, Stored update 5300 of 50000
-2 * Loglike = 1016.565454
Actual update 5350 of 50000, Stored update 5350 of 50000
-2 * Loglike = 1014.025173
Actual update 5400 of 50000, Stored update 5400 of 50000
-2 * Loglike = 1017.817837
Actual update 5450 of 50000, Stored update 5450 of 50000
-2 * Loglike = 1014.595849
Actual update 5500 of 50000, Stored update 5500 of 50000
-2 * Loglike = 1030.471663
Actual update 5550 of 50000, Stored update 5550 of 50000
-2 * Loglike = 1022.483837
Actual update 5600 of 50000, Stored update 5600 of 50000
-2 * Loglike = 1018.042110
Actual update 5650 of 50000, Stored update 5650 of 50000
-2 * Loglike = 1012.885574
Actual update 5700 of 50000, Stored update 5700 of 50000
-2 * Loglike = 1016.189808
Actual update 5750 of 50000, Stored update 5750 of 50000
-2 * Loglike = 1011.106721
Actual update 5800 of 50000, Stored update 5800 of 50000
-2 * Loglike = 1012.598813
Actual update 5850 of 50000, Stored update 5850 of 50000
-2 * Loglike = 1022.565045
Actual update 5900 of 50000, Stored update 5900 of 50000
-2 * Loglike = 1004.418480
Actual update 5950 of 50000, Stored update 5950 of 50000
-2 * Loglike = 1016.348989
Actual update 6000 of 50000, Stored update 6000 of 50000
-2 * Loglike = 1013.595084
Actual update 6050 of 50000, Stored update 6050 of 50000
-2 * Loglike = 1011.264494
Actual update 6100 of 50000, Stored update 6100 of 50000
-2 * Loglike = 1015.217719
Actual update 6150 of 50000, Stored update 6150 of 50000
-2 * Loglike = 1010.845722
Actual update 6200 of 50000, Stored update 6200 of 50000
-2 * Loglike = 1021.520700
Actual update 6250 of 50000, Stored update 6250 of 50000
-2 * Loglike = 1016.725922
Actual update 6300 of 50000, Stored update 6300 of 50000
-2 * Loglike = 1016.120195
Actual update 6350 of 50000, Stored update 6350 of 50000
-2 * Loglike = 1014.078924
Actual update 6400 of 50000, Stored update 6400 of 50000
-2 * Loglike = 1018.960654
Actual update 6450 of 50000, Stored update 6450 of 50000
-2 * Loglike = 1008.901435
Actual update 6500 of 50000, Stored update 6500 of 50000
-2 * Loglike = 1011.444744
Actual update 6550 of 50000, Stored update 6550 of 50000
-2 * Loglike = 1013.420129
Actual update 6600 of 50000, Stored update 6600 of 50000
-2 * Loglike = 1007.864058
Actual update 6650 of 50000, Stored update 6650 of 50000
-2 * Loglike = 1018.205973
Actual update 6700 of 50000, Stored update 6700 of 50000
-2 * Loglike = 1013.274486
Actual update 6750 of 50000, Stored update 6750 of 50000
-2 * Loglike = 1016.110789
Actual update 6800 of 50000, Stored update 6800 of 50000
-2 * Loglike = 1017.608655
Actual update 6850 of 50000, Stored update 6850 of 50000
-2 * Loglike = 1013.546396
Actual update 6900 of 50000, Stored update 6900 of 50000
-2 * Loglike = 1025.164005
Actual update 6950 of 50000, Stored update 6950 of 50000
-2 * Loglike = 1013.537496
Actual update 7000 of 50000, Stored update 7000 of 50000
-2 * Loglike = 1016.381952
Actual update 7050 of 50000, Stored update 7050 of 50000
-2 * Loglike = 1017.875625
Actual update 7100 of 50000, Stored update 7100 of 50000
-2 * Loglike = 1008.188043
Actual update 7150 of 50000, Stored update 7150 of 50000
-2 * Loglike = 1019.416455
Actual update 7200 of 50000, Stored update 7200 of 50000
-2 * Loglike = 1012.508957
Actual update 7250 of 50000, Stored update 7250 of 50000
-2 * Loglike = 1008.972442
Actual update 7300 of 50000, Stored update 7300 of 50000
-2 * Loglike = 1014.152011
Actual update 7350 of 50000, Stored update 7350 of 50000
-2 * Loglike = 1013.915934
Actual update 7400 of 50000, Stored update 7400 of 50000
-2 * Loglike = 1006.275735
Actual update 7450 of 50000, Stored update 7450 of 50000
-2 * Loglike = 1015.176184
Actual update 7500 of 50000, Stored update 7500 of 50000
-2 * Loglike = 1013.728997
Actual update 7550 of 50000, Stored update 7550 of 50000
-2 * Loglike = 1020.678858
Actual update 7600 of 50000, Stored update 7600 of 50000
-2 * Loglike = 1012.035842
Actual update 7650 of 50000, Stored update 7650 of 50000
-2 * Loglike = 1006.766038
Actual update 7700 of 50000, Stored update 7700 of 50000
-2 * Loglike = 1016.564075
Actual update 7750 of 50000, Stored update 7750 of 50000
-2 * Loglike = 1014.050578
Actual update 7800 of 50000, Stored update 7800 of 50000
-2 * Loglike = 1013.929140
Actual update 7850 of 50000, Stored update 7850 of 50000
-2 * Loglike = 1008.341894
Actual update 7900 of 50000, Stored update 7900 of 50000
-2 * Loglike = 1023.470895
Actual update 7950 of 50000, Stored update 7950 of 50000
-2 * Loglike = 1012.510404
Actual update 8000 of 50000, Stored update 8000 of 50000
-2 * Loglike = 1017.547756
Actual update 8050 of 50000, Stored update 8050 of 50000
-2 * Loglike = 1017.961700
Actual update 8100 of 50000, Stored update 8100 of 50000
-2 * Loglike = 1018.208763
Actual update 8150 of 50000, Stored update 8150 of 50000
-2 * Loglike = 1023.720373
Actual update 8200 of 50000, Stored update 8200 of 50000
-2 * Loglike = 1022.182625
Actual update 8250 of 50000, Stored update 8250 of 50000
-2 * Loglike = 1026.144048
Actual update 8300 of 50000, Stored update 8300 of 50000
-2 * Loglike = 1027.865113
Actual update 8350 of 50000, Stored update 8350 of 50000
-2 * Loglike = 1025.751672
Actual update 8400 of 50000, Stored update 8400 of 50000
-2 * Loglike = 1021.474865
Actual update 8450 of 50000, Stored update 8450 of 50000
-2 * Loglike = 1013.114197
Actual update 8500 of 50000, Stored update 8500 of 50000
-2 * Loglike = 1013.144770
Actual update 8550 of 50000, Stored update 8550 of 50000
-2 * Loglike = 1023.445052
Actual update 8600 of 50000, Stored update 8600 of 50000
-2 * Loglike = 1034.558138
Actual update 8650 of 50000, Stored update 8650 of 50000
-2 * Loglike = 1022.642475
Actual update 8700 of 50000, Stored update 8700 of 50000
-2 * Loglike = 1010.628490
Actual update 8750 of 50000, Stored update 8750 of 50000
-2 * Loglike = 1012.931440
Actual update 8800 of 50000, Stored update 8800 of 50000
-2 * Loglike = 1010.421994
Actual update 8850 of 50000, Stored update 8850 of 50000
-2 * Loglike = 1020.767053
Actual update 8900 of 50000, Stored update 8900 of 50000
-2 * Loglike = 1012.968031
Actual update 8950 of 50000, Stored update 8950 of 50000
-2 * Loglike = 1012.515236
Actual update 9000 of 50000, Stored update 9000 of 50000
-2 * Loglike = 1016.661845
Actual update 9050 of 50000, Stored update 9050 of 50000
-2 * Loglike = 1022.831293
Actual update 9100 of 50000, Stored update 9100 of 50000
-2 * Loglike = 1012.099324
Actual update 9150 of 50000, Stored update 9150 of 50000
-2 * Loglike = 1011.757689
Actual update 9200 of 50000, Stored update 9200 of 50000
-2 * Loglike = 1016.779354
Actual update 9250 of 50000, Stored update 9250 of 50000
-2 * Loglike = 1006.889424
Actual update 9300 of 50000, Stored update 9300 of 50000
-2 * Loglike = 1006.889134
Actual update 9350 of 50000, Stored update 9350 of 50000
-2 * Loglike = 1013.323594
Actual update 9400 of 50000, Stored update 9400 of 50000
-2 * Loglike = 1011.090975
Actual update 9450 of 50000, Stored update 9450 of 50000
-2 * Loglike = 1014.260062
Actual update 9500 of 50000, Stored update 9500 of 50000
-2 * Loglike = 1013.368589
Actual update 9550 of 50000, Stored update 9550 of 50000
-2 * Loglike = 1007.431861
Actual update 9600 of 50000, Stored update 9600 of 50000
-2 * Loglike = 1020.973556
Actual update 9650 of 50000, Stored update 9650 of 50000
-2 * Loglike = 1008.025566
Actual update 9700 of 50000, Stored update 9700 of 50000
-2 * Loglike = 1010.059710
Actual update 9750 of 50000, Stored update 9750 of 50000
-2 * Loglike = 1018.419980
Actual update 9800 of 50000, Stored update 9800 of 50000
-2 * Loglike = 1018.019092
Actual update 9850 of 50000, Stored update 9850 of 50000
-2 * Loglike = 1008.379689
Actual update 9900 of 50000, Stored update 9900 of 50000
-2 * Loglike = 1011.359852
Actual update 9950 of 50000, Stored update 9950 of 50000
-2 * Loglike = 1021.091746
Actual update 10000 of 50000, Stored update 10000 of 50000
-2 * Loglike = 1021.518916
Actual update 10050 of 50000, Stored update 10050 of 50000
-2 * Loglike = 1016.078643
Actual update 10100 of 50000, Stored update 10100 of 50000
-2 * Loglike = 1012.363887
Actual update 10150 of 50000, Stored update 10150 of 50000
-2 * Loglike = 1017.441452
Actual update 10200 of 50000, Stored update 10200 of 50000
-2 * Loglike = 1030.613342
Actual update 10250 of 50000, Stored update 10250 of 50000
-2 * Loglike = 1012.345129
Actual update 10300 of 50000, Stored update 10300 of 50000
-2 * Loglike = 1018.585096
Actual update 10350 of 50000, Stored update 10350 of 50000
-2 * Loglike = 1015.761809
Actual update 10400 of 50000, Stored update 10400 of 50000
-2 * Loglike = 1025.300878
Actual update 10450 of 50000, Stored update 10450 of 50000
-2 * Loglike = 1014.265210
Actual update 10500 of 50000, Stored update 10500 of 50000
-2 * Loglike = 1009.709390
Actual update 10550 of 50000, Stored update 10550 of 50000
-2 * Loglike = 1006.547863
Actual update 10600 of 50000, Stored update 10600 of 50000
-2 * Loglike = 1015.950591
Actual update 10650 of 50000, Stored update 10650 of 50000
-2 * Loglike = 1012.360603
Actual update 10700 of 50000, Stored update 10700 of 50000
-2 * Loglike = 1028.352130
Actual update 10750 of 50000, Stored update 10750 of 50000
-2 * Loglike = 1014.240675
Actual update 10800 of 50000, Stored update 10800 of 50000
-2 * Loglike = 1017.647386
Actual update 10850 of 50000, Stored update 10850 of 50000
-2 * Loglike = 1012.507521
Actual update 10900 of 50000, Stored update 10900 of 50000
-2 * Loglike = 1017.436609
Actual update 10950 of 50000, Stored update 10950 of 50000
-2 * Loglike = 1020.401706
Actual update 11000 of 50000, Stored update 11000 of 50000
-2 * Loglike = 1017.892344
Actual update 11050 of 50000, Stored update 11050 of 50000
-2 * Loglike = 1010.564051
Actual update 11100 of 50000, Stored update 11100 of 50000
-2 * Loglike = 1013.367025
Actual update 11150 of 50000, Stored update 11150 of 50000
-2 * Loglike = 1021.587847
Actual update 11200 of 50000, Stored update 11200 of 50000
-2 * Loglike = 1018.930712
Actual update 11250 of 50000, Stored update 11250 of 50000
-2 * Loglike = 1031.074578
Actual update 11300 of 50000, Stored update 11300 of 50000
-2 * Loglike = 1015.678243
Actual update 11350 of 50000, Stored update 11350 of 50000
-2 * Loglike = 1011.960416
Actual update 11400 of 50000, Stored update 11400 of 50000
-2 * Loglike = 1017.905070
Actual update 11450 of 50000, Stored update 11450 of 50000
-2 * Loglike = 1011.313374
Actual update 11500 of 50000, Stored update 11500 of 50000
-2 * Loglike = 1017.383281
Actual update 11550 of 50000, Stored update 11550 of 50000
-2 * Loglike = 1016.598992
Actual update 11600 of 50000, Stored update 11600 of 50000
-2 * Loglike = 1025.809351
Actual update 11650 of 50000, Stored update 11650 of 50000
-2 * Loglike = 1004.367094
Actual update 11700 of 50000, Stored update 11700 of 50000
-2 * Loglike = 1010.269951
Actual update 11750 of 50000, Stored update 11750 of 50000
-2 * Loglike = 1014.574390
Actual update 11800 of 50000, Stored update 11800 of 50000
-2 * Loglike = 1016.814490
Actual update 11850 of 50000, Stored update 11850 of 50000
-2 * Loglike = 1029.963367
Actual update 11900 of 50000, Stored update 11900 of 50000
-2 * Loglike = 1022.895894
Actual update 11950 of 50000, Stored update 11950 of 50000
-2 * Loglike = 1031.103133
Actual update 12000 of 50000, Stored update 12000 of 50000
-2 * Loglike = 1013.210749
Actual update 12050 of 50000, Stored update 12050 of 50000
-2 * Loglike = 1024.602991
Actual update 12100 of 50000, Stored update 12100 of 50000
-2 * Loglike = 1014.755734
Actual update 12150 of 50000, Stored update 12150 of 50000
-2 * Loglike = 1020.777593
Actual update 12200 of 50000, Stored update 12200 of 50000
-2 * Loglike = 1014.768563
Actual update 12250 of 50000, Stored update 12250 of 50000
-2 * Loglike = 1025.956509
Actual update 12300 of 50000, Stored update 12300 of 50000
-2 * Loglike = 1020.813694
Actual update 12350 of 50000, Stored update 12350 of 50000
-2 * Loglike = 1017.562069
Actual update 12400 of 50000, Stored update 12400 of 50000
-2 * Loglike = 1019.747766
Actual update 12450 of 50000, Stored update 12450 of 50000
-2 * Loglike = 1022.155495
Actual update 12500 of 50000, Stored update 12500 of 50000
-2 * Loglike = 1029.228346
Actual update 12550 of 50000, Stored update 12550 of 50000
-2 * Loglike = 1037.810849
Actual update 12600 of 50000, Stored update 12600 of 50000
-2 * Loglike = 1025.184707
Actual update 12650 of 50000, Stored update 12650 of 50000
-2 * Loglike = 1016.696596
Actual update 12700 of 50000, Stored update 12700 of 50000
-2 * Loglike = 1015.431775
Actual update 12750 of 50000, Stored update 12750 of 50000
-2 * Loglike = 1020.044487
Actual update 12800 of 50000, Stored update 12800 of 50000
-2 * Loglike = 1012.377032
Actual update 12850 of 50000, Stored update 12850 of 50000
-2 * Loglike = 1010.798001
Actual update 12900 of 50000, Stored update 12900 of 50000
-2 * Loglike = 1029.641312
Actual update 12950 of 50000, Stored update 12950 of 50000
-2 * Loglike = 1023.494558
Actual update 13000 of 50000, Stored update 13000 of 50000
-2 * Loglike = 1019.669746
Actual update 13050 of 50000, Stored update 13050 of 50000
-2 * Loglike = 1012.291378
Actual update 13100 of 50000, Stored update 13100 of 50000
-2 * Loglike = 1012.453680
Actual update 13150 of 50000, Stored update 13150 of 50000
-2 * Loglike = 1018.641583
Actual update 13200 of 50000, Stored update 13200 of 50000
-2 * Loglike = 1007.142214
Actual update 13250 of 50000, Stored update 13250 of 50000
-2 * Loglike = 1011.633253
Actual update 13300 of 50000, Stored update 13300 of 50000
-2 * Loglike = 1020.050497
Actual update 13350 of 50000, Stored update 13350 of 50000
-2 * Loglike = 1014.234335
Actual update 13400 of 50000, Stored update 13400 of 50000
-2 * Loglike = 1007.830551
Actual update 13450 of 50000, Stored update 13450 of 50000
-2 * Loglike = 1027.487521
Actual update 13500 of 50000, Stored update 13500 of 50000
-2 * Loglike = 1019.412450
Actual update 13550 of 50000, Stored update 13550 of 50000
-2 * Loglike = 1018.049675
Actual update 13600 of 50000, Stored update 13600 of 50000
-2 * Loglike = 1021.036879
Actual update 13650 of 50000, Stored update 13650 of 50000
-2 * Loglike = 1023.330570
Actual update 13700 of 50000, Stored update 13700 of 50000
-2 * Loglike = 1009.096437
Actual update 13750 of 50000, Stored update 13750 of 50000
-2 * Loglike = 1016.202675
Actual update 13800 of 50000, Stored update 13800 of 50000
-2 * Loglike = 1019.543385
Actual update 13850 of 50000, Stored update 13850 of 50000
-2 * Loglike = 1010.583577
Actual update 13900 of 50000, Stored update 13900 of 50000
-2 * Loglike = 1019.759116
Actual update 13950 of 50000, Stored update 13950 of 50000
-2 * Loglike = 1010.772336
Actual update 14000 of 50000, Stored update 14000 of 50000
-2 * Loglike = 1013.678018
Actual update 14050 of 50000, Stored update 14050 of 50000
-2 * Loglike = 1011.119690
Actual update 14100 of 50000, Stored update 14100 of 50000
-2 * Loglike = 1012.057994
Actual update 14150 of 50000, Stored update 14150 of 50000
-2 * Loglike = 1011.637830
Actual update 14200 of 50000, Stored update 14200 of 50000
-2 * Loglike = 1013.807278
Actual update 14250 of 50000, Stored update 14250 of 50000
-2 * Loglike = 1017.441221
Actual update 14300 of 50000, Stored update 14300 of 50000
-2 * Loglike = 1013.613454
Actual update 14350 of 50000, Stored update 14350 of 50000
-2 * Loglike = 1014.288373
Actual update 14400 of 50000, Stored update 14400 of 50000
-2 * Loglike = 1014.150426
Actual update 14450 of 50000, Stored update 14450 of 50000
-2 * Loglike = 1015.373787
Actual update 14500 of 50000, Stored update 14500 of 50000
-2 * Loglike = 1015.458492
Actual update 14550 of 50000, Stored update 14550 of 50000
-2 * Loglike = 1019.723960
Actual update 14600 of 50000, Stored update 14600 of 50000
-2 * Loglike = 1019.057411
Actual update 14650 of 50000, Stored update 14650 of 50000
-2 * Loglike = 1016.766330
Actual update 14700 of 50000, Stored update 14700 of 50000
-2 * Loglike = 1013.115791
Actual update 14750 of 50000, Stored update 14750 of 50000
-2 * Loglike = 1017.357341
Actual update 14800 of 50000, Stored update 14800 of 50000
-2 * Loglike = 1019.144170
Actual update 14850 of 50000, Stored update 14850 of 50000
-2 * Loglike = 1015.512496
Actual update 14900 of 50000, Stored update 14900 of 50000
-2 * Loglike = 1017.338403
Actual update 14950 of 50000, Stored update 14950 of 50000
-2 * Loglike = 1010.210537
Actual update 15000 of 50000, Stored update 15000 of 50000
-2 * Loglike = 1011.540402
Actual update 15050 of 50000, Stored update 15050 of 50000
-2 * Loglike = 1012.599622
Actual update 15100 of 50000, Stored update 15100 of 50000
-2 * Loglike = 1012.264768
Actual update 15150 of 50000, Stored update 15150 of 50000
-2 * Loglike = 1008.869309
Actual update 15200 of 50000, Stored update 15200 of 50000
-2 * Loglike = 1017.630672
Actual update 15250 of 50000, Stored update 15250 of 50000
-2 * Loglike = 1012.223249
Actual update 15300 of 50000, Stored update 15300 of 50000
-2 * Loglike = 1009.892588
Actual update 15350 of 50000, Stored update 15350 of 50000
-2 * Loglike = 1016.227911
Actual update 15400 of 50000, Stored update 15400 of 50000
-2 * Loglike = 1027.478128
Actual update 15450 of 50000, Stored update 15450 of 50000
-2 * Loglike = 1007.285343
Actual update 15500 of 50000, Stored update 15500 of 50000
-2 * Loglike = 1015.928894
Actual update 15550 of 50000, Stored update 15550 of 50000
-2 * Loglike = 1005.200543
Actual update 15600 of 50000, Stored update 15600 of 50000
-2 * Loglike = 1009.954585
Actual update 15650 of 50000, Stored update 15650 of 50000
-2 * Loglike = 1015.249704
Actual update 15700 of 50000, Stored update 15700 of 50000
-2 * Loglike = 1010.882554
Actual update 15750 of 50000, Stored update 15750 of 50000
-2 * Loglike = 1014.196848
Actual update 15800 of 50000, Stored update 15800 of 50000
-2 * Loglike = 1010.499864
Actual update 15850 of 50000, Stored update 15850 of 50000
-2 * Loglike = 1019.649682
Actual update 15900 of 50000, Stored update 15900 of 50000
-2 * Loglike = 1011.057320
Actual update 15950 of 50000, Stored update 15950 of 50000
-2 * Loglike = 1021.499947
Actual update 16000 of 50000, Stored update 16000 of 50000
-2 * Loglike = 1011.276774
Actual update 16050 of 50000, Stored update 16050 of 50000
-2 * Loglike = 1005.445654
Actual update 16100 of 50000, Stored update 16100 of 50000
-2 * Loglike = 1010.780604
Actual update 16150 of 50000, Stored update 16150 of 50000
-2 * Loglike = 1015.842631
Actual update 16200 of 50000, Stored update 16200 of 50000
-2 * Loglike = 1016.736483
Actual update 16250 of 50000, Stored update 16250 of 50000
-2 * Loglike = 1014.200249
Actual update 16300 of 50000, Stored update 16300 of 50000
-2 * Loglike = 1017.642637
Actual update 16350 of 50000, Stored update 16350 of 50000
-2 * Loglike = 1010.785966
Actual update 16400 of 50000, Stored update 16400 of 50000
-2 * Loglike = 1006.556109
Actual update 16450 of 50000, Stored update 16450 of 50000
-2 * Loglike = 1012.408513
Actual update 16500 of 50000, Stored update 16500 of 50000
-2 * Loglike = 1009.497224
Actual update 16550 of 50000, Stored update 16550 of 50000
-2 * Loglike = 1019.820925
Actual update 16600 of 50000, Stored update 16600 of 50000
-2 * Loglike = 1012.318203
Actual update 16650 of 50000, Stored update 16650 of 50000
-2 * Loglike = 1014.889575
Actual update 16700 of 50000, Stored update 16700 of 50000
-2 * Loglike = 1015.471913
Actual update 16750 of 50000, Stored update 16750 of 50000
-2 * Loglike = 1006.010054
Actual update 16800 of 50000, Stored update 16800 of 50000
-2 * Loglike = 1012.154717
Actual update 16850 of 50000, Stored update 16850 of 50000
-2 * Loglike = 1016.868583
Actual update 16900 of 50000, Stored update 16900 of 50000
-2 * Loglike = 1014.738077
Actual update 16950 of 50000, Stored update 16950 of 50000
-2 * Loglike = 1016.035198
Actual update 17000 of 50000, Stored update 17000 of 50000
-2 * Loglike = 1011.179329
Actual update 17050 of 50000, Stored update 17050 of 50000
-2 * Loglike = 1025.645854
Actual update 17100 of 50000, Stored update 17100 of 50000
-2 * Loglike = 1007.695052
Actual update 17150 of 50000, Stored update 17150 of 50000
-2 * Loglike = 1007.495788
Actual update 17200 of 50000, Stored update 17200 of 50000
-2 * Loglike = 1014.123291
Actual update 17250 of 50000, Stored update 17250 of 50000
-2 * Loglike = 1010.850138
Actual update 17300 of 50000, Stored update 17300 of 50000
-2 * Loglike = 1014.257317
Actual update 17350 of 50000, Stored update 17350 of 50000
-2 * Loglike = 1010.357786
Actual update 17400 of 50000, Stored update 17400 of 50000
-2 * Loglike = 1010.343306
Actual update 17450 of 50000, Stored update 17450 of 50000
-2 * Loglike = 1018.387430
Actual update 17500 of 50000, Stored update 17500 of 50000
-2 * Loglike = 1004.950759
Actual update 17550 of 50000, Stored update 17550 of 50000
-2 * Loglike = 1023.750674
Actual update 17600 of 50000, Stored update 17600 of 50000
-2 * Loglike = 1021.147627
Actual update 17650 of 50000, Stored update 17650 of 50000
-2 * Loglike = 1012.289048
Actual update 17700 of 50000, Stored update 17700 of 50000
-2 * Loglike = 1020.990506
Actual update 17750 of 50000, Stored update 17750 of 50000
-2 * Loglike = 1015.888403
Actual update 17800 of 50000, Stored update 17800 of 50000
-2 * Loglike = 1019.785225
Actual update 17850 of 50000, Stored update 17850 of 50000
-2 * Loglike = 1025.669394
Actual update 17900 of 50000, Stored update 17900 of 50000
-2 * Loglike = 1025.831908
Actual update 17950 of 50000, Stored update 17950 of 50000
-2 * Loglike = 1026.393815
Actual update 18000 of 50000, Stored update 18000 of 50000
-2 * Loglike = 1030.660983
Actual update 18050 of 50000, Stored update 18050 of 50000
-2 * Loglike = 1020.287160
Actual update 18100 of 50000, Stored update 18100 of 50000
-2 * Loglike = 1037.137904
Actual update 18150 of 50000, Stored update 18150 of 50000
-2 * Loglike = 1016.264596
Actual update 18200 of 50000, Stored update 18200 of 50000
-2 * Loglike = 1012.032685
Actual update 18250 of 50000, Stored update 18250 of 50000
-2 * Loglike = 1030.012754
Actual update 18300 of 50000, Stored update 18300 of 50000
-2 * Loglike = 1017.941574
Actual update 18350 of 50000, Stored update 18350 of 50000
-2 * Loglike = 1015.634703
Actual update 18400 of 50000, Stored update 18400 of 50000
-2 * Loglike = 1015.822441
Actual update 18450 of 50000, Stored update 18450 of 50000
-2 * Loglike = 1014.423459
Actual update 18500 of 50000, Stored update 18500 of 50000
-2 * Loglike = 1013.883768
Actual update 18550 of 50000, Stored update 18550 of 50000
-2 * Loglike = 1012.864958
Actual update 18600 of 50000, Stored update 18600 of 50000
-2 * Loglike = 1020.069001
Actual update 18650 of 50000, Stored update 18650 of 50000
-2 * Loglike = 1017.835995
Actual update 18700 of 50000, Stored update 18700 of 50000
-2 * Loglike = 1021.724054
Actual update 18750 of 50000, Stored update 18750 of 50000
-2 * Loglike = 1008.407774
Actual update 18800 of 50000, Stored update 18800 of 50000
-2 * Loglike = 1019.157667
Actual update 18850 of 50000, Stored update 18850 of 50000
-2 * Loglike = 1016.032794
Actual update 18900 of 50000, Stored update 18900 of 50000
-2 * Loglike = 1017.849092
Actual update 18950 of 50000, Stored update 18950 of 50000
-2 * Loglike = 1012.267856
Actual update 19000 of 50000, Stored update 19000 of 50000
-2 * Loglike = 1018.368819
Actual update 19050 of 50000, Stored update 19050 of 50000
-2 * Loglike = 1013.444017
Actual update 19100 of 50000, Stored update 19100 of 50000
-2 * Loglike = 1011.169727
Actual update 19150 of 50000, Stored update 19150 of 50000
-2 * Loglike = 1012.538362
Actual update 19200 of 50000, Stored update 19200 of 50000
-2 * Loglike = 1021.094641
Actual update 19250 of 50000, Stored update 19250 of 50000
-2 * Loglike = 1018.970402
Actual update 19300 of 50000, Stored update 19300 of 50000
-2 * Loglike = 1019.508837
Actual update 19350 of 50000, Stored update 19350 of 50000
-2 * Loglike = 1013.984086
Actual update 19400 of 50000, Stored update 19400 of 50000
-2 * Loglike = 1009.607082
Actual update 19450 of 50000, Stored update 19450 of 50000
-2 * Loglike = 1020.395255
Actual update 19500 of 50000, Stored update 19500 of 50000
-2 * Loglike = 1018.349668
Actual update 19550 of 50000, Stored update 19550 of 50000
-2 * Loglike = 1009.024769
Actual update 19600 of 50000, Stored update 19600 of 50000
-2 * Loglike = 1014.276184
Actual update 19650 of 50000, Stored update 19650 of 50000
-2 * Loglike = 1019.012843
Actual update 19700 of 50000, Stored update 19700 of 50000
-2 * Loglike = 1017.132933
Actual update 19750 of 50000, Stored update 19750 of 50000
-2 * Loglike = 1013.572575
Actual update 19800 of 50000, Stored update 19800 of 50000
-2 * Loglike = 1011.797495
Actual update 19850 of 50000, Stored update 19850 of 50000
-2 * Loglike = 1017.233783
Actual update 19900 of 50000, Stored update 19900 of 50000
-2 * Loglike = 1014.327298
Actual update 19950 of 50000, Stored update 19950 of 50000
-2 * Loglike = 1015.542208
Actual update 20000 of 50000, Stored update 20000 of 50000
-2 * Loglike = 1008.116090
Actual update 20050 of 50000, Stored update 20050 of 50000
-2 * Loglike = 1014.914085
Actual update 20100 of 50000, Stored update 20100 of 50000
-2 * Loglike = 1012.834730
Actual update 20150 of 50000, Stored update 20150 of 50000
-2 * Loglike = 1009.549109
Actual update 20200 of 50000, Stored update 20200 of 50000
-2 * Loglike = 1013.028455
Actual update 20250 of 50000, Stored update 20250 of 50000
-2 * Loglike = 1015.499438
Actual update 20300 of 50000, Stored update 20300 of 50000
-2 * Loglike = 1009.818919
Actual update 20350 of 50000, Stored update 20350 of 50000
-2 * Loglike = 1015.150664
Actual update 20400 of 50000, Stored update 20400 of 50000
-2 * Loglike = 1015.032915
Actual update 20450 of 50000, Stored update 20450 of 50000
-2 * Loglike = 1017.024621
Actual update 20500 of 50000, Stored update 20500 of 50000
-2 * Loglike = 1007.462265
Actual update 20550 of 50000, Stored update 20550 of 50000
-2 * Loglike = 1016.620455
Actual update 20600 of 50000, Stored update 20600 of 50000
-2 * Loglike = 1013.586755
Actual update 20650 of 50000, Stored update 20650 of 50000
-2 * Loglike = 1006.485915
Actual update 20700 of 50000, Stored update 20700 of 50000
-2 * Loglike = 1007.539889
Actual update 20750 of 50000, Stored update 20750 of 50000
-2 * Loglike = 1010.329178
Actual update 20800 of 50000, Stored update 20800 of 50000
-2 * Loglike = 1008.442124
Actual update 20850 of 50000, Stored update 20850 of 50000
-2 * Loglike = 1024.599475
Actual update 20900 of 50000, Stored update 20900 of 50000
-2 * Loglike = 1017.229863
Actual update 20950 of 50000, Stored update 20950 of 50000
-2 * Loglike = 1013.321109
Actual update 21000 of 50000, Stored update 21000 of 50000
-2 * Loglike = 1023.862171
Actual update 21050 of 50000, Stored update 21050 of 50000
-2 * Loglike = 1008.407557
Actual update 21100 of 50000, Stored update 21100 of 50000
-2 * Loglike = 1019.710814
Actual update 21150 of 50000, Stored update 21150 of 50000
-2 * Loglike = 1008.380720
Actual update 21200 of 50000, Stored update 21200 of 50000
-2 * Loglike = 1012.487123
Actual update 21250 of 50000, Stored update 21250 of 50000
-2 * Loglike = 1010.139225
Actual update 21300 of 50000, Stored update 21300 of 50000
-2 * Loglike = 1036.094170
Actual update 21350 of 50000, Stored update 21350 of 50000
-2 * Loglike = 1008.983449
Actual update 21400 of 50000, Stored update 21400 of 50000
-2 * Loglike = 1007.819138
Actual update 21450 of 50000, Stored update 21450 of 50000
-2 * Loglike = 1024.796560
Actual update 21500 of 50000, Stored update 21500 of 50000
-2 * Loglike = 1012.264836
Actual update 21550 of 50000, Stored update 21550 of 50000
-2 * Loglike = 1019.780883
Actual update 21600 of 50000, Stored update 21600 of 50000
-2 * Loglike = 1009.844781
Actual update 21650 of 50000, Stored update 21650 of 50000
-2 * Loglike = 1019.589689
Actual update 21700 of 50000, Stored update 21700 of 50000
-2 * Loglike = 1017.872601
Actual update 21750 of 50000, Stored update 21750 of 50000
-2 * Loglike = 1012.065225
Actual update 21800 of 50000, Stored update 21800 of 50000
-2 * Loglike = 1013.900490
Actual update 21850 of 50000, Stored update 21850 of 50000
-2 * Loglike = 1018.926663
Actual update 21900 of 50000, Stored update 21900 of 50000
-2 * Loglike = 1017.417450
Actual update 21950 of 50000, Stored update 21950 of 50000
-2 * Loglike = 1010.973600
Actual update 22000 of 50000, Stored update 22000 of 50000
-2 * Loglike = 1011.283403
Actual update 22050 of 50000, Stored update 22050 of 50000
-2 * Loglike = 1026.154802
Actual update 22100 of 50000, Stored update 22100 of 50000
-2 * Loglike = 1007.712712
Actual update 22150 of 50000, Stored update 22150 of 50000
-2 * Loglike = 1015.853452
Actual update 22200 of 50000, Stored update 22200 of 50000
-2 * Loglike = 1011.661123
Actual update 22250 of 50000, Stored update 22250 of 50000
-2 * Loglike = 1022.476893
Actual update 22300 of 50000, Stored update 22300 of 50000
-2 * Loglike = 1016.345279
Actual update 22350 of 50000, Stored update 22350 of 50000
-2 * Loglike = 1019.022622
Actual update 22400 of 50000, Stored update 22400 of 50000
-2 * Loglike = 1009.927381
Actual update 22450 of 50000, Stored update 22450 of 50000
-2 * Loglike = 1019.391549
Actual update 22500 of 50000, Stored update 22500 of 50000
-2 * Loglike = 1011.258039
Actual update 22550 of 50000, Stored update 22550 of 50000
-2 * Loglike = 1011.922250
Actual update 22600 of 50000, Stored update 22600 of 50000
-2 * Loglike = 1011.506125
Actual update 22650 of 50000, Stored update 22650 of 50000
-2 * Loglike = 1014.258347
Actual update 22700 of 50000, Stored update 22700 of 50000
-2 * Loglike = 1009.573482
Actual update 22750 of 50000, Stored update 22750 of 50000
-2 * Loglike = 1008.085839
Actual update 22800 of 50000, Stored update 22800 of 50000
-2 * Loglike = 1016.461796
Actual update 22850 of 50000, Stored update 22850 of 50000
-2 * Loglike = 1012.640404
Actual update 22900 of 50000, Stored update 22900 of 50000
-2 * Loglike = 1015.479433
Actual update 22950 of 50000, Stored update 22950 of 50000
-2 * Loglike = 1027.991799
Actual update 23000 of 50000, Stored update 23000 of 50000
-2 * Loglike = 1022.470886
Actual update 23050 of 50000, Stored update 23050 of 50000
-2 * Loglike = 1022.293930
Actual update 23100 of 50000, Stored update 23100 of 50000
-2 * Loglike = 1012.815088
Actual update 23150 of 50000, Stored update 23150 of 50000
-2 * Loglike = 1016.091044
Actual update 23200 of 50000, Stored update 23200 of 50000
-2 * Loglike = 1017.577212
Actual update 23250 of 50000, Stored update 23250 of 50000
-2 * Loglike = 1012.858533
Actual update 23300 of 50000, Stored update 23300 of 50000
-2 * Loglike = 1009.311580
Actual update 23350 of 50000, Stored update 23350 of 50000
-2 * Loglike = 1018.190200
Actual update 23400 of 50000, Stored update 23400 of 50000
-2 * Loglike = 1024.010302
Actual update 23450 of 50000, Stored update 23450 of 50000
-2 * Loglike = 1007.711802
Actual update 23500 of 50000, Stored update 23500 of 50000
-2 * Loglike = 1017.284549
Actual update 23550 of 50000, Stored update 23550 of 50000
-2 * Loglike = 1009.118059
Actual update 23600 of 50000, Stored update 23600 of 50000
-2 * Loglike = 1010.148871
Actual update 23650 of 50000, Stored update 23650 of 50000
-2 * Loglike = 1011.883918
Actual update 23700 of 50000, Stored update 23700 of 50000
-2 * Loglike = 1013.954376
Actual update 23750 of 50000, Stored update 23750 of 50000
-2 * Loglike = 1027.504648
Actual update 23800 of 50000, Stored update 23800 of 50000
-2 * Loglike = 1016.285125
Actual update 23850 of 50000, Stored update 23850 of 50000
-2 * Loglike = 1013.885297
Actual update 23900 of 50000, Stored update 23900 of 50000
-2 * Loglike = 1009.496975
Actual update 23950 of 50000, Stored update 23950 of 50000
-2 * Loglike = 1010.159219
Actual update 24000 of 50000, Stored update 24000 of 50000
-2 * Loglike = 1015.208346
Actual update 24050 of 50000, Stored update 24050 of 50000
-2 * Loglike = 1006.575928
Actual update 24100 of 50000, Stored update 24100 of 50000
-2 * Loglike = 1016.359051
Actual update 24150 of 50000, Stored update 24150 of 50000
-2 * Loglike = 1015.047199
Actual update 24200 of 50000, Stored update 24200 of 50000
-2 * Loglike = 1014.277479
Actual update 24250 of 50000, Stored update 24250 of 50000
-2 * Loglike = 1009.933575
Actual update 24300 of 50000, Stored update 24300 of 50000
-2 * Loglike = 1018.023459
Actual update 24350 of 50000, Stored update 24350 of 50000
-2 * Loglike = 1024.662163
Actual update 24400 of 50000, Stored update 24400 of 50000
-2 * Loglike = 1015.149503
Actual update 24450 of 50000, Stored update 24450 of 50000
-2 * Loglike = 1016.083246
Actual update 24500 of 50000, Stored update 24500 of 50000
-2 * Loglike = 1007.938614
Actual update 24550 of 50000, Stored update 24550 of 50000
-2 * Loglike = 1009.436042
Actual update 24600 of 50000, Stored update 24600 of 50000
-2 * Loglike = 1015.870148
Actual update 24650 of 50000, Stored update 24650 of 50000
-2 * Loglike = 1013.974038
Actual update 24700 of 50000, Stored update 24700 of 50000
-2 * Loglike = 1014.578963
Actual update 24750 of 50000, Stored update 24750 of 50000
-2 * Loglike = 1018.437836
Actual update 24800 of 50000, Stored update 24800 of 50000
-2 * Loglike = 1020.471033
Actual update 24850 of 50000, Stored update 24850 of 50000
-2 * Loglike = 1026.170236
Actual update 24900 of 50000, Stored update 24900 of 50000
-2 * Loglike = 1012.251847
Actual update 24950 of 50000, Stored update 24950 of 50000
-2 * Loglike = 1014.197894
Actual update 25000 of 50000, Stored update 25000 of 50000
-2 * Loglike = 1015.683607
Actual update 25050 of 50000, Stored update 25050 of 50000
-2 * Loglike = 1024.076611
Actual update 25100 of 50000, Stored update 25100 of 50000
-2 * Loglike = 1018.810024
Actual update 25150 of 50000, Stored update 25150 of 50000
-2 * Loglike = 1019.427145
Actual update 25200 of 50000, Stored update 25200 of 50000
-2 * Loglike = 1019.123683
Actual update 25250 of 50000, Stored update 25250 of 50000
-2 * Loglike = 1017.905714
Actual update 25300 of 50000, Stored update 25300 of 50000
-2 * Loglike = 1021.801265
Actual update 25350 of 50000, Stored update 25350 of 50000
-2 * Loglike = 1017.205694
Actual update 25400 of 50000, Stored update 25400 of 50000
-2 * Loglike = 1034.984244
Actual update 25450 of 50000, Stored update 25450 of 50000
-2 * Loglike = 1028.986578
Actual update 25500 of 50000, Stored update 25500 of 50000
-2 * Loglike = 1022.258114
Actual update 25550 of 50000, Stored update 25550 of 50000
-2 * Loglike = 1016.331365
Actual update 25600 of 50000, Stored update 25600 of 50000
-2 * Loglike = 1019.412121
Actual update 25650 of 50000, Stored update 25650 of 50000
-2 * Loglike = 1013.917172
Actual update 25700 of 50000, Stored update 25700 of 50000
-2 * Loglike = 1016.047232
Actual update 25750 of 50000, Stored update 25750 of 50000
-2 * Loglike = 1025.497118
Actual update 25800 of 50000, Stored update 25800 of 50000
-2 * Loglike = 1006.337525
Actual update 25850 of 50000, Stored update 25850 of 50000
-2 * Loglike = 1014.478307
Actual update 25900 of 50000, Stored update 25900 of 50000
-2 * Loglike = 1011.359115
Actual update 25950 of 50000, Stored update 25950 of 50000
-2 * Loglike = 1006.365110
Actual update 26000 of 50000, Stored update 26000 of 50000
-2 * Loglike = 1015.962959
Actual update 26050 of 50000, Stored update 26050 of 50000
-2 * Loglike = 1016.674367
Actual update 26100 of 50000, Stored update 26100 of 50000
-2 * Loglike = 1022.780432
Actual update 26150 of 50000, Stored update 26150 of 50000
-2 * Loglike = 1017.202858
Actual update 26200 of 50000, Stored update 26200 of 50000
-2 * Loglike = 1009.130041
Actual update 26250 of 50000, Stored update 26250 of 50000
-2 * Loglike = 1006.740061
Actual update 26300 of 50000, Stored update 26300 of 50000
-2 * Loglike = 1011.636715
Actual update 26350 of 50000, Stored update 26350 of 50000
-2 * Loglike = 1013.597405
Actual update 26400 of 50000, Stored update 26400 of 50000
-2 * Loglike = 1010.239446
Actual update 26450 of 50000, Stored update 26450 of 50000
-2 * Loglike = 1011.104915
Actual update 26500 of 50000, Stored update 26500 of 50000
-2 * Loglike = 1008.104440
Actual update 26550 of 50000, Stored update 26550 of 50000
-2 * Loglike = 1015.799914
Actual update 26600 of 50000, Stored update 26600 of 50000
-2 * Loglike = 1003.848428
Actual update 26650 of 50000, Stored update 26650 of 50000
-2 * Loglike = 1013.888753
Actual update 26700 of 50000, Stored update 26700 of 50000
-2 * Loglike = 1014.513030
Actual update 26750 of 50000, Stored update 26750 of 50000
-2 * Loglike = 1021.671935
Actual update 26800 of 50000, Stored update 26800 of 50000
-2 * Loglike = 1024.967841
Actual update 26850 of 50000, Stored update 26850 of 50000
-2 * Loglike = 1025.698001
Actual update 26900 of 50000, Stored update 26900 of 50000
-2 * Loglike = 1018.576035
Actual update 26950 of 50000, Stored update 26950 of 50000
-2 * Loglike = 1021.518675
Actual update 27000 of 50000, Stored update 27000 of 50000
-2 * Loglike = 1019.957116
Actual update 27050 of 50000, Stored update 27050 of 50000
-2 * Loglike = 1012.178108
Actual update 27100 of 50000, Stored update 27100 of 50000
-2 * Loglike = 1024.896844
Actual update 27150 of 50000, Stored update 27150 of 50000
-2 * Loglike = 1017.839483
Actual update 27200 of 50000, Stored update 27200 of 50000
-2 * Loglike = 1030.867788
Actual update 27250 of 50000, Stored update 27250 of 50000
-2 * Loglike = 1019.810063
Actual update 27300 of 50000, Stored update 27300 of 50000
-2 * Loglike = 1017.959577
Actual update 27350 of 50000, Stored update 27350 of 50000
-2 * Loglike = 1017.828613
Actual update 27400 of 50000, Stored update 27400 of 50000
-2 * Loglike = 1009.238477
Actual update 27450 of 50000, Stored update 27450 of 50000
-2 * Loglike = 1010.569009
Actual update 27500 of 50000, Stored update 27500 of 50000
-2 * Loglike = 1006.831235
Actual update 27550 of 50000, Stored update 27550 of 50000
-2 * Loglike = 1017.999358
Actual update 27600 of 50000, Stored update 27600 of 50000
-2 * Loglike = 1019.127802
Actual update 27650 of 50000, Stored update 27650 of 50000
-2 * Loglike = 1016.251828
Actual update 27700 of 50000, Stored update 27700 of 50000
-2 * Loglike = 1016.217064
Actual update 27750 of 50000, Stored update 27750 of 50000
-2 * Loglike = 1027.971835
Actual update 27800 of 50000, Stored update 27800 of 50000
-2 * Loglike = 1027.486288
Actual update 27850 of 50000, Stored update 27850 of 50000
-2 * Loglike = 1028.153723
Actual update 27900 of 50000, Stored update 27900 of 50000
-2 * Loglike = 1016.470015
Actual update 27950 of 50000, Stored update 27950 of 50000
-2 * Loglike = 1018.313147
Actual update 28000 of 50000, Stored update 28000 of 50000
-2 * Loglike = 1026.553348
Actual update 28050 of 50000, Stored update 28050 of 50000
-2 * Loglike = 1009.381410
Actual update 28100 of 50000, Stored update 28100 of 50000
-2 * Loglike = 1009.652690
Actual update 28150 of 50000, Stored update 28150 of 50000
-2 * Loglike = 1011.109354
Actual update 28200 of 50000, Stored update 28200 of 50000
-2 * Loglike = 1015.281891
Actual update 28250 of 50000, Stored update 28250 of 50000
-2 * Loglike = 1023.550608
Actual update 28300 of 50000, Stored update 28300 of 50000
-2 * Loglike = 1017.778718
Actual update 28350 of 50000, Stored update 28350 of 50000
-2 * Loglike = 1016.051687
Actual update 28400 of 50000, Stored update 28400 of 50000
-2 * Loglike = 1027.813728
Actual update 28450 of 50000, Stored update 28450 of 50000
-2 * Loglike = 1019.695040
Actual update 28500 of 50000, Stored update 28500 of 50000
-2 * Loglike = 1011.982923
Actual update 28550 of 50000, Stored update 28550 of 50000
-2 * Loglike = 1015.188889
Actual update 28600 of 50000, Stored update 28600 of 50000
-2 * Loglike = 1012.322938
Actual update 28650 of 50000, Stored update 28650 of 50000
-2 * Loglike = 1007.622548
Actual update 28700 of 50000, Stored update 28700 of 50000
-2 * Loglike = 1008.562231
Actual update 28750 of 50000, Stored update 28750 of 50000
-2 * Loglike = 1024.625866
Actual update 28800 of 50000, Stored update 28800 of 50000
-2 * Loglike = 1019.351011
Actual update 28850 of 50000, Stored update 28850 of 50000
-2 * Loglike = 1019.129922
Actual update 28900 of 50000, Stored update 28900 of 50000
-2 * Loglike = 1016.072980
Actual update 28950 of 50000, Stored update 28950 of 50000
-2 * Loglike = 1013.525024
Actual update 29000 of 50000, Stored update 29000 of 50000
-2 * Loglike = 1010.359947
Actual update 29050 of 50000, Stored update 29050 of 50000
-2 * Loglike = 1017.285945
Actual update 29100 of 50000, Stored update 29100 of 50000
-2 * Loglike = 1011.272531
Actual update 29150 of 50000, Stored update 29150 of 50000
-2 * Loglike = 1007.555593
Actual update 29200 of 50000, Stored update 29200 of 50000
-2 * Loglike = 1018.459172
Actual update 29250 of 50000, Stored update 29250 of 50000
-2 * Loglike = 1018.409756
Actual update 29300 of 50000, Stored update 29300 of 50000
-2 * Loglike = 1009.104247
Actual update 29350 of 50000, Stored update 29350 of 50000
-2 * Loglike = 1010.540608
Actual update 29400 of 50000, Stored update 29400 of 50000
-2 * Loglike = 1008.644610
Actual update 29450 of 50000, Stored update 29450 of 50000
-2 * Loglike = 1010.165215
Actual update 29500 of 50000, Stored update 29500 of 50000
-2 * Loglike = 1011.637166
Actual update 29550 of 50000, Stored update 29550 of 50000
-2 * Loglike = 1018.944187
Actual update 29600 of 50000, Stored update 29600 of 50000
-2 * Loglike = 1010.738152
Actual update 29650 of 50000, Stored update 29650 of 50000
-2 * Loglike = 1014.169370
Actual update 29700 of 50000, Stored update 29700 of 50000
-2 * Loglike = 1025.067415
Actual update 29750 of 50000, Stored update 29750 of 50000
-2 * Loglike = 1027.997261
Actual update 29800 of 50000, Stored update 29800 of 50000
-2 * Loglike = 1028.244813
Actual update 29850 of 50000, Stored update 29850 of 50000
-2 * Loglike = 1023.139731
Actual update 29900 of 50000, Stored update 29900 of 50000
-2 * Loglike = 1019.390160
Actual update 29950 of 50000, Stored update 29950 of 50000
-2 * Loglike = 1015.549473
Actual update 30000 of 50000, Stored update 30000 of 50000
-2 * Loglike = 1012.559741
Actual update 30050 of 50000, Stored update 30050 of 50000
-2 * Loglike = 1016.445101
Actual update 30100 of 50000, Stored update 30100 of 50000
-2 * Loglike = 1009.298084
Actual update 30150 of 50000, Stored update 30150 of 50000
-2 * Loglike = 1019.888292
Actual update 30200 of 50000, Stored update 30200 of 50000
-2 * Loglike = 1016.419789
Actual update 30250 of 50000, Stored update 30250 of 50000
-2 * Loglike = 1007.438551
Actual update 30300 of 50000, Stored update 30300 of 50000
-2 * Loglike = 1020.539933
Actual update 30350 of 50000, Stored update 30350 of 50000
-2 * Loglike = 1010.305648
Actual update 30400 of 50000, Stored update 30400 of 50000
-2 * Loglike = 1010.449264
Actual update 30450 of 50000, Stored update 30450 of 50000
-2 * Loglike = 1009.854220
Actual update 30500 of 50000, Stored update 30500 of 50000
-2 * Loglike = 1012.158465
Actual update 30550 of 50000, Stored update 30550 of 50000
-2 * Loglike = 1012.921442
Actual update 30600 of 50000, Stored update 30600 of 50000
-2 * Loglike = 1008.520352
Actual update 30650 of 50000, Stored update 30650 of 50000
-2 * Loglike = 1009.828564
Actual update 30700 of 50000, Stored update 30700 of 50000
-2 * Loglike = 1010.665656
Actual update 30750 of 50000, Stored update 30750 of 50000
-2 * Loglike = 1008.827683
Actual update 30800 of 50000, Stored update 30800 of 50000
-2 * Loglike = 1013.748627
Actual update 30850 of 50000, Stored update 30850 of 50000
-2 * Loglike = 1014.016626
Actual update 30900 of 50000, Stored update 30900 of 50000
-2 * Loglike = 1009.488027
Actual update 30950 of 50000, Stored update 30950 of 50000
-2 * Loglike = 1018.841304
Actual update 31000 of 50000, Stored update 31000 of 50000
-2 * Loglike = 1010.640491
Actual update 31050 of 50000, Stored update 31050 of 50000
-2 * Loglike = 1010.569249
Actual update 31100 of 50000, Stored update 31100 of 50000
-2 * Loglike = 1015.587757
Actual update 31150 of 50000, Stored update 31150 of 50000
-2 * Loglike = 1013.065733
Actual update 31200 of 50000, Stored update 31200 of 50000
-2 * Loglike = 1015.544068
Actual update 31250 of 50000, Stored update 31250 of 50000
-2 * Loglike = 1019.720311
Actual update 31300 of 50000, Stored update 31300 of 50000
-2 * Loglike = 1015.961348
Actual update 31350 of 50000, Stored update 31350 of 50000
-2 * Loglike = 1014.750244
Actual update 31400 of 50000, Stored update 31400 of 50000
-2 * Loglike = 1012.929770
Actual update 31450 of 50000, Stored update 31450 of 50000
-2 * Loglike = 1013.852600
Actual update 31500 of 50000, Stored update 31500 of 50000
-2 * Loglike = 1015.369091
Actual update 31550 of 50000, Stored update 31550 of 50000
-2 * Loglike = 1015.497320
Actual update 31600 of 50000, Stored update 31600 of 50000
-2 * Loglike = 1015.071488
Actual update 31650 of 50000, Stored update 31650 of 50000
-2 * Loglike = 1016.338807
Actual update 31700 of 50000, Stored update 31700 of 50000
-2 * Loglike = 1012.218343
Actual update 31750 of 50000, Stored update 31750 of 50000
-2 * Loglike = 1008.066460
Actual update 31800 of 50000, Stored update 31800 of 50000
-2 * Loglike = 1011.421021
Actual update 31850 of 50000, Stored update 31850 of 50000
-2 * Loglike = 1014.135595
Actual update 31900 of 50000, Stored update 31900 of 50000
-2 * Loglike = 1013.763067
Actual update 31950 of 50000, Stored update 31950 of 50000
-2 * Loglike = 1010.865429
Actual update 32000 of 50000, Stored update 32000 of 50000
-2 * Loglike = 1009.491903
Actual update 32050 of 50000, Stored update 32050 of 50000
-2 * Loglike = 1019.843158
Actual update 32100 of 50000, Stored update 32100 of 50000
-2 * Loglike = 1009.841154
Actual update 32150 of 50000, Stored update 32150 of 50000
-2 * Loglike = 1010.354794
Actual update 32200 of 50000, Stored update 32200 of 50000
-2 * Loglike = 1009.453904
Actual update 32250 of 50000, Stored update 32250 of 50000
-2 * Loglike = 1011.909461
Actual update 32300 of 50000, Stored update 32300 of 50000
-2 * Loglike = 1012.370847
Actual update 32350 of 50000, Stored update 32350 of 50000
-2 * Loglike = 1021.754810
Actual update 32400 of 50000, Stored update 32400 of 50000
-2 * Loglike = 1013.972030
Actual update 32450 of 50000, Stored update 32450 of 50000
-2 * Loglike = 1017.908593
Actual update 32500 of 50000, Stored update 32500 of 50000
-2 * Loglike = 1035.470390
Actual update 32550 of 50000, Stored update 32550 of 50000
-2 * Loglike = 1020.341232
Actual update 32600 of 50000, Stored update 32600 of 50000
-2 * Loglike = 1015.267760
Actual update 32650 of 50000, Stored update 32650 of 50000
-2 * Loglike = 1024.877665
Actual update 32700 of 50000, Stored update 32700 of 50000
-2 * Loglike = 1015.458002
Actual update 32750 of 50000, Stored update 32750 of 50000
-2 * Loglike = 1027.818467
Actual update 32800 of 50000, Stored update 32800 of 50000
-2 * Loglike = 1020.135009
Actual update 32850 of 50000, Stored update 32850 of 50000
-2 * Loglike = 1012.826095
Actual update 32900 of 50000, Stored update 32900 of 50000
-2 * Loglike = 1009.663431
Actual update 32950 of 50000, Stored update 32950 of 50000
-2 * Loglike = 1015.366448
Actual update 33000 of 50000, Stored update 33000 of 50000
-2 * Loglike = 1008.169912
Actual update 33050 of 50000, Stored update 33050 of 50000
-2 * Loglike = 1020.652011
Actual update 33100 of 50000, Stored update 33100 of 50000
-2 * Loglike = 1009.426871
Actual update 33150 of 50000, Stored update 33150 of 50000
-2 * Loglike = 1019.629770
Actual update 33200 of 50000, Stored update 33200 of 50000
-2 * Loglike = 1013.531113
Actual update 33250 of 50000, Stored update 33250 of 50000
-2 * Loglike = 1005.268348
Actual update 33300 of 50000, Stored update 33300 of 50000
-2 * Loglike = 1007.902460
Actual update 33350 of 50000, Stored update 33350 of 50000
-2 * Loglike = 1004.988357
Actual update 33400 of 50000, Stored update 33400 of 50000
-2 * Loglike = 1013.838627
Actual update 33450 of 50000, Stored update 33450 of 50000
-2 * Loglike = 1012.397271
Actual update 33500 of 50000, Stored update 33500 of 50000
-2 * Loglike = 1013.227700
Actual update 33550 of 50000, Stored update 33550 of 50000
-2 * Loglike = 1020.590010
Actual update 33600 of 50000, Stored update 33600 of 50000
-2 * Loglike = 1012.941044
Actual update 33650 of 50000, Stored update 33650 of 50000
-2 * Loglike = 1011.052135
Actual update 33700 of 50000, Stored update 33700 of 50000
-2 * Loglike = 1026.050789
Actual update 33750 of 50000, Stored update 33750 of 50000
-2 * Loglike = 1012.151825
Actual update 33800 of 50000, Stored update 33800 of 50000
-2 * Loglike = 1022.222604
Actual update 33850 of 50000, Stored update 33850 of 50000
-2 * Loglike = 1027.136631
Actual update 33900 of 50000, Stored update 33900 of 50000
-2 * Loglike = 1021.946383
Actual update 33950 of 50000, Stored update 33950 of 50000
-2 * Loglike = 1020.448561
Actual update 34000 of 50000, Stored update 34000 of 50000
-2 * Loglike = 1013.133517
Actual update 34050 of 50000, Stored update 34050 of 50000
-2 * Loglike = 1015.814014
Actual update 34100 of 50000, Stored update 34100 of 50000
-2 * Loglike = 1012.744238
Actual update 34150 of 50000, Stored update 34150 of 50000
-2 * Loglike = 1020.464689
Actual update 34200 of 50000, Stored update 34200 of 50000
-2 * Loglike = 1026.260038
Actual update 34250 of 50000, Stored update 34250 of 50000
-2 * Loglike = 1019.960607
Actual update 34300 of 50000, Stored update 34300 of 50000
-2 * Loglike = 1016.353921
Actual update 34350 of 50000, Stored update 34350 of 50000
-2 * Loglike = 1011.551472
Actual update 34400 of 50000, Stored update 34400 of 50000
-2 * Loglike = 1016.657709
Actual update 34450 of 50000, Stored update 34450 of 50000
-2 * Loglike = 1010.933789
Actual update 34500 of 50000, Stored update 34500 of 50000
-2 * Loglike = 1008.900530
Actual update 34550 of 50000, Stored update 34550 of 50000
-2 * Loglike = 1015.680262
Actual update 34600 of 50000, Stored update 34600 of 50000
-2 * Loglike = 1017.629415
Actual update 34650 of 50000, Stored update 34650 of 50000
-2 * Loglike = 1015.354174
Actual update 34700 of 50000, Stored update 34700 of 50000
-2 * Loglike = 1012.156396
Actual update 34750 of 50000, Stored update 34750 of 50000
-2 * Loglike = 1020.760882
Actual update 34800 of 50000, Stored update 34800 of 50000
-2 * Loglike = 1020.982360
Actual update 34850 of 50000, Stored update 34850 of 50000
-2 * Loglike = 1024.377973
Actual update 34900 of 50000, Stored update 34900 of 50000
-2 * Loglike = 1018.592918
Actual update 34950 of 50000, Stored update 34950 of 50000
-2 * Loglike = 1018.727559
Actual update 35000 of 50000, Stored update 35000 of 50000
-2 * Loglike = 1023.397213
Actual update 35050 of 50000, Stored update 35050 of 50000
-2 * Loglike = 1015.548981
Actual update 35100 of 50000, Stored update 35100 of 50000
-2 * Loglike = 1011.803388
Actual update 35150 of 50000, Stored update 35150 of 50000
-2 * Loglike = 1013.825515
Actual update 35200 of 50000, Stored update 35200 of 50000
-2 * Loglike = 1013.229506
Actual update 35250 of 50000, Stored update 35250 of 50000
-2 * Loglike = 1012.129489
Actual update 35300 of 50000, Stored update 35300 of 50000
-2 * Loglike = 1013.919177
Actual update 35350 of 50000, Stored update 35350 of 50000
-2 * Loglike = 1013.545770
Actual update 35400 of 50000, Stored update 35400 of 50000
-2 * Loglike = 1018.590127
Actual update 35450 of 50000, Stored update 35450 of 50000
-2 * Loglike = 1012.696902
Actual update 35500 of 50000, Stored update 35500 of 50000
-2 * Loglike = 1019.111370
Actual update 35550 of 50000, Stored update 35550 of 50000
-2 * Loglike = 1009.397282
Actual update 35600 of 50000, Stored update 35600 of 50000
-2 * Loglike = 1014.033750
Actual update 35650 of 50000, Stored update 35650 of 50000
-2 * Loglike = 1015.681981
Actual update 35700 of 50000, Stored update 35700 of 50000
-2 * Loglike = 1019.474222
Actual update 35750 of 50000, Stored update 35750 of 50000
-2 * Loglike = 1020.951001
Actual update 35800 of 50000, Stored update 35800 of 50000
-2 * Loglike = 1015.564042
Actual update 35850 of 50000, Stored update 35850 of 50000
-2 * Loglike = 1023.304159
Actual update 35900 of 50000, Stored update 35900 of 50000
-2 * Loglike = 1012.054709
Actual update 35950 of 50000, Stored update 35950 of 50000
-2 * Loglike = 1018.019549
Actual update 36000 of 50000, Stored update 36000 of 50000
-2 * Loglike = 1006.783736
Actual update 36050 of 50000, Stored update 36050 of 50000
-2 * Loglike = 1017.440424
Actual update 36100 of 50000, Stored update 36100 of 50000
-2 * Loglike = 1016.330838
Actual update 36150 of 50000, Stored update 36150 of 50000
-2 * Loglike = 1014.159914
Actual update 36200 of 50000, Stored update 36200 of 50000
-2 * Loglike = 1012.568194
Actual update 36250 of 50000, Stored update 36250 of 50000
-2 * Loglike = 1015.315274
Actual update 36300 of 50000, Stored update 36300 of 50000
-2 * Loglike = 1017.131358
Actual update 36350 of 50000, Stored update 36350 of 50000
-2 * Loglike = 1013.262969
Actual update 36400 of 50000, Stored update 36400 of 50000
-2 * Loglike = 1015.374368
Actual update 36450 of 50000, Stored update 36450 of 50000
-2 * Loglike = 1018.705718
Actual update 36500 of 50000, Stored update 36500 of 50000
-2 * Loglike = 1025.043830
Actual update 36550 of 50000, Stored update 36550 of 50000
-2 * Loglike = 1022.762267
Actual update 36600 of 50000, Stored update 36600 of 50000
-2 * Loglike = 1025.127996
Actual update 36650 of 50000, Stored update 36650 of 50000
-2 * Loglike = 1029.206830
Actual update 36700 of 50000, Stored update 36700 of 50000
-2 * Loglike = 1022.386266
Actual update 36750 of 50000, Stored update 36750 of 50000
-2 * Loglike = 1027.754723
Actual update 36800 of 50000, Stored update 36800 of 50000
-2 * Loglike = 1032.647243
Actual update 36850 of 50000, Stored update 36850 of 50000
-2 * Loglike = 1021.307476
Actual update 36900 of 50000, Stored update 36900 of 50000
-2 * Loglike = 1010.026086
Actual update 36950 of 50000, Stored update 36950 of 50000
-2 * Loglike = 1006.994246
Actual update 37000 of 50000, Stored update 37000 of 50000
-2 * Loglike = 1016.317403
Actual update 37050 of 50000, Stored update 37050 of 50000
-2 * Loglike = 1012.847418
Actual update 37100 of 50000, Stored update 37100 of 50000
-2 * Loglike = 1016.037799
Actual update 37150 of 50000, Stored update 37150 of 50000
-2 * Loglike = 1013.059868
Actual update 37200 of 50000, Stored update 37200 of 50000
-2 * Loglike = 1008.529177
Actual update 37250 of 50000, Stored update 37250 of 50000
-2 * Loglike = 1009.061621
Actual update 37300 of 50000, Stored update 37300 of 50000
-2 * Loglike = 1019.354116
Actual update 37350 of 50000, Stored update 37350 of 50000
-2 * Loglike = 1014.197516
Actual update 37400 of 50000, Stored update 37400 of 50000
-2 * Loglike = 1026.139144
Actual update 37450 of 50000, Stored update 37450 of 50000
-2 * Loglike = 1011.405860
Actual update 37500 of 50000, Stored update 37500 of 50000
-2 * Loglike = 1010.975653
Actual update 37550 of 50000, Stored update 37550 of 50000
-2 * Loglike = 1012.891119
Actual update 37600 of 50000, Stored update 37600 of 50000
-2 * Loglike = 1011.632942
Actual update 37650 of 50000, Stored update 37650 of 50000
-2 * Loglike = 1018.769010
Actual update 37700 of 50000, Stored update 37700 of 50000
-2 * Loglike = 1011.807313
Actual update 37750 of 50000, Stored update 37750 of 50000
-2 * Loglike = 1008.242646
Actual update 37800 of 50000, Stored update 37800 of 50000
-2 * Loglike = 1011.905329
Actual update 37850 of 50000, Stored update 37850 of 50000
-2 * Loglike = 1009.383226
Actual update 37900 of 50000, Stored update 37900 of 50000
-2 * Loglike = 1009.079231
Actual update 37950 of 50000, Stored update 37950 of 50000
-2 * Loglike = 1011.196232
Actual update 38000 of 50000, Stored update 38000 of 50000
-2 * Loglike = 1009.596920
Actual update 38050 of 50000, Stored update 38050 of 50000
-2 * Loglike = 1009.357757
Actual update 38100 of 50000, Stored update 38100 of 50000
-2 * Loglike = 1012.354520
Actual update 38150 of 50000, Stored update 38150 of 50000
-2 * Loglike = 1008.633654
Actual update 38200 of 50000, Stored update 38200 of 50000
-2 * Loglike = 1007.787091
Actual update 38250 of 50000, Stored update 38250 of 50000
-2 * Loglike = 1019.915410
Actual update 38300 of 50000, Stored update 38300 of 50000
-2 * Loglike = 1018.277126
Actual update 38350 of 50000, Stored update 38350 of 50000
-2 * Loglike = 1012.971383
Actual update 38400 of 50000, Stored update 38400 of 50000
-2 * Loglike = 1018.313935
Actual update 38450 of 50000, Stored update 38450 of 50000
-2 * Loglike = 1014.647677
Actual update 38500 of 50000, Stored update 38500 of 50000
-2 * Loglike = 1011.328150
Actual update 38550 of 50000, Stored update 38550 of 50000
-2 * Loglike = 1015.963671
Actual update 38600 of 50000, Stored update 38600 of 50000
-2 * Loglike = 1026.570854
Actual update 38650 of 50000, Stored update 38650 of 50000
-2 * Loglike = 1015.712596
Actual update 38700 of 50000, Stored update 38700 of 50000
-2 * Loglike = 1011.918340
Actual update 38750 of 50000, Stored update 38750 of 50000
-2 * Loglike = 1031.302700
Actual update 38800 of 50000, Stored update 38800 of 50000
-2 * Loglike = 1012.284909
Actual update 38850 of 50000, Stored update 38850 of 50000
-2 * Loglike = 1007.900828
Actual update 38900 of 50000, Stored update 38900 of 50000
-2 * Loglike = 1021.741261
Actual update 38950 of 50000, Stored update 38950 of 50000
-2 * Loglike = 1013.899148
Actual update 39000 of 50000, Stored update 39000 of 50000
-2 * Loglike = 1027.280571
Actual update 39050 of 50000, Stored update 39050 of 50000
-2 * Loglike = 1020.622980
Actual update 39100 of 50000, Stored update 39100 of 50000
-2 * Loglike = 1020.591263
Actual update 39150 of 50000, Stored update 39150 of 50000
-2 * Loglike = 1023.001835
Actual update 39200 of 50000, Stored update 39200 of 50000
-2 * Loglike = 1018.867199
Actual update 39250 of 50000, Stored update 39250 of 50000
-2 * Loglike = 1021.601699
Actual update 39300 of 50000, Stored update 39300 of 50000
-2 * Loglike = 1030.030185
Actual update 39350 of 50000, Stored update 39350 of 50000
-2 * Loglike = 1011.796341
Actual update 39400 of 50000, Stored update 39400 of 50000
-2 * Loglike = 1014.129014
Actual update 39450 of 50000, Stored update 39450 of 50000
-2 * Loglike = 1025.046845
Actual update 39500 of 50000, Stored update 39500 of 50000
-2 * Loglike = 1014.857856
Actual update 39550 of 50000, Stored update 39550 of 50000
-2 * Loglike = 1024.664074
Actual update 39600 of 50000, Stored update 39600 of 50000
-2 * Loglike = 1013.072882
Actual update 39650 of 50000, Stored update 39650 of 50000
-2 * Loglike = 1013.469086
Actual update 39700 of 50000, Stored update 39700 of 50000
-2 * Loglike = 1010.017115
Actual update 39750 of 50000, Stored update 39750 of 50000
-2 * Loglike = 1015.126914
Actual update 39800 of 50000, Stored update 39800 of 50000
-2 * Loglike = 1013.384563
Actual update 39850 of 50000, Stored update 39850 of 50000
-2 * Loglike = 1011.412561
Actual update 39900 of 50000, Stored update 39900 of 50000
-2 * Loglike = 1018.904279
Actual update 39950 of 50000, Stored update 39950 of 50000
-2 * Loglike = 1010.714926
Actual update 40000 of 50000, Stored update 40000 of 50000
-2 * Loglike = 1017.231939
Actual update 40050 of 50000, Stored update 40050 of 50000
-2 * Loglike = 1013.026520
Actual update 40100 of 50000, Stored update 40100 of 50000
-2 * Loglike = 1019.188418
Actual update 40150 of 50000, Stored update 40150 of 50000
-2 * Loglike = 1018.659504
Actual update 40200 of 50000, Stored update 40200 of 50000
-2 * Loglike = 1013.097888
Actual update 40250 of 50000, Stored update 40250 of 50000
-2 * Loglike = 1014.779167
Actual update 40300 of 50000, Stored update 40300 of 50000
-2 * Loglike = 1013.355276
Actual update 40350 of 50000, Stored update 40350 of 50000
-2 * Loglike = 1008.748162
Actual update 40400 of 50000, Stored update 40400 of 50000
-2 * Loglike = 1008.588878
Actual update 40450 of 50000, Stored update 40450 of 50000
-2 * Loglike = 1019.625305
Actual update 40500 of 50000, Stored update 40500 of 50000
-2 * Loglike = 1017.510828
Actual update 40550 of 50000, Stored update 40550 of 50000
-2 * Loglike = 1017.392080
Actual update 40600 of 50000, Stored update 40600 of 50000
-2 * Loglike = 1015.891251
Actual update 40650 of 50000, Stored update 40650 of 50000
-2 * Loglike = 1012.005119
Actual update 40700 of 50000, Stored update 40700 of 50000
-2 * Loglike = 1005.015218
Actual update 40750 of 50000, Stored update 40750 of 50000
-2 * Loglike = 1014.367835
Actual update 40800 of 50000, Stored update 40800 of 50000
-2 * Loglike = 1012.979552
Actual update 40850 of 50000, Stored update 40850 of 50000
-2 * Loglike = 1014.770270
Actual update 40900 of 50000, Stored update 40900 of 50000
-2 * Loglike = 1020.709008
Actual update 40950 of 50000, Stored update 40950 of 50000
-2 * Loglike = 1035.122600
Actual update 41000 of 50000, Stored update 41000 of 50000
-2 * Loglike = 1014.023050
Actual update 41050 of 50000, Stored update 41050 of 50000
-2 * Loglike = 1014.165334
Actual update 41100 of 50000, Stored update 41100 of 50000
-2 * Loglike = 1018.502552
Actual update 41150 of 50000, Stored update 41150 of 50000
-2 * Loglike = 1016.269849
Actual update 41200 of 50000, Stored update 41200 of 50000
-2 * Loglike = 1008.727931
Actual update 41250 of 50000, Stored update 41250 of 50000
-2 * Loglike = 1010.209905
Actual update 41300 of 50000, Stored update 41300 of 50000
-2 * Loglike = 1013.716790
Actual update 41350 of 50000, Stored update 41350 of 50000
-2 * Loglike = 1008.854623
Actual update 41400 of 50000, Stored update 41400 of 50000
-2 * Loglike = 1007.867554
Actual update 41450 of 50000, Stored update 41450 of 50000
-2 * Loglike = 1019.844819
Actual update 41500 of 50000, Stored update 41500 of 50000
-2 * Loglike = 1012.865395
Actual update 41550 of 50000, Stored update 41550 of 50000
-2 * Loglike = 1010.862695
Actual update 41600 of 50000, Stored update 41600 of 50000
-2 * Loglike = 1013.184797
Actual update 41650 of 50000, Stored update 41650 of 50000
-2 * Loglike = 1013.370214
Actual update 41700 of 50000, Stored update 41700 of 50000
-2 * Loglike = 1017.001634
Actual update 41750 of 50000, Stored update 41750 of 50000
-2 * Loglike = 1011.451369
Actual update 41800 of 50000, Stored update 41800 of 50000
-2 * Loglike = 1020.339165
Actual update 41850 of 50000, Stored update 41850 of 50000
-2 * Loglike = 1021.522763
Actual update 41900 of 50000, Stored update 41900 of 50000
-2 * Loglike = 1020.968648
Actual update 41950 of 50000, Stored update 41950 of 50000
-2 * Loglike = 1009.717406
Actual update 42000 of 50000, Stored update 42000 of 50000
-2 * Loglike = 1030.639046
Actual update 42050 of 50000, Stored update 42050 of 50000
-2 * Loglike = 1018.397010
Actual update 42100 of 50000, Stored update 42100 of 50000
-2 * Loglike = 1013.547798
Actual update 42150 of 50000, Stored update 42150 of 50000
-2 * Loglike = 1010.420071
Actual update 42200 of 50000, Stored update 42200 of 50000
-2 * Loglike = 1011.766066
Actual update 42250 of 50000, Stored update 42250 of 50000
-2 * Loglike = 1010.380988
Actual update 42300 of 50000, Stored update 42300 of 50000
-2 * Loglike = 1011.507273
Actual update 42350 of 50000, Stored update 42350 of 50000
-2 * Loglike = 1016.504284
Actual update 42400 of 50000, Stored update 42400 of 50000
-2 * Loglike = 1016.155978
Actual update 42450 of 50000, Stored update 42450 of 50000
-2 * Loglike = 1016.787153
Actual update 42500 of 50000, Stored update 42500 of 50000
-2 * Loglike = 1008.216854
Actual update 42550 of 50000, Stored update 42550 of 50000
-2 * Loglike = 1012.760015
Actual update 42600 of 50000, Stored update 42600 of 50000
-2 * Loglike = 1012.037720
Actual update 42650 of 50000, Stored update 42650 of 50000
-2 * Loglike = 1022.118399
Actual update 42700 of 50000, Stored update 42700 of 50000
-2 * Loglike = 1011.166316
Actual update 42750 of 50000, Stored update 42750 of 50000
-2 * Loglike = 1019.940299
Actual update 42800 of 50000, Stored update 42800 of 50000
-2 * Loglike = 1017.190702
Actual update 42850 of 50000, Stored update 42850 of 50000
-2 * Loglike = 1013.278190
Actual update 42900 of 50000, Stored update 42900 of 50000
-2 * Loglike = 1010.200481
Actual update 42950 of 50000, Stored update 42950 of 50000
-2 * Loglike = 1021.507479
Actual update 43000 of 50000, Stored update 43000 of 50000
-2 * Loglike = 1017.927708
Actual update 43050 of 50000, Stored update 43050 of 50000
-2 * Loglike = 1009.208135
Actual update 43100 of 50000, Stored update 43100 of 50000
-2 * Loglike = 1020.376719
Actual update 43150 of 50000, Stored update 43150 of 50000
-2 * Loglike = 1020.535095
Actual update 43200 of 50000, Stored update 43200 of 50000
-2 * Loglike = 1016.009378
Actual update 43250 of 50000, Stored update 43250 of 50000
-2 * Loglike = 1008.689612
Actual update 43300 of 50000, Stored update 43300 of 50000
-2 * Loglike = 1007.654226
Actual update 43350 of 50000, Stored update 43350 of 50000
-2 * Loglike = 1008.987489
Actual update 43400 of 50000, Stored update 43400 of 50000
-2 * Loglike = 1015.178445
Actual update 43450 of 50000, Stored update 43450 of 50000
-2 * Loglike = 1014.867810
Actual update 43500 of 50000, Stored update 43500 of 50000
-2 * Loglike = 1020.030957
Actual update 43550 of 50000, Stored update 43550 of 50000
-2 * Loglike = 1007.610400
Actual update 43600 of 50000, Stored update 43600 of 50000
-2 * Loglike = 1010.148931
Actual update 43650 of 50000, Stored update 43650 of 50000
-2 * Loglike = 1008.806970
Actual update 43700 of 50000, Stored update 43700 of 50000
-2 * Loglike = 1016.064025
Actual update 43750 of 50000, Stored update 43750 of 50000
-2 * Loglike = 1009.738028
Actual update 43800 of 50000, Stored update 43800 of 50000
-2 * Loglike = 1013.415551
Actual update 43850 of 50000, Stored update 43850 of 50000
-2 * Loglike = 1016.603614
Actual update 43900 of 50000, Stored update 43900 of 50000
-2 * Loglike = 1018.773351
Actual update 43950 of 50000, Stored update 43950 of 50000
-2 * Loglike = 1029.086015
Actual update 44000 of 50000, Stored update 44000 of 50000
-2 * Loglike = 1025.305952
Actual update 44050 of 50000, Stored update 44050 of 50000
-2 * Loglike = 1020.996233
Actual update 44100 of 50000, Stored update 44100 of 50000
-2 * Loglike = 1010.781738
Actual update 44150 of 50000, Stored update 44150 of 50000
-2 * Loglike = 1018.208929
Actual update 44200 of 50000, Stored update 44200 of 50000
-2 * Loglike = 1008.377905
Actual update 44250 of 50000, Stored update 44250 of 50000
-2 * Loglike = 1012.540591
Actual update 44300 of 50000, Stored update 44300 of 50000
-2 * Loglike = 1024.381462
Actual update 44350 of 50000, Stored update 44350 of 50000
-2 * Loglike = 1008.168419
Actual update 44400 of 50000, Stored update 44400 of 50000
-2 * Loglike = 1005.106823
Actual update 44450 of 50000, Stored update 44450 of 50000
-2 * Loglike = 1023.840554
Actual update 44500 of 50000, Stored update 44500 of 50000
-2 * Loglike = 1011.912229
Actual update 44550 of 50000, Stored update 44550 of 50000
-2 * Loglike = 1012.566833
Actual update 44600 of 50000, Stored update 44600 of 50000
-2 * Loglike = 1017.629798
Actual update 44650 of 50000, Stored update 44650 of 50000
-2 * Loglike = 1019.934424
Actual update 44700 of 50000, Stored update 44700 of 50000
-2 * Loglike = 1015.956862
Actual update 44750 of 50000, Stored update 44750 of 50000
-2 * Loglike = 1014.499485
Actual update 44800 of 50000, Stored update 44800 of 50000
-2 * Loglike = 1021.009293
Actual update 44850 of 50000, Stored update 44850 of 50000
-2 * Loglike = 1010.520371
Actual update 44900 of 50000, Stored update 44900 of 50000
-2 * Loglike = 1017.373001
Actual update 44950 of 50000, Stored update 44950 of 50000
-2 * Loglike = 1013.058972
Actual update 45000 of 50000, Stored update 45000 of 50000
-2 * Loglike = 1026.855867
Actual update 45050 of 50000, Stored update 45050 of 50000
-2 * Loglike = 1013.111981
Actual update 45100 of 50000, Stored update 45100 of 50000
-2 * Loglike = 1017.537453
Actual update 45150 of 50000, Stored update 45150 of 50000
-2 * Loglike = 1012.799952
Actual update 45200 of 50000, Stored update 45200 of 50000
-2 * Loglike = 1011.114277
Actual update 45250 of 50000, Stored update 45250 of 50000
-2 * Loglike = 1013.771763
Actual update 45300 of 50000, Stored update 45300 of 50000
-2 * Loglike = 1018.792987
Actual update 45350 of 50000, Stored update 45350 of 50000
-2 * Loglike = 1013.366007
Actual update 45400 of 50000, Stored update 45400 of 50000
-2 * Loglike = 1024.310437
Actual update 45450 of 50000, Stored update 45450 of 50000
-2 * Loglike = 1014.381273
Actual update 45500 of 50000, Stored update 45500 of 50000
-2 * Loglike = 1014.674063
Actual update 45550 of 50000, Stored update 45550 of 50000
-2 * Loglike = 1019.336134
Actual update 45600 of 50000, Stored update 45600 of 50000
-2 * Loglike = 1019.248715
Actual update 45650 of 50000, Stored update 45650 of 50000
-2 * Loglike = 1007.541375
Actual update 45700 of 50000, Stored update 45700 of 50000
-2 * Loglike = 1012.254447
Actual update 45750 of 50000, Stored update 45750 of 50000
-2 * Loglike = 1014.266055
Actual update 45800 of 50000, Stored update 45800 of 50000
-2 * Loglike = 1016.259406
Actual update 45850 of 50000, Stored update 45850 of 50000
-2 * Loglike = 1021.590632
Actual update 45900 of 50000, Stored update 45900 of 50000
-2 * Loglike = 1014.222469
Actual update 45950 of 50000, Stored update 45950 of 50000
-2 * Loglike = 1010.156540
Actual update 46000 of 50000, Stored update 46000 of 50000
-2 * Loglike = 1012.608806
Actual update 46050 of 50000, Stored update 46050 of 50000
-2 * Loglike = 1024.177928
Actual update 46100 of 50000, Stored update 46100 of 50000
-2 * Loglike = 1012.818571
Actual update 46150 of 50000, Stored update 46150 of 50000
-2 * Loglike = 1013.945776
Actual update 46200 of 50000, Stored update 46200 of 50000
-2 * Loglike = 1009.852223
Actual update 46250 of 50000, Stored update 46250 of 50000
-2 * Loglike = 1014.777561
Actual update 46300 of 50000, Stored update 46300 of 50000
-2 * Loglike = 1011.414043
Actual update 46350 of 50000, Stored update 46350 of 50000
-2 * Loglike = 1011.371580
Actual update 46400 of 50000, Stored update 46400 of 50000
-2 * Loglike = 1009.434461
Actual update 46450 of 50000, Stored update 46450 of 50000
-2 * Loglike = 1007.559185
Actual update 46500 of 50000, Stored update 46500 of 50000
-2 * Loglike = 1018.883699
Actual update 46550 of 50000, Stored update 46550 of 50000
-2 * Loglike = 1009.948922
Actual update 46600 of 50000, Stored update 46600 of 50000
-2 * Loglike = 1013.289432
Actual update 46650 of 50000, Stored update 46650 of 50000
-2 * Loglike = 1014.074470
Actual update 46700 of 50000, Stored update 46700 of 50000
-2 * Loglike = 1018.217341
Actual update 46750 of 50000, Stored update 46750 of 50000
-2 * Loglike = 1017.515874
Actual update 46800 of 50000, Stored update 46800 of 50000
-2 * Loglike = 1015.569662
Actual update 46850 of 50000, Stored update 46850 of 50000
-2 * Loglike = 1012.528388
Actual update 46900 of 50000, Stored update 46900 of 50000
-2 * Loglike = 1018.810858
Actual update 46950 of 50000, Stored update 46950 of 50000
-2 * Loglike = 1012.827777
Actual update 47000 of 50000, Stored update 47000 of 50000
-2 * Loglike = 1012.915862
Actual update 47050 of 50000, Stored update 47050 of 50000
-2 * Loglike = 1009.873152
Actual update 47100 of 50000, Stored update 47100 of 50000
-2 * Loglike = 1017.736243
Actual update 47150 of 50000, Stored update 47150 of 50000
-2 * Loglike = 1015.345371
Actual update 47200 of 50000, Stored update 47200 of 50000
-2 * Loglike = 1021.933511
Actual update 47250 of 50000, Stored update 47250 of 50000
-2 * Loglike = 1019.088429
Actual update 47300 of 50000, Stored update 47300 of 50000
-2 * Loglike = 1022.587543
Actual update 47350 of 50000, Stored update 47350 of 50000
-2 * Loglike = 1022.362336
Actual update 47400 of 50000, Stored update 47400 of 50000
-2 * Loglike = 1018.435267
Actual update 47450 of 50000, Stored update 47450 of 50000
-2 * Loglike = 1011.434498
Actual update 47500 of 50000, Stored update 47500 of 50000
-2 * Loglike = 1018.836900
Actual update 47550 of 50000, Stored update 47550 of 50000
-2 * Loglike = 1015.627159
Actual update 47600 of 50000, Stored update 47600 of 50000
-2 * Loglike = 1026.712738
Actual update 47650 of 50000, Stored update 47650 of 50000
-2 * Loglike = 1020.561584
Actual update 47700 of 50000, Stored update 47700 of 50000
-2 * Loglike = 1011.664197
Actual update 47750 of 50000, Stored update 47750 of 50000
-2 * Loglike = 1010.175849
Actual update 47800 of 50000, Stored update 47800 of 50000
-2 * Loglike = 1015.636093
Actual update 47850 of 50000, Stored update 47850 of 50000
-2 * Loglike = 1009.957975
Actual update 47900 of 50000, Stored update 47900 of 50000
-2 * Loglike = 1014.648497
Actual update 47950 of 50000, Stored update 47950 of 50000
-2 * Loglike = 1008.800804
Actual update 48000 of 50000, Stored update 48000 of 50000
-2 * Loglike = 1008.559100
Actual update 48050 of 50000, Stored update 48050 of 50000
-2 * Loglike = 1017.132258
Actual update 48100 of 50000, Stored update 48100 of 50000
-2 * Loglike = 1011.514625
Actual update 48150 of 50000, Stored update 48150 of 50000
-2 * Loglike = 1012.955873
Actual update 48200 of 50000, Stored update 48200 of 50000
-2 * Loglike = 1013.574219
Actual update 48250 of 50000, Stored update 48250 of 50000
-2 * Loglike = 1007.334774
Actual update 48300 of 50000, Stored update 48300 of 50000
-2 * Loglike = 1011.256255
Actual update 48350 of 50000, Stored update 48350 of 50000
-2 * Loglike = 1016.131127
Actual update 48400 of 50000, Stored update 48400 of 50000
-2 * Loglike = 1013.767156
Actual update 48450 of 50000, Stored update 48450 of 50000
-2 * Loglike = 1012.470402
Actual update 48500 of 50000, Stored update 48500 of 50000
-2 * Loglike = 1015.086505
Actual update 48550 of 50000, Stored update 48550 of 50000
-2 * Loglike = 1010.612767
Actual update 48600 of 50000, Stored update 48600 of 50000
-2 * Loglike = 1016.878634
Actual update 48650 of 50000, Stored update 48650 of 50000
-2 * Loglike = 1012.380493
Actual update 48700 of 50000, Stored update 48700 of 50000
-2 * Loglike = 1016.985610
Actual update 48750 of 50000, Stored update 48750 of 50000
-2 * Loglike = 1013.807346
Actual update 48800 of 50000, Stored update 48800 of 50000
-2 * Loglike = 1012.994436
Actual update 48850 of 50000, Stored update 48850 of 50000
-2 * Loglike = 1012.485581
Actual update 48900 of 50000, Stored update 48900 of 50000
-2 * Loglike = 1007.228853
Actual update 48950 of 50000, Stored update 48950 of 50000
-2 * Loglike = 1009.238441
Actual update 49000 of 50000, Stored update 49000 of 50000
-2 * Loglike = 1015.673783
Actual update 49050 of 50000, Stored update 49050 of 50000
-2 * Loglike = 1013.867770
Actual update 49100 of 50000, Stored update 49100 of 50000
-2 * Loglike = 1014.033355
Actual update 49150 of 50000, Stored update 49150 of 50000
-2 * Loglike = 1009.060836
Actual update 49200 of 50000, Stored update 49200 of 50000
-2 * Loglike = 1028.496445
Actual update 49250 of 50000, Stored update 49250 of 50000
-2 * Loglike = 1012.580715
Actual update 49300 of 50000, Stored update 49300 of 50000
-2 * Loglike = 1015.614574
Actual update 49350 of 50000, Stored update 49350 of 50000
-2 * Loglike = 1013.408860
Actual update 49400 of 50000, Stored update 49400 of 50000
-2 * Loglike = 1016.566013
Actual update 49450 of 50000, Stored update 49450 of 50000
-2 * Loglike = 1011.656152
Actual update 49500 of 50000, Stored update 49500 of 50000
-2 * Loglike = 1016.559243
Actual update 49550 of 50000, Stored update 49550 of 50000
-2 * Loglike = 1013.854577
Actual update 49600 of 50000, Stored update 49600 of 50000
-2 * Loglike = 1015.676116
Actual update 49650 of 50000, Stored update 49650 of 50000
-2 * Loglike = 1015.930715
Actual update 49700 of 50000, Stored update 49700 of 50000
-2 * Loglike = 1011.734354
Actual update 49750 of 50000, Stored update 49750 of 50000
-2 * Loglike = 1022.297470
Actual update 49800 of 50000, Stored update 49800 of 50000
-2 * Loglike = 1018.200043
Actual update 49850 of 50000, Stored update 49850 of 50000
-2 * Loglike = 1012.411804
Actual update 49900 of 50000, Stored update 49900 of 50000
-2 * Loglike = 1020.945468
Actual update 49950 of 50000, Stored update 49950 of 50000
-2 * Loglike = 1013.121629
Actual update 50000 of 50000, Stored update 50000 of 50000
-2 * Loglike = 1016.110930
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =          50000
Average =          1015.5
S.D.    =          5.8399
S.E.M.  =        0.026117

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  multilevel model (Multivariate Normal) 
Estimation algorithm:  MCMC      Elapsed time : 3.5s 
Number of obs:  30 (from total 30)          Number of iter.: 50000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
1015.492   1006.066   9.426      1024.918   
--------------------------------------------------------------------------------------------------- 
The model formula:
c(y8, y15, y22, y29, y36) ~ 1 + (1 | rat)
Level 1: rat      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
                    Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]     ESS 
Intercept_y8    152.17942     2.48879  61.15          0  ***   147.26931   157.10175   54056 
Intercept_y15   201.78200     2.82935  71.32          0  ***   196.20846   207.39190   50000 
Intercept_y22   245.04770     3.37835  72.53          0  ***   238.38848   251.74274   50000 
Intercept_y29   289.52491     4.19432  69.03          0  ***   281.27384   297.85980   50000 
Intercept_y36   324.82494     4.27035  76.07          0  ***   316.37428   333.28894   50000 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the rat level: 
                                      Coef.   Std. Err.   [95% Cred.   Interval]   ESS 
var_Intercept_y8                  184.31374    59.43683    100.84944   334.80092   239 
cov_Intercept_y8_Intercept_y15    197.37939    65.54327    105.97165   364.16411   209 
var_Intercept_y15                 237.99477    77.86234    129.72295   436.48498   223 
cov_Intercept_y8_Intercept_y22    222.45205    77.19883    115.25896   420.62690   180 
cov_Intercept_y15_Intercept_y22   268.00309    91.00943    142.85305   501.32784   181 
var_Intercept_y22                 339.38254   113.14830    184.20810   633.02593   196 
cov_Intercept_y8_Intercept_y29    260.95261    94.45845    130.83844   507.71541   172 
cov_Intercept_y15_Intercept_y29   314.22007   110.97162    162.83473   605.80377   165 
cov_Intercept_y22_Intercept_y29   397.51232   136.76851    211.48635   753.85170   179 
var_Intercept_y29                 523.46746   175.75299    285.63715   975.19910   205 
cov_Intercept_y8_Intercept_y36    250.40602    92.92772    122.29886   491.81793   183 
cov_Intercept_y15_Intercept_y36   301.46519   109.19109    151.94555   589.57541   182 
cov_Intercept_y22_Intercept_y36   381.18483   134.04153    198.21371   732.98502   193 
cov_Intercept_y29_Intercept_y36   501.51490   170.68462    268.88069   952.66896   205 
var_Intercept_y36                 540.75002   177.81520    297.56802   999.60314   271 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> covM2 <- matrix(, 5, 5)
> colnames(covM2) <- rownames(covM2) <- c("cons.y8", "cons.y15", "cons.y22", "cons.y29", "cons.y36")
> covM2[upper.tri(covM2, diag = TRUE)] <- mymodel@RP
> # covM2[lower.tri(covM2)] <- t(covM2)[lower.tri(covM2)]
> round(t(covM2), 3)
         cons.y8 cons.y15 cons.y22 cons.y29 cons.y36
cons.y8  184.314       NA       NA       NA       NA
cons.y15 197.379  237.995       NA       NA       NA
cons.y22 222.452  268.003  339.383       NA       NA
cons.y29 260.953  314.220  397.512  523.467       NA
cons.y36 250.406  301.465  381.185  501.515   540.75
> round(cov2cor(t(covM2)), 3)
         cons.y8 cons.y15 cons.y22 cons.y29 cons.y36
cons.y8    1.000       NA       NA       NA       NA
cons.y15   0.942    1.000       NA       NA       NA
cons.y22   0.889    0.943    1.000       NA       NA
cons.y29   0.840    0.890    0.943    1.000       NA
cons.y36   0.793    0.840    0.890    0.943        1
> 
> # Chapter learning outcomes . . . . . . . . . . . . . . . . . . . . . . .128
> 
> 
> 
> 
> 
> ############################################################################
> 
> proc.time()
   user  system elapsed 
   9.59    1.56   23.93