R version 4.4.1 (2024-06-14 ucrt) -- "Race for Your Life"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

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.

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
> #
> # 5   Prior Distributions, Starting Values and Random Number Seeds . . . .61
> #
> #     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/
> ############################################################################
> 
> # 5.1 Prior distributions . . . . . . . . . . . . . . . . . . . . . . . . 61
> 
> # 5.2 Uniform on variance scale priors . . . . . . . . . . . . . . . . . .61
> 
> 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.11/
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 tutorial data
> data(tutorial, package = "R2MLwiN")
> 
> ## IGLS
> (mymodel1 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), data = tutorial))
/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
NEXT
iteration 2
iteration 3

Convergence achieved
ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  IGLS        Elapsed time : 0.08s 
Number of obs:  4059 (from total 4059)        The model converged after 4 iterations.
Log likelihood:      -4678.6 
Deviance statistic:  9357.2 
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
              Coef.   Std. Err.       z   Pr(>|z|)         [95% Conf.   Interval] 
Intercept   0.00239     0.04002    0.06     0.9524           -0.07605     0.08083 
standlrt    0.56337     0.01247   45.19          0   ***      0.53894     0.58780 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the school level: 
                  Coef.   Std. Err. 
var_Intercept   0.09213     0.01815 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err. 
var_Intercept   0.56573     0.01266 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ## Diffuse priors (Gamma priors)
> (mymodel2 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1), data = tutorial))
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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500
-2 * Loglike = 9206.193266
Burning in for 100 iterations out of 500
-2 * Loglike = 9216.262176
Burning in for 150 iterations out of 500
-2 * Loglike = 9241.057582
Burning in for 200 iterations out of 500
-2 * Loglike = 9214.675885
Burning in for 250 iterations out of 500
-2 * Loglike = 9198.255145
Burning in for 300 iterations out of 500
-2 * Loglike = 9227.186873
Burning in for 350 iterations out of 500
-2 * Loglike = 9222.091977
Burning in for 400 iterations out of 500
-2 * Loglike = 9214.630238
Burning in for 450 iterations out of 500
-2 * Loglike = 9216.522216
Burning in for 500 iterations out of 500
-2 * Loglike = 9201.148603

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9212.190791
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9208.377297
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9207.886882
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9195.385976
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9238.483062
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9215.621742
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9205.716499
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9207.158850
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9210.248427
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9228.495732
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9221.683303
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9202.166694
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9237.981658
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9201.878383
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9198.726654
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9214.104227
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9227.156460
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9197.279082
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9191.816267
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9219.178731
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9214.163387
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9218.308723
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9188.214048
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9188.278770
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9208.819913
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9193.939811
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9193.651680
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9215.606479
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9195.925780
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9218.422125
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9202.244183
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9204.907738
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9194.883902
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9214.140517
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9211.397609
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9207.941176
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9216.158993
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9216.821254
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9209.099646
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9205.356036
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9226.820798
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9215.889068
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9204.485567
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9239.773102
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9204.613020
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9205.832435
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9197.787494
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9202.635017
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9202.243753
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9199.686027
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9213.435702
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9215.862131
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9203.440604
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9197.917390
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9225.824248
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9222.623770
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9226.897380
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9197.383677
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9211.390320
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9207.286579
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9217.534200
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9227.872396
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9198.935844
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9203.578259
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9193.257434
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9204.815469
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9200.291358
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9219.578031
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9214.657834
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9213.144417
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9217.997575
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9225.109666
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9212.063094
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9226.746381
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9214.159422
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9204.966122
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9190.735716
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9206.324394
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9205.512200
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9201.017143
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9211.234780
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9211.765768
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9200.968089
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9229.231239
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9207.636298
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9197.507235
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9208.445006
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9231.667774
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9201.470719
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9208.425209
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9216.996959
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9225.405793
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9201.746288
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9212.656436
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9202.976127
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9218.846042
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9200.637551
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9219.146880
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9215.841845
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9212.109834
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9209.1
S.D.    =          11.916
S.E.M.  =         0.16852

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 1.56s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
9209.078   9149.069   60.010     9269.088   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
              Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   0.00318     0.04420   0.07     0.9427         -0.08330     0.09189    223 
