R version 4.3.2 (2023-10-31 ucrt) -- "Eye Holes"
Copyright (C) 2023 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-2024 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.09/
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
RUNNING ADAPTIVE PROCEDURE AND BURNING IN...
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 12693.142081
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 12633.435052
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 12626.287424
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 12711.568924
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 12691.888707
Adapting finished and took 600 iterations
Adapting took 600 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 12706.820313
Burning in for 100 iterations out of 500
-2 * Loglike = 12680.003708
Burning in for 150 iterations out of 500
-2 * Loglike = 12754.493736
Burning in for 200 iterations out of 500
-2 * Loglike = 12718.890569
Burning in for 250 iterations out of 500
-2 * Loglike = 12736.891399
Burning in for 300 iterations out of 500
-2 * Loglike = 12663.804243
Burning in for 350 iterations out of 500
-2 * Loglike = 12632.308989
Burning in for 400 iterations out of 500
-2 * Loglike = 12635.602976
Burning in for 450 iterations out of 500
-2 * Loglike = 12671.092005
Burning in for 500 iterations out of 500
-2 * Loglike = 12740.891388

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 = 12807.791975
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 12641.661566
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 12718.713584
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 12715.324608
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 12680.595855
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 12677.733868
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 12695.519896
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 12743.622529
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 12718.062384
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 12665.445471
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 12781.619363
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 12683.858425
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 12688.692164
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 12631.853448
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 12658.904941
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 12668.259689
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 12636.797785
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 12708.512875
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 12667.228870
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 12687.499871
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 12640.034386
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 12645.339793
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 12744.428788
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 12703.238361
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 12704.164123
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 12620.568286
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 12711.893573
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 12648.088361
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 12695.702446
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 12764.239416
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 12675.957253
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 12682.495876
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 12727.032363
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 12708.676448
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 12677.362947
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 12766.828212
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 12694.700132
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 12639.169744
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 12791.219243
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 12641.909301
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 12711.189131
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 12676.165963
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 12710.767976
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 12661.804844
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 12697.409744
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 12732.211293
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 12709.786976
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 12650.650718
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 12699.336069
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 12804.009949
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 12719.314396
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 12728.369301
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 12660.104862
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 12678.043912
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 12763.983726
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 12731.328537
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 12649.566710
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 12625.613313
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 12734.831035
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 12676.376062
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 12705.316787
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 12676.704256
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 12670.752686
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 12733.095099
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 12661.996692
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 12709.952387
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 12745.800296
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 12649.630027
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 12789.117110
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 12722.070779
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 12764.778102
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 12701.145819
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 12585.218159
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 12757.649844
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 12618.066050
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 12748.897366
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 12682.125963
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 12724.444280
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 12612.248601
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 12742.397192
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 12558.128414
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 12801.000705
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 12632.308041
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 12665.810970
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 12772.048231
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 12697.634394
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 12673.779381
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 12692.138168
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 12653.712942
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 12743.388607
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 12631.854068
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 12732.461995
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 12753.118908
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 12686.854974
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 12790.498686
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 12722.348651
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 12734.608580
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 12736.613513
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 12783.885027
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 12667.758698
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          12696.
S.D.    =          54.126
S.E.M.  =         0.76546

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  multilevel model (Mixed) 
Estimation algorithm:  MCMC      Elapsed time : 3.48s 
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.97689     1.09830  -6.35   2.119e-10  ***    -9.17039    -4.86454   5000 
sexmale_behaviour    -0.40019     0.08605  -4.65   3.308e-06  ***    -0.56150    -0.22784   1405 
ravens_english        1.75430     0.09099  19.28   7.831e-83  ***     1.57991     1.92902   5000 
ravens_behaviour      0.05643     0.00736   7.67   1.686e-14  ***     0.04173     0.07035   1489 
fluent_1              6.74508     1.31762   5.12   3.069e-07  ***     4.20866     9.33295   4451 
Intercept_1         -11.75602     3.18476  -3.69   0.0002231  ***   -18.00297    -5.44557   5000 
Intercept_2          -0.38972     0.18250  -2.14     0.03272  *      -0.73348    -0.03014   1782 
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.29710     0.73252      4.92041     7.74367    336 
var_Intercept_1           330.94600    14.23168    304.24223   360.59835    766 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 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
RUNNING ADAPTIVE PROCEDURE AND BURNING IN...
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 12536.149788
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 12585.717124
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 12444.391777
Adapting finished and took 400 iterations
Adapting took 400 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 12579.176377
Burning in for 100 iterations out of 500
-2 * Loglike = 12515.131882
Burning in for 150 iterations out of 500
-2 * Loglike = 12614.607752
Burning in for 200 iterations out of 500
-2 * Loglike = 12568.693381
Burning in for 250 iterations out of 500
-2 * Loglike = 12559.752140
Burning in for 300 iterations out of 500
-2 * Loglike = 12527.326924
Burning in for 350 iterations out of 500
-2 * Loglike = 12484.072449
Burning in for 400 iterations out of 500
-2 * Loglike = 12435.580925
Burning in for 450 iterations out of 500
-2 * Loglike = 12475.570542
Burning in for 500 iterations out of 500
-2 * Loglike = 12535.032570

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 = 12597.843704
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 12457.186261
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 12536.387587
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 12548.616782
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 12527.964332
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 12568.974422
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 12454.888900
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 12522.023997
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 12577.549771
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 12527.106968
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 12515.163443
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 12546.165344
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 12586.919012
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 12587.113640
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 12499.474443
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 12459.947737
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 12509.725762
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 12583.676748
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 12631.213135
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 12581.829947
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 12547.670681
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 12511.365531
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 12599.246844
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 12444.543881
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 12515.845942
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 12617.981897
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 12492.843713
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 12674.902359
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 12501.247046
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 12510.327603
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 12565.566224
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 12714.733601
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 12581.643906
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 12571.585204
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 12557.771879
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 12532.084909
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 12471.369978
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 12528.455471
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 12582.441419
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 12643.449847
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 12568.796011
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 12575.083616
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 12566.911963
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 12507.069972
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 12469.157205
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 12581.803732
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 12446.506281
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 12590.779143
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 12570.539315
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 12631.619064
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 12579.010521
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 12689.975253
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 12496.906558
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 12516.566176
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 12575.320996
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 12529.789908
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 12535.470033
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 12561.352425
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 12671.230786
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 12582.579881
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 12533.147470
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 12567.410857
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 12535.664330
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 12665.307915
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 12472.272621
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 12394.924296
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 12534.587981
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 12570.372898
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 12595.534320
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 12538.276503
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 12520.938389
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 12513.768691
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 12481.694870
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 12472.266231
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 12538.775647
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 12592.002896
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 12610.406139
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 12575.957195
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 12505.315855
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 12599.470864
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 12470.124532
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 12587.448029
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 12530.642285
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 12466.136063
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 12575.855786
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 12497.280694
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 12468.822350
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 12547.808457
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 12480.278603
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 12477.648313
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 12591.087538
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 12527.113687
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 12505.882637
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 12564.345371
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 12471.043401
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 12545.962210
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 12513.499562
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 12546.274233
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 12604.814956
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 12580.475354
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          12542.
S.D.    =          56.731
S.E.M.  =         0.80230

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  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.77s 
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.25243     1.04997  -5.95   2.603e-09  ***    -8.30831    -4.17068   4674 
sexmale_behaviour   -0.41671     0.09151  -4.55   5.274e-06  ***    -0.59736    -0.23591   1086 
ravens_english       1.65817     0.09137  18.15   1.313e-73  ***     1.47724     1.84114   3706 
ravens_behaviour     0.05620     0.00759   7.41   1.268e-13  ***     0.04139     0.07098   1292 
fluent_1             6.37506     1.27140   5.01   5.325e-07  ***     3.86689     8.89825   3749 
Intercept_1         -9.11965     3.30985  -2.76    0.005864  **    -15.48075    -2.69822   3040 
Intercept_2         -0.33781     0.19233  -1.76     0.07902  .      -0.71015     0.03909   1316 
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.53026    11.40708     22.79914    66.37680   1767 
cov_Intercept_1_Intercept_2    0.05692     0.42241     -0.74792     0.93309    647 
var_Intercept_2                0.06424     0.02931      0.02268     0.13419    451 
--------------------------------------------------------------------------------------------------- 
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.17684     0.74534      4.67787     7.59995    315 
var_Intercept_1           291.67313    12.91537    266.64852   319.77538    637 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 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...
Burning in for 50 iterations out of 500
-2 * Loglike = 1012.661040
Burning in for 100 iterations out of 500
-2 * Loglike = 1014.011632
Burning in for 150 iterations out of 500
-2 * Loglike = 1005.748207
Burning in for 200 iterations out of 500
-2 * Loglike = 1005.517134
Burning in for 250 iterations out of 500
-2 * Loglike = 1018.903915
Burning in for 300 iterations out of 500
-2 * Loglike = 1004.354927
Burning in for 350 iterations out of 500
-2 * Loglike = 1007.551571
Burning in for 400 iterations out of 500
-2 * Loglike = 1009.914258
Burning in for 450 iterations out of 500
-2 * Loglike = 1014.321248
Burning in for 500 iterations out of 500
-2 * Loglike = 1013.925423

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 = 1007.454750
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 1004.212412
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 1001.233755
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 1009.195712
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 1014.962130
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 1009.118809
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 1006.963438
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 1004.987457
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 1008.616276
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 1012.917743
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 1006.069992
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 1012.595147
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 1008.851846
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 1003.735234
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 1013.086054
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 1007.016049
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 1007.807141
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 1010.790419
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 1009.957117
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 1017.100269
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 1022.759183
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 1006.662868
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 1005.273793
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 1017.292523
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 1009.117765
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 1003.474246
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 1020.678615
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 1006.909563
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 1002.116169
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 1016.051596
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 1012.979385
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 1008.673159
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 1013.422769
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 1010.989110
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 999.080702
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 1006.276316
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 1020.435805
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 1006.087253
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 1025.540698
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 1006.513193
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 1014.622884
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 1012.733423
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 998.431612
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 1006.917249
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 1013.541946
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 1012.892504
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 1021.003381
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 1003.297699
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 1010.150906
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 1007.738569
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 1003.926733
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 998.536786
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 1003.194447
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 1010.113962
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 1024.328744
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 1012.574128
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 1008.569130
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 1010.415554
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 1001.502422
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 1004.547392
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 1005.786131
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 1000.451027
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 1013.746751
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 1011.172871
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 1015.763066
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 1003.384757
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 1011.833740
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 1007.068176
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 1018.982099
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 1006.693649
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 998.314613
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 1013.739038
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 1014.187737
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 1005.693900
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 1011.961285
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 1009.089766
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 1006.625336
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 1020.493370
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 1002.472077
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 1019.024436
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 1005.563823
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 1009.886860
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 999.580343
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 1007.157973
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 1007.747610
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 1011.431200
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 1014.293789
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 1004.763867
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 1007.668730
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 1004.240551
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 1008.826466
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 1013.227764
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 1007.053812
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 1011.279993
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 1010.476324
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 998.339476
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 1005.835322
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 1009.020815
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 999.091833
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 1007.955096
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          1008.7
S.D.    =          5.9352
S.E.M.  =        0.083936

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  multilevel model (Multivariate Normal) 
Estimation algorithm:  MCMC      Elapsed time : 0.36s 
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.744   993.331    15.413     1024.158   
--------------------------------------------------------------------------------------------------- 
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.17643     2.20463  69.03          0  ***   147.76010   156.44198   5000 
Intercept_y15   201.77279     2.49569  80.85          0  ***   196.91725   206.71000   4732 
Intercept_y22   245.04163     3.05538  80.20          0  ***   239.06703   251.04587   5000 
Intercept_y29   289.52110     3.82770  75.64          0  ***   282.06752   297.08368   5000 
Intercept_y36   324.79665     3.92172  82.82          0  ***   317.02291   332.47919   5000 
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.44783    41.02490     89.57631   246.19454   4759 
cov_Intercept_y8_Intercept_y15    156.56749    44.80222     90.96581   263.54542   4788 
var_Intercept_y15                 193.73990    53.02276    115.85750   320.58344   4939 
cov_Intercept_y8_Intercept_y22    171.02243    52.19967     94.71162   299.05885   4837 
cov_Intercept_y15_Intercept_y22   217.18174    62.30167    127.68802   365.71242   4983 
var_Intercept_y22                 286.17880    79.10952    173.88062   475.92816   5002 
cov_Intercept_y8_Intercept_y29    187.92012    62.36200     96.60231   339.43036   4823 
cov_Intercept_y15_Intercept_y29   247.85443    74.96997    139.41152   425.18577   4985 
cov_Intercept_y22_Intercept_y29   338.20100    96.27904    202.11554   561.68398   5057 
var_Intercept_y29                 444.56787   123.71311    271.56103   730.93164   4845 
cov_Intercept_y8_Intercept_y36    161.93346    59.95086     71.26897   304.33830   4934 
cov_Intercept_y15_Intercept_y36   209.72304    70.90976    106.04118   380.19339   4985 
cov_Intercept_y22_Intercept_y36   311.58165    93.18058    180.63540   535.38993   4988 
cov_Intercept_y29_Intercept_y36   415.67384   120.19575    246.44677   701.61143   4911 
var_Intercept_y36                 459.09514   126.67347    277.19021   761.02348   4898 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> 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.448       NA       NA       NA       NA
cons.y15 156.567  193.740       NA       NA       NA
cons.y22 171.022  217.182  286.179       NA       NA
cons.y29 187.920  247.854  338.201  444.568       NA
cons.y36 161.933  209.723  311.582  415.674  459.095
> 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.727    0.845    0.948     1.00       NA
cons.y36   0.616    0.703    0.860     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
RUNNING ADAPTIVE PROCEDURE AND BURNING IN...
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 1015.281601
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 1006.685972
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 1016.244204
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 1012.153839
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 1013.765262
Adapting for 600 iterations (Maximum 5000)
-2 * Loglike = 1012.095955
Adapting for 700 iterations (Maximum 5000)
-2 * Loglike = 1006.700647
Adapting for 800 iterations (Maximum 5000)
-2 * Loglike = 1011.952056
Adapting for 900 iterations (Maximum 5000)
-2 * Loglike = 1011.831544
Adapting for 1000 iterations (Maximum 5000)
-2 * Loglike = 1020.987309
Adapting finished and took 1100 iterations
Adapting took 1100 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 1009.419250
Burning in for 100 iterations out of 500
-2 * Loglike = 1016.907075
Burning in for 150 iterations out of 500
-2 * Loglike = 1018.645999
Burning in for 200 iterations out of 500
-2 * Loglike = 1010.513030
Burning in for 250 iterations out of 500
-2 * Loglike = 1013.499704
Burning in for 300 iterations out of 500
-2 * Loglike = 1018.167524
Burning in for 350 iterations out of 500
-2 * Loglike = 1017.586584
Burning in for 400 iterations out of 500
-2 * Loglike = 1020.316302
Burning in for 450 iterations out of 500
-2 * Loglike = 1029.491529
Burning in for 500 iterations out of 500
-2 * Loglike = 1023.550239

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 = 1023.017979
Actual update 100 of 50000, Stored update 100 of 50000
-2 * Loglike = 1011.405505
Actual update 150 of 50000, Stored update 150 of 50000
-2 * Loglike = 1025.551505
Actual update 200 of 50000, Stored update 200 of 50000
-2 * Loglike = 1014.207688
Actual update 250 of 50000, Stored update 250 of 50000
-2 * Loglike = 1014.299858
Actual update 300 of 50000, Stored update 300 of 50000
-2 * Loglike = 1013.430934
Actual update 350 of 50000, Stored update 350 of 50000
-2 * Loglike = 1018.533518
Actual update 400 of 50000, Stored update 400 of 50000
-2 * Loglike = 1010.569799
Actual update 450 of 50000, Stored update 450 of 50000
-2 * Loglike = 1018.231424
Actual update 500 of 50000, Stored update 500 of 50000
-2 * Loglike = 1014.435392
Actual update 550 of 50000, Stored update 550 of 50000
-2 * Loglike = 1009.624736
Actual update 600 of 50000, Stored update 600 of 50000
-2 * Loglike = 1009.904976
Actual update 650 of 50000, Stored update 650 of 50000
-2 * Loglike = 1007.881980
Actual update 700 of 50000, Stored update 700 of 50000
-2 * Loglike = 1013.648611
Actual update 750 of 50000, Stored update 750 of 50000
-2 * Loglike = 1013.726368
Actual update 800 of 50000, Stored update 800 of 50000
-2 * Loglike = 1013.581074
Actual update 850 of 50000, Stored update 850 of 50000
-2 * Loglike = 1013.550843
Actual update 900 of 50000, Stored update 900 of 50000
-2 * Loglike = 1009.328128
Actual update 950 of 50000, Stored update 950 of 50000
-2 * Loglike = 1014.972894
Actual update 1000 of 50000, Stored update 1000 of 50000
-2 * Loglike = 1007.838605
Actual update 1050 of 50000, Stored update 1050 of 50000
-2 * Loglike = 1032.372451
Actual update 1100 of 50000, Stored update 1100 of 50000
-2 * Loglike = 1019.954796
Actual update 1150 of 50000, Stored update 1150 of 50000
-2 * Loglike = 1007.320474
Actual update 1200 of 50000, Stored update 1200 of 50000
-2 * Loglike = 1016.306548
Actual update 1250 of 50000, Stored update 1250 of 50000
-2 * Loglike = 1014.906122
Actual update 1300 of 50000, Stored update 1300 of 50000
-2 * Loglike = 1012.641473
Actual update 1350 of 50000, Stored update 1350 of 50000
-2 * Loglike = 1014.162851
Actual update 1400 of 50000, Stored update 1400 of 50000
-2 * Loglike = 1017.447613
Actual update 1450 of 50000, Stored update 1450 of 50000
-2 * Loglike = 1010.988470
Actual update 1500 of 50000, Stored update 1500 of 50000
-2 * Loglike = 1025.002481
Actual update 1550 of 50000, Stored update 1550 of 50000
-2 * Loglike = 1010.359791
Actual update 1600 of 50000, Stored update 1600 of 50000
-2 * Loglike = 1017.300048
Actual update 1650 of 50000, Stored update 1650 of 50000
-2 * Loglike = 1013.774648
Actual update 1700 of 50000, Stored update 1700 of 50000
-2 * Loglike = 1011.357388
Actual update 1750 of 50000, Stored update 1750 of 50000
-2 * Loglike = 1009.323816
Actual update 1800 of 50000, Stored update 1800 of 50000
-2 * Loglike = 1008.830230
Actual update 1850 of 50000, Stored update 1850 of 50000
-2 * Loglike = 1018.730288
Actual update 1900 of 50000, Stored update 1900 of 50000
-2 * Loglike = 1008.799317
Actual update 1950 of 50000, Stored update 1950 of 50000
-2 * Loglike = 1016.317208
Actual update 2000 of 50000, Stored update 2000 of 50000
-2 * Loglike = 1020.547804
Actual update 2050 of 50000, Stored update 2050 of 50000
-2 * Loglike = 1014.279919
Actual update 2100 of 50000, Stored update 2100 of 50000
-2 * Loglike = 1013.214540
Actual update 2150 of 50000, Stored update 2150 of 50000
-2 * Loglike = 1013.145049
Actual update 2200 of 50000, Stored update 2200 of 50000
-2 * Loglike = 1014.673402
Actual update 2250 of 50000, Stored update 2250 of 50000
-2 * Loglike = 1010.454817
Actual update 2300 of 50000, Stored update 2300 of 50000
-2 * Loglike = 1010.053827
Actual update 2350 of 50000, Stored update 2350 of 50000
-2 * Loglike = 1009.219991
Actual update 2400 of 50000, Stored update 2400 of 50000
-2 * Loglike = 1009.731448
Actual update 2450 of 50000, Stored update 2450 of 50000
-2 * Loglike = 1007.589706
Actual update 2500 of 50000, Stored update 2500 of 50000
-2 * Loglike = 1012.943066
Actual update 2550 of 50000, Stored update 2550 of 50000
-2 * Loglike = 1010.694472
Actual update 2600 of 50000, Stored update 2600 of 50000
-2 * Loglike = 1011.854459
Actual update 2650 of 50000, Stored update 2650 of 50000
-2 * Loglike = 1018.578500
Actual update 2700 of 50000, Stored update 2700 of 50000
-2 * Loglike = 1018.894767
Actual update 2750 of 50000, Stored update 2750 of 50000
-2 * Loglike = 1029.315143
Actual update 2800 of 50000, Stored update 2800 of 50000
-2 * Loglike = 1015.342452
Actual update 2850 of 50000, Stored update 2850 of 50000
-2 * Loglike = 1009.665000
Actual update 2900 of 50000, Stored update 2900 of 50000
-2 * Loglike = 1014.098452
Actual update 2950 of 50000, Stored update 2950 of 50000
-2 * Loglike = 1015.744375
Actual update 3000 of 50000, Stored update 3000 of 50000
-2 * Loglike = 1019.845718
Actual update 3050 of 50000, Stored update 3050 of 50000
-2 * Loglike = 1012.434663
Actual update 3100 of 50000, Stored update 3100 of 50000
-2 * Loglike = 1019.341844
Actual update 3150 of 50000, Stored update 3150 of 50000
-2 * Loglike = 1015.224016
Actual update 3200 of 50000, Stored update 3200 of 50000
-2 * Loglike = 1012.410147
Actual update 3250 of 50000, Stored update 3250 of 50000
-2 * Loglike = 1010.255334
Actual update 3300 of 50000, Stored update 3300 of 50000
-2 * Loglike = 1016.854196
Actual update 3350 of 50000, Stored update 3350 of 50000
-2 * Loglike = 1021.091161
Actual update 3400 of 50000, Stored update 3400 of 50000
-2 * Loglike = 1008.940692
Actual update 3450 of 50000, Stored update 3450 of 50000
-2 * Loglike = 1011.761772
Actual update 3500 of 50000, Stored update 3500 of 50000
-2 * Loglike = 1021.377219
Actual update 3550 of 50000, Stored update 3550 of 50000
-2 * Loglike = 1016.027490
Actual update 3600 of 50000, Stored update 3600 of 50000
-2 * Loglike = 1018.301234
Actual update 3650 of 50000, Stored update 3650 of 50000
-2 * Loglike = 1025.897432
Actual update 3700 of 50000, Stored update 3700 of 50000
-2 * Loglike = 1012.987783
Actual update 3750 of 50000, Stored update 3750 of 50000
-2 * Loglike = 1010.271775
Actual update 3800 of 50000, Stored update 3800 of 50000
-2 * Loglike = 1015.465728
Actual update 3850 of 50000, Stored update 3850 of 50000
-2 * Loglike = 1011.749796
Actual update 3900 of 50000, Stored update 3900 of 50000
-2 * Loglike = 1007.495827
Actual update 3950 of 50000, Stored update 3950 of 50000
-2 * Loglike = 1010.326046
Actual update 4000 of 50000, Stored update 4000 of 50000
-2 * Loglike = 1021.361475
Actual update 4050 of 50000, Stored update 4050 of 50000
-2 * Loglike = 1014.167444
Actual update 4100 of 50000, Stored update 4100 of 50000
-2 * Loglike = 1018.993019
Actual update 4150 of 50000, Stored update 4150 of 50000
-2 * Loglike = 1019.165654
Actual update 4200 of 50000, Stored update 4200 of 50000
-2 * Loglike = 1013.822779
Actual update 4250 of 50000, Stored update 4250 of 50000
-2 * Loglike = 1016.693355
Actual update 4300 of 50000, Stored update 4300 of 50000
-2 * Loglike = 1010.683507
Actual update 4350 of 50000, Stored update 4350 of 50000
-2 * Loglike = 1023.717417
Actual update 4400 of 50000, Stored update 4400 of 50000
-2 * Loglike = 1009.211183
Actual update 4450 of 50000, Stored update 4450 of 50000
-2 * Loglike = 1016.241309
Actual update 4500 of 50000, Stored update 4500 of 50000
-2 * Loglike = 1018.898505
Actual update 4550 of 50000, Stored update 4550 of 50000
-2 * Loglike = 1010.556895
Actual update 4600 of 50000, Stored update 4600 of 50000
-2 * Loglike = 1010.748989
Actual update 4650 of 50000, Stored update 4650 of 50000
-2 * Loglike = 1011.905371
Actual update 4700 of 50000, Stored update 4700 of 50000
-2 * Loglike = 1012.108001
Actual update 4750 of 50000, Stored update 4750 of 50000
-2 * Loglike = 1019.151655
Actual update 4800 of 50000, Stored update 4800 of 50000
-2 * Loglike = 1004.957745
Actual update 4850 of 50000, Stored update 4850 of 50000
-2 * Loglike = 1013.728555
Actual update 4900 of 50000, Stored update 4900 of 50000
-2 * Loglike = 1013.104159
Actual update 4950 of 50000, Stored update 4950 of 50000
-2 * Loglike = 1031.642758
Actual update 5000 of 50000, Stored update 5000 of 50000
-2 * Loglike = 1021.036845
Actual update 5050 of 50000, Stored update 5050 of 50000
-2 * Loglike = 1018.340942
Actual update 5100 of 50000, Stored update 5100 of 50000
-2 * Loglike = 1016.082404
Actual update 5150 of 50000, Stored update 5150 of 50000
-2 * Loglike = 1017.450996
Actual update 5200 of 50000, Stored update 5200 of 50000
-2 * Loglike = 1006.799921
Actual update 5250 of 50000, Stored update 5250 of 50000
-2 * Loglike = 1003.501577
Actual update 5300 of 50000, Stored update 5300 of 50000
-2 * Loglike = 1020.805505
Actual update 5350 of 50000, Stored update 5350 of 50000
-2 * Loglike = 1010.782525
Actual update 5400 of 50000, Stored update 5400 of 50000
-2 * Loglike = 1028.301373
Actual update 5450 of 50000, Stored update 5450 of 50000
-2 * Loglike = 1026.593729
Actual update 5500 of 50000, Stored update 5500 of 50000
-2 * Loglike = 1024.979449
Actual update 5550 of 50000, Stored update 5550 of 50000
-2 * Loglike = 1030.376446
Actual update 5600 of 50000, Stored update 5600 of 50000
-2 * Loglike = 1029.210005
Actual update 5650 of 50000, Stored update 5650 of 50000
-2 * Loglike = 1014.771338
Actual update 5700 of 50000, Stored update 5700 of 50000
-2 * Loglike = 1022.132661
Actual update 5750 of 50000, Stored update 5750 of 50000
-2 * Loglike = 1015.744661
Actual update 5800 of 50000, Stored update 5800 of 50000
-2 * Loglike = 1013.424261
Actual update 5850 of 50000, Stored update 5850 of 50000
-2 * Loglike = 1005.856308
Actual update 5900 of 50000, Stored update 5900 of 50000
-2 * Loglike = 1016.125989
Actual update 5950 of 50000, Stored update 5950 of 50000
-2 * Loglike = 1013.578927
Actual update 6000 of 50000, Stored update 6000 of 50000
-2 * Loglike = 1013.483425
Actual update 6050 of 50000, Stored update 6050 of 50000
-2 * Loglike = 1014.326422
Actual update 6100 of 50000, Stored update 6100 of 50000
-2 * Loglike = 1008.915200
Actual update 6150 of 50000, Stored update 6150 of 50000
-2 * Loglike = 1011.301687
Actual update 6200 of 50000, Stored update 6200 of 50000
-2 * Loglike = 1014.566714
Actual update 6250 of 50000, Stored update 6250 of 50000
-2 * Loglike = 1008.962502
Actual update 6300 of 50000, Stored update 6300 of 50000
-2 * Loglike = 1018.211058
Actual update 6350 of 50000, Stored update 6350 of 50000
-2 * Loglike = 1009.010984
Actual update 6400 of 50000, Stored update 6400 of 50000
-2 * Loglike = 1017.209165
Actual update 6450 of 50000, Stored update 6450 of 50000
-2 * Loglike = 1017.227923
Actual update 6500 of 50000, Stored update 6500 of 50000
-2 * Loglike = 1012.291742
Actual update 6550 of 50000, Stored update 6550 of 50000
-2 * Loglike = 1014.000802
Actual update 6600 of 50000, Stored update 6600 of 50000
-2 * Loglike = 1020.953722
Actual update 6650 of 50000, Stored update 6650 of 50000
-2 * Loglike = 1014.885246
Actual update 6700 of 50000, Stored update 6700 of 50000
-2 * Loglike = 1012.854997
Actual update 6750 of 50000, Stored update 6750 of 50000
-2 * Loglike = 1010.460369
Actual update 6800 of 50000, Stored update 6800 of 50000
-2 * Loglike = 1017.019503
Actual update 6850 of 50000, Stored update 6850 of 50000
-2 * Loglike = 1020.658404
Actual update 6900 of 50000, Stored update 6900 of 50000
-2 * Loglike = 1019.994900
Actual update 6950 of 50000, Stored update 6950 of 50000
-2 * Loglike = 1017.169788
Actual update 7000 of 50000, Stored update 7000 of 50000
-2 * Loglike = 1010.100263
Actual update 7050 of 50000, Stored update 7050 of 50000
-2 * Loglike = 1006.632735
Actual update 7100 of 50000, Stored update 7100 of 50000
-2 * Loglike = 1008.333765
Actual update 7150 of 50000, Stored update 7150 of 50000
-2 * Loglike = 1012.183680
Actual update 7200 of 50000, Stored update 7200 of 50000
-2 * Loglike = 1014.348219
Actual update 7250 of 50000, Stored update 7250 of 50000
-2 * Loglike = 1015.331828
Actual update 7300 of 50000, Stored update 7300 of 50000
-2 * Loglike = 1022.008010
Actual update 7350 of 50000, Stored update 7350 of 50000
-2 * Loglike = 1020.572149
Actual update 7400 of 50000, Stored update 7400 of 50000
-2 * Loglike = 1007.515463
Actual update 7450 of 50000, Stored update 7450 of 50000
-2 * Loglike = 1009.790463
Actual update 7500 of 50000, Stored update 7500 of 50000
-2 * Loglike = 1009.733920
Actual update 7550 of 50000, Stored update 7550 of 50000
-2 * Loglike = 1016.548845
Actual update 7600 of 50000, Stored update 7600 of 50000
-2 * Loglike = 1019.690678
Actual update 7650 of 50000, Stored update 7650 of 50000
-2 * Loglike = 1013.176140
Actual update 7700 of 50000, Stored update 7700 of 50000
-2 * Loglike = 1013.177816
Actual update 7750 of 50000, Stored update 7750 of 50000
-2 * Loglike = 1022.324018
Actual update 7800 of 50000, Stored update 7800 of 50000
-2 * Loglike = 1018.889491
Actual update 7850 of 50000, Stored update 7850 of 50000
-2 * Loglike = 1013.331551
Actual update 7900 of 50000, Stored update 7900 of 50000
-2 * Loglike = 1014.257197
Actual update 7950 of 50000, Stored update 7950 of 50000
-2 * Loglike = 1010.443750
Actual update 8000 of 50000, Stored update 8000 of 50000
-2 * Loglike = 1019.756750
Actual update 8050 of 50000, Stored update 8050 of 50000
-2 * Loglike = 1023.733017
Actual update 8100 of 50000, Stored update 8100 of 50000
-2 * Loglike = 1013.940806
Actual update 8150 of 50000, Stored update 8150 of 50000
-2 * Loglike = 1025.855444
Actual update 8200 of 50000, Stored update 8200 of 50000
-2 * Loglike = 1009.752750
Actual update 8250 of 50000, Stored update 8250 of 50000
-2 * Loglike = 1008.702513
Actual update 8300 of 50000, Stored update 8300 of 50000
-2 * Loglike = 1007.264576
Actual update 8350 of 50000, Stored update 8350 of 50000
-2 * Loglike = 1011.286483
Actual update 8400 of 50000, Stored update 8400 of 50000
-2 * Loglike = 1015.798855
Actual update 8450 of 50000, Stored update 8450 of 50000
-2 * Loglike = 1015.861931
Actual update 8500 of 50000, Stored update 8500 of 50000
-2 * Loglike = 1015.615287
Actual update 8550 of 50000, Stored update 8550 of 50000
-2 * Loglike = 1013.897330
Actual update 8600 of 50000, Stored update 8600 of 50000
-2 * Loglike = 1019.669509
Actual update 8650 of 50000, Stored update 8650 of 50000
-2 * Loglike = 1015.876941
Actual update 8700 of 50000, Stored update 8700 of 50000
-2 * Loglike = 1021.470787
Actual update 8750 of 50000, Stored update 8750 of 50000
-2 * Loglike = 1017.497280
Actual update 8800 of 50000, Stored update 8800 of 50000
-2 * Loglike = 1015.356243
Actual update 8850 of 50000, Stored update 8850 of 50000
-2 * Loglike = 1016.480388
Actual update 8900 of 50000, Stored update 8900 of 50000
-2 * Loglike = 1011.280918
Actual update 8950 of 50000, Stored update 8950 of 50000
-2 * Loglike = 1017.439254
Actual update 9000 of 50000, Stored update 9000 of 50000
-2 * Loglike = 1010.044420
Actual update 9050 of 50000, Stored update 9050 of 50000
-2 * Loglike = 1010.947687
Actual update 9100 of 50000, Stored update 9100 of 50000
-2 * Loglike = 1011.588180
Actual update 9150 of 50000, Stored update 9150 of 50000
-2 * Loglike = 1020.931749
Actual update 9200 of 50000, Stored update 9200 of 50000
-2 * Loglike = 1017.118680
Actual update 9250 of 50000, Stored update 9250 of 50000
-2 * Loglike = 1015.577160
Actual update 9300 of 50000, Stored update 9300 of 50000
-2 * Loglike = 1016.202255
Actual update 9350 of 50000, Stored update 9350 of 50000
-2 * Loglike = 1014.839827
Actual update 9400 of 50000, Stored update 9400 of 50000
-2 * Loglike = 1009.764047
Actual update 9450 of 50000, Stored update 9450 of 50000
-2 * Loglike = 1013.385215
Actual update 9500 of 50000, Stored update 9500 of 50000
-2 * Loglike = 1031.644243
Actual update 9550 of 50000, Stored update 9550 of 50000
-2 * Loglike = 1025.877350
Actual update 9600 of 50000, Stored update 9600 of 50000
-2 * Loglike = 1038.321311
Actual update 9650 of 50000, Stored update 9650 of 50000
-2 * Loglike = 1024.265715
Actual update 9700 of 50000, Stored update 9700 of 50000
-2 * Loglike = 1020.261734
Actual update 9750 of 50000, Stored update 9750 of 50000
-2 * Loglike = 1036.486776
Actual update 9800 of 50000, Stored update 9800 of 50000
-2 * Loglike = 1037.600996
Actual update 9850 of 50000, Stored update 9850 of 50000
-2 * Loglike = 1014.656906
Actual update 9900 of 50000, Stored update 9900 of 50000
-2 * Loglike = 1016.380565
Actual update 9950 of 50000, Stored update 9950 of 50000
-2 * Loglike = 1016.416364
Actual update 10000 of 50000, Stored update 10000 of 50000
-2 * Loglike = 1019.463400
Actual update 10050 of 50000, Stored update 10050 of 50000
-2 * Loglike = 1019.762377
Actual update 10100 of 50000, Stored update 10100 of 50000
-2 * Loglike = 1016.638067
Actual update 10150 of 50000, Stored update 10150 of 50000
-2 * Loglike = 1020.441978
Actual update 10200 of 50000, Stored update 10200 of 50000
-2 * Loglike = 1015.430948
Actual update 10250 of 50000, Stored update 10250 of 50000
-2 * Loglike = 1011.006190
Actual update 10300 of 50000, Stored update 10300 of 50000
-2 * Loglike = 1010.924053
Actual update 10350 of 50000, Stored update 10350 of 50000
-2 * Loglike = 1020.233901
Actual update 10400 of 50000, Stored update 10400 of 50000
-2 * Loglike = 1010.527962
Actual update 10450 of 50000, Stored update 10450 of 50000
-2 * Loglike = 1012.576783
Actual update 10500 of 50000, Stored update 10500 of 50000
-2 * Loglike = 1016.559291
Actual update 10550 of 50000, Stored update 10550 of 50000
-2 * Loglike = 1010.739885
Actual update 10600 of 50000, Stored update 10600 of 50000
-2 * Loglike = 1015.173917
Actual update 10650 of 50000, Stored update 10650 of 50000
-2 * Loglike = 1012.174022
Actual update 10700 of 50000, Stored update 10700 of 50000
-2 * Loglike = 1006.156335
Actual update 10750 of 50000, Stored update 10750 of 50000
-2 * Loglike = 1019.351777
Actual update 10800 of 50000, Stored update 10800 of 50000
-2 * Loglike = 1010.305549
Actual update 10850 of 50000, Stored update 10850 of 50000
-2 * Loglike = 1011.716021
Actual update 10900 of 50000, Stored update 10900 of 50000
-2 * Loglike = 1035.803743
Actual update 10950 of 50000, Stored update 10950 of 50000
-2 * Loglike = 1010.733000
Actual update 11000 of 50000, Stored update 11000 of 50000
-2 * Loglike = 1012.306113
Actual update 11050 of 50000, Stored update 11050 of 50000
-2 * Loglike = 1009.760896
Actual update 11100 of 50000, Stored update 11100 of 50000
-2 * Loglike = 1021.255062
Actual update 11150 of 50000, Stored update 11150 of 50000
-2 * Loglike = 1009.539631
Actual update 11200 of 50000, Stored update 11200 of 50000
-2 * Loglike = 1020.259200
Actual update 11250 of 50000, Stored update 11250 of 50000
-2 * Loglike = 1019.215856
Actual update 11300 of 50000, Stored update 11300 of 50000
-2 * Loglike = 1015.959355
Actual update 11350 of 50000, Stored update 11350 of 50000
-2 * Loglike = 1023.857218
Actual update 11400 of 50000, Stored update 11400 of 50000
-2 * Loglike = 1025.197945
Actual update 11450 of 50000, Stored update 11450 of 50000
-2 * Loglike = 1019.246089
Actual update 11500 of 50000, Stored update 11500 of 50000
-2 * Loglike = 1032.633896
Actual update 11550 of 50000, Stored update 11550 of 50000
-2 * Loglike = 1012.828143
Actual update 11600 of 50000, Stored update 11600 of 50000
-2 * Loglike = 1019.601907
Actual update 11650 of 50000, Stored update 11650 of 50000
-2 * Loglike = 1011.429172
Actual update 11700 of 50000, Stored update 11700 of 50000
-2 * Loglike = 1020.703058
Actual update 11750 of 50000, Stored update 11750 of 50000
-2 * Loglike = 1019.139498
Actual update 11800 of 50000, Stored update 11800 of 50000
-2 * Loglike = 1019.786802
Actual update 11850 of 50000, Stored update 11850 of 50000
-2 * Loglike = 1017.896996
Actual update 11900 of 50000, Stored update 11900 of 50000
-2 * Loglike = 1011.019026
Actual update 11950 of 50000, Stored update 11950 of 50000
-2 * Loglike = 1010.728086
Actual update 12000 of 50000, Stored update 12000 of 50000
-2 * Loglike = 1013.457610
Actual update 12050 of 50000, Stored update 12050 of 50000
-2 * Loglike = 1010.454703
Actual update 12100 of 50000, Stored update 12100 of 50000
-2 * Loglike = 1014.271822
Actual update 12150 of 50000, Stored update 12150 of 50000
-2 * Loglike = 1020.457160
Actual update 12200 of 50000, Stored update 12200 of 50000
-2 * Loglike = 1010.294347
Actual update 12250 of 50000, Stored update 12250 of 50000
-2 * Loglike = 1016.004537
Actual update 12300 of 50000, Stored update 12300 of 50000
-2 * Loglike = 1024.194039
Actual update 12350 of 50000, Stored update 12350 of 50000
-2 * Loglike = 1025.417387
Actual update 12400 of 50000, Stored update 12400 of 50000
-2 * Loglike = 1021.021878
Actual update 12450 of 50000, Stored update 12450 of 50000
-2 * Loglike = 1018.249356
Actual update 12500 of 50000, Stored update 12500 of 50000
-2 * Loglike = 1016.896367
Actual update 12550 of 50000, Stored update 12550 of 50000
-2 * Loglike = 1021.870389
Actual update 12600 of 50000, Stored update 12600 of 50000
-2 * Loglike = 1021.664685
Actual update 12650 of 50000, Stored update 12650 of 50000
-2 * Loglike = 1013.659373
Actual update 12700 of 50000, Stored update 12700 of 50000
-2 * Loglike = 1015.937538
Actual update 12750 of 50000, Stored update 12750 of 50000
-2 * Loglike = 1015.868564
Actual update 12800 of 50000, Stored update 12800 of 50000
-2 * Loglike = 1015.808367
Actual update 12850 of 50000, Stored update 12850 of 50000
-2 * Loglike = 1013.533447
Actual update 12900 of 50000, Stored update 12900 of 50000
-2 * Loglike = 1016.999083
Actual update 12950 of 50000, Stored update 12950 of 50000
-2 * Loglike = 1018.719864
Actual update 13000 of 50000, Stored update 13000 of 50000
-2 * Loglike = 1007.938979
Actual update 13050 of 50000, Stored update 13050 of 50000
-2 * Loglike = 1018.360823
Actual update 13100 of 50000, Stored update 13100 of 50000
-2 * Loglike = 1022.042710
Actual update 13150 of 50000, Stored update 13150 of 50000
-2 * Loglike = 1030.466064
Actual update 13200 of 50000, Stored update 13200 of 50000
-2 * Loglike = 1031.544969
Actual update 13250 of 50000, Stored update 13250 of 50000
-2 * Loglike = 1035.210082
Actual update 13300 of 50000, Stored update 13300 of 50000
-2 * Loglike = 1020.880845
Actual update 13350 of 50000, Stored update 13350 of 50000
-2 * Loglike = 1041.956144
Actual update 13400 of 50000, Stored update 13400 of 50000
-2 * Loglike = 1021.304282
Actual update 13450 of 50000, Stored update 13450 of 50000
-2 * Loglike = 1022.337329
Actual update 13500 of 50000, Stored update 13500 of 50000
-2 * Loglike = 1020.473528
Actual update 13550 of 50000, Stored update 13550 of 50000
-2 * Loglike = 1028.368702
Actual update 13600 of 50000, Stored update 13600 of 50000
-2 * Loglike = 1013.481833
Actual update 13650 of 50000, Stored update 13650 of 50000
-2 * Loglike = 1019.017810
Actual update 13700 of 50000, Stored update 13700 of 50000
-2 * Loglike = 1012.081614
Actual update 13750 of 50000, Stored update 13750 of 50000
-2 * Loglike = 1015.224128
Actual update 13800 of 50000, Stored update 13800 of 50000
-2 * Loglike = 1009.019261
Actual update 13850 of 50000, Stored update 13850 of 50000
-2 * Loglike = 1014.333502
Actual update 13900 of 50000, Stored update 13900 of 50000
-2 * Loglike = 1017.308416
Actual update 13950 of 50000, Stored update 13950 of 50000
-2 * Loglike = 1011.325308
Actual update 14000 of 50000, Stored update 14000 of 50000
-2 * Loglike = 1010.784101
Actual update 14050 of 50000, Stored update 14050 of 50000
-2 * Loglike = 1010.597782
Actual update 14100 of 50000, Stored update 14100 of 50000
-2 * Loglike = 1014.932823
Actual update 14150 of 50000, Stored update 14150 of 50000
-2 * Loglike = 1011.972253
Actual update 14200 of 50000, Stored update 14200 of 50000
-2 * Loglike = 1018.413086
Actual update 14250 of 50000, Stored update 14250 of 50000
-2 * Loglike = 1012.066502
Actual update 14300 of 50000, Stored update 14300 of 50000
-2 * Loglike = 1015.326537
Actual update 14350 of 50000, Stored update 14350 of 50000
-2 * Loglike = 1016.413550
Actual update 14400 of 50000, Stored update 14400 of 50000
-2 * Loglike = 1009.011141
Actual update 14450 of 50000, Stored update 14450 of 50000
-2 * Loglike = 1016.541651
Actual update 14500 of 50000, Stored update 14500 of 50000
-2 * Loglike = 1012.283804
Actual update 14550 of 50000, Stored update 14550 of 50000
-2 * Loglike = 1015.605436
Actual update 14600 of 50000, Stored update 14600 of 50000
-2 * Loglike = 1014.041233
Actual update 14650 of 50000, Stored update 14650 of 50000
-2 * Loglike = 1021.796929
Actual update 14700 of 50000, Stored update 14700 of 50000
-2 * Loglike = 1014.481514
Actual update 14750 of 50000, Stored update 14750 of 50000
-2 * Loglike = 1010.836707
Actual update 14800 of 50000, Stored update 14800 of 50000
-2 * Loglike = 1020.825998
Actual update 14850 of 50000, Stored update 14850 of 50000
-2 * Loglike = 1006.995931
Actual update 14900 of 50000, Stored update 14900 of 50000
-2 * Loglike = 1018.466918
Actual update 14950 of 50000, Stored update 14950 of 50000
-2 * Loglike = 1017.393839
Actual update 15000 of 50000, Stored update 15000 of 50000
-2 * Loglike = 1009.580628
Actual update 15050 of 50000, Stored update 15050 of 50000
-2 * Loglike = 1016.124882
Actual update 15100 of 50000, Stored update 15100 of 50000
-2 * Loglike = 1016.313563
Actual update 15150 of 50000, Stored update 15150 of 50000
-2 * Loglike = 1014.914667
Actual update 15200 of 50000, Stored update 15200 of 50000
-2 * Loglike = 1018.751501
Actual update 15250 of 50000, Stored update 15250 of 50000
-2 * Loglike = 1017.500081
Actual update 15300 of 50000, Stored update 15300 of 50000
-2 * Loglike = 1016.608002
Actual update 15350 of 50000, Stored update 15350 of 50000
-2 * Loglike = 1024.986276
Actual update 15400 of 50000, Stored update 15400 of 50000
-2 * Loglike = 1019.008891
Actual update 15450 of 50000, Stored update 15450 of 50000
-2 * Loglike = 1027.301937
Actual update 15500 of 50000, Stored update 15500 of 50000
-2 * Loglike = 1029.406345
Actual update 15550 of 50000, Stored update 15550 of 50000
-2 * Loglike = 1021.635595
Actual update 15600 of 50000, Stored update 15600 of 50000
-2 * Loglike = 1013.168982
Actual update 15650 of 50000, Stored update 15650 of 50000
-2 * Loglike = 1015.021378
Actual update 15700 of 50000, Stored update 15700 of 50000
-2 * Loglike = 1022.964211
Actual update 15750 of 50000, Stored update 15750 of 50000
-2 * Loglike = 1007.489503
Actual update 15800 of 50000, Stored update 15800 of 50000
-2 * Loglike = 1018.835986
Actual update 15850 of 50000, Stored update 15850 of 50000
-2 * Loglike = 1013.093773
Actual update 15900 of 50000, Stored update 15900 of 50000
-2 * Loglike = 1005.214963
Actual update 15950 of 50000, Stored update 15950 of 50000
-2 * Loglike = 1013.764749
Actual update 16000 of 50000, Stored update 16000 of 50000
-2 * Loglike = 1023.883619
Actual update 16050 of 50000, Stored update 16050 of 50000
-2 * Loglike = 1010.223663
Actual update 16100 of 50000, Stored update 16100 of 50000
-2 * Loglike = 1008.480388
Actual update 16150 of 50000, Stored update 16150 of 50000
-2 * Loglike = 1016.240743
Actual update 16200 of 50000, Stored update 16200 of 50000
-2 * Loglike = 1013.689226
Actual update 16250 of 50000, Stored update 16250 of 50000
-2 * Loglike = 1021.936797
Actual update 16300 of 50000, Stored update 16300 of 50000
-2 * Loglike = 1011.654139
Actual update 16350 of 50000, Stored update 16350 of 50000
-2 * Loglike = 1013.624027
Actual update 16400 of 50000, Stored update 16400 of 50000
-2 * Loglike = 1022.347933
Actual update 16450 of 50000, Stored update 16450 of 50000
-2 * Loglike = 1008.239716
Actual update 16500 of 50000, Stored update 16500 of 50000
-2 * Loglike = 1012.635765
Actual update 16550 of 50000, Stored update 16550 of 50000
-2 * Loglike = 1014.155260
Actual update 16600 of 50000, Stored update 16600 of 50000
-2 * Loglike = 1015.046112
Actual update 16650 of 50000, Stored update 16650 of 50000
-2 * Loglike = 1008.492752
Actual update 16700 of 50000, Stored update 16700 of 50000
-2 * Loglike = 1015.236819
Actual update 16750 of 50000, Stored update 16750 of 50000
-2 * Loglike = 1009.712367
Actual update 16800 of 50000, Stored update 16800 of 50000
-2 * Loglike = 1012.396161
Actual update 16850 of 50000, Stored update 16850 of 50000
-2 * Loglike = 1016.554581
Actual update 16900 of 50000, Stored update 16900 of 50000
-2 * Loglike = 1012.187885
Actual update 16950 of 50000, Stored update 16950 of 50000
-2 * Loglike = 1012.177173
Actual update 17000 of 50000, Stored update 17000 of 50000
-2 * Loglike = 1006.192803
Actual update 17050 of 50000, Stored update 17050 of 50000
-2 * Loglike = 1016.568261
Actual update 17100 of 50000, Stored update 17100 of 50000
-2 * Loglike = 1011.628739
Actual update 17150 of 50000, Stored update 17150 of 50000
-2 * Loglike = 1016.936143
Actual update 17200 of 50000, Stored update 17200 of 50000
-2 * Loglike = 1010.220929
Actual update 17250 of 50000, Stored update 17250 of 50000
-2 * Loglike = 1018.977661
Actual update 17300 of 50000, Stored update 17300 of 50000
-2 * Loglike = 1007.678359
Actual update 17350 of 50000, Stored update 17350 of 50000
-2 * Loglike = 1009.453960
Actual update 17400 of 50000, Stored update 17400 of 50000
-2 * Loglike = 1016.844453
Actual update 17450 of 50000, Stored update 17450 of 50000
-2 * Loglike = 1011.789970
Actual update 17500 of 50000, Stored update 17500 of 50000
-2 * Loglike = 1010.389456
Actual update 17550 of 50000, Stored update 17550 of 50000
-2 * Loglike = 1012.023642
Actual update 17600 of 50000, Stored update 17600 of 50000
-2 * Loglike = 1017.445736
Actual update 17650 of 50000, Stored update 17650 of 50000
-2 * Loglike = 1016.219416
Actual update 17700 of 50000, Stored update 17700 of 50000
-2 * Loglike = 1021.602534
Actual update 17750 of 50000, Stored update 17750 of 50000
-2 * Loglike = 1007.586293
Actual update 17800 of 50000, Stored update 17800 of 50000
-2 * Loglike = 1006.882112
Actual update 17850 of 50000, Stored update 17850 of 50000
-2 * Loglike = 1016.134379
Actual update 17900 of 50000, Stored update 17900 of 50000
-2 * Loglike = 1022.028735
Actual update 17950 of 50000, Stored update 17950 of 50000
-2 * Loglike = 1013.624661
Actual update 18000 of 50000, Stored update 18000 of 50000
-2 * Loglike = 1008.887709
Actual update 18050 of 50000, Stored update 18050 of 50000
-2 * Loglike = 1012.345714
Actual update 18100 of 50000, Stored update 18100 of 50000
-2 * Loglike = 1012.625410
Actual update 18150 of 50000, Stored update 18150 of 50000
-2 * Loglike = 1012.581397
Actual update 18200 of 50000, Stored update 18200 of 50000
-2 * Loglike = 1017.492760
Actual update 18250 of 50000, Stored update 18250 of 50000
-2 * Loglike = 1021.451531
Actual update 18300 of 50000, Stored update 18300 of 50000
-2 * Loglike = 1019.530274
Actual update 18350 of 50000, Stored update 18350 of 50000
-2 * Loglike = 1017.346354
Actual update 18400 of 50000, Stored update 18400 of 50000
-2 * Loglike = 1018.865134
Actual update 18450 of 50000, Stored update 18450 of 50000
-2 * Loglike = 1011.516889
Actual update 18500 of 50000, Stored update 18500 of 50000
-2 * Loglike = 1015.789904
Actual update 18550 of 50000, Stored update 18550 of 50000
-2 * Loglike = 1012.657122
Actual update 18600 of 50000, Stored update 18600 of 50000
-2 * Loglike = 1021.131996
Actual update 18650 of 50000, Stored update 18650 of 50000
-2 * Loglike = 1010.291226
Actual update 18700 of 50000, Stored update 18700 of 50000
-2 * Loglike = 1012.225802
Actual update 18750 of 50000, Stored update 18750 of 50000
-2 * Loglike = 1018.569506
Actual update 18800 of 50000, Stored update 18800 of 50000
-2 * Loglike = 1012.088039
Actual update 18850 of 50000, Stored update 18850 of 50000
-2 * Loglike = 1017.904340
Actual update 18900 of 50000, Stored update 18900 of 50000
-2 * Loglike = 1016.350178
Actual update 18950 of 50000, Stored update 18950 of 50000
-2 * Loglike = 1011.614848
Actual update 19000 of 50000, Stored update 19000 of 50000
-2 * Loglike = 1008.327090
Actual update 19050 of 50000, Stored update 19050 of 50000
-2 * Loglike = 1012.347279
Actual update 19100 of 50000, Stored update 19100 of 50000
-2 * Loglike = 1014.097134
Actual update 19150 of 50000, Stored update 19150 of 50000
-2 * Loglike = 1019.775042
Actual update 19200 of 50000, Stored update 19200 of 50000
-2 * Loglike = 1007.828581
Actual update 19250 of 50000, Stored update 19250 of 50000
-2 * Loglike = 1009.592257
Actual update 19300 of 50000, Stored update 19300 of 50000
-2 * Loglike = 1013.883629
Actual update 19350 of 50000, Stored update 19350 of 50000
-2 * Loglike = 1016.039881
Actual update 19400 of 50000, Stored update 19400 of 50000
-2 * Loglike = 1010.220647
Actual update 19450 of 50000, Stored update 19450 of 50000
-2 * Loglike = 1012.569985
Actual update 19500 of 50000, Stored update 19500 of 50000
-2 * Loglike = 1016.282431
Actual update 19550 of 50000, Stored update 19550 of 50000
-2 * Loglike = 1017.707805
Actual update 19600 of 50000, Stored update 19600 of 50000
-2 * Loglike = 1010.852535
Actual update 19650 of 50000, Stored update 19650 of 50000
-2 * Loglike = 1019.742741
Actual update 19700 of 50000, Stored update 19700 of 50000
-2 * Loglike = 1013.957142
Actual update 19750 of 50000, Stored update 19750 of 50000
-2 * Loglike = 1016.175068
Actual update 19800 of 50000, Stored update 19800 of 50000
-2 * Loglike = 1006.571496
Actual update 19850 of 50000, Stored update 19850 of 50000
-2 * Loglike = 1018.482082
Actual update 19900 of 50000, Stored update 19900 of 50000
-2 * Loglike = 1015.258720
Actual update 19950 of 50000, Stored update 19950 of 50000
-2 * Loglike = 1016.414805
Actual update 20000 of 50000, Stored update 20000 of 50000
-2 * Loglike = 1012.654273
Actual update 20050 of 50000, Stored update 20050 of 50000
-2 * Loglike = 1017.400031
Actual update 20100 of 50000, Stored update 20100 of 50000
-2 * Loglike = 1009.953099
Actual update 20150 of 50000, Stored update 20150 of 50000
-2 * Loglike = 1010.887775
Actual update 20200 of 50000, Stored update 20200 of 50000
-2 * Loglike = 1019.522479
Actual update 20250 of 50000, Stored update 20250 of 50000
-2 * Loglike = 1013.284626
Actual update 20300 of 50000, Stored update 20300 of 50000
-2 * Loglike = 1010.450056
Actual update 20350 of 50000, Stored update 20350 of 50000
-2 * Loglike = 1017.072792
Actual update 20400 of 50000, Stored update 20400 of 50000
-2 * Loglike = 1016.085792
Actual update 20450 of 50000, Stored update 20450 of 50000
-2 * Loglike = 1008.117317
Actual update 20500 of 50000, Stored update 20500 of 50000
-2 * Loglike = 1020.257413
Actual update 20550 of 50000, Stored update 20550 of 50000
-2 * Loglike = 1009.432710
Actual update 20600 of 50000, Stored update 20600 of 50000
-2 * Loglike = 1017.336392
Actual update 20650 of 50000, Stored update 20650 of 50000
-2 * Loglike = 1011.492222
Actual update 20700 of 50000, Stored update 20700 of 50000
-2 * Loglike = 1021.838976
Actual update 20750 of 50000, Stored update 20750 of 50000
-2 * Loglike = 1009.851031
Actual update 20800 of 50000, Stored update 20800 of 50000
-2 * Loglike = 1014.858660
Actual update 20850 of 50000, Stored update 20850 of 50000
-2 * Loglike = 1018.964814
Actual update 20900 of 50000, Stored update 20900 of 50000
-2 * Loglike = 1012.476454
Actual update 20950 of 50000, Stored update 20950 of 50000
-2 * Loglike = 1012.907973
Actual update 21000 of 50000, Stored update 21000 of 50000
-2 * Loglike = 1010.653028
Actual update 21050 of 50000, Stored update 21050 of 50000
-2 * Loglike = 1009.902065
Actual update 21100 of 50000, Stored update 21100 of 50000
-2 * Loglike = 1011.162404
Actual update 21150 of 50000, Stored update 21150 of 50000
-2 * Loglike = 1006.788870
Actual update 21200 of 50000, Stored update 21200 of 50000
-2 * Loglike = 1013.842707
Actual update 21250 of 50000, Stored update 21250 of 50000
-2 * Loglike = 1023.882709
Actual update 21300 of 50000, Stored update 21300 of 50000
-2 * Loglike = 1019.241748
Actual update 21350 of 50000, Stored update 21350 of 50000
-2 * Loglike = 1020.112523
Actual update 21400 of 50000, Stored update 21400 of 50000
-2 * Loglike = 1038.029423
Actual update 21450 of 50000, Stored update 21450 of 50000
-2 * Loglike = 1017.387483
Actual update 21500 of 50000, Stored update 21500 of 50000
-2 * Loglike = 1012.815004
Actual update 21550 of 50000, Stored update 21550 of 50000
-2 * Loglike = 1015.672150
Actual update 21600 of 50000, Stored update 21600 of 50000
-2 * Loglike = 1008.312414
Actual update 21650 of 50000, Stored update 21650 of 50000
-2 * Loglike = 1010.838035
Actual update 21700 of 50000, Stored update 21700 of 50000
-2 * Loglike = 1011.921770
Actual update 21750 of 50000, Stored update 21750 of 50000
-2 * Loglike = 1016.646219
Actual update 21800 of 50000, Stored update 21800 of 50000
-2 * Loglike = 1013.532396
Actual update 21850 of 50000, Stored update 21850 of 50000
-2 * Loglike = 1016.770767
Actual update 21900 of 50000, Stored update 21900 of 50000
-2 * Loglike = 1015.394716
Actual update 21950 of 50000, Stored update 21950 of 50000
-2 * Loglike = 1010.430592
Actual update 22000 of 50000, Stored update 22000 of 50000
-2 * Loglike = 1008.196574
Actual update 22050 of 50000, Stored update 22050 of 50000
-2 * Loglike = 1010.315885
Actual update 22100 of 50000, Stored update 22100 of 50000
-2 * Loglike = 1007.848493
Actual update 22150 of 50000, Stored update 22150 of 50000
-2 * Loglike = 1009.773631
Actual update 22200 of 50000, Stored update 22200 of 50000
-2 * Loglike = 1019.037696
Actual update 22250 of 50000, Stored update 22250 of 50000
-2 * Loglike = 1004.259763
Actual update 22300 of 50000, Stored update 22300 of 50000
-2 * Loglike = 1017.942394
Actual update 22350 of 50000, Stored update 22350 of 50000
-2 * Loglike = 1011.305120
Actual update 22400 of 50000, Stored update 22400 of 50000
-2 * Loglike = 1039.508833
Actual update 22450 of 50000, Stored update 22450 of 50000
-2 * Loglike = 1018.198684
Actual update 22500 of 50000, Stored update 22500 of 50000
-2 * Loglike = 1012.285124
Actual update 22550 of 50000, Stored update 22550 of 50000
-2 * Loglike = 1019.700527
Actual update 22600 of 50000, Stored update 22600 of 50000
-2 * Loglike = 1010.957452
Actual update 22650 of 50000, Stored update 22650 of 50000
-2 * Loglike = 1017.281714
Actual update 22700 of 50000, Stored update 22700 of 50000
-2 * Loglike = 1018.221782
Actual update 22750 of 50000, Stored update 22750 of 50000
-2 * Loglike = 1028.292097
Actual update 22800 of 50000, Stored update 22800 of 50000
-2 * Loglike = 1028.656488
Actual update 22850 of 50000, Stored update 22850 of 50000
-2 * Loglike = 1026.209942
Actual update 22900 of 50000, Stored update 22900 of 50000
-2 * Loglike = 1018.498988
Actual update 22950 of 50000, Stored update 22950 of 50000
-2 * Loglike = 1016.274403
Actual update 23000 of 50000, Stored update 23000 of 50000
-2 * Loglike = 1017.934734
Actual update 23050 of 50000, Stored update 23050 of 50000
-2 * Loglike = 1010.699630
Actual update 23100 of 50000, Stored update 23100 of 50000
-2 * Loglike = 1018.979646
Actual update 23150 of 50000, Stored update 23150 of 50000
-2 * Loglike = 1010.482821
Actual update 23200 of 50000, Stored update 23200 of 50000
-2 * Loglike = 1014.858286
Actual update 23250 of 50000, Stored update 23250 of 50000
-2 * Loglike = 1010.816911
Actual update 23300 of 50000, Stored update 23300 of 50000
-2 * Loglike = 1011.139480
Actual update 23350 of 50000, Stored update 23350 of 50000
-2 * Loglike = 1020.160572
Actual update 23400 of 50000, Stored update 23400 of 50000
-2 * Loglike = 1025.359837
Actual update 23450 of 50000, Stored update 23450 of 50000
-2 * Loglike = 1018.556225
Actual update 23500 of 50000, Stored update 23500 of 50000
-2 * Loglike = 1012.255923
Actual update 23550 of 50000, Stored update 23550 of 50000
-2 * Loglike = 1014.427148
Actual update 23600 of 50000, Stored update 23600 of 50000
-2 * Loglike = 1009.062289
Actual update 23650 of 50000, Stored update 23650 of 50000
-2 * Loglike = 1010.851088
Actual update 23700 of 50000, Stored update 23700 of 50000
-2 * Loglike = 1013.441039
Actual update 23750 of 50000, Stored update 23750 of 50000
-2 * Loglike = 1009.227298
Actual update 23800 of 50000, Stored update 23800 of 50000
-2 * Loglike = 1010.534207
Actual update 23850 of 50000, Stored update 23850 of 50000
-2 * Loglike = 1006.654405
Actual update 23900 of 50000, Stored update 23900 of 50000
-2 * Loglike = 1015.165317
Actual update 23950 of 50000, Stored update 23950 of 50000
-2 * Loglike = 1013.255872
Actual update 24000 of 50000, Stored update 24000 of 50000
-2 * Loglike = 1020.139116
Actual update 24050 of 50000, Stored update 24050 of 50000
-2 * Loglike = 1011.155759
Actual update 24100 of 50000, Stored update 24100 of 50000
-2 * Loglike = 1023.783073
Actual update 24150 of 50000, Stored update 24150 of 50000
-2 * Loglike = 1013.796028
Actual update 24200 of 50000, Stored update 24200 of 50000
-2 * Loglike = 1018.448922
Actual update 24250 of 50000, Stored update 24250 of 50000
-2 * Loglike = 1013.876324
Actual update 24300 of 50000, Stored update 24300 of 50000
-2 * Loglike = 1014.046071
Actual update 24350 of 50000, Stored update 24350 of 50000
-2 * Loglike = 1019.398219
Actual update 24400 of 50000, Stored update 24400 of 50000
-2 * Loglike = 1026.895831
Actual update 24450 of 50000, Stored update 24450 of 50000
-2 * Loglike = 1012.355177
Actual update 24500 of 50000, Stored update 24500 of 50000
-2 * Loglike = 1016.131289
Actual update 24550 of 50000, Stored update 24550 of 50000
-2 * Loglike = 1012.436682
Actual update 24600 of 50000, Stored update 24600 of 50000
-2 * Loglike = 1010.493411
Actual update 24650 of 50000, Stored update 24650 of 50000
-2 * Loglike = 1011.643799
Actual update 24700 of 50000, Stored update 24700 of 50000
-2 * Loglike = 1009.203911
Actual update 24750 of 50000, Stored update 24750 of 50000
-2 * Loglike = 1011.646917
Actual update 24800 of 50000, Stored update 24800 of 50000
-2 * Loglike = 1015.431406
Actual update 24850 of 50000, Stored update 24850 of 50000
-2 * Loglike = 1013.936755
Actual update 24900 of 50000, Stored update 24900 of 50000
-2 * Loglike = 1011.982734
Actual update 24950 of 50000, Stored update 24950 of 50000
-2 * Loglike = 1010.948474
Actual update 25000 of 50000, Stored update 25000 of 50000
-2 * Loglike = 1009.009952
Actual update 25050 of 50000, Stored update 25050 of 50000
-2 * Loglike = 1013.136461
Actual update 25100 of 50000, Stored update 25100 of 50000
-2 * Loglike = 1016.732398
Actual update 25150 of 50000, Stored update 25150 of 50000
-2 * Loglike = 1012.213446
Actual update 25200 of 50000, Stored update 25200 of 50000
-2 * Loglike = 1032.365096
Actual update 25250 of 50000, Stored update 25250 of 50000
-2 * Loglike = 1022.182635
Actual update 25300 of 50000, Stored update 25300 of 50000
-2 * Loglike = 1022.767844
Actual update 25350 of 50000, Stored update 25350 of 50000
-2 * Loglike = 1014.175344
Actual update 25400 of 50000, Stored update 25400 of 50000
-2 * Loglike = 1017.053689
Actual update 25450 of 50000, Stored update 25450 of 50000
-2 * Loglike = 1019.958157
Actual update 25500 of 50000, Stored update 25500 of 50000
-2 * Loglike = 1013.650877
Actual update 25550 of 50000, Stored update 25550 of 50000
-2 * Loglike = 1011.417913
Actual update 25600 of 50000, Stored update 25600 of 50000
-2 * Loglike = 1022.682661
Actual update 25650 of 50000, Stored update 25650 of 50000
-2 * Loglike = 1020.528314
Actual update 25700 of 50000, Stored update 25700 of 50000
-2 * Loglike = 1016.229365
Actual update 25750 of 50000, Stored update 25750 of 50000
-2 * Loglike = 1017.059486
Actual update 25800 of 50000, Stored update 25800 of 50000
-2 * Loglike = 1012.561260
Actual update 25850 of 50000, Stored update 25850 of 50000
-2 * Loglike = 1009.138870
Actual update 25900 of 50000, Stored update 25900 of 50000
-2 * Loglike = 1010.611791
Actual update 25950 of 50000, Stored update 25950 of 50000
-2 * Loglike = 1012.694357
Actual update 26000 of 50000, Stored update 26000 of 50000
-2 * Loglike = 1009.596683
Actual update 26050 of 50000, Stored update 26050 of 50000
-2 * Loglike = 1007.960633
Actual update 26100 of 50000, Stored update 26100 of 50000
-2 * Loglike = 1010.697691
Actual update 26150 of 50000, Stored update 26150 of 50000
-2 * Loglike = 1025.989317
Actual update 26200 of 50000, Stored update 26200 of 50000
-2 * Loglike = 1012.456987
Actual update 26250 of 50000, Stored update 26250 of 50000
-2 * Loglike = 1012.335325
Actual update 26300 of 50000, Stored update 26300 of 50000
-2 * Loglike = 1009.639279
Actual update 26350 of 50000, Stored update 26350 of 50000
-2 * Loglike = 1011.356000
Actual update 26400 of 50000, Stored update 26400 of 50000
-2 * Loglike = 1009.473359
Actual update 26450 of 50000, Stored update 26450 of 50000
-2 * Loglike = 1009.992389
Actual update 26500 of 50000, Stored update 26500 of 50000
-2 * Loglike = 1009.024505
Actual update 26550 of 50000, Stored update 26550 of 50000
-2 * Loglike = 1019.090683
Actual update 26600 of 50000, Stored update 26600 of 50000
-2 * Loglike = 1012.034493
Actual update 26650 of 50000, Stored update 26650 of 50000
-2 * Loglike = 1020.382374
Actual update 26700 of 50000, Stored update 26700 of 50000
-2 * Loglike = 1008.376960
Actual update 26750 of 50000, Stored update 26750 of 50000
-2 * Loglike = 1011.747451
Actual update 26800 of 50000, Stored update 26800 of 50000
-2 * Loglike = 1018.421122
Actual update 26850 of 50000, Stored update 26850 of 50000
-2 * Loglike = 1012.321074
Actual update 26900 of 50000, Stored update 26900 of 50000
-2 * Loglike = 1018.502465
Actual update 26950 of 50000, Stored update 26950 of 50000
-2 * Loglike = 1009.032835
Actual update 27000 of 50000, Stored update 27000 of 50000
-2 * Loglike = 1008.854592
Actual update 27050 of 50000, Stored update 27050 of 50000
-2 * Loglike = 1007.556362
Actual update 27100 of 50000, Stored update 27100 of 50000
-2 * Loglike = 1015.106514
Actual update 27150 of 50000, Stored update 27150 of 50000
-2 * Loglike = 1011.311125
Actual update 27200 of 50000, Stored update 27200 of 50000
-2 * Loglike = 1013.050611
Actual update 27250 of 50000, Stored update 27250 of 50000
-2 * Loglike = 1016.884845
Actual update 27300 of 50000, Stored update 27300 of 50000
-2 * Loglike = 1012.028100
Actual update 27350 of 50000, Stored update 27350 of 50000
-2 * Loglike = 1014.690877
Actual update 27400 of 50000, Stored update 27400 of 50000
-2 * Loglike = 1011.121396
Actual update 27450 of 50000, Stored update 27450 of 50000
-2 * Loglike = 1013.158024
Actual update 27500 of 50000, Stored update 27500 of 50000
-2 * Loglike = 1012.907078
Actual update 27550 of 50000, Stored update 27550 of 50000
-2 * Loglike = 1010.318697
Actual update 27600 of 50000, Stored update 27600 of 50000
-2 * Loglike = 1010.922786
Actual update 27650 of 50000, Stored update 27650 of 50000
-2 * Loglike = 1024.845790
Actual update 27700 of 50000, Stored update 27700 of 50000
-2 * Loglike = 1011.106695
Actual update 27750 of 50000, Stored update 27750 of 50000
-2 * Loglike = 1015.091502
Actual update 27800 of 50000, Stored update 27800 of 50000
-2 * Loglike = 1008.886231
Actual update 27850 of 50000, Stored update 27850 of 50000
-2 * Loglike = 1010.078843
Actual update 27900 of 50000, Stored update 27900 of 50000
-2 * Loglike = 1011.118544
Actual update 27950 of 50000, Stored update 27950 of 50000
-2 * Loglike = 1015.696325
Actual update 28000 of 50000, Stored update 28000 of 50000
-2 * Loglike = 1013.462786
Actual update 28050 of 50000, Stored update 28050 of 50000
-2 * Loglike = 1015.872271
Actual update 28100 of 50000, Stored update 28100 of 50000
-2 * Loglike = 1014.351654
Actual update 28150 of 50000, Stored update 28150 of 50000
-2 * Loglike = 1019.282124
Actual update 28200 of 50000, Stored update 28200 of 50000
-2 * Loglike = 1017.146642
Actual update 28250 of 50000, Stored update 28250 of 50000
-2 * Loglike = 1012.808539
Actual update 28300 of 50000, Stored update 28300 of 50000
-2 * Loglike = 1016.767534
Actual update 28350 of 50000, Stored update 28350 of 50000
-2 * Loglike = 1014.398837
Actual update 28400 of 50000, Stored update 28400 of 50000
-2 * Loglike = 1016.253059
Actual update 28450 of 50000, Stored update 28450 of 50000
-2 * Loglike = 1014.784052
Actual update 28500 of 50000, Stored update 28500 of 50000
-2 * Loglike = 1013.001700
Actual update 28550 of 50000, Stored update 28550 of 50000
-2 * Loglike = 1016.786639
Actual update 28600 of 50000, Stored update 28600 of 50000
-2 * Loglike = 1011.791597
Actual update 28650 of 50000, Stored update 28650 of 50000
-2 * Loglike = 1007.439351
Actual update 28700 of 50000, Stored update 28700 of 50000
-2 * Loglike = 1016.439929
Actual update 28750 of 50000, Stored update 28750 of 50000
-2 * Loglike = 1018.399758
Actual update 28800 of 50000, Stored update 28800 of 50000
-2 * Loglike = 1013.839641
Actual update 28850 of 50000, Stored update 28850 of 50000
-2 * Loglike = 1011.812521
Actual update 28900 of 50000, Stored update 28900 of 50000
-2 * Loglike = 1008.408231
Actual update 28950 of 50000, Stored update 28950 of 50000
-2 * Loglike = 1006.035459
Actual update 29000 of 50000, Stored update 29000 of 50000
-2 * Loglike = 1013.999699
Actual update 29050 of 50000, Stored update 29050 of 50000
-2 * Loglike = 1020.890350
Actual update 29100 of 50000, Stored update 29100 of 50000
-2 * Loglike = 1018.005710
Actual update 29150 of 50000, Stored update 29150 of 50000
-2 * Loglike = 1006.944528
Actual update 29200 of 50000, Stored update 29200 of 50000
-2 * Loglike = 1029.935852
Actual update 29250 of 50000, Stored update 29250 of 50000
-2 * Loglike = 1016.175165
Actual update 29300 of 50000, Stored update 29300 of 50000
-2 * Loglike = 1018.862350
Actual update 29350 of 50000, Stored update 29350 of 50000
-2 * Loglike = 1006.851386
Actual update 29400 of 50000, Stored update 29400 of 50000
-2 * Loglike = 1013.370776
Actual update 29450 of 50000, Stored update 29450 of 50000
-2 * Loglike = 1007.837351
Actual update 29500 of 50000, Stored update 29500 of 50000
-2 * Loglike = 1008.808932
Actual update 29550 of 50000, Stored update 29550 of 50000
-2 * Loglike = 1015.940819
Actual update 29600 of 50000, Stored update 29600 of 50000
-2 * Loglike = 1016.026242
Actual update 29650 of 50000, Stored update 29650 of 50000
-2 * Loglike = 1013.039253
Actual update 29700 of 50000, Stored update 29700 of 50000
-2 * Loglike = 1011.413323
Actual update 29750 of 50000, Stored update 29750 of 50000
-2 * Loglike = 1013.885066
Actual update 29800 of 50000, Stored update 29800 of 50000
-2 * Loglike = 1016.252476
Actual update 29850 of 50000, Stored update 29850 of 50000
-2 * Loglike = 1015.352376
Actual update 29900 of 50000, Stored update 29900 of 50000
-2 * Loglike = 1023.607692
Actual update 29950 of 50000, Stored update 29950 of 50000
-2 * Loglike = 1012.680063
Actual update 30000 of 50000, Stored update 30000 of 50000
-2 * Loglike = 1011.506823
Actual update 30050 of 50000, Stored update 30050 of 50000
-2 * Loglike = 1014.483104
Actual update 30100 of 50000, Stored update 30100 of 50000
-2 * Loglike = 1017.950036
Actual update 30150 of 50000, Stored update 30150 of 50000
-2 * Loglike = 1014.665623
Actual update 30200 of 50000, Stored update 30200 of 50000
-2 * Loglike = 1019.371571
Actual update 30250 of 50000, Stored update 30250 of 50000
-2 * Loglike = 1012.194696
Actual update 30300 of 50000, Stored update 30300 of 50000
-2 * Loglike = 1011.844851
Actual update 30350 of 50000, Stored update 30350 of 50000
-2 * Loglike = 1015.147301
Actual update 30400 of 50000, Stored update 30400 of 50000
-2 * Loglike = 1015.229233
Actual update 30450 of 50000, Stored update 30450 of 50000
-2 * Loglike = 1013.825092
Actual update 30500 of 50000, Stored update 30500 of 50000
-2 * Loglike = 1015.221662
Actual update 30550 of 50000, Stored update 30550 of 50000
-2 * Loglike = 1009.497076
Actual update 30600 of 50000, Stored update 30600 of 50000
-2 * Loglike = 1013.313032
Actual update 30650 of 50000, Stored update 30650 of 50000
-2 * Loglike = 1018.670330
Actual update 30700 of 50000, Stored update 30700 of 50000
-2 * Loglike = 1017.045338
Actual update 30750 of 50000, Stored update 30750 of 50000
-2 * Loglike = 1008.593064
Actual update 30800 of 50000, Stored update 30800 of 50000
-2 * Loglike = 1010.829865
Actual update 30850 of 50000, Stored update 30850 of 50000
-2 * Loglike = 1010.465625
Actual update 30900 of 50000, Stored update 30900 of 50000
-2 * Loglike = 1013.073947
Actual update 30950 of 50000, Stored update 30950 of 50000
-2 * Loglike = 1017.124056
Actual update 31000 of 50000, Stored update 31000 of 50000
-2 * Loglike = 1007.805153
Actual update 31050 of 50000, Stored update 31050 of 50000
-2 * Loglike = 1035.996048
Actual update 31100 of 50000, Stored update 31100 of 50000
-2 * Loglike = 1011.017615
Actual update 31150 of 50000, Stored update 31150 of 50000
-2 * Loglike = 1017.343489
Actual update 31200 of 50000, Stored update 31200 of 50000
-2 * Loglike = 1024.031001
Actual update 31250 of 50000, Stored update 31250 of 50000
-2 * Loglike = 1014.609349
Actual update 31300 of 50000, Stored update 31300 of 50000
-2 * Loglike = 1015.205120
Actual update 31350 of 50000, Stored update 31350 of 50000
-2 * Loglike = 1012.071746
Actual update 31400 of 50000, Stored update 31400 of 50000
-2 * Loglike = 1009.994638
Actual update 31450 of 50000, Stored update 31450 of 50000
-2 * Loglike = 1028.765995
Actual update 31500 of 50000, Stored update 31500 of 50000
-2 * Loglike = 1024.742962
Actual update 31550 of 50000, Stored update 31550 of 50000
-2 * Loglike = 1009.606143
Actual update 31600 of 50000, Stored update 31600 of 50000
-2 * Loglike = 1011.312699
Actual update 31650 of 50000, Stored update 31650 of 50000
-2 * Loglike = 1013.446782
Actual update 31700 of 50000, Stored update 31700 of 50000
-2 * Loglike = 1012.984742
Actual update 31750 of 50000, Stored update 31750 of 50000
-2 * Loglike = 1020.262826
Actual update 31800 of 50000, Stored update 31800 of 50000
-2 * Loglike = 1020.428723
Actual update 31850 of 50000, Stored update 31850 of 50000
-2 * Loglike = 1019.044503
Actual update 31900 of 50000, Stored update 31900 of 50000
-2 * Loglike = 1018.627481
Actual update 31950 of 50000, Stored update 31950 of 50000
-2 * Loglike = 1011.454141
Actual update 32000 of 50000, Stored update 32000 of 50000
-2 * Loglike = 1012.475811
Actual update 32050 of 50000, Stored update 32050 of 50000
-2 * Loglike = 1017.569499
Actual update 32100 of 50000, Stored update 32100 of 50000
-2 * Loglike = 1011.813123
Actual update 32150 of 50000, Stored update 32150 of 50000
-2 * Loglike = 1011.879180
Actual update 32200 of 50000, Stored update 32200 of 50000
-2 * Loglike = 1010.976230
Actual update 32250 of 50000, Stored update 32250 of 50000
-2 * Loglike = 1013.425863
Actual update 32300 of 50000, Stored update 32300 of 50000
-2 * Loglike = 1011.628873
Actual update 32350 of 50000, Stored update 32350 of 50000
-2 * Loglike = 1012.742792
Actual update 32400 of 50000, Stored update 32400 of 50000
-2 * Loglike = 1020.022273
Actual update 32450 of 50000, Stored update 32450 of 50000
-2 * Loglike = 1025.788723
Actual update 32500 of 50000, Stored update 32500 of 50000
-2 * Loglike = 1020.644352
Actual update 32550 of 50000, Stored update 32550 of 50000
-2 * Loglike = 1020.933693
Actual update 32600 of 50000, Stored update 32600 of 50000
-2 * Loglike = 1015.228959
Actual update 32650 of 50000, Stored update 32650 of 50000
-2 * Loglike = 1014.891305
Actual update 32700 of 50000, Stored update 32700 of 50000
-2 * Loglike = 1018.609216
Actual update 32750 of 50000, Stored update 32750 of 50000
-2 * Loglike = 1024.384338
Actual update 32800 of 50000, Stored update 32800 of 50000
-2 * Loglike = 1009.055806
Actual update 32850 of 50000, Stored update 32850 of 50000
-2 * Loglike = 1011.949861
Actual update 32900 of 50000, Stored update 32900 of 50000
-2 * Loglike = 1007.399662
Actual update 32950 of 50000, Stored update 32950 of 50000
-2 * Loglike = 1014.238873
Actual update 33000 of 50000, Stored update 33000 of 50000
-2 * Loglike = 1010.218718
Actual update 33050 of 50000, Stored update 33050 of 50000
-2 * Loglike = 1011.519601
Actual update 33100 of 50000, Stored update 33100 of 50000
-2 * Loglike = 1011.285200
Actual update 33150 of 50000, Stored update 33150 of 50000
-2 * Loglike = 1008.241595
Actual update 33200 of 50000, Stored update 33200 of 50000
-2 * Loglike = 1022.726646
Actual update 33250 of 50000, Stored update 33250 of 50000
-2 * Loglike = 1024.251233
Actual update 33300 of 50000, Stored update 33300 of 50000
-2 * Loglike = 1016.728459
Actual update 33350 of 50000, Stored update 33350 of 50000
-2 * Loglike = 1009.741257
Actual update 33400 of 50000, Stored update 33400 of 50000
-2 * Loglike = 1012.933561
Actual update 33450 of 50000, Stored update 33450 of 50000
-2 * Loglike = 1014.531110
Actual update 33500 of 50000, Stored update 33500 of 50000
-2 * Loglike = 1018.790830
Actual update 33550 of 50000, Stored update 33550 of 50000
-2 * Loglike = 1016.742208
Actual update 33600 of 50000, Stored update 33600 of 50000
-2 * Loglike = 1011.565638
Actual update 33650 of 50000, Stored update 33650 of 50000
-2 * Loglike = 1010.209983
Actual update 33700 of 50000, Stored update 33700 of 50000
-2 * Loglike = 1017.773415
Actual update 33750 of 50000, Stored update 33750 of 50000
-2 * Loglike = 1017.985726
Actual update 33800 of 50000, Stored update 33800 of 50000
-2 * Loglike = 1026.344923
Actual update 33850 of 50000, Stored update 33850 of 50000
-2 * Loglike = 1016.865428
Actual update 33900 of 50000, Stored update 33900 of 50000
-2 * Loglike = 1015.628954
Actual update 33950 of 50000, Stored update 33950 of 50000
-2 * Loglike = 1014.440038
Actual update 34000 of 50000, Stored update 34000 of 50000
-2 * Loglike = 1019.383510
Actual update 34050 of 50000, Stored update 34050 of 50000
-2 * Loglike = 1012.906030
Actual update 34100 of 50000, Stored update 34100 of 50000
-2 * Loglike = 1008.115850
Actual update 34150 of 50000, Stored update 34150 of 50000
-2 * Loglike = 1011.197575
Actual update 34200 of 50000, Stored update 34200 of 50000
-2 * Loglike = 1015.155783
Actual update 34250 of 50000, Stored update 34250 of 50000
-2 * Loglike = 1018.570905
Actual update 34300 of 50000, Stored update 34300 of 50000
-2 * Loglike = 1016.655962
Actual update 34350 of 50000, Stored update 34350 of 50000
-2 * Loglike = 1013.687611
Actual update 34400 of 50000, Stored update 34400 of 50000
-2 * Loglike = 1015.475632
Actual update 34450 of 50000, Stored update 34450 of 50000
-2 * Loglike = 1011.535526
Actual update 34500 of 50000, Stored update 34500 of 50000
-2 * Loglike = 1017.538951
Actual update 34550 of 50000, Stored update 34550 of 50000
-2 * Loglike = 1009.529608
Actual update 34600 of 50000, Stored update 34600 of 50000
-2 * Loglike = 1012.660714
Actual update 34650 of 50000, Stored update 34650 of 50000
-2 * Loglike = 1015.524140
Actual update 34700 of 50000, Stored update 34700 of 50000
-2 * Loglike = 1017.643296
Actual update 34750 of 50000, Stored update 34750 of 50000
-2 * Loglike = 1011.184569
Actual update 34800 of 50000, Stored update 34800 of 50000
-2 * Loglike = 1025.510893
Actual update 34850 of 50000, Stored update 34850 of 50000
-2 * Loglike = 1012.934228
Actual update 34900 of 50000, Stored update 34900 of 50000
-2 * Loglike = 1008.051407
Actual update 34950 of 50000, Stored update 34950 of 50000
-2 * Loglike = 1013.867223
Actual update 35000 of 50000, Stored update 35000 of 50000
-2 * Loglike = 1015.610339
Actual update 35050 of 50000, Stored update 35050 of 50000
-2 * Loglike = 1013.268010
Actual update 35100 of 50000, Stored update 35100 of 50000
-2 * Loglike = 1016.234616
Actual update 35150 of 50000, Stored update 35150 of 50000
-2 * Loglike = 1019.223870
Actual update 35200 of 50000, Stored update 35200 of 50000
-2 * Loglike = 1017.620650
Actual update 35250 of 50000, Stored update 35250 of 50000
-2 * Loglike = 1019.090994
Actual update 35300 of 50000, Stored update 35300 of 50000
-2 * Loglike = 1024.579189
Actual update 35350 of 50000, Stored update 35350 of 50000
-2 * Loglike = 1016.953481
Actual update 35400 of 50000, Stored update 35400 of 50000
-2 * Loglike = 1013.727411
Actual update 35450 of 50000, Stored update 35450 of 50000
-2 * Loglike = 1022.168272
Actual update 35500 of 50000, Stored update 35500 of 50000
-2 * Loglike = 1026.369643
Actual update 35550 of 50000, Stored update 35550 of 50000
-2 * Loglike = 1017.765494
Actual update 35600 of 50000, Stored update 35600 of 50000
-2 * Loglike = 1020.834589
Actual update 35650 of 50000, Stored update 35650 of 50000
-2 * Loglike = 1023.012842
Actual update 35700 of 50000, Stored update 35700 of 50000
-2 * Loglike = 1013.868587
Actual update 35750 of 50000, Stored update 35750 of 50000
-2 * Loglike = 1015.191781
Actual update 35800 of 50000, Stored update 35800 of 50000
-2 * Loglike = 1009.579886
Actual update 35850 of 50000, Stored update 35850 of 50000
-2 * Loglike = 1021.669174
Actual update 35900 of 50000, Stored update 35900 of 50000
-2 * Loglike = 1015.532884
Actual update 35950 of 50000, Stored update 35950 of 50000
-2 * Loglike = 1021.175026
Actual update 36000 of 50000, Stored update 36000 of 50000
-2 * Loglike = 1017.729193
Actual update 36050 of 50000, Stored update 36050 of 50000
-2 * Loglike = 1023.237692
Actual update 36100 of 50000, Stored update 36100 of 50000
-2 * Loglike = 1028.603063
Actual update 36150 of 50000, Stored update 36150 of 50000
-2 * Loglike = 1018.668304
Actual update 36200 of 50000, Stored update 36200 of 50000
-2 * Loglike = 1029.241625
Actual update 36250 of 50000, Stored update 36250 of 50000
-2 * Loglike = 1025.949979
Actual update 36300 of 50000, Stored update 36300 of 50000
-2 * Loglike = 1034.853912
Actual update 36350 of 50000, Stored update 36350 of 50000
-2 * Loglike = 1020.249395
Actual update 36400 of 50000, Stored update 36400 of 50000
-2 * Loglike = 1019.540404
Actual update 36450 of 50000, Stored update 36450 of 50000
-2 * Loglike = 1019.907318
Actual update 36500 of 50000, Stored update 36500 of 50000
-2 * Loglike = 1014.228262
Actual update 36550 of 50000, Stored update 36550 of 50000
-2 * Loglike = 1010.615286
Actual update 36600 of 50000, Stored update 36600 of 50000
-2 * Loglike = 1027.263674
Actual update 36650 of 50000, Stored update 36650 of 50000
-2 * Loglike = 1009.580459
Actual update 36700 of 50000, Stored update 36700 of 50000
-2 * Loglike = 1009.630336
Actual update 36750 of 50000, Stored update 36750 of 50000
-2 * Loglike = 1013.244616
Actual update 36800 of 50000, Stored update 36800 of 50000
-2 * Loglike = 1012.795628
Actual update 36850 of 50000, Stored update 36850 of 50000
-2 * Loglike = 1015.670309
Actual update 36900 of 50000, Stored update 36900 of 50000
-2 * Loglike = 1007.633600
Actual update 36950 of 50000, Stored update 36950 of 50000
-2 * Loglike = 1012.153551
Actual update 37000 of 50000, Stored update 37000 of 50000
-2 * Loglike = 1013.846125
Actual update 37050 of 50000, Stored update 37050 of 50000
-2 * Loglike = 1013.426087
Actual update 37100 of 50000, Stored update 37100 of 50000
-2 * Loglike = 1014.822932
Actual update 37150 of 50000, Stored update 37150 of 50000
-2 * Loglike = 1012.857011
Actual update 37200 of 50000, Stored update 37200 of 50000
-2 * Loglike = 1014.525586
Actual update 37250 of 50000, Stored update 37250 of 50000
-2 * Loglike = 1012.826960
Actual update 37300 of 50000, Stored update 37300 of 50000
-2 * Loglike = 1025.154903
Actual update 37350 of 50000, Stored update 37350 of 50000
-2 * Loglike = 1013.267033
Actual update 37400 of 50000, Stored update 37400 of 50000
-2 * Loglike = 1027.439747
Actual update 37450 of 50000, Stored update 37450 of 50000
-2 * Loglike = 1009.967629
Actual update 37500 of 50000, Stored update 37500 of 50000
-2 * Loglike = 1018.426390
Actual update 37550 of 50000, Stored update 37550 of 50000
-2 * Loglike = 1018.930702
Actual update 37600 of 50000, Stored update 37600 of 50000
-2 * Loglike = 1011.933841
Actual update 37650 of 50000, Stored update 37650 of 50000
-2 * Loglike = 1018.379476
Actual update 37700 of 50000, Stored update 37700 of 50000
-2 * Loglike = 1012.213997
Actual update 37750 of 50000, Stored update 37750 of 50000
-2 * Loglike = 1014.349920
Actual update 37800 of 50000, Stored update 37800 of 50000
-2 * Loglike = 1017.414453
Actual update 37850 of 50000, Stored update 37850 of 50000
-2 * Loglike = 1007.442575
Actual update 37900 of 50000, Stored update 37900 of 50000
-2 * Loglike = 1011.015450
Actual update 37950 of 50000, Stored update 37950 of 50000
-2 * Loglike = 1027.317097
Actual update 38000 of 50000, Stored update 38000 of 50000
-2 * Loglike = 1019.949236
Actual update 38050 of 50000, Stored update 38050 of 50000
-2 * Loglike = 1031.025497
Actual update 38100 of 50000, Stored update 38100 of 50000
-2 * Loglike = 1019.268610
Actual update 38150 of 50000, Stored update 38150 of 50000
-2 * Loglike = 1015.351948
Actual update 38200 of 50000, Stored update 38200 of 50000
-2 * Loglike = 1016.777414
Actual update 38250 of 50000, Stored update 38250 of 50000
-2 * Loglike = 1019.482160
Actual update 38300 of 50000, Stored update 38300 of 50000
-2 * Loglike = 1017.113572
Actual update 38350 of 50000, Stored update 38350 of 50000
-2 * Loglike = 1013.686688
Actual update 38400 of 50000, Stored update 38400 of 50000
-2 * Loglike = 1016.682589
Actual update 38450 of 50000, Stored update 38450 of 50000
-2 * Loglike = 1026.631735
Actual update 38500 of 50000, Stored update 38500 of 50000
-2 * Loglike = 1014.259554
Actual update 38550 of 50000, Stored update 38550 of 50000
-2 * Loglike = 1020.046944
Actual update 38600 of 50000, Stored update 38600 of 50000
-2 * Loglike = 1015.901544
Actual update 38650 of 50000, Stored update 38650 of 50000
-2 * Loglike = 1019.109739
Actual update 38700 of 50000, Stored update 38700 of 50000
-2 * Loglike = 1013.790301
Actual update 38750 of 50000, Stored update 38750 of 50000
-2 * Loglike = 1024.274104
Actual update 38800 of 50000, Stored update 38800 of 50000
-2 * Loglike = 1009.308985
Actual update 38850 of 50000, Stored update 38850 of 50000
-2 * Loglike = 1017.721623
Actual update 38900 of 50000, Stored update 38900 of 50000
-2 * Loglike = 1018.364468
Actual update 38950 of 50000, Stored update 38950 of 50000
-2 * Loglike = 1022.589230
Actual update 39000 of 50000, Stored update 39000 of 50000
-2 * Loglike = 1023.218167
Actual update 39050 of 50000, Stored update 39050 of 50000
-2 * Loglike = 1015.464636
Actual update 39100 of 50000, Stored update 39100 of 50000
-2 * Loglike = 1013.365776
Actual update 39150 of 50000, Stored update 39150 of 50000
-2 * Loglike = 1013.834657
Actual update 39200 of 50000, Stored update 39200 of 50000
-2 * Loglike = 1013.276261
Actual update 39250 of 50000, Stored update 39250 of 50000
-2 * Loglike = 1027.136349
Actual update 39300 of 50000, Stored update 39300 of 50000
-2 * Loglike = 1011.953578
Actual update 39350 of 50000, Stored update 39350 of 50000
-2 * Loglike = 1021.253360
Actual update 39400 of 50000, Stored update 39400 of 50000
-2 * Loglike = 1017.246615
Actual update 39450 of 50000, Stored update 39450 of 50000
-2 * Loglike = 1013.089311
Actual update 39500 of 50000, Stored update 39500 of 50000
-2 * Loglike = 1022.580497
Actual update 39550 of 50000, Stored update 39550 of 50000
-2 * Loglike = 1009.595976
Actual update 39600 of 50000, Stored update 39600 of 50000
-2 * Loglike = 1009.652332
Actual update 39650 of 50000, Stored update 39650 of 50000
-2 * Loglike = 1026.346520
Actual update 39700 of 50000, Stored update 39700 of 50000
-2 * Loglike = 1023.709867
Actual update 39750 of 50000, Stored update 39750 of 50000
-2 * Loglike = 1012.929684
Actual update 39800 of 50000, Stored update 39800 of 50000
-2 * Loglike = 1017.935648
Actual update 39850 of 50000, Stored update 39850 of 50000
-2 * Loglike = 1021.273383
Actual update 39900 of 50000, Stored update 39900 of 50000
-2 * Loglike = 1008.621794
Actual update 39950 of 50000, Stored update 39950 of 50000
-2 * Loglike = 1012.923481
Actual update 40000 of 50000, Stored update 40000 of 50000
-2 * Loglike = 1018.720435
Actual update 40050 of 50000, Stored update 40050 of 50000
-2 * Loglike = 1019.307354
Actual update 40100 of 50000, Stored update 40100 of 50000
-2 * Loglike = 1015.234294
Actual update 40150 of 50000, Stored update 40150 of 50000
-2 * Loglike = 1016.766313
Actual update 40200 of 50000, Stored update 40200 of 50000
-2 * Loglike = 1011.409181
Actual update 40250 of 50000, Stored update 40250 of 50000
-2 * Loglike = 1015.696719
Actual update 40300 of 50000, Stored update 40300 of 50000
-2 * Loglike = 1011.626855
Actual update 40350 of 50000, Stored update 40350 of 50000
-2 * Loglike = 1011.609198
Actual update 40400 of 50000, Stored update 40400 of 50000
-2 * Loglike = 1013.067471
Actual update 40450 of 50000, Stored update 40450 of 50000
-2 * Loglike = 1012.760804
Actual update 40500 of 50000, Stored update 40500 of 50000
-2 * Loglike = 1009.619017
Actual update 40550 of 50000, Stored update 40550 of 50000
-2 * Loglike = 1007.737983
Actual update 40600 of 50000, Stored update 40600 of 50000
-2 * Loglike = 1026.428751
Actual update 40650 of 50000, Stored update 40650 of 50000
-2 * Loglike = 1009.307452
Actual update 40700 of 50000, Stored update 40700 of 50000
-2 * Loglike = 1010.652319
Actual update 40750 of 50000, Stored update 40750 of 50000
-2 * Loglike = 1007.765919
Actual update 40800 of 50000, Stored update 40800 of 50000
-2 * Loglike = 1017.944424
Actual update 40850 of 50000, Stored update 40850 of 50000
-2 * Loglike = 1017.419814
Actual update 40900 of 50000, Stored update 40900 of 50000
-2 * Loglike = 1012.361062
Actual update 40950 of 50000, Stored update 40950 of 50000
-2 * Loglike = 1010.130151
Actual update 41000 of 50000, Stored update 41000 of 50000
-2 * Loglike = 1015.939877
Actual update 41050 of 50000, Stored update 41050 of 50000
-2 * Loglike = 1010.179255
Actual update 41100 of 50000, Stored update 41100 of 50000
-2 * Loglike = 1009.835432
Actual update 41150 of 50000, Stored update 41150 of 50000
-2 * Loglike = 1016.697811
Actual update 41200 of 50000, Stored update 41200 of 50000
-2 * Loglike = 1015.579601
Actual update 41250 of 50000, Stored update 41250 of 50000
-2 * Loglike = 1013.644970
Actual update 41300 of 50000, Stored update 41300 of 50000
-2 * Loglike = 1015.415643
Actual update 41350 of 50000, Stored update 41350 of 50000
-2 * Loglike = 1016.451752
Actual update 41400 of 50000, Stored update 41400 of 50000
-2 * Loglike = 1016.217666
Actual update 41450 of 50000, Stored update 41450 of 50000
-2 * Loglike = 1018.065986
Actual update 41500 of 50000, Stored update 41500 of 50000
-2 * Loglike = 1016.763248
Actual update 41550 of 50000, Stored update 41550 of 50000
-2 * Loglike = 1011.008421
Actual update 41600 of 50000, Stored update 41600 of 50000
-2 * Loglike = 1007.201419
Actual update 41650 of 50000, Stored update 41650 of 50000
-2 * Loglike = 1013.574334
Actual update 41700 of 50000, Stored update 41700 of 50000
-2 * Loglike = 1018.654852
Actual update 41750 of 50000, Stored update 41750 of 50000
-2 * Loglike = 1006.347415
Actual update 41800 of 50000, Stored update 41800 of 50000
-2 * Loglike = 1017.660107
Actual update 41850 of 50000, Stored update 41850 of 50000
-2 * Loglike = 1006.964707
Actual update 41900 of 50000, Stored update 41900 of 50000
-2 * Loglike = 1021.683434
Actual update 41950 of 50000, Stored update 41950 of 50000
-2 * Loglike = 1014.572015
Actual update 42000 of 50000, Stored update 42000 of 50000
-2 * Loglike = 1017.215517
Actual update 42050 of 50000, Stored update 42050 of 50000
-2 * Loglike = 1013.038388
Actual update 42100 of 50000, Stored update 42100 of 50000
-2 * Loglike = 1013.659412
Actual update 42150 of 50000, Stored update 42150 of 50000
-2 * Loglike = 1015.861877
Actual update 42200 of 50000, Stored update 42200 of 50000
-2 * Loglike = 1012.470942
Actual update 42250 of 50000, Stored update 42250 of 50000
-2 * Loglike = 1017.931465
Actual update 42300 of 50000, Stored update 42300 of 50000
-2 * Loglike = 1011.471872
Actual update 42350 of 50000, Stored update 42350 of 50000
-2 * Loglike = 1005.559930
Actual update 42400 of 50000, Stored update 42400 of 50000
-2 * Loglike = 1009.286208
Actual update 42450 of 50000, Stored update 42450 of 50000
-2 * Loglike = 1006.248901
Actual update 42500 of 50000, Stored update 42500 of 50000
-2 * Loglike = 1011.377708
Actual update 42550 of 50000, Stored update 42550 of 50000
-2 * Loglike = 1014.588306
Actual update 42600 of 50000, Stored update 42600 of 50000
-2 * Loglike = 1007.667555
Actual update 42650 of 50000, Stored update 42650 of 50000
-2 * Loglike = 1017.149706
Actual update 42700 of 50000, Stored update 42700 of 50000
-2 * Loglike = 1015.749133
Actual update 42750 of 50000, Stored update 42750 of 50000
-2 * Loglike = 1023.165612
Actual update 42800 of 50000, Stored update 42800 of 50000
-2 * Loglike = 1011.073359
Actual update 42850 of 50000, Stored update 42850 of 50000
-2 * Loglike = 1026.110836
Actual update 42900 of 50000, Stored update 42900 of 50000
-2 * Loglike = 1020.029657
Actual update 42950 of 50000, Stored update 42950 of 50000
-2 * Loglike = 1007.179671
Actual update 43000 of 50000, Stored update 43000 of 50000
-2 * Loglike = 1013.832312
Actual update 43050 of 50000, Stored update 43050 of 50000
-2 * Loglike = 1013.849582
Actual update 43100 of 50000, Stored update 43100 of 50000
-2 * Loglike = 1019.648720
Actual update 43150 of 50000, Stored update 43150 of 50000
-2 * Loglike = 1012.336025
Actual update 43200 of 50000, Stored update 43200 of 50000
-2 * Loglike = 1017.298653
Actual update 43250 of 50000, Stored update 43250 of 50000
-2 * Loglike = 1012.632440
Actual update 43300 of 50000, Stored update 43300 of 50000
-2 * Loglike = 1014.479484
Actual update 43350 of 50000, Stored update 43350 of 50000
-2 * Loglike = 1021.559446
Actual update 43400 of 50000, Stored update 43400 of 50000
-2 * Loglike = 1018.008961
Actual update 43450 of 50000, Stored update 43450 of 50000
-2 * Loglike = 1018.119187
Actual update 43500 of 50000, Stored update 43500 of 50000
-2 * Loglike = 1028.136304
Actual update 43550 of 50000, Stored update 43550 of 50000
-2 * Loglike = 1023.718373
Actual update 43600 of 50000, Stored update 43600 of 50000
-2 * Loglike = 1013.610340
Actual update 43650 of 50000, Stored update 43650 of 50000
-2 * Loglike = 1022.560843
Actual update 43700 of 50000, Stored update 43700 of 50000
-2 * Loglike = 1021.661887
Actual update 43750 of 50000, Stored update 43750 of 50000
-2 * Loglike = 1023.872263
Actual update 43800 of 50000, Stored update 43800 of 50000
-2 * Loglike = 1019.592850
Actual update 43850 of 50000, Stored update 43850 of 50000
-2 * Loglike = 1033.023415
Actual update 43900 of 50000, Stored update 43900 of 50000
-2 * Loglike = 1016.068031
Actual update 43950 of 50000, Stored update 43950 of 50000
-2 * Loglike = 1021.791338
Actual update 44000 of 50000, Stored update 44000 of 50000
-2 * Loglike = 1013.375598
Actual update 44050 of 50000, Stored update 44050 of 50000
-2 * Loglike = 1020.736708
Actual update 44100 of 50000, Stored update 44100 of 50000
-2 * Loglike = 1020.809018
Actual update 44150 of 50000, Stored update 44150 of 50000
-2 * Loglike = 1016.917299
Actual update 44200 of 50000, Stored update 44200 of 50000
-2 * Loglike = 1016.008858
Actual update 44250 of 50000, Stored update 44250 of 50000
-2 * Loglike = 1013.933323
Actual update 44300 of 50000, Stored update 44300 of 50000
-2 * Loglike = 1012.416308
Actual update 44350 of 50000, Stored update 44350 of 50000
-2 * Loglike = 1014.453843
Actual update 44400 of 50000, Stored update 44400 of 50000
-2 * Loglike = 1021.424604
Actual update 44450 of 50000, Stored update 44450 of 50000
-2 * Loglike = 1014.929860
Actual update 44500 of 50000, Stored update 44500 of 50000
-2 * Loglike = 1011.058740
Actual update 44550 of 50000, Stored update 44550 of 50000
-2 * Loglike = 1018.685427
Actual update 44600 of 50000, Stored update 44600 of 50000
-2 * Loglike = 1014.438178
Actual update 44650 of 50000, Stored update 44650 of 50000
-2 * Loglike = 1009.099571
Actual update 44700 of 50000, Stored update 44700 of 50000
-2 * Loglike = 1019.649136
Actual update 44750 of 50000, Stored update 44750 of 50000
-2 * Loglike = 1011.329398
Actual update 44800 of 50000, Stored update 44800 of 50000
-2 * Loglike = 1014.330768
Actual update 44850 of 50000, Stored update 44850 of 50000
-2 * Loglike = 1018.096415
Actual update 44900 of 50000, Stored update 44900 of 50000
-2 * Loglike = 1016.692859
Actual update 44950 of 50000, Stored update 44950 of 50000
-2 * Loglike = 1034.898666
Actual update 45000 of 50000, Stored update 45000 of 50000
-2 * Loglike = 1009.498898
Actual update 45050 of 50000, Stored update 45050 of 50000
-2 * Loglike = 1016.413502
Actual update 45100 of 50000, Stored update 45100 of 50000
-2 * Loglike = 1015.745939
Actual update 45150 of 50000, Stored update 45150 of 50000
-2 * Loglike = 1018.742556
Actual update 45200 of 50000, Stored update 45200 of 50000
-2 * Loglike = 1019.306511
Actual update 45250 of 50000, Stored update 45250 of 50000
-2 * Loglike = 1014.452897
Actual update 45300 of 50000, Stored update 45300 of 50000
-2 * Loglike = 1009.849451
Actual update 45350 of 50000, Stored update 45350 of 50000
-2 * Loglike = 1007.813453
Actual update 45400 of 50000, Stored update 45400 of 50000
-2 * Loglike = 1022.318190
Actual update 45450 of 50000, Stored update 45450 of 50000
-2 * Loglike = 1016.420654
Actual update 45500 of 50000, Stored update 45500 of 50000
-2 * Loglike = 1026.353085
Actual update 45550 of 50000, Stored update 45550 of 50000
-2 * Loglike = 1007.877553
Actual update 45600 of 50000, Stored update 45600 of 50000
-2 * Loglike = 1011.016672
Actual update 45650 of 50000, Stored update 45650 of 50000
-2 * Loglike = 1010.121166
Actual update 45700 of 50000, Stored update 45700 of 50000
-2 * Loglike = 1013.949134
Actual update 45750 of 50000, Stored update 45750 of 50000
-2 * Loglike = 1010.597545
Actual update 45800 of 50000, Stored update 45800 of 50000
-2 * Loglike = 1017.941617
Actual update 45850 of 50000, Stored update 45850 of 50000
-2 * Loglike = 1006.185291
Actual update 45900 of 50000, Stored update 45900 of 50000
-2 * Loglike = 1017.665340
Actual update 45950 of 50000, Stored update 45950 of 50000
-2 * Loglike = 1017.400967
Actual update 46000 of 50000, Stored update 46000 of 50000
-2 * Loglike = 1009.746855
Actual update 46050 of 50000, Stored update 46050 of 50000
-2 * Loglike = 1012.986622
Actual update 46100 of 50000, Stored update 46100 of 50000
-2 * Loglike = 1010.183813
Actual update 46150 of 50000, Stored update 46150 of 50000
-2 * Loglike = 1007.758249
Actual update 46200 of 50000, Stored update 46200 of 50000
-2 * Loglike = 1014.166081
Actual update 46250 of 50000, Stored update 46250 of 50000
-2 * Loglike = 1019.943380
Actual update 46300 of 50000, Stored update 46300 of 50000
-2 * Loglike = 1010.882182
Actual update 46350 of 50000, Stored update 46350 of 50000
-2 * Loglike = 1016.344244
Actual update 46400 of 50000, Stored update 46400 of 50000
-2 * Loglike = 1009.166636
Actual update 46450 of 50000, Stored update 46450 of 50000
-2 * Loglike = 1013.033375
Actual update 46500 of 50000, Stored update 46500 of 50000
-2 * Loglike = 1011.619360
Actual update 46550 of 50000, Stored update 46550 of 50000
-2 * Loglike = 1026.279838
Actual update 46600 of 50000, Stored update 46600 of 50000
-2 * Loglike = 1016.779558
Actual update 46650 of 50000, Stored update 46650 of 50000
-2 * Loglike = 1017.482440
Actual update 46700 of 50000, Stored update 46700 of 50000
-2 * Loglike = 1035.164280
Actual update 46750 of 50000, Stored update 46750 of 50000
-2 * Loglike = 1022.019867
Actual update 46800 of 50000, Stored update 46800 of 50000
-2 * Loglike = 1015.398784
Actual update 46850 of 50000, Stored update 46850 of 50000
-2 * Loglike = 1015.143969
Actual update 46900 of 50000, Stored update 46900 of 50000
-2 * Loglike = 1008.607687
Actual update 46950 of 50000, Stored update 46950 of 50000
-2 * Loglike = 1016.864691
Actual update 47000 of 50000, Stored update 47000 of 50000
-2 * Loglike = 1020.149287
Actual update 47050 of 50000, Stored update 47050 of 50000
-2 * Loglike = 1014.384425
Actual update 47100 of 50000, Stored update 47100 of 50000
-2 * Loglike = 1009.033449
Actual update 47150 of 50000, Stored update 47150 of 50000
-2 * Loglike = 1013.944668
Actual update 47200 of 50000, Stored update 47200 of 50000
-2 * Loglike = 1011.215152
Actual update 47250 of 50000, Stored update 47250 of 50000
-2 * Loglike = 1022.830737
Actual update 47300 of 50000, Stored update 47300 of 50000
-2 * Loglike = 1009.989252
Actual update 47350 of 50000, Stored update 47350 of 50000
-2 * Loglike = 1009.750907
Actual update 47400 of 50000, Stored update 47400 of 50000
-2 * Loglike = 1016.511961
Actual update 47450 of 50000, Stored update 47450 of 50000
-2 * Loglike = 1015.576075
Actual update 47500 of 50000, Stored update 47500 of 50000
-2 * Loglike = 1013.078344
Actual update 47550 of 50000, Stored update 47550 of 50000
-2 * Loglike = 1017.021667
Actual update 47600 of 50000, Stored update 47600 of 50000
-2 * Loglike = 1009.400686
Actual update 47650 of 50000, Stored update 47650 of 50000
-2 * Loglike = 1012.030414
Actual update 47700 of 50000, Stored update 47700 of 50000
-2 * Loglike = 1013.202048
Actual update 47750 of 50000, Stored update 47750 of 50000
-2 * Loglike = 1008.563831
Actual update 47800 of 50000, Stored update 47800 of 50000
-2 * Loglike = 1010.145806
Actual update 47850 of 50000, Stored update 47850 of 50000
-2 * Loglike = 1015.663682
Actual update 47900 of 50000, Stored update 47900 of 50000
-2 * Loglike = 1012.917417
Actual update 47950 of 50000, Stored update 47950 of 50000
-2 * Loglike = 1013.074049
Actual update 48000 of 50000, Stored update 48000 of 50000
-2 * Loglike = 1019.192202
Actual update 48050 of 50000, Stored update 48050 of 50000
-2 * Loglike = 1010.451228
Actual update 48100 of 50000, Stored update 48100 of 50000
-2 * Loglike = 1013.146405
Actual update 48150 of 50000, Stored update 48150 of 50000
-2 * Loglike = 1011.272163
Actual update 48200 of 50000, Stored update 48200 of 50000
-2 * Loglike = 1026.175229
Actual update 48250 of 50000, Stored update 48250 of 50000
-2 * Loglike = 1015.247527
Actual update 48300 of 50000, Stored update 48300 of 50000
-2 * Loglike = 1008.312349
Actual update 48350 of 50000, Stored update 48350 of 50000
-2 * Loglike = 1012.424343
Actual update 48400 of 50000, Stored update 48400 of 50000
-2 * Loglike = 1008.094262
Actual update 48450 of 50000, Stored update 48450 of 50000
-2 * Loglike = 1017.146038
Actual update 48500 of 50000, Stored update 48500 of 50000
-2 * Loglike = 1009.376353
Actual update 48550 of 50000, Stored update 48550 of 50000
-2 * Loglike = 1012.845103
Actual update 48600 of 50000, Stored update 48600 of 50000
-2 * Loglike = 1018.827319
Actual update 48650 of 50000, Stored update 48650 of 50000
-2 * Loglike = 1007.456953
Actual update 48700 of 50000, Stored update 48700 of 50000
-2 * Loglike = 1008.678272
Actual update 48750 of 50000, Stored update 48750 of 50000
-2 * Loglike = 1012.165366
Actual update 48800 of 50000, Stored update 48800 of 50000
-2 * Loglike = 1015.348380
Actual update 48850 of 50000, Stored update 48850 of 50000
-2 * Loglike = 1004.712810
Actual update 48900 of 50000, Stored update 48900 of 50000
-2 * Loglike = 1008.726124
Actual update 48950 of 50000, Stored update 48950 of 50000
-2 * Loglike = 1008.605697
Actual update 49000 of 50000, Stored update 49000 of 50000
-2 * Loglike = 1014.806138
Actual update 49050 of 50000, Stored update 49050 of 50000
-2 * Loglike = 1011.869344
Actual update 49100 of 50000, Stored update 49100 of 50000
-2 * Loglike = 1008.158289
Actual update 49150 of 50000, Stored update 49150 of 50000
-2 * Loglike = 1013.971231
Actual update 49200 of 50000, Stored update 49200 of 50000
-2 * Loglike = 1013.664295
Actual update 49250 of 50000, Stored update 49250 of 50000
-2 * Loglike = 1019.209775
Actual update 49300 of 50000, Stored update 49300 of 50000
-2 * Loglike = 1023.973890
Actual update 49350 of 50000, Stored update 49350 of 50000
-2 * Loglike = 1012.283011
Actual update 49400 of 50000, Stored update 49400 of 50000
-2 * Loglike = 1009.184983
Actual update 49450 of 50000, Stored update 49450 of 50000
-2 * Loglike = 1012.493149
Actual update 49500 of 50000, Stored update 49500 of 50000
-2 * Loglike = 1011.023499
Actual update 49550 of 50000, Stored update 49550 of 50000
-2 * Loglike = 1014.696888
Actual update 49600 of 50000, Stored update 49600 of 50000
-2 * Loglike = 1007.361718
Actual update 49650 of 50000, Stored update 49650 of 50000
-2 * Loglike = 1012.970215
Actual update 49700 of 50000, Stored update 49700 of 50000
-2 * Loglike = 1017.844700
Actual update 49750 of 50000, Stored update 49750 of 50000
-2 * Loglike = 1014.576744
Actual update 49800 of 50000, Stored update 49800 of 50000
-2 * Loglike = 1008.832715
Actual update 49850 of 50000, Stored update 49850 of 50000
-2 * Loglike = 1016.122609
Actual update 49900 of 50000, Stored update 49900 of 50000
-2 * Loglike = 1009.735502
Actual update 49950 of 50000, Stored update 49950 of 50000
-2 * Loglike = 1011.263113
Actual update 50000 of 50000, Stored update 50000 of 50000
-2 * Loglike = 1010.892600
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =          50000
Average =          1015.4
S.D.    =          5.8158
S.E.M.  =        0.026009

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  multilevel model (Multivariate Normal) 
Estimation algorithm:  MCMC      Elapsed time : 3.84s 
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.423   1006.005   9.418      1024.840   
--------------------------------------------------------------------------------------------------- 
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.16955     2.48435  61.25          0  ***   147.25693   157.09464   50000 
Intercept_y15   201.76484     2.81909  71.57          0  ***   196.19856   207.31011   50000 
Intercept_y22   245.03782     3.35730  72.99          0  ***   238.42476   251.67358   50000 
Intercept_y29   289.50561     4.16949  69.43          0  ***   281.27449   297.73565   50000 
Intercept_y36   324.80671     4.23349  76.72          0  ***   316.44088   333.18680   49056 
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.99986    64.63109     99.27080    349.14376   266 
cov_Intercept_y8_Intercept_y15    197.96303    70.90982    104.18641    384.19912   229 
var_Intercept_y15                 238.45038    83.02111    128.03163    456.41404   257 
cov_Intercept_y8_Intercept_y22    222.66138    82.27239    113.50299    441.85902   200 
cov_Intercept_y15_Intercept_y22   267.96912    95.66329    139.64078    524.43019   221 
var_Intercept_y22                 338.66755   116.70093    180.02655    648.20730   229 
cov_Intercept_y8_Intercept_y29    260.67292   100.06371    127.81774    534.59883   190 
cov_Intercept_y15_Intercept_y29   313.58451   116.24887    157.67839    628.60641   205 
cov_Intercept_y22_Intercept_y29   395.91211   140.85835    204.83120    778.70819   207 
var_Intercept_y29                 520.39377   179.25479    277.65158   1008.49029   225 
cov_Intercept_y8_Intercept_y36    250.00869    98.12261    119.43459    517.81668   204 
cov_Intercept_y15_Intercept_y36   300.69907   113.98774    148.09726    612.21703   200 
cov_Intercept_y22_Intercept_y36   379.48084   137.84916    193.14231    755.83215   206 
cov_Intercept_y29_Intercept_y36   498.41114   174.27542    262.93141    977.45230   236 
var_Intercept_y36                 537.40121   180.64997    292.34950   1017.60914   241 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> 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  185.000       NA       NA       NA       NA
cons.y15 197.963  238.450       NA       NA       NA
cons.y22 222.661  267.969  338.668       NA       NA
cons.y29 260.673  313.585  395.912  520.394       NA
cons.y36 250.009  300.699  379.481  498.411  537.401
> 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.943    1.000       NA       NA       NA
cons.y22   0.890    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.942        1
> 
> # Chapter learning outcomes . . . . . . . . . . . . . . . . . . . . . . .128
> 
> 
> 
> 
> 
> ############################################################################
> 
> proc.time()
   user  system elapsed 
   9.40    1.59   24.39