R version 4.3.2 (2023-10-31 ucrt) -- "Eye Holes"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (64-bit)

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

  Natural language support but running in an English locale

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

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

> ############################################################################
> #     MLwiN MCMC Manual
> #
> # 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.09/
To change this use: options(MLwiN_path="<path to MLwiN>")
> # MLwiN folder
> mlwin <- getOption("MLwiN_path")
> while (!file.access(mlwin, mode = 1) == 0) {
+   cat("Please specify the root MLwiN folder or the full path to the MLwiN executable:\n")
+   mlwin <- scan(what = character(0), sep = "\n")
+   mlwin <- gsub("\\", "/", mlwin, fixed = TRUE)
+ }
> options(MLwiN_path = mlwin)
> 
> # User's input if necessary
> 
> ## Read 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.09)  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.09)  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.44s 
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.09)  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.52s 
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.09)  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.09)  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.52s 
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.09)  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.49s 
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.09)  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.24s 
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=8680 on localhost:11440 at 11:28:33.433
starting worker pid=7616 on localhost:11440 at 11:28:33.425
starting worker pid=11708 on localhost:11440 at 11:28:33.402
starting worker pid=6932 on localhost:11440 at 11:28:33.421
Loading required package: R2MLwiN
Loading required package: R2MLwiN
Loading required package: R2MLwiN
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)

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)

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},
  }


To cite R2MLwiN in publications use:

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

A BibTeX entry for LaTeX users is

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

The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.09/
The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.09/
To change this use: options(MLwiN_path="<path to MLwiN>")
To change this use: options(MLwiN_path="<path to MLwiN>")
loaded R2MLwiN and set parent environment
loaded R2MLwiN and set parent environment

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

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

To cite R2MLwiN in publications use:

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

A BibTeX entry for LaTeX users is

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

The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.09/
To change this use: options(MLwiN_path="<path to MLwiN>")
loaded R2MLwiN and set parent environment

To cite R2MLwiN in publications use:

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

A BibTeX entry for LaTeX users is

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

The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.09/
To change this use: options(MLwiN_path="<path to MLwiN>")
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...

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

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 = 9200.609099
Burning in for 100 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 = 9216.262176
Burning in for 150 iterations out of 500
-2 * Loglike = 9200.366607
Burning in for 150 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 = 9241.057582
Burning in for 200 iterations out of 500
-2 * Loglike = 9209.601793
Burning in for 200 iterations out of 500
-2 * Loglike = 9209.632813
Burning in for 200 iterations out of 500
-2 * Loglike = 9200.368178
Burning in for 200 iterations out of 500
-2 * Loglike = 9214.675885
Burning in for 250 iterations out of 500
-2 * Loglike = 9191.581836
Burning in for 250 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 = 9198.255145
Burning in for 300 iterations out of 500
-2 * Loglike = 9212.840404
Burning in for 300 iterations out of 500
-2 * Loglike = 9209.709202
Burning in for 300 iterations out of 500
-2 * Loglike = 9222.120004
Burning in for 300 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 = 9189.314825
Burning in for 350 iterations out of 500
-2 * Loglike = 9202.902457
Burning in for 350 iterations out of 500
-2 * Loglike = 9222.091977
Burning in for 400 iterations out of 500
-2 * Loglike = 9226.679111
Burning in for 400 iterations out of 500
-2 * Loglike = 9210.065393
Burning in for 400 iterations out of 500
-2 * Loglike = 9211.953489
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 = 9201.240831
Burning in for 450 iterations out of 500
-2 * Loglike = 9213.270899
Burning in for 450 iterations out of 500
-2 * Loglike = 9216.522216
Burning in for 500 iterations out of 500
-2 * Loglike = 9200.733163
Burning in for 500 iterations out of 500
-2 * Loglike = 9204.605317
Burning in for 500 iterations out of 500
-2 * Loglike = 9215.653768
Burning in for 500 iterations out of 500
-2 * Loglike = 9201.148603

-2 * Loglike = 9215.837665

-2 * Loglike = 9205.650695