standlrt    0.56322     0.01258  44.77          0  ***     0.53844     0.58780   3821 
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   0.09742     0.01988      0.06472     0.14144   3126 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.56617     0.01249      0.54212     0.59090   4751 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ## Diffuse priors (Uniform priors)
> (mymodel3 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcMeth = list(priorcode = 0)), 
+   data = tutorial))
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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 0 1
BURNING IN...
Burning in for 50 iterations out of 500
-2 * Loglike = 9205.810806
Burning in for 100 iterations out of 500
-2 * Loglike = 9215.875072
Burning in for 150 iterations out of 500
-2 * Loglike = 9241.027909
Burning in for 200 iterations out of 500
-2 * Loglike = 9214.690189
Burning in for 250 iterations out of 500
-2 * Loglike = 9197.915505
Burning in for 300 iterations out of 500
-2 * Loglike = 9226.691192
Burning in for 350 iterations out of 500
-2 * Loglike = 9221.939812
Burning in for 400 iterations out of 500
-2 * Loglike = 9214.648403
Burning in for 450 iterations out of 500
-2 * Loglike = 9216.436187
Burning in for 500 iterations out of 500
-2 * Loglike = 9200.962722

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9211.829316
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9208.252254
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9207.474658
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9195.303044
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9237.875695
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9215.279544
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9205.216055
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9207.194053
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9209.670491
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9228.147390
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9221.142218
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9202.025591
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9237.310950
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9201.609653
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9198.651720
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9213.676697
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9226.654472
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9197.091761
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9191.332457
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9203.168020
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9219.930342
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9190.572166
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9187.820769
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9191.319309
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9208.424962
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9193.578324
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9193.362651
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9215.800487
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9195.697241
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9217.649827
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9202.221875
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9204.702977
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9194.734325
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9214.075925
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9210.850079
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9207.673997
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9215.689888
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9216.562637
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9208.939484
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9204.992174
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9226.559863
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9215.628259
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9204.220816
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9239.698267
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9204.146545
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9205.754234
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9197.376984
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9202.261581
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9201.776501
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9199.483581
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9212.886749
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9215.503650
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9203.140920
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9197.880593
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9225.594123
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9222.127125
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9226.428207
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9197.088472
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9211.288577
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9206.916470
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9217.312858
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9227.433842
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9198.775966
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9203.263025
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9192.770836
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9204.577115
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9199.755428
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9223.824791
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9214.510921
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9212.806462
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9217.731049
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9224.753294
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9211.983406
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9226.618849
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9212.927802
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9194.951292
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9202.101969
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9203.534025
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9201.808440
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9203.452422
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9192.064532
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9206.754094
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9192.398448
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9197.709294
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9235.276558
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9213.693064
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9215.554988
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9227.538264
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9212.019772
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9215.221028
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9229.428466
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9227.351341
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9215.531289
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9208.716057
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9217.974852
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9206.761836
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9195.291703
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9203.115266
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9202.826876
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9211.036201
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9208.8
S.D.    =          11.875
S.E.M.  =         0.16793

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 1.55s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
9208.775   9148.626   60.149     9268.924   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
              Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   0.00351     0.04393   0.08     0.9364         -0.08436     0.09203    228 
standlrt    0.56291     0.01262  44.61          0  ***     0.53789     0.58746   4151 
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   0.10125     0.02083      0.06713     0.14856   3105 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.56647     0.01264      0.54210     0.59173   4686 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> if (!require(texreg)) {
+   warning("texreg package required to use screenreg() function")
+ } else {
+   screenreg(list(mymodel1, mymodel2, mymodel3), custom.model.names=c("IGLS", "Gibbs1", "Gibbs2"), groups = list("Fixed Part" = 1:2, "Level-2" = 3:3, "Level-1" = 4:4),
+    stars = numeric(0), include.nobs=FALSE, include.loglik=FALSE, include.deviance=FALSE, include.dbar=FALSE, include.dthetabar=FALSE, include.pd=FALSE, include.dic=FALSE)
+ }
Loading required package: texreg
Version:  1.39.3
Date:     2023-11-09
Author:   Philip Leifeld (University of Essex)

Consider submitting praise using the praise or praise_interactive functions.
Please cite the JSS article in your publications -- see citation("texreg").

===========================================================
                       IGLS    Gibbs1         Gibbs2       
-----------------------------------------------------------
Fixed Part                                                 
                                                           
    FP_Intercept        0.00     0.00           0.00       
                       (0.04)  [-0.08; 0.09]  [-0.08; 0.09]
    FP_standlrt         0.56     0.56 *         0.56 *     
                       (0.01)  [ 0.54; 0.59]  [ 0.54; 0.59]
