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

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

  Natural language support but running in an English locale

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

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

> ############################################################################
> #     MLwiN MCMC Manual
> #
> # 22  Using the Structured MVN framework for models . . . . . . . . . . .341
> #
> #     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/
> ############################################################################
> 
> # 22.1 MCMC theory for Structured MVN models . . . . . . . . . . . . . . 341
> 
> # 22.2 Using the SMVN framework in practice . . . . . . . . . . . . . . .344
> 
> library(R2MLwiN)

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

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

To cite R2MLwiN in publications use:

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

A BibTeX entry for LaTeX users is

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

The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.05/
To change this use: options(MLwiN_path="<path to MLwiN>")
> # MLwiN folder
> mlwin <- getOption("MLwiN_path")
> while (!file.access(mlwin, mode = 1) == 0) {
+   cat("Please specify the root MLwiN folder or the full path to the MLwiN executable:\n")
+   mlwin <- scan(what = character(0), sep = "\n")
+   mlwin <- gsub("\\", "/", mlwin, fixed = TRUE)
+ }
> options(MLwiN_path = mlwin)
> 
> # User's input if necessary
> 
> ## Read tutorial data
> data(tutorial, package = "R2MLwiN")
> 
> ## Define the model
> 
> (mymodel <- runMLwiN(normexam ~ 1 + (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

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 for 50 iterations out of 500
-2 * Loglike = 10847.820270
Burning in for 100 iterations out of 500
-2 * Loglike = 10860.697552
Burning in for 150 iterations out of 500
-2 * Loglike = 10849.534371
Burning in for 200 iterations out of 500
-2 * Loglike = 10850.670040
Burning in for 250 iterations out of 500
-2 * Loglike = 10830.922843
Burning in for 300 iterations out of 500
-2 * Loglike = 10851.202135
Burning in for 350 iterations out of 500
-2 * Loglike = 10831.186056
Burning in for 400 iterations out of 500
-2 * Loglike = 10838.507313
Burning in for 450 iterations out of 500
-2 * Loglike = 10841.701201
Burning in for 500 iterations out of 500
-2 * Loglike = 10848.918096

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 = 10822.788614
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 10865.181513
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 10837.691070
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 10880.285861
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 10833.110430
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 10844.265345
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 10858.607350
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 10831.891795
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 10856.861249
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 10824.758762
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 10845.858542
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 10844.826830
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 10862.143643
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 10842.581836
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 10834.259843
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 10851.536666
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 10852.274906
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 10859.691363
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 10833.246570
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 10833.952931
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 10857.442861
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 10860.615023
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 10852.796682
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 10867.313991
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 10846.393868
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 10863.772440
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 10851.851137
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 10849.506033
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 10865.054268
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 10841.493048
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 10837.389773
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 10842.953887
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 10844.488177
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 10844.711989
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 10844.804713
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 10839.594112
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 10849.154572
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 10860.437248
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 10868.448213
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 10835.531579
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 10854.971270
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 10864.444900
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 10850.386079
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 10837.019473
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 10838.210265
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 10890.708656
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 10835.721477
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 10860.325960
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 10882.421342
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 10873.131199
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 10832.176706
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 10858.781475
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 10842.603451
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 10845.443216
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 10846.696437
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 10852.980567
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 10870.685600
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 10865.558939
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 10839.289750
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 10856.743050
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 10846.803279
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 10845.939962
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 10859.518286
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 10844.083664
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 10853.951421
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 10827.939256
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 10855.188332
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 10848.802296
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 10857.682950
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 10847.743737
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 10854.692642
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 10845.608599
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 10853.097920
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 10846.784688
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 10847.376478
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 10851.908312
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 10841.962641
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 10844.944032
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 10866.722823
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 10838.317977
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 10826.617080
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 10844.415737
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 10846.888156
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 10859.615806
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 10834.763065
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 10852.353281
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 10856.812512
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 10852.181962
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 10848.798576
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 10864.858735
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 10844.724386
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 10840.243623
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 10834.740697
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 10848.424253
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 10854.432505
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 10829.687734
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 10840.477795
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 10873.693521
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 10849.245570
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 10856.432891
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          10850.
S.D.    =          11.751
S.E.M.  =         0.16618

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.37s 
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
10849.609  10789.729  59.879     10909.488  
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + (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.01474     0.05650  -0.26     0.7942         -0.12039     0.09603   172 
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.17775     0.03642      0.11920     0.26236   3278 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.84847     0.01900      0.81147     0.88711   5000 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> 
> ## Structured MVN
> 
> (mymodel <- runMLwiN(normexam ~ 1 + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcOptions = list(smvn = 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

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
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 11012.863626
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 11012.663456
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 11012.578448
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 11014.280881
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 11012.328840
Adapting for 600 iterations (Maximum 5000)
-2 * Loglike = 11012.580639
Adapting for 700 iterations (Maximum 5000)
-2 * Loglike = 11013.041361
Adapting for 800 iterations (Maximum 5000)
-2 * Loglike = 11013.849121
Adapting for 900 iterations (Maximum 5000)
-2 * Loglike = 11012.693450
Adapting for 1000 iterations (Maximum 5000)
-2 * Loglike = 11014.530885
Adapting finished and took 1100 iterations
Adapting took 1100 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 11012.148070
Burning in for 100 iterations out of 500
-2 * Loglike = 11014.851665
Burning in for 150 iterations out of 500
-2 * Loglike = 11012.170059
Burning in for 200 iterations out of 500
-2 * Loglike = 11012.821938
Burning in for 250 iterations out of 500
-2 * Loglike = 11012.146179
Burning in for 300 iterations out of 500
-2 * Loglike = 11010.688811
Burning in for 350 iterations out of 500
-2 * Loglike = 11013.354383
Burning in for 400 iterations out of 500
-2 * Loglike = 11023.641698
Burning in for 450 iterations out of 500
-2 * Loglike = 11017.075559
Burning in for 500 iterations out of 500
-2 * Loglike = 11017.567066

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 = 11020.689846
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 11012.974385
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 11014.963344
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 11012.083844
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 11012.343608
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 11010.925270
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 11013.441563
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 11017.649880
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 11018.579112
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 11011.730889
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 11010.947546
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 11013.493241
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 11014.037641
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 11011.975148
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 11014.024250
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 11016.316016
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 11014.257367
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 11011.006754
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 11016.284807
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 11014.241245
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 11011.247204
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 11014.376159
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 11013.413291
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 11015.471457
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 11020.668214
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 11012.976457
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 11014.201102
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 11011.683040
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 11014.977871
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 11014.027952
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 11013.750331
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 11014.221547
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 11012.383737
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 11011.745350
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 11017.110765
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 11011.876114
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 11012.017375
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 11010.963100
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 11011.568996
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 11016.520281
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 11011.996574
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 11014.440404
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 11010.912783
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 11014.053047
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 11014.942963
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 11013.048878
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 11012.498652
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 11014.501208
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 11013.223714
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 11012.704304
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 11016.789300
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 11013.953733
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 11013.186484
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 11013.471564
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 11014.987223
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 11018.847044
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 11013.566918
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 11012.179517
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 11013.702834
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 11015.560406
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 11010.850516
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 11021.633081
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 11011.349557
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 11013.858916
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 11012.510326
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 11013.956546
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 11013.774374
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 11015.397246
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 11012.921941
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 11012.587118
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 11015.589959
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 11011.366234
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 11012.350562
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 11013.608037
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 11018.625548
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 11012.947028
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 11011.452294
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 11012.833282
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 11013.199609
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 11017.776365
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 11016.668728
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 11012.336122
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 11013.916875
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 11012.945658
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 11018.289892
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 11013.128028
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 11011.683053
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 11011.040994
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 11011.578062
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 11014.806784
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 11013.981855
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 11015.870830
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 11012.581314
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 11013.127389
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 11015.877685
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 11012.720185
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 11012.931875
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 11016.195554
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 11014.582923
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 11012.080539
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          11014.
S.D.    =          2.3388
S.E.M.  =        0.033076

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.33s 
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
11013.694  11010.848  2.847      11016.541  
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + (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.01415     0.05499  -0.26      0.797         -0.12492     0.09344   1048 
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.18392     0.03572      0.12385     0.26388   1009 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.84898     0.01910      0.81177     0.88752   1108 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 22.3 Model Comparison and structured MVN models . . . . . . . . . . . .349
> 
> ## Define the model
> 
> ## Gibbs
> 
> (mymodel <- 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 for 50 iterations out of 500
-2 * Loglike = 9204.262556
Burning in for 100 iterations out of 500
-2 * Loglike = 9225.063419
Burning in for 150 iterations out of 500
-2 * Loglike = 9205.620423
Burning in for 200 iterations out of 500
-2 * Loglike = 9207.286175
Burning in for 250 iterations out of 500
-2 * Loglike = 9189.882541
Burning in for 300 iterations out of 500
-2 * Loglike = 9213.131781
Burning in for 350 iterations out of 500
-2 * Loglike = 9191.967839
Burning in for 400 iterations out of 500
-2 * Loglike = 9200.072789
Burning in for 450 iterations out of 500
-2 * Loglike = 9197.978878
Burning in for 500 iterations out of 500
-2 * Loglike = 9206.312775

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 = 9182.514808
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9224.799767
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9192.533264
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9235.302728
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9192.588497
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9200.798346
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9215.971258
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9191.471053
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9216.792403
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9185.644474
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9203.752807
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9201.223665
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9215.804872
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9204.519349
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9189.609771
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9213.416013
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9211.107554
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9219.303837
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9195.273065
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9191.690604
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9215.236489
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9218.185496
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9210.517749
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9225.032838
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9202.460499
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9225.333397
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9208.977623
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9210.427472
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9225.575111
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9202.428957
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9200.886427
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9198.521334
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9207.777166
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9198.900296
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9208.265053
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9203.987061
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9207.776243
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9215.497296
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9224.467108
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9194.146736
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9212.845807
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9224.944089
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9213.612843
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9198.265355
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9195.222358
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9250.643212
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9194.598687
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9214.421500
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9240.315838
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9235.012216
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9188.825586
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9220.910444
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9198.293627
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9205.222052
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9209.264188
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9211.039764
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9230.996469
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9224.580436
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9199.499799
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9215.292584
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9207.040587
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9202.604250
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9216.541549
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9205.256474
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9211.029403
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9186.699312
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9210.823399
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9207.867998
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9219.418322
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9206.509889
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9211.940598
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9209.426437
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9217.104187
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9206.433860
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9208.530610
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9212.678566
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9201.101902
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9202.943453
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9226.646620
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9199.272707
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9185.496662
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9204.693995
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9205.338782
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9218.098365
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9192.878156
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9211.471250
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9217.216161
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9214.835575
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9208.972779
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9223.102763
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9203.956845
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9201.960121
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9194.101999
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9207.476225
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9215.003009
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9186.582930
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9202.744286
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9228.521835
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9210.720494
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9214.075214
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9208.8
S.D.    =          12.005
S.E.M.  =         0.16978

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.38s 
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.817   9148.971   59.846     9268.664   
--------------------------------------------------------------------------------------------------- 
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.00133     0.04216   0.03     0.9748         -0.07762     0.08519    208 
standlrt    0.56331     0.01254  44.93          0  ***     0.53892     0.58799   4017 
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.09732     0.02049      0.06441     0.14454   2981 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.56634     0.01269      0.54175     0.59197   5000 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ## SMCMC
> 
> (mymodel <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcOptions = list(smcm = 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 for 50 iterations out of 500
-2 * Loglike = 9200.709700
Burning in for 100 iterations out of 500
-2 * Loglike = 9201.435737
Burning in for 150 iterations out of 500
-2 * Loglike = 9209.262379
Burning in for 200 iterations out of 500
-2 * Loglike = 9209.226410
Burning in for 250 iterations out of 500
-2 * Loglike = 9206.491463
Burning in for 300 iterations out of 500
-2 * Loglike = 9208.492028
Burning in for 350 iterations out of 500
-2 * Loglike = 9209.635411
Burning in for 400 iterations out of 500
-2 * Loglike = 9210.443709
Burning in for 450 iterations out of 500
-2 * Loglike = 9235.721253
Burning in for 500 iterations out of 500
-2 * Loglike = 9203.839599

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 = 9230.332549
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9202.341831
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9232.953634
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9215.606536
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9208.176270
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9218.484728
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9219.124827
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9198.457795
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9209.506110
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9192.683097
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9209.708961
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9224.596985
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9249.224240
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9188.090515
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9208.281629
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9217.868590
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9224.766120
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9204.701012
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9223.380787
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9206.882543
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9213.963221
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9211.473028
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9214.147432
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9195.994075
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9190.030516
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9220.520423
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9215.829728
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9225.875679
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9206.247760
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9233.564768
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9202.626280
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9207.437502
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9213.058495
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9197.621676
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9186.292236
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9197.477940
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9209.525959
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9186.879900
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9203.214335
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9191.904673
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9216.848175
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9220.738003
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9194.424145
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9242.598093
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9229.582570
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9204.395010
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9219.687213
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9205.258652
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9232.778476
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9227.551018
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9219.618475
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9226.394056
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9204.207297
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9199.645637
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9195.691836
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9202.292002
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9225.527835
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9201.874058
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9224.988804
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9240.038427
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9220.056620
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9213.839309
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9220.803399
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9202.668798
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9196.566125
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9225.443709
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9186.522076
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9217.861308
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9202.261338
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9205.562007
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9219.727788
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9198.388846
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9246.979029
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9194.679872
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9216.034785
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9231.797241
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9217.960398
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9227.059506
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9211.872387
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9230.500703
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9228.345856
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9216.840453
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9224.178411
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9219.965704
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9209.807301
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9206.488689
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9208.714671
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9216.697776
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9205.877082
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9193.365786
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9189.222370
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9212.293205
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9209.165921
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9216.354484
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9216.972597
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9206.307658
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9194.518261
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9219.047699
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9207.300055
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9207.704573
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9208.9
S.D.    =          11.773
S.E.M.  =         0.16650

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.4s 
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.861   9149.139   59.722     9268.583   
--------------------------------------------------------------------------------------------------- 
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.00344     0.04122   0.08     0.9336         -0.07884     0.08547   4919 
standlrt    0.56373     0.01246  45.25          0  ***     0.53935     0.58816   5000 
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.09722     0.02047      0.06398     0.14410   3191 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.56607     0.01268      0.54237     0.59158   5000 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ## Structured MVN
> 
> (mymodel <- runMLwiN(normexam ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcOptions = list(smvn = 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
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 9362.825198
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 9360.594862
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 9363.073927
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 9361.421867
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 9362.391617
Adapting for 600 iterations (Maximum 5000)
-2 * Loglike = 9358.085197
Adapting for 700 iterations (Maximum 5000)
-2 * Loglike = 9363.357474
Adapting for 800 iterations (Maximum 5000)
-2 * Loglike = 9360.990998
Adapting for 900 iterations (Maximum 5000)
-2 * Loglike = 9360.296437
Adapting for 1000 iterations (Maximum 5000)
-2 * Loglike = 9360.691919
Adapting finished and took 1100 iterations
Adapting took 1100 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 9362.146224
Burning in for 100 iterations out of 500
-2 * Loglike = 9364.698399
Burning in for 150 iterations out of 500
-2 * Loglike = 9362.534714
Burning in for 200 iterations out of 500
-2 * Loglike = 9360.489043
Burning in for 250 iterations out of 500
-2 * Loglike = 9364.454890
Burning in for 300 iterations out of 500
-2 * Loglike = 9359.088457
Burning in for 350 iterations out of 500
-2 * Loglike = 9361.805831
Burning in for 400 iterations out of 500
-2 * Loglike = 9359.550780
Burning in for 450 iterations out of 500
-2 * Loglike = 9360.612186
Burning in for 500 iterations out of 500
-2 * Loglike = 9365.283089

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 = 9365.452212
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9360.937117
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9364.026290
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9363.752298
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9361.746827
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9361.582523
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9363.956382
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9365.064894
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9361.450031
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9360.778354
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9369.173199
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9361.851720
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9360.442694
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9364.875244
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9360.947518
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9366.934844
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9362.788994
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9359.700733
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9357.490608
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9360.521812
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9364.754975
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9368.279920
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9359.206277
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9358.277126
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9360.024906
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9364.341239
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9364.475794
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9363.318276
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9359.682949
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9359.390200
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9362.989099
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9367.513539
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9362.708470
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9362.276539
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9367.856516
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9361.213239
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9362.600631
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9364.007244
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9372.857313
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9358.906690
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9360.824365
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9364.229995
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9360.619689
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9362.585844
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9358.967858
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9357.537609
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9364.019432
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9362.159584
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9366.148311
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9359.647616
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9362.092035
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9357.978092
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9361.152792
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9359.286245
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9362.073225
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9362.110430
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9359.755414
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9359.492089
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9358.678900
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9361.569683
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9364.792387
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9365.415614
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9357.986652
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9363.320770
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9357.529539
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9357.625892
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9362.425055
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9362.224139
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9359.137827
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9367.467738
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9359.788856
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9361.258638
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9363.233216
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9358.832153
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9360.702278
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9366.031663
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9359.558954
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9359.620431
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9362.757603
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9363.212602
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9364.196227
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9362.466529
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9362.645632
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9359.730576
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9357.399877
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9359.575070
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9359.681231
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9358.145070
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9363.513302
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9358.807638
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9359.566248
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9358.357916
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9360.937777
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9358.934388
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9358.493370
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9358.336089
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9361.954951
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9361.371640
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9364.677774
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9360.380043
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9361.5
S.D.    =          2.9802
S.E.M.  =        0.042147

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.36s 
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
9361.457   9357.472   3.985      9365.443   
--------------------------------------------------------------------------------------------------- 
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.04390   0.08     0.9362         -0.08178     0.09048   1093 
standlrt    0.56262     0.01252  44.92          0  ***     0.53745     0.58656   1120 
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.10146     0.02121      0.06727     0.15104   1003 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.56610     0.01250      0.54177     0.59062   1062 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 22.4 Assessing the need for the level 2 variance . . . . . . . . . . . 350
> 
> sixway(mymodel@chains[, "RP2_var_Intercept", drop = FALSE], "sigma2u0")
> 
> set.seed(1)
> tutorial$temp <- rnorm(4059)
> 
> ## Define the model
> 
> ## IGLS
> 
> (mymodel <- runMLwiN(temp ~ 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 achieved
TOLE 2
MAXI 20
NEXT
iteration 2

Convergence achieved
ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.03s 
Number of obs:  4059 (from total 4059)        The model converged after 3 iterations.
Log likelihood:      -5904.1 
Deviance statistic:  11808.3 
--------------------------------------------------------------------------------------------------- 
The model formula:
temp ~ 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.00085     0.01627    0.05     0.9584           -0.03103     0.03273 
standlrt    -0.00606     0.01638   -0.37     0.7114           -0.03816     0.02604 
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.00000     0.00000 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err. 
var_Intercept   1.07389     0.02384 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ## Gibbs
> 
> (mymodel <- runMLwiN(temp ~ 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 achieved
TOLE 2
MAXI 20
BATC 1

Batch mode is ON
NEXT
iteration 2

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
Burning in for 50 iterations out of 500
-2 * Loglike = 11808.142902
Burning in for 100 iterations out of 500
-2 * Loglike = 11818.701337
Burning in for 150 iterations out of 500
-2 * Loglike = 11809.750288
Burning in for 200 iterations out of 500
-2 * Loglike = 11803.809296
Burning in for 250 iterations out of 500
-2 * Loglike = 11808.102753
Burning in for 300 iterations out of 500
-2 * Loglike = 11809.850726
Burning in for 350 iterations out of 500
-2 * Loglike = 11811.591596
Burning in for 400 iterations out of 500
-2 * Loglike = 11806.665959
Burning in for 450 iterations out of 500
-2 * Loglike = 11806.971022
Burning in for 500 iterations out of 500
-2 * Loglike = 11816.329114

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 = 11805.528416
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 11816.054648
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 11809.518147
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 11805.670977
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 11803.508696
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 11808.110473
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 11810.791396
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 11801.107681
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 11809.786714
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 11804.398000
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 11804.711673
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 11803.687937
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 11816.068385
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 11808.573758
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 11812.092009
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 11819.578714
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 11814.102862
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 11810.174823
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 11811.488875
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 11802.564815
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 11812.683232
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 11809.293671
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 11802.724650
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 11809.351976
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 11803.580429
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 11808.815671
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 11813.152128
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 11810.060530
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 11809.315957
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 11800.512852
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 11805.290150
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 11803.306609
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 11804.534866
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 11808.279503
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 11813.408176
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 11803.773026
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 11806.037478
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 11808.975355
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 11822.650844
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 11803.943054
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 11811.556640
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 11813.369826
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 11810.158132
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 11814.643543
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 11805.837990
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 11826.820069
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 11808.403943
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 11812.380541
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 11812.779757
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 11811.994570
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 11813.078072
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 11817.477923
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 11810.696625
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 11809.642024
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 11807.014420
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 11806.392512
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 11824.199638
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 11807.055519
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 11802.828644
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 11813.683393
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 11809.722197
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 11808.093174
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 11811.935271
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 11806.210866
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 11809.267239
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 11807.885339
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 11814.307281
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 11801.563219
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 11810.477934
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 11809.253014
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 11812.326950
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 11810.904155
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 11811.408584
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 11809.151897
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 11804.079441
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 11803.878928
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 11805.467359
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 11808.979023
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 11808.821776
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 11815.916841
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 11806.834459
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 11807.454242
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 11814.916535
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 11806.855922
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 11797.062941
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 11806.547215
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 11811.997772
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 11809.444423
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 11795.749420
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 11806.470851
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 11810.190948
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 11815.330083
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 11803.388315
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 11806.867092
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 11806.992070
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 11805.282597
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 11808.955478
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 11812.066714
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 11815.255229
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 11807.715520
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          11809.
S.D.    =          4.4412
S.E.M.  =        0.062808

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.36s 
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
11809.329  11800.997  8.332      11817.661  
--------------------------------------------------------------------------------------------------- 
The model formula:
temp ~ 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.00072     0.01729   0.04     0.9668         -0.03358     0.03422   3971 
standlrt    -0.00595     0.01645  -0.36     0.7178         -0.03764     0.02611   4787 
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.00163     0.00120      0.00031     0.00463   248 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   1.07478     0.02398      1.02899     1.12291   5000 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ## Structured MVN
> 
> (mymodel <- runMLwiN(temp ~ 1 + standlrt + (1 | school) + (1 | student), estoptions = list(EstM = 1, mcmcOptions = list(smvn = 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 achieved
TOLE 2
MAXI 20
BATC 1

Batch mode is ON
NEXT
iteration 2

Convergence achieved
ECHO 0

Echoing is ON
MCMC   0 500 1 5.8 50 10 G30[1] G30[2]  1 1 1 1 1 1
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 11891.495848
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 11808.073415
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 11807.555947
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 11807.137368
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 11820.094150
Adapting for 600 iterations (Maximum 5000)
-2 * Loglike = 11808.183755
Adapting for 700 iterations (Maximum 5000)
-2 * Loglike = 11807.459363
Adapting for 800 iterations (Maximum 5000)
-2 * Loglike = 11810.792895
Adapting for 900 iterations (Maximum 5000)
-2 * Loglike = 11819.698293
Adapting for 1000 iterations (Maximum 5000)
-2 * Loglike = 11809.330693
Adapting for 1100 iterations (Maximum 5000)
-2 * Loglike = 11810.206487
Adapting for 1200 iterations (Maximum 5000)
-2 * Loglike = 11807.492050
Adapting finished and took 1300 iterations
Adapting took 1300 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 11806.991986
Burning in for 100 iterations out of 500
-2 * Loglike = 11812.747221
Burning in for 150 iterations out of 500
-2 * Loglike = 11811.361448
Burning in for 200 iterations out of 500
-2 * Loglike = 11812.047488
Burning in for 250 iterations out of 500
-2 * Loglike = 11812.814481
Burning in for 300 iterations out of 500
-2 * Loglike = 11808.151226
Burning in for 350 iterations out of 500
-2 * Loglike = 11809.433431
Burning in for 400 iterations out of 500
-2 * Loglike = 11809.185816
Burning in for 450 iterations out of 500
-2 * Loglike = 11811.378355
Burning in for 500 iterations out of 500
-2 * Loglike = 11806.372698

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 = 11808.690641
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 11808.201529
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 11815.633524
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 11816.225637
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 11811.113585
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 11807.857988
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 11811.580508
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 11811.618279
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 11807.809019
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 11806.641639
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 11807.221276
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 11807.864041
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 11807.028415
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 11813.632625
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 11807.110016
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 11811.761273
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 11807.177802
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 11815.398973
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 11810.746410
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 11809.227083
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 11810.596040
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 11810.739626
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 11814.029041
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 11809.346063
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 11809.694301
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 11808.315813
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 11812.904789
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 11814.085218
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 11809.169998
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 11813.293497
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 11810.271966
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 11810.189730
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 11810.318328
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 11807.888556
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 11811.307024
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 11807.525222
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 11808.689747
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 11808.468179
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 11808.900187
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 11809.069073
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 11810.618065
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 11809.031899
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 11808.947728
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 11808.088541
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 11810.415588
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 11807.809598
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 11808.328051
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 11806.763987
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 11810.127480
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 11809.727930
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 11810.399947
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 11807.469276
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 11807.330412
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 11807.925481
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 11810.362979
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 11807.945746
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 11808.818566
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 11808.468987
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 11807.797501
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 11810.008999
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 11810.019026
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 11808.489968
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 11809.509858
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 11814.215594
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 11806.934486
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 11806.975037
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 11810.328335
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 11809.411828
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 11813.057990
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 11814.733279
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 11810.276050
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 11814.346464
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 11812.028001
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 11814.518320
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 11807.956478
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 11810.741747
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 11807.810339
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 11808.860877
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 11807.550677
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 11811.193736
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 11808.211707
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 11809.878520
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 11813.323606
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 11809.815090
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 11810.822129
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 11806.790801
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 11809.239931
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 11807.953527
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 11813.527933
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 11806.866374
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 11808.727554
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 11809.457511
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 11809.728917
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 11811.604603
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 11815.448445
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 11812.243321
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 11817.200234
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 11810.575611
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 11809.024420
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 11808.433868
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          11810.
S.D.    =          2.8177
S.E.M.  =        0.039848

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.05)  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.32s 
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
11810.121  11806.585  3.537      11813.658  
--------------------------------------------------------------------------------------------------- 
The model formula:
temp ~ 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.00025     0.01453   0.02     0.9861         -0.02771     0.02897   1082 
standlrt    -0.00624     0.01620  -0.39     0.7001         -0.03750     0.02568   1036 
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.00217     0.00216     -0.00517     0.00306   739 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   1.07763     0.02292      1.03313     1.12278   1157 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> summary(mymodel@chains[, "RP2_var_Intercept"])

Iterations = 1:5000
Thinning interval = 1 
Number of chains = 1 
Sample size per chain = 5000 

1. Empirical mean and standard deviation for each variable,
   plus standard error of the mean:

          Mean             SD       Naive SE Time-series SE 
    -2.175e-03      2.159e-03      3.053e-05      7.945e-05 

2. Quantiles for each variable:

     2.5%       25%       50%       75%     97.5% 
-0.005174 -0.003776 -0.002543 -0.000952  0.003058 

> sixway(mymodel@chains[, "RP2_var_Intercept", drop = FALSE], "sigma2u0")
> 
> # Chapter learning outcomes . . . . . . . . . . . . . . . . . . . . . . .355
> 
> 
> 
> 
> 
> ############################################################################
> 
> proc.time()
   user  system elapsed 
   4.67    0.65   13.28