-2 * Loglike = 9204.115454

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
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 = 9211.426484
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 = 9196.808263
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9199.987992
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 = 9193.251090
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9192.092122
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 = 9189.156727
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 = 9225.593562
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 = 9208.899478
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9217.924542
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9221.025874
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 = 9191.635017
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 = 9205.716499
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9204.652958
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9201.279205
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 = 9207.158850
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9220.281051
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9208.585127
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9210.248427
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9212.428999
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9215.304042
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9223.741372
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9222.186272
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9228.495732
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9214.890801
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9202.604628
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9213.319899
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9221.683303
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9216.279711
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9231.834103
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9224.335795
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9203.028947
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 = 9195.828276
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9224.268752
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9196.046922
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 = 9202.177193
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9237.981658
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9200.012106
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9204.102958
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9197.193870
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 = 9214.794631
Actual update 750 of 5000, Stored update 750 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 = 9221.116965
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9198.177466
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9213.648448
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9192.249867
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9199.150201
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9214.104227
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9198.579514
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9203.221257
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9204.619839
Actual update 900 of 5000, Stored update 900 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 = 9203.345896
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9188.863509
Actual update 850 of 5000, Stored update 850 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 = 9212.994917
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9205.396874
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9215.281319
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 = 9192.114800
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9216.935902
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9208.641266
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 = 9194.305711
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9206.033213
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 = 9201.489802
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9211.297370
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9213.127568
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 = 9196.692674
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9212.707453
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9211.035687
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 = 9217.019911
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9209.784406
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9215.285002
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.780187
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9213.478807
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9208.160610
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9208.819913
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9214.394137
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9192.727608
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9205.291307
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9193.939811
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9207.058767
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9228.308587
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9213.582699
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9193.651680
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9226.133690
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9225.084081
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9216.004168
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9215.606479
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9211.051554
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9224.493065
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9222.536523
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9195.925780
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9211.914868
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9194.758430
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 = 9222.573230
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9198.414226
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9228.144914
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 = 9218.156790
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9211.806227
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 = 9210.119962
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9205.528701
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 = 9213.070127
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9206.138262
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9194.883902
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9217.912243
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9217.234648
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9219.249197
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9231.465362
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 = 9204.566824
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9198.318630
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 = 9195.204182
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9203.448633
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9197.216009
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 = 9212.704366
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9216.644577
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9195.266954
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9216.158993
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9224.027647
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9209.768881
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9222.603976
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9216.821254
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9210.092157
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9198.810339
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9210.117979
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9209.099646
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9195.423341
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9208.635545
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9205.356036
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9219.725147
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9215.286472
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9225.161893
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9226.820798
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9211.861240
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9186.129983
Actual update 2050 of 5000, Stored update 2050 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 = 9211.904328
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9207.716548
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9212.555828
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9204.485567
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9223.058748
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9206.248565
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9198.974980
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9198.054775
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9239.773102
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9221.103720
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9214.691676
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9191.361316
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9204.613020
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9209.288321
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9220.424840
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9185.020836
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9205.832435
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9212.267938
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9195.326250
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9221.708398
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9197.787494
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9194.648832
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9201.412595
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9210.376876
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9202.635017
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9208.042101
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9205.740712
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9216.264655
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9202.243753
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9228.516251
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9224.934371
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9201.098674
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9199.686027
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9207.086769
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9192.633871
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9212.969987
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9213.435702
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9199.001420
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9203.808438
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9214.543462
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9215.862131
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9187.587677
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9203.440604
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9214.099955
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9228.244909
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9207.384781
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9197.917390
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9188.001440
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9222.213115
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9191.264804
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9225.824248
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9208.527365
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9222.623770
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9228.129572
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9217.040436
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9214.700301
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9226.897380
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9231.720113
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9197.383677
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9210.342063
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9208.811215
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9213.497604
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9211.390320
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9212.367623
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9195.549581
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9208.504279
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9207.286579
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9200.358770
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9217.534200
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9222.209846
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9205.264187
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9203.040085
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9227.872396
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9208.216150
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9206.060326
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9198.935844
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9206.919224
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9217.250158
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9203.578259
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9213.603963
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9231.068743
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9205.845134
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9193.257434
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9202.725373
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9204.815469
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9197.638023
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9198.221864
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9219.997563
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9200.291358
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9201.561854
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9192.182015
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9202.379355
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9219.578031
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9204.661285
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9216.063459
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9217.921427
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9202.001472
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9210.166699
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9214.657834
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9212.138475
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9205.701131
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 = 9205.467869
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9213.144417
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9257.194444
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9193.584281
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9213.028813
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9217.997575
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9211.360067
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9214.807017
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9225.109666
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9199.191002
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9192.738284
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9190.821237
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9205.260076
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9212.063094
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9219.171539
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9195.831784
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9203.526438
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9226.746381
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9199.425383
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9201.925078
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9186.474187
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9214.159422
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9194.223636
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9201.784637
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9216.844055
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9204.966122
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9210.617870
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9203.962725
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9190.735716
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9208.496500
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9223.620521
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9190.546147
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9206.324394
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9198.847084
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9226.140355
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9212.724490
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9205.512200
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9209.556744
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9221.365302
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9210.761516
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9201.017143
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9207.447260
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9206.893648
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9195.339542
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9211.234780
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9223.111260
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9195.049840
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9211.765768
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9193.782401
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9204.745394
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9203.875559
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9200.968089
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9225.827063
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9192.696656
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9199.652712
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9229.231239
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9208.704966
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9216.944129
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9203.701447
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9207.636298
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9211.972280
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9214.847221
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9192.747742
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9197.507235
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9226.847535
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9228.292624
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9208.750673
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9208.445006
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9218.677910
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9209.774569
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9242.850710
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9231.667774
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9208.428869
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9197.541074
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9195.406259
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9201.470719
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9200.985839
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9216.500091
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9208.397974
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9208.425209
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9216.507349
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9211.461829
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9228.250155
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9216.996959
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9185.466396
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9200.327046
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9203.297389
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9225.405793
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9212.550413
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9207.836886
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9209.754549
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9218.017798
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9216.204867
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9216.845179
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9201.746288
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9209.041996
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9216.152979
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9219.880589
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9212.656436
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9219.625871
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9190.602337
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9198.137286
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9202.976127
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9220.467826
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9205.869858
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9206.964668
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9218.846042
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9206.364287
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9205.874076
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9200.637551
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9208.178207
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9219.038505
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9219.146880
Actual update 4950 of 5000, Stored update 4950 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
-2 * Loglike = 9191.717666
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9215.841845
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9193.691572
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9195.952868
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9217.473815
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9204.239398
Actual update 4700 of 5000, Stored update 4700 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 = 9216.478442
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9212.453348
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9203.793999
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9207.805347
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9200.803309
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9241.694771
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9199.768306
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9216.028873
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9230.761641
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9208.494945
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9211.206432
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9205.899950
Actual update 5000 of 5000, Stored update 5000 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 = 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


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  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 : 7.29s 
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 
   5.51    1.25   22.62