Level-2                                                    
                                                           
    RP2_var_Intercept   0.09     0.10 *         0.10 *     
                       (0.02)  [ 0.06; 0.14]  [ 0.07; 0.15]
Level-1                                                    
                                                           
    RP1_var_Intercept   0.57     0.57 *         0.57 *     
                       (0.01)  [ 0.54; 0.59]  [ 0.54; 0.59]
===========================================================
* Null hypothesis value outside the confidence interval.
> rm(list = c("mymodel1", "mymodel2", "mymodel3"))
> 
> # 5.3 Using informative priors . . . . . . . . . . . . . . . . . . . . . .62
> 
> ## Informative normal prior for beta_1 Fit the model
> (mymodel4 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcMeth = list(priorParam = list(fixe = list(standlrt = c(1, 
+   0.01))))), data = tutorial))
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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2] c1092 1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500
-2 * Loglike = 9706.897984
Burning in for 100 iterations out of 500
-2 * Loglike = 9692.204734
Burning in for 150 iterations out of 500
-2 * Loglike = 9763.640621
Burning in for 200 iterations out of 500
-2 * Loglike = 9711.797029
Burning in for 250 iterations out of 500
-2 * Loglike = 9686.897247
Burning in for 300 iterations out of 500
-2 * Loglike = 9737.126921
Burning in for 350 iterations out of 500
-2 * Loglike = 9704.681847
Burning in for 400 iterations out of 500
-2 * Loglike = 9743.766532
Burning in for 450 iterations out of 500
-2 * Loglike = 9732.145050
Burning in for 500 iterations out of 500
-2 * Loglike = 9688.598828

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9663.749821
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9694.728849
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9656.802947
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9717.991675
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9668.526644
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9726.106157
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9721.929674
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9639.131284
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9650.353119
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9708.131249
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9744.590154
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9676.893208
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9770.852251
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9683.471463
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9663.269389
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9730.620727
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9739.250161
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9710.446081
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9706.962764
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9772.993395
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9675.097206
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9718.700524
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9630.493486
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9682.242045
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9694.439944
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9658.361906
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9684.709040
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9698.544988
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9689.057655
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9723.925780
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9692.764087
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9685.469009
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9696.294834
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9701.814457
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9692.436443
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9659.608527
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9668.086384
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9688.381724
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9715.856864
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9712.438810
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9690.653990
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9687.137126
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9695.711341
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9681.837627
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9692.423154
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9710.025959
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9673.628870
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9670.984715
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9695.745300
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9652.224320
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9750.598002
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9701.024944
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9639.931450
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9643.622179
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9643.021300
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9729.152943
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9738.401912
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9692.696238
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9697.932981
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9705.171569
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9675.544005
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9712.743829
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9720.615734
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9691.682590
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9672.517163
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9672.209919
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9667.287666
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9722.144632
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9705.917305
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9719.438042
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9645.434029
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9665.169280
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9708.598819
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9761.714940
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9729.405106
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9702.728690
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9694.532237
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9719.237985
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9644.240477
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9653.359411
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9710.737794
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9700.668476
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9711.283281
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9695.412798
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9661.960027
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9597.190179
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9664.232782
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9729.116340
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9659.267146
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9665.682121
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9723.944958
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9669.154975
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9707.298068
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9762.556156
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9736.773417
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9673.146156
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9650.841687
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9720.659811
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9765.776432
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9669.513096
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9693.6
S.D.    =          29.477
S.E.M.  =         0.41687

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 1.51s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
9693.590   9636.339   57.251     9750.841   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
              Coef.   Std. Err.       z   Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   0.01107     0.04037    0.27      0.784         -0.06769     0.09277    296 
standlrt    0.84138     0.00836  100.68          0  ***     0.82527     0.85794   3907 
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   0.07966     0.01695      0.05182     0.11744   2784 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.63797     0.01463      0.61068     0.66676   4059 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> sixway(mymodel4@chains[, "FP_standlrt", drop = FALSE], "beta_1")
> 
> ## Informative normal prior for beta_1
> (mymodel5 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcMeth = list(priorParam = list(fixe = list(standlrt = c(1, 
+   0.1))))), data = tutorial))
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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2] c1092 1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500
-2 * Loglike = 9207.985216
Burning in for 100 iterations out of 500
-2 * Loglike = 9216.102662
Burning in for 150 iterations out of 500
-2 * Loglike = 9243.571422
Burning in for 200 iterations out of 500
-2 * Loglike = 9215.690174
Burning in for 250 iterations out of 500
-2 * Loglike = 9198.775689
Burning in for 300 iterations out of 500
-2 * Loglike = 9228.659484
Burning in for 350 iterations out of 500
-2 * Loglike = 9222.722984
Burning in for 400 iterations out of 500
-2 * Loglike = 9217.078317
Burning in for 450 iterations out of 500
-2 * Loglike = 9217.631078
Burning in for 500 iterations out of 500
-2 * Loglike = 9201.764797

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9211.696867
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9209.129098
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9207.357335
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9197.845491
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9237.241753
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9217.281889
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9207.791478
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9205.998727
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9209.784194
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9228.940042
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9224.415860
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9202.348184
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9241.028627
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9202.402208
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9198.495657
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9215.963780
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9228.807716
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9198.929861
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9193.682076
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9222.747857
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9213.741806
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9219.166484
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9186.863587
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9189.394709
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9209.553019
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9193.626728
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9193.969821
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9215.228990
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9196.527283
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9220.350380
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9202.866631
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9205.142903
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9195.973112
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9214.431721
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9211.950249
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9208.156255
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9215.788040
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9216.395300
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9210.668706
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9207.000238
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9226.619605
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9215.860598
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9205.092949
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9237.822492
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9205.112151
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9207.222986
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9198.231301
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9202.622871
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9203.176608
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9198.846688
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9215.870073
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9216.738036
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9202.367064
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9197.129153
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9224.019985
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9224.641693
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9228.591495
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9198.709548
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9211.568953
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9208.823713
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9217.227464
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9228.497068
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9201.141041
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9204.506055
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9193.694823
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9204.619453
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9200.189500
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9220.632169
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9215.657935
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9214.690498
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9216.348847
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9224.153216
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9213.344602
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9229.784443
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9216.149315
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9206.085336
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9192.403630
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9208.110191
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9204.746366
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9200.713896
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9212.165574
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9212.740803
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9202.220604
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9229.238261
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9206.983202
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9195.685974
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9208.183131
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9232.357468
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9201.258067
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9208.538861
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9218.194152
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9224.126889
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9202.649450
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9216.105846
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9205.914714
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9218.390524
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9200.528545
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9220.807429
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9219.297955
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9211.516391
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9209.7
S.D.    =          12.007
S.E.M.  =         0.16981

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 1.58s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
9209.703   9149.746   59.958     9269.661   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
              Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   0.00337     0.04407   0.08     0.9391         -0.08278     0.09181    224 
standlrt    0.56993     0.01248  45.65          0  ***     0.54536     0.59425   3833 
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   0.09682     0.01977      0.06426     0.14053   3120 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.56626     0.01249      0.54221     0.59094   4749 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> sixway(mymodel5@chains[, "FP_standlrt", drop = FALSE], "beta_1")
> 
> # 5.4 Specifying an informative prior for a random parameter . . . . . . .65
> 
> ## Specifies an informative prior for sigma2u
> (mymodel6 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcMeth = list(priorParam = list(rp2 = list(estimate = 0.2, 
+   size = 100)))), data = tutorial))
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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2] c1092 1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500
-2 * Loglike = 9200.975178
Burning in for 100 iterations out of 500
-2 * Loglike = 9211.642830
Burning in for 150 iterations out of 500
-2 * Loglike = 9225.925396
Burning in for 200 iterations out of 500
-2 * Loglike = 9208.675167
Burning in for 250 iterations out of 500
-2 * Loglike = 9194.596250
Burning in for 300 iterations out of 500
-2 * Loglike = 9221.816208
Burning in for 350 iterations out of 500
-2 * Loglike = 9192.758664
Burning in for 400 iterations out of 500
-2 * Loglike = 9227.020649
Burning in for 450 iterations out of 500
-2 * Loglike = 9213.221638
Burning in for 500 iterations out of 500
-2 * Loglike = 9188.837295

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9208.066750
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9203.796364
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9228.740811
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9209.595852
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9231.189494
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9234.944884
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9201.655433
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9189.857434
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9199.006124
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9214.990026
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9201.639794
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9206.960474
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9212.982401
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9199.895259
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9200.381874
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9214.080645
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9221.560540
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9196.072222
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9184.763615
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9199.447618
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9209.644626
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9209.123179
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9197.197996
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9214.218123
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9210.314687
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9207.494891
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9191.292131
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9205.562004
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9201.761205
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9209.508691
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9202.325080
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9220.680576
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9211.893144
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9198.670044
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9199.001567
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9198.338025
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9193.416439
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9222.008733
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9206.954432
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9241.075478
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9224.942279
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9214.051172
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9202.108798
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9220.245636
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9213.662462
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9216.085188
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9205.581201
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9234.970242
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9203.459254
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9234.625322
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9197.310536
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9222.465143
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9196.876344
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9214.202155
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9208.363738
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9187.398774
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9199.548630
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9199.385393
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9208.001648
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9192.965523
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9187.360690
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9192.991865
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9198.193132
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9209.102644
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9216.296302
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9219.762555
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9198.225017
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9201.045935
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9206.216911
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9220.358976
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9210.605130
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9193.728209
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9196.798648
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9199.080693
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9190.899116
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9209.740332
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9190.024250
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9209.121341
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9206.355193
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9221.732453
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9197.507346
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9208.494837
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9200.569805
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9195.967277
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9215.909703
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9203.182960
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9205.478887
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9231.021055
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9225.252870
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9204.459594
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9205.334930
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9208.844135
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9205.264687
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9200.525936
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9224.002779
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9206.100319
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9210.750466
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9194.457806
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9212.002155
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9234.213885
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9206.4
S.D.    =          11.659
S.E.M.  =         0.16489

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 1.55s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
9206.419   9144.020   62.400     9268.819   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   -0.00226     0.05170  -0.04     0.9651         -0.09927     0.10185    134 
standlrt     0.56178     0.01255  44.78          0  ***     0.53741     0.58645   4366 
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   0.16292     0.01872      0.13042     0.20240   4562 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.56588     0.01287      0.54100     0.59181   4735 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> sixway(mymodel6@chains[, "RP2_var_Intercept", drop = FALSE], "sigma^2_u0")
> 
> # 5.5 Changing the random number seed and the parameter starting values  .66
> 
> ## Set starting values for random and fixed parameter estimates
> FP.b <- c(-2, 5)
> names(FP.b) <- c("FP_Intercept", "FP_standlrt")
> RP.b <- c(2, 4)
> names(RP.b) <- c("RP2_var_Intercept", "RP1_var_Intercept")
> startval <- list(FP.b = FP.b, RP.b = RP.b)
> 
> (mymodel7 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcMeth = list(burnin = 0, 
+   iterations = 500), startval = startval), data = tutorial))
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
JOIN -2 5 '_FP_b'
JOIN 2 4 '_RP_b'
ECHO 0

Echoing is ON
MCMC   0 0 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
BURNING IN...

ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 500 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 500, Stored update 50 of 500
-2 * Loglike = 9203.151396
Actual update 100 of 500, Stored update 100 of 500
-2 * Loglike = 9191.883805
Actual update 150 of 500, Stored update 150 of 500
-2 * Loglike = 9182.198346
Actual update 200 of 500, Stored update 200 of 500
-2 * Loglike = 9203.293638
Actual update 250 of 500, Stored update 250 of 500
-2 * Loglike = 9212.362161
Actual update 300 of 500, Stored update 300 of 500
-2 * Loglike = 9207.851148
Actual update 350 of 500, Stored update 350 of 500
-2 * Loglike = 9197.330971
Actual update 400 of 500, Stored update 400 of 500
-2 * Loglike = 9217.243681
Actual update 450 of 500, Stored update 450 of 500
-2 * Loglike = 9209.551458
Actual update 500 of 500, Stored update 500 of 500
-2 * Loglike = 9199.247938
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =            500
Average =          9235.0
S.D.    =          487.65
S.E.M.  =          21.809

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 0.26s 
Number of obs:  4059 (from total 4059)          Number of iter.: 500  Chains: 1  Burn-in: 0 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
9234.998   9157.571   77.427     9312.425   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.      z    Pr(>|z|)       [95% Cred.   Interval]   ESS 
Intercept   -0.41657     0.62787  -0.66       0.507         -1.77741     0.08865     1 
standlrt     0.56303     0.02685  20.97   1.181e-97  ***     0.53948     0.58684   405 
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   0.70599     1.07994      0.07201     3.53361     3 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]   ESS 
var_Intercept   0.60511     0.78702      0.54187     0.59526   406 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> rm(list = c("mymodel4", "mymodel5", "mymodel6", "mymodel7"))
> 
> ## Use different seeds
> (mymodel8 <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcMeth = list(nchains=4, seed = 1:4)),
+   data = tutorial))
MLwiN is running, please wait......
starting worker pid=7984 on localhost:11856 at 14:23:11.602
starting worker pid=1808 on localhost:11856 at 14:23:11.623
starting worker pid=6680 on localhost:11856 at 14:23:11.626
starting worker pid=11420 on localhost:11856 at 14:23:11.626
Loading required package: 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.11/
To change this use: options(MLwiN_path="<path to MLwiN>")
loaded R2MLwiN and set parent environment
Loading required package: 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.11/
To change this use: options(MLwiN_path="<path to MLwiN>")
loaded R2MLwiN and set parent environment
Loading required package: 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.11/
To change this use: options(MLwiN_path="<path to MLwiN>")
loaded R2MLwiN and set parent environment
Loading required package: 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.11/
To change this use: options(MLwiN_path="<path to MLwiN>")
loaded R2MLwiN and set parent environment
/nogui option ignored
/nogui option ignored
/nogui option ignored
/nogui option ignored
ECHO    0

ECHO    0

ECHO    0

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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500

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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500

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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
BURNING IN...

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

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500
Burning in for 50 iterations out of 500
-2 * Loglike = 9206.193266
Burning in for 100 iterations out of 500
-2 * Loglike = 9200.609099
Burning in for 100 iterations out of 500
-2 * Loglike = 9216.262176
Burning in for 150 iterations out of 500
-2 * Loglike = 9203.175563
Burning in for 100 iterations out of 500
-2 * Loglike = 9215.293382
Burning in for 100 iterations out of 500
-2 * Loglike = 9200.366607
Burning in for 150 iterations out of 500
-2 * Loglike = 9241.057582
Burning in for 200 iterations out of 500
-2 * Loglike = 9209.601793
Burning in for 200 iterations out of 500
-2 * Loglike = 9214.675885
Burning in for 250 iterations out of 500
-2 * Loglike = 9210.592682
Burning in for 150 iterations out of 500
-2 * Loglike = 9204.601702
Burning in for 150 iterations out of 500
-2 * Loglike = 9191.581836
Burning in for 250 iterations out of 500
-2 * Loglike = 9198.255145
Burning in for 300 iterations out of 500
-2 * Loglike = 9209.632813
Burning in for 200 iterations out of 500
-2 * Loglike = 9212.840404
Burning in for 300 iterations out of 500
-2 * Loglike = 9200.368178
Burning in for 200 iterations out of 500
-2 * Loglike = 9227.186873
Burning in for 350 iterations out of 500
-2 * Loglike = 9209.689105
Burning in for 350 iterations out of 500
-2 * Loglike = 9222.091977
Burning in for 400 iterations out of 500
-2 * Loglike = 9209.324694
Burning in for 250 iterations out of 500
-2 * Loglike = 9211.614600
Burning in for 250 iterations out of 500
-2 * Loglike = 9226.679111
Burning in for 400 iterations out of 500
-2 * Loglike = 9214.630238
Burning in for 450 iterations out of 500
-2 * Loglike = 9191.232999
Burning in for 450 iterations out of 500
-2 * Loglike = 9216.522216
Burning in for 500 iterations out of 500
-2 * Loglike = 9222.120004
Burning in for 300 iterations out of 500
-2 * Loglike = 9209.709202
Burning in for 300 iterations out of 500
-2 * Loglike = 9200.733163
Burning in for 500 iterations out of 500
-2 * Loglike = 9201.148603

-2 * Loglike = 9215.837665

-2 * Loglike = 9202.902457
Burning in for 350 iterations out of 500
-2 * Loglike = 9189.314825
Burning in for 350 iterations out of 500
ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9212.190791
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9211.953489
Burning in for 400 iterations out of 500
-2 * Loglike = 9210.065393
Burning in for 400 iterations out of 500
-2 * Loglike = 9196.808263
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9208.377297
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9189.156727
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9207.886882
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9213.270899
Burning in for 450 iterations out of 500
-2 * Loglike = 9201.240831
Burning in for 450 iterations out of 500
-2 * Loglike = 9208.899478
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9195.385976
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9215.653768
Burning in for 500 iterations out of 500
-2 * Loglike = 9204.605317
Burning in for 500 iterations out of 500
-2 * Loglike = 9191.635017
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9238.483062
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9201.279205
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9215.621742
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9204.115454

-2 * Loglike = 9205.650695

-2 * Loglike = 9205.716499
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9220.281051
Actual update 350 of 5000, Stored update 350 of 5000
ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9207.158850
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9215.304042
Actual update 400 of 5000, Stored update 400 of 5000
ERAS G30
LINK 0 G30
ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 9210.248427
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9214.890801
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9211.426484
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9193.251090
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9228.495732
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9224.335795
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9221.683303
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9224.268752
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9199.987992
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9198.799074
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9202.166694
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9202.589985
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9192.092122
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9237.981658
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9217.924542
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9204.102958
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9201.878383
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9225.593562
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9208.559151
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9208.542998
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9198.726654
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9213.648448
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9214.104227
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9221.025874
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9204.413553
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9198.579514
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9227.156460
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9203.928314
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9204.652958
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9212.428999
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9197.279082
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9208.752577
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9191.816267
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9222.186272
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9208.585127
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9205.396874
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9219.178731
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9213.319899
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9216.935902
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9214.163387
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9231.834103
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9223.741372
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9206.033213
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9218.308723
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9195.828276
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9213.127568
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9188.214048
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9202.604628
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9202.177193
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9211.035687
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9188.278770
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9216.279711
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9197.193870
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9215.285002
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9208.819913
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9213.478807
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9193.939811
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9203.028947
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9198.177466
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9192.727608
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9193.651680
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9228.308587
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9196.046922
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9215.606479
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9192.249867
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9225.084081
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9195.925780
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9200.012106
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9203.221257
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9224.493065
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9218.422125
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9194.758430
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9202.244183
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9214.794631
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9188.863509
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9228.144914
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9204.907738
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9203.041639
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9221.116965
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9194.883902
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9215.281319
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9217.912243
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9214.140517
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9199.150201
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9208.641266
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9231.465362
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9211.397609
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9204.619839
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9201.489802
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9195.204182
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9207.941176
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9216.158993
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9195.266954
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9203.345896
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9196.692674
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9216.821254
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9222.603976
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9212.994917
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9209.099646
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9210.117979
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9217.019911
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9205.356036
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9215.286472
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9192.114800
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9226.820798
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9208.780187
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9199.573660
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9215.889068
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9212.555828
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9194.305711
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9214.394137
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9204.485567
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9198.974980
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9239.773102
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9211.297370
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9214.691676
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9207.058767
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9204.613020
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9220.424840
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9212.707453
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9226.133690
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9205.832435
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9195.326250
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9197.787494
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9201.412595
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9209.784406
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9211.051554
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9202.635017
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9205.740712
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9202.243753
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9208.160610
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9224.934371
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9211.914868
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9199.686027
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9192.633871
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9205.291307
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9213.435702
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9222.573230
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9199.001420
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9215.862131
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9214.099955
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9213.582699
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9218.156790
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9203.440604
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9222.213115
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9197.917390
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9216.004168
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9228.129572
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9210.119962
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9225.824248
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9210.342063
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9222.536523
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9222.623770
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9213.070127
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9195.549581
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9226.897380
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9222.209846
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9198.414226
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9217.234648
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9197.383677
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9206.060326
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9211.390320
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9211.806227
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9198.318630
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9213.603963
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9207.286579
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9197.638023
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9205.528701
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9217.534200
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9197.216009
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9202.379355
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9227.872396
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9206.138262
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9217.921427
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9216.644577
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9198.935844
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9189.073511
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9219.249197
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9203.578259
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9213.028813
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9209.768881
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9193.257434
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9204.566824
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9192.738284
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9198.810339
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9204.815469
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9203.448633
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9205.260076
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9208.635545
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9200.291358
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9212.704366
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9203.526438
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9225.161893
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9219.578031
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9224.027647
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9186.129983
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9186.474187
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9210.092157
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9214.657834
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9207.716548
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9201.784637
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9195.423341
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9213.144417
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9206.248565
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9203.962725
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9217.997575
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9219.725147
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9221.103720
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9190.546147
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9225.109666
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9211.861240
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9209.288321
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9212.724490
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9212.063094
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9211.904328
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9212.267938
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9210.761516
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9226.746381
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9223.058748
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9194.648832
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9206.893648
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9214.159422
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9198.054775
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9208.042101
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9195.049840
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9191.361316
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9204.966122
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9228.516251
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9203.875559
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9185.020836
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9190.735716
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9207.086769
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9199.652712
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9221.708398
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9206.324394
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9203.808438
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9203.701447
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9210.376876
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9205.512200
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9228.244909
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9192.747742
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9216.264655
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9201.017143
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9191.264804
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9208.750673
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9201.098674
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9211.234780
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9217.040436
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9209.774569
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9212.969987
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9211.765768
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9208.811215
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9197.541074
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9214.543462
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9200.968089
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9208.504279
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9216.500091
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9187.587677
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9229.231239
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9205.264187
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9211.461829
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9207.384781
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9207.636298
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9206.919224
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9200.327046
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9188.001440
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9197.507235
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9231.068743
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9207.836886
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9208.527365
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9208.445006
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9198.221864
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9214.700301
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9218.017798
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9231.667774
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9192.182015
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9231.720113
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9209.041996
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9201.470719
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9216.063459
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9213.497604
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9219.625871
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9208.425209
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9210.166699
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9212.367623
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9220.467826
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9216.996959
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9205.701131
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9200.358770
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9206.364287
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9225.405793
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9257.194444
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9203.040085
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9208.178207
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9201.746288
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9211.360067
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9208.216150
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9191.717666
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9212.656436
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9199.191002
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9217.250158
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9195.952868
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9202.976127
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9195.831784
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9205.845134
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9204.239398
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9218.846042
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9201.925078
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9212.453348
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9200.637551
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9194.223636
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9202.725373
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9207.805347
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9210.617870
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9219.146880
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9241.694771
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9219.997563
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9208.496500
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9215.841845
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9216.028873
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9198.847084
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9201.561854
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9212.109834
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9209.1
S.D.    =          11.916
S.E.M.  =         0.16852

ECHO 0
Execution completed
-2 * Loglike = 9204.661285
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9208.494945
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9209.556744
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9202.001472
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9205.899950
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9195.339542
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9212.138475
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9193.782401
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9205.467869
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9225.827063
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9208.704966
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9193.584281
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9231.051021
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9208.9
S.D.    =          11.965
S.E.M.  =         0.16921

ECHO 0
Execution completed
-2 * Loglike = 9211.972280
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9214.807017
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9226.847535
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9190.821237
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9218.677910
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9219.171539
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9208.428869
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9199.425383
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9200.985839
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9216.844055
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9216.507349
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9223.620521
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9185.466396
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9226.140355
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9212.550413
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9221.365302
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9216.204867
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9207.447260
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9216.152979
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9223.111260
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9190.602337
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9204.745394
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9205.869858
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9192.696656
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9205.874076
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9216.944129
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9219.038505
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9214.847221
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9193.691572
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9228.292624
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9217.473815
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9242.850710
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9216.478442
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9195.406259
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9203.793999
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9208.397974
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9228.250155
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9200.803309
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9203.297389
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9199.768306
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9209.754549
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9230.761641
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9216.845179
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9211.206432
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9219.880589
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9198.137286
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9209.635081
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9208.9
S.D.    =          11.754
S.E.M.  =         0.16623

ECHO 0
Execution completed
-2 * Loglike = 9206.964668
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9208.598434
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9209.3
S.D.    =          11.903
S.E.M.  =         0.16833

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 14.97s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 4  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
9209.034   9149.120   59.915     9268.949   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | school) + (1 | student)
Level 2: school     Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
              Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]     ESS 
Intercept   0.00407     0.04216   0.10      0.923         -0.07969     0.08818     825 
standlrt    0.56320     0.01259  44.72          0  ***     0.53849     0.58803   16207 
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   0.09690     0.02013      0.06434     0.14226   11456 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]     ESS 
var_Intercept   0.56618     0.01250      0.54232     0.59091   19472 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ctable <- round(as.data.frame(lapply(mymodel8@chains[, c(names(mymodel8@FP), names(mymodel8@RP))], colMeans)), 3)
> colnames(ctable) <- c("Seed1", "Seed2", "Seed3", "Seed4")
> print(ctable)
                  Seed1 Seed2 Seed3 Seed4
FP_Intercept      0.003 0.005 0.001 0.007
FP_standlrt       0.563 0.563 0.563 0.563
RP2_var_Intercept 0.097 0.096 0.097 0.097
RP1_var_Intercept 0.566 0.566 0.566 0.566
> rm(list = "mymodel8")
> 
> # 5.6 Improving the speed of MCMC Estimation . . . . . . . . . . . . . . .69
> 
> # Chapter learning outcomes . . . . . . . . . . . . . . . . . . . . . . . 70
> 
> 
> 
> 
> 
> ############################################################################
> 
> proc.time()
   user  system elapsed 
   7.43    1.29   31.81