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
> #
> # 14  Adjusting for Measurement Errors in Predictor Variables . . . . . .199
> #
> #     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/
> ############################################################################
> 
> 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)
> 
> ## Read tutorial data
> data(tutorial, package = "R2MLwiN")
> 
> # 14.1 Effects of measurement error on predictors . . . . . . . . . . . .200
> set.seed(1)
> error <- rnorm(length(tutorial$standlrt), 0, sqrt(0.2))
> obslrt <- tutorial$standlrt + error
> tutorial <- cbind(tutorial, error, obslrt)
> 
> (mymodel <- runMLwiN(normexam ~ 1 + standlrt + (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.09)  multilevel model (Normal) 
Estimation algorithm:  IGLS        Elapsed time : 0.06s 
Number of obs:  4059 (from total 4059)        The model converged after 3 iterations.
Log likelihood:      -4880.3 
Deviance statistic:  9760.5 
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 | student)
Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.       z   Pr(>|z|)         [95% Conf.   Interval] 
Intercept   -0.00119     0.01264   -0.09     0.9249           -0.02596     0.02358 
standlrt     0.59506     0.01273   46.76          0   ***      0.57011     0.62000 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err. 
var_Intercept   0.64842     0.01439 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> (mymodel <- runMLwiN(normexam ~ 1 + error + (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.09)  multilevel model (Normal) 
Estimation algorithm:  IGLS        Elapsed time : 0.07s 
Number of obs:  4059 (from total 4059)        The model converged after 3 iterations.
Log likelihood:      -5754.2 
Deviance statistic:  11508.4 
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + error + (1 | student)
Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.       z   Pr(>|z|)         [95% Conf.   Interval] 
Intercept   -0.00010     0.01568   -0.01     0.9948           -0.03083     0.03062 
error       -0.03305     0.03382   -0.98     0.3285           -0.09935     0.03324 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err. 
var_Intercept   0.99741     0.02214 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> (mymodel <- runMLwiN(normexam ~ 1 + obslrt + (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.09)  multilevel model (Normal) 
Estimation algorithm:  IGLS        Elapsed time : 0.06s 
Number of obs:  4059 (from total 4059)        The model converged after 3 iterations.
Log likelihood:      -5082.9 
Deviance statistic:  10165.9 
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + obslrt + (1 | student)
Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.       z   Pr(>|z|)         [95% Conf.   Interval] 
Intercept   -0.00117     0.01329   -0.09     0.9296           -0.02721     0.02487 
obslrt       0.48489     0.01215   39.91          0   ***      0.46107     0.50870 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err. 
var_Intercept   0.71652     0.01590 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> (mymodel <- runMLwiN(normexam ~ 1 + obslrt + (1 | student), estoptions = list(EstM = 1, merr = c(N = 1, "obslrt",
+   0.2)), 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   1 1 1 1 1 1
BURNING IN...
Burning in for 50 iterations out of 500
-2 * Loglike = 10030.845657
Burning in for 100 iterations out of 500
-2 * Loglike = 9851.851496
Burning in for 150 iterations out of 500
-2 * Loglike = 9843.264679
Burning in for 200 iterations out of 500
-2 * Loglike = 9858.702314
Burning in for 250 iterations out of 500
-2 * Loglike = 9888.231519
Burning in for 300 iterations out of 500
-2 * Loglike = 9813.536891
Burning in for 350 iterations out of 500
-2 * Loglike = 9843.921642
Burning in for 400 iterations out of 500
-2 * Loglike = 9800.714375
Burning in for 450 iterations out of 500
-2 * Loglike = 9840.312040
Burning in for 500 iterations out of 500
-2 * Loglike = 9788.990109

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 = 9847.304628
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9819.677891
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9848.074774
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9791.465112
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9921.508411
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9778.157787
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9842.020229
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9810.932119
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9906.691621
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9848.405205
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9808.855434
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9839.063777
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9869.949887
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9846.491754
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9854.445409
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9775.309901
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9846.537354
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9906.111436
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9831.622497
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9849.783603
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9925.350094
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9879.287368
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9902.331039
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9857.586986
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9869.943413
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9851.253545
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9849.119034
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9823.024482
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9890.138041
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9824.142453
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9868.086119
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9792.192429
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9867.982866
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9739.912920
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9860.880063
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9861.782844
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9854.273092
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9807.031052
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9767.102476
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9859.825103
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9806.644858
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9828.001232
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9810.584953
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9815.807291
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9795.376865
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9794.356658
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9850.553847
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9766.698580
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9881.067990
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9863.191879
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9847.504459
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9812.147786
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9819.640761
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9807.046431
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9803.070480
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9793.072528
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9823.526435
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9866.810586
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9836.905477
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9870.657468
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9786.882076
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9811.108759
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9875.196027
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9837.623099
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9870.371265
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9833.170599
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9852.286089
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9869.407387
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9885.701904
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9774.055522
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9853.793481
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9805.597142
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9902.754888
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9880.965631
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9802.835661
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9790.205208
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9872.982077
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9872.558436
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9796.176164
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9860.425046
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9815.898081
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9881.146968
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9845.893271
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9739.407179
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9884.568957
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9864.781319
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9837.937826
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9790.670662
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9771.457137
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9804.427734
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9873.549246
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9826.286372
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9878.335481
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9803.820622
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9860.198670
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9890.513414
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9875.162411
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9788.761250
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9814.122794
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9849.955140
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9834.0
S.D.    =          41.190
S.E.M.  =         0.58251

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  multilevel model (Normal) 
Estimation algorithm:  MCMC      Elapsed time : 3.81s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Deviance statistic:  10165.9 
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + obslrt + (1 | student)
Level 1: student      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.      z   Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   -0.00154     0.01341  -0.11     0.9087         -0.02766     0.02476   4782 
obslrt       0.58279     0.01447  40.29          0  ***     0.55500     0.61126   2850 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.66053     0.01623      0.62956     0.69350   3486 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 14.2 Measurement error modelling in multilevel models . . . . . . . . .205
> 
> (mymodel <- runMLwiN(normexam ~ 1 + standlrt + (1 + standlrt | 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 = 9136.615416
Burning in for 100 iterations out of 500
-2 * Loglike = 9129.472590
Burning in for 150 iterations out of 500
-2 * Loglike = 9137.752388
Burning in for 200 iterations out of 500
-2 * Loglike = 9149.641457
Burning in for 250 iterations out of 500
-2 * Loglike = 9113.737768
Burning in for 300 iterations out of 500
-2 * Loglike = 9127.854908
Burning in for 350 iterations out of 500
-2 * Loglike = 9099.670195
Burning in for 400 iterations out of 500
-2 * Loglike = 9107.133017
Burning in for 450 iterations out of 500
-2 * Loglike = 9108.686421
Burning in for 500 iterations out of 500
-2 * Loglike = 9142.796685

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 = 9114.789855
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9128.936360
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9120.626663
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9140.221260
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9133.366931
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9112.633267
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9110.267312
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9134.271816
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9146.903828
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9118.578020
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9132.630258
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9105.115711
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9126.654041
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9138.820724
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9119.239813
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9145.879653
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9132.483550
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9094.657080
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9118.502105
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9094.902922
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9135.645832
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9128.926707
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9104.509733
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9108.989181
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9111.759420
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9103.538840
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9120.142978
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9117.533977
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9114.717113
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9118.929584
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9146.845851
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9126.654978
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9127.854549
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9114.945398
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9156.130168
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9134.732656
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9142.867734
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9137.471434
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9127.702153
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9133.939751
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9107.640871
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9137.519015
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9147.472132
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9104.767457
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9118.396071
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9143.494727
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9130.468641
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9137.527730
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9106.747459
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9154.348069
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9142.340919
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9100.831693
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9131.653716
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9151.916582
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9165.855130
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9125.508254
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9146.519009
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9122.524343
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9128.808395
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9125.880481
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9130.006670
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9116.446060
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9105.444348
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9141.749060
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9122.819747
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9140.592905
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9119.521430
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9137.983211
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9128.855413
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9130.638160
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9119.619021
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9124.980114
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9135.552922
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9127.477138
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9114.227379
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9110.953194
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9118.839082
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9118.763419
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9117.590737
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9106.053640
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9112.274103
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9115.311806
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9122.276612
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9120.226668
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9140.403968
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9153.535294
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9117.253885
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9116.415887
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9136.007855
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9119.650550
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9108.406619
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9104.027872
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9104.517145
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9102.506754
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9125.791725
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9121.949129
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9123.201133
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9116.574595
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9120.080554
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9120.628553
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9123.2
S.D.    =          16.792
S.E.M.  =         0.23747

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 : 2.08s 
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
9123.206   9031.714   91.492     9214.698   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + standlrt + (1 + standlrt | 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.01257     0.04089  -0.31       0.7585         -0.09055     0.06865   209 
standlrt     0.55645     0.02084  26.70   4.362e-157  ***     0.51501     0.59701   596 
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.09711     0.02001      0.06465     0.14261   2549 
cov_Intercept_standlrt   0.01933     0.00739      0.00660     0.03537   1740 
var_standlrt             0.01526     0.00485      0.00772     0.02662    846 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.55440     0.01242      0.53052     0.57904   4750 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> (mymodel <- runMLwiN(normexam ~ 1 + obslrt + (1 + obslrt | 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 = 9557.668282
Burning in for 100 iterations out of 500
-2 * Loglike = 9550.917606
Burning in for 150 iterations out of 500
-2 * Loglike = 9549.118836
Burning in for 200 iterations out of 500
-2 * Loglike = 9565.853907
Burning in for 250 iterations out of 500
-2 * Loglike = 9535.755436
Burning in for 300 iterations out of 500
-2 * Loglike = 9551.604339
Burning in for 350 iterations out of 500
-2 * Loglike = 9521.217086
Burning in for 400 iterations out of 500
-2 * Loglike = 9525.186785
Burning in for 450 iterations out of 500
-2 * Loglike = 9527.220279
Burning in for 500 iterations out of 500
-2 * Loglike = 9562.996780

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 = 9533.942692
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9546.830416
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9538.011865
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9562.720763
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9549.270526
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9525.594138
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9532.863206
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9550.859701
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9559.755437
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9539.984608
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9548.853954
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9528.929272
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9547.606762
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9565.177846
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9535.605809
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9559.737804
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9554.399230
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9510.102497
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9542.541735
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9514.007889
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9558.415358
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9551.683094
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9532.865730
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9530.380909
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9534.912076
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9519.856443
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9539.607552
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9541.111941
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9524.165669
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9539.106284
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9563.877231
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9551.022654
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9541.087252
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9531.774728
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9573.923755
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9557.885911
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9562.648458
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9561.604055
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9545.540060
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9554.142479
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9525.332602
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9557.005210
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9566.388943
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9526.952193
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9535.427972
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9559.943565
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9555.797602
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9558.778354
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9528.578565
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9569.218993
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9561.648845
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9521.225257
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9555.237439
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9570.927445
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9591.336281
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9547.903420
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9567.588973
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9534.928934
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9549.500399
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9542.623386
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9549.179479
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9545.705284
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9524.902854
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9564.493064
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9543.741719
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9560.419950
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9538.132705
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9554.308780
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9551.226625
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9548.161279
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9543.936913
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9545.821637
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9552.829729
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9549.571849
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9530.718690
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9534.998514
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9540.139851
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9544.511048
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9538.782197
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9528.921241
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9525.715522
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9537.007117
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9535.882353
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9556.707698
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9563.200049
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9569.650403
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9539.258262
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9529.417533
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9557.965434
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9540.101298
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9530.353551
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9526.039907
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9525.832007
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9518.499086
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9551.095047
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9541.579475
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9545.171747
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9541.329254
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9539.143612
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9536.765172
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9543.3
S.D.    =          16.788
S.E.M.  =         0.23742

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 : 2.05s 
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
9543.340   9454.185   89.155     9632.495   
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + obslrt + (1 + obslrt | 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.01368     0.04347  -0.31       0.7529         -0.09629     0.07262   205 
obslrt       0.44666     0.01849  24.16   6.015e-129  ***     0.40942     0.48240   677 
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.11021     0.02258      0.07358     0.16132   2592 
cov_Intercept_obslrt   0.01666     0.00691      0.00432     0.03161   1549 
var_obslrt             0.01111     0.00385      0.00525     0.02022    697 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.61486     0.01376      0.58831     0.64217   4739 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> (mymodel <- runMLwiN(normexam ~ 1 + obslrt + (1 + obslrt | school) + (1 | student), estoptions = list(EstM = 1, merr = c(N = 1,
+   "obslrt", 0.2)), 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 = 9320.351125
Burning in for 100 iterations out of 500
-2 * Loglike = 9260.818285
Burning in for 150 iterations out of 500
-2 * Loglike = 9141.510893
Burning in for 200 iterations out of 500
-2 * Loglike = 9145.593816
Burning in for 250 iterations out of 500
-2 * Loglike = 9194.006232
Burning in for 300 iterations out of 500
-2 * Loglike = 9137.387311
Burning in for 350 iterations out of 500
-2 * Loglike = 9066.721982
Burning in for 400 iterations out of 500
-2 * Loglike = 9132.959682
Burning in for 450 iterations out of 500
-2 * Loglike = 9115.614642
Burning in for 500 iterations out of 500
-2 * Loglike = 9165.227555

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 = 9169.822790
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 9125.315605
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 9180.133121
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 9139.562343
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 9108.393027
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 9174.636434
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 9082.428470
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 9166.382884
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 9122.760679
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 9146.073877
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 9171.913143
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 9162.283468
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 9182.275881
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 9027.633098
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 9187.103140
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 9206.396795
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 9106.301202
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 9122.710869
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 9102.116353
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 9068.202328
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 9165.731705
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 9235.685589
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 9160.124956
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 9061.767625
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 9096.306304
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 9079.234017
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 9165.745558
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 9135.493488
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 9149.678508
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 9162.716973
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 9219.169513
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 9159.935102
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 9113.345853
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 9199.043005
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 9063.196993
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 9188.132645
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 9121.092220
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 9052.393725
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 9198.934878
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 9183.629757
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 9104.895322
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 9164.201968
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 9182.605933
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 9199.759322
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 9173.826678
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 9090.962133
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 9197.778300
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 9157.893644
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 9139.076110
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 9148.674366
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 9196.099735
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 9064.372964
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 9108.716338
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 9047.232875
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 9159.336633
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 9101.936371
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 9127.966133
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 9114.384562
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 9179.478007
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 9140.143147
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 9224.961460
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 9179.675745
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 9142.885800
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 9110.116998
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 9052.127550
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 9008.559547
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 9206.072783
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 9129.905029
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 9030.933147
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 9219.652280
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 9188.472382
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 9130.086829
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 9005.777033
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 9091.525305
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 9181.997443
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 9134.847759
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 9111.981805
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 9113.637146
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 9185.335924
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 9128.918247
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 9082.486609
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 9135.323714
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 9032.296682
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 9201.895934
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 9107.381195
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 9119.637739
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 9086.643888
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 9129.511940
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 9190.244905
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 9237.996677
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 9148.937698
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 9149.452312
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 9124.560876
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 9097.198341
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 9140.730696
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 9223.723116
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 9142.459031
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 9107.239612
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 9160.518346
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 9174.181651
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          9143.0
S.D.    =          49.641
S.E.M.  =         0.70203

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 : 6.86s 
Number of obs:  4059 (from total 4059)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Deviance statistic:  9732.6 
--------------------------------------------------------------------------------------------------- 
The model formula:
normexam ~ 1 + obslrt + (1 + obslrt | 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.01346     0.04326  -0.31       0.7557         -0.09414     0.07359   194 
obslrt       0.53729     0.02265  23.72   2.156e-124  ***     0.49286     0.58115   659 
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.10976     0.02265      0.07298     0.15995   2359 
cov_Intercept_obslrt   0.02210     0.00903      0.00619     0.04150   1568 
var_obslrt             0.02114     0.00666      0.01061     0.03644    749 
--------------------------------------------------------------------------------------------------- 
The random part estimates at the student level: 
                  Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_Intercept   0.55723     0.01437      0.52972     0.58654   2703 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 14.3 Measurement errors in binomial models . . . . . . . . . . . . . . 208
> 
> 
> ## Read bang1 data
> data(bang1, package = "R2MLwiN")
> 
> set.seed(1)
> bang1$obsage <- bang1$age + rnorm(length(bang1$age), 0, 5)
> 
> (mymodel <- runMLwiN(logit(use) ~ 1 + age, D = "Binomial", estoptions = list(EstM = 1), data = bang1))
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   2 2 2 1 1 2
RUNNING ADAPTIVE PROCEDURE AND BURNING IN...
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 2589.951603
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 2590.870573
Adapting finished and took 300 iterations
Adapting took 300 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 2591.646899
Burning in for 100 iterations out of 500
-2 * Loglike = 2589.408544
Burning in for 150 iterations out of 500
-2 * Loglike = 2591.703242
Burning in for 200 iterations out of 500
-2 * Loglike = 2589.578707
Burning in for 250 iterations out of 500
-2 * Loglike = 2591.143776
Burning in for 300 iterations out of 500
-2 * Loglike = 2589.497722
Burning in for 350 iterations out of 500
-2 * Loglike = 2589.866468
Burning in for 400 iterations out of 500
-2 * Loglike = 2591.366495
Burning in for 450 iterations out of 500
-2 * Loglike = 2601.081625
Burning in for 500 iterations out of 500
-2 * Loglike = 2590.826384

ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 2
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 2593.768209
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 2594.974846
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 2592.681292
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 2595.029755
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 2589.956166
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 2589.700843
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 2590.249050
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 2595.764631
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 2594.215259
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 2593.583255
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 2589.738290
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 2599.665863
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 2590.851598
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 2594.064439
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 2590.602086
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 2592.016924
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 2590.673363
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 2592.401756
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 2592.806336
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 2589.859640
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 2589.674262
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 2591.332550
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 2590.929091
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 2590.027123
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 2590.651481
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 2593.305020
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 2590.219008
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 2592.793939
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 2590.448211
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 2591.191911
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 2591.086373
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 2594.539469
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 2591.927492
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 2592.508995
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 2589.905070
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 2591.755927
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 2591.259256
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 2590.056615
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 2591.229959
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 2590.161462
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 2589.898548
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 2590.782298
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 2591.749425
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 2589.438829
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 2589.692805
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 2589.387252
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 2591.241413
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 2591.300298
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 2594.469366
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 2594.291484
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 2591.115828
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 2591.134031
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 2590.254941
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 2593.757296
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 2589.491965
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 2590.248874
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 2592.171248
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 2589.329180
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 2590.887602
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 2589.460348
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 2590.804734
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 2589.476827
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 2589.524551
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 2590.356007
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 2590.321576
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 2590.285680
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 2591.197588
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 2595.992028
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 2589.900265
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 2590.613041
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 2592.566051
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 2591.248283
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 2590.675361
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 2589.528497
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 2593.472607
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 2589.815054
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 2591.423243
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 2589.429422
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 2590.500183
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 2592.090979
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 2592.064172
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 2589.588351
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 2591.532194
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 2600.825764
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 2592.877217
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 2590.619297
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 2592.853876
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 2593.211902
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 2589.607546
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 2591.407865
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 2590.268036
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 2590.484512
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 2589.891584
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 2589.687159
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 2591.627120
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 2591.106576
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 2591.048963
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 2589.781838
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 2591.584549
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 2591.501194
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          2591.3
S.D.    =          1.9941
S.E.M.  =        0.028200

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  multilevel model (Binomial) 
Estimation algorithm:  MCMC      Elapsed time : 2.58s 
Number of obs:  1934 (from total 1934)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
2591.294   2589.293   2.001      2593.294   
--------------------------------------------------------------------------------------------------- 
The model formula:
logit(use) ~ 1 + age
Level 1: l1id      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.      z    Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   -0.43858     0.04731  -9.27   1.863e-20  ***    -0.53074    -0.34677   1100 
age          0.00681     0.00509   1.34      0.1803         -0.00333     0.01649   1145 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the l1id level: 
                Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_bcons_1   1.00000       1e-05      1.00000     1.00000   5000 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> (mymodel <- runMLwiN(logit(use) ~ 1 + obsage, D = "Binomial", estoptions = list(EstM = 1), data = bang1))
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   2 2 2 1 1 2
RUNNING ADAPTIVE PROCEDURE AND BURNING IN...
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 2595.199073
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 2589.902225
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 2589.864198
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 2589.867866
Adapting finished and took 500 iterations
Adapting took 500 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 2591.696711
Burning in for 100 iterations out of 500
-2 * Loglike = 2589.680114
Burning in for 150 iterations out of 500
-2 * Loglike = 2592.340952
Burning in for 200 iterations out of 500
-2 * Loglike = 2589.744447
Burning in for 250 iterations out of 500
-2 * Loglike = 2592.048609
Burning in for 300 iterations out of 500
-2 * Loglike = 2590.549333
Burning in for 350 iterations out of 500
-2 * Loglike = 2589.886888
Burning in for 400 iterations out of 500
-2 * Loglike = 2589.990956
Burning in for 450 iterations out of 500
-2 * Loglike = 2590.173068
Burning in for 500 iterations out of 500
-2 * Loglike = 2592.398467

ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 2
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 2589.817573
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 2590.623709
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 2594.023814
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 2593.019489
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 2592.433573
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 2590.196465
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 2592.624468
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 2592.679076
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 2592.581293
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 2591.826219
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 2591.729285
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 2592.411348
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 2590.485122
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 2589.950383
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 2590.657924
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 2590.648427
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 2592.723986
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 2590.801623
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 2589.923736
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 2597.092208
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 2592.668731
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 2589.903505
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 2591.721020
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 2591.208793
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 2590.656182
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 2592.717453
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 2591.016605
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 2591.492345
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 2590.918941
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 2590.742885
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 2592.586752
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 2589.857299
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 2590.341719
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 2590.702930
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 2589.933571
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 2589.929880
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 2589.597299
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 2590.533599
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 2590.202207
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 2592.445180
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 2589.996578
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 2590.419777
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 2591.436971
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 2590.467230
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 2590.675127
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 2596.279303
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 2592.571497
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 2592.423897
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 2591.213799
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 2590.297146
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 2589.591340
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 2591.205876
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 2591.674572
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 2592.545471
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 2590.538674
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 2594.789342
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 2590.733650
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 2593.186276
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 2591.251757
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 2591.329570
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 2590.566061
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 2597.085175
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 2595.773919
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 2589.998117
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 2591.815930
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 2591.496575
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 2590.334875
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 2595.209172
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 2591.518678
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 2591.211781
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 2590.951866
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 2591.897413
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 2590.477940
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 2589.951612
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 2590.679531
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 2590.851501
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 2590.412308
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 2589.985272
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 2596.651283
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 2590.312295
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 2592.998750
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 2591.372354
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 2589.705697
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 2590.894771
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 2592.051769
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 2589.574243
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 2591.654868
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 2591.490101
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 2591.570147
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 2592.357774
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 2591.596151
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 2589.629862
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 2590.503448
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 2590.528930
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 2592.032511
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 2590.244396
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 2589.769054
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 2590.714765
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 2589.781194
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 2591.713482
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          2591.5
S.D.    =          1.9662
S.E.M.  =        0.027806

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  multilevel model (Binomial) 
Estimation algorithm:  MCMC      Elapsed time : 2.69s 
Number of obs:  1934 (from total 1934)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Bayesian Deviance Information Criterion (DIC)
Dbar      D(thetabar)    pD      DIC
2591.522   2589.572   1.950      2593.472   
--------------------------------------------------------------------------------------------------- 
The model formula:
logit(use) ~ 1 + obsage
Level 1: l1id      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.      z    Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   -0.43710     0.04510  -9.69   3.281e-22  ***    -0.52199    -0.34799   1111 
obsage       0.00490     0.00450   1.09      0.2762         -0.00414     0.01355   1175 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the l1id level: 
                Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_bcons_1   1.00000       1e-05      1.00000     1.00000   5000 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> ## Adjust for the measurement errors
> (mymodel <- runMLwiN(logit(use) ~ 1 + obsage, D = "Binomial", estoptions = list(EstM = 1, merr = c(N = 1,
+   "obsage", 25)), data = bang1))
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   2 2 2 1 1 2
RUNNING ADAPTIVE PROCEDURE AND BURNING IN...
Adapting for 100 iterations (Maximum 5000)
-2 * Loglike = 2587.296697
Adapting for 200 iterations (Maximum 5000)
-2 * Loglike = 2590.917904
Adapting for 300 iterations (Maximum 5000)
-2 * Loglike = 2589.351124
Adapting for 400 iterations (Maximum 5000)
-2 * Loglike = 2591.896449
Adapting for 500 iterations (Maximum 5000)
-2 * Loglike = 2591.706105
Adapting for 600 iterations (Maximum 5000)
-2 * Loglike = 2591.515255
Adapting for 700 iterations (Maximum 5000)
-2 * Loglike = 2591.949838
Adapting for 800 iterations (Maximum 5000)
-2 * Loglike = 2594.326103
Adapting for 900 iterations (Maximum 5000)
-2 * Loglike = 2588.150069
Adapting for 1000 iterations (Maximum 5000)
-2 * Loglike = 2589.923420
Adapting for 1100 iterations (Maximum 5000)
-2 * Loglike = 2594.209432
Adapting for 1200 iterations (Maximum 5000)
-2 * Loglike = 2589.788355
Adapting for 1300 iterations (Maximum 5000)
-2 * Loglike = 2594.735322
Adapting for 1400 iterations (Maximum 5000)
-2 * Loglike = 2589.281035
Adapting for 1500 iterations (Maximum 5000)
-2 * Loglike = 2591.736834
Adapting for 1600 iterations (Maximum 5000)
-2 * Loglike = 2586.839461
Adapting for 1700 iterations (Maximum 5000)
-2 * Loglike = 2588.688721
Adapting finished and took 1800 iterations
Adapting took 1800 iterations
Burning in for 50 iterations out of 500
-2 * Loglike = 2590.576466
Burning in for 100 iterations out of 500
-2 * Loglike = 2591.915352
Burning in for 150 iterations out of 500
-2 * Loglike = 2592.777747
Burning in for 200 iterations out of 500
-2 * Loglike = 2590.017124
Burning in for 250 iterations out of 500
-2 * Loglike = 2588.691006
Burning in for 300 iterations out of 500
-2 * Loglike = 2590.087763
Burning in for 350 iterations out of 500
-2 * Loglike = 2589.218394
Burning in for 400 iterations out of 500
-2 * Loglike = 2591.784592
Burning in for 450 iterations out of 500
-2 * Loglike = 2589.934634
Burning in for 500 iterations out of 500
-2 * Loglike = 2591.717224

ERAS  c1090 c1091

MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 2
Actual update 50 of 5000, Stored update 50 of 5000
-2 * Loglike = 2593.801457
Actual update 100 of 5000, Stored update 100 of 5000
-2 * Loglike = 2593.732311
Actual update 150 of 5000, Stored update 150 of 5000
-2 * Loglike = 2591.501635
Actual update 200 of 5000, Stored update 200 of 5000
-2 * Loglike = 2589.968686
Actual update 250 of 5000, Stored update 250 of 5000
-2 * Loglike = 2589.932230
Actual update 300 of 5000, Stored update 300 of 5000
-2 * Loglike = 2596.655870
Actual update 350 of 5000, Stored update 350 of 5000
-2 * Loglike = 2590.127903
Actual update 400 of 5000, Stored update 400 of 5000
-2 * Loglike = 2590.067537
Actual update 450 of 5000, Stored update 450 of 5000
-2 * Loglike = 2588.771925
Actual update 500 of 5000, Stored update 500 of 5000
-2 * Loglike = 2591.800016
Actual update 550 of 5000, Stored update 550 of 5000
-2 * Loglike = 2589.154095
Actual update 600 of 5000, Stored update 600 of 5000
-2 * Loglike = 2589.218438
Actual update 650 of 5000, Stored update 650 of 5000
-2 * Loglike = 2587.029599
Actual update 700 of 5000, Stored update 700 of 5000
-2 * Loglike = 2587.352441
Actual update 750 of 5000, Stored update 750 of 5000
-2 * Loglike = 2591.050325
Actual update 800 of 5000, Stored update 800 of 5000
-2 * Loglike = 2590.050677
Actual update 850 of 5000, Stored update 850 of 5000
-2 * Loglike = 2592.606194
Actual update 900 of 5000, Stored update 900 of 5000
-2 * Loglike = 2589.928948
Actual update 950 of 5000, Stored update 950 of 5000
-2 * Loglike = 2591.498326
Actual update 1000 of 5000, Stored update 1000 of 5000
-2 * Loglike = 2595.436023
Actual update 1050 of 5000, Stored update 1050 of 5000
-2 * Loglike = 2589.268145
Actual update 1100 of 5000, Stored update 1100 of 5000
-2 * Loglike = 2594.293472
Actual update 1150 of 5000, Stored update 1150 of 5000
-2 * Loglike = 2593.869591
Actual update 1200 of 5000, Stored update 1200 of 5000
-2 * Loglike = 2589.049693
Actual update 1250 of 5000, Stored update 1250 of 5000
-2 * Loglike = 2592.859466
Actual update 1300 of 5000, Stored update 1300 of 5000
-2 * Loglike = 2591.997995
Actual update 1350 of 5000, Stored update 1350 of 5000
-2 * Loglike = 2592.253875
Actual update 1400 of 5000, Stored update 1400 of 5000
-2 * Loglike = 2588.742808
Actual update 1450 of 5000, Stored update 1450 of 5000
-2 * Loglike = 2591.575817
Actual update 1500 of 5000, Stored update 1500 of 5000
-2 * Loglike = 2592.124739
Actual update 1550 of 5000, Stored update 1550 of 5000
-2 * Loglike = 2588.054433
Actual update 1600 of 5000, Stored update 1600 of 5000
-2 * Loglike = 2591.351245
Actual update 1650 of 5000, Stored update 1650 of 5000
-2 * Loglike = 2591.128259
Actual update 1700 of 5000, Stored update 1700 of 5000
-2 * Loglike = 2591.167754
Actual update 1750 of 5000, Stored update 1750 of 5000
-2 * Loglike = 2592.132583
Actual update 1800 of 5000, Stored update 1800 of 5000
-2 * Loglike = 2590.972674
Actual update 1850 of 5000, Stored update 1850 of 5000
-2 * Loglike = 2592.087644
Actual update 1900 of 5000, Stored update 1900 of 5000
-2 * Loglike = 2591.336572
Actual update 1950 of 5000, Stored update 1950 of 5000
-2 * Loglike = 2592.243831
Actual update 2000 of 5000, Stored update 2000 of 5000
-2 * Loglike = 2590.513777
Actual update 2050 of 5000, Stored update 2050 of 5000
-2 * Loglike = 2592.682514
Actual update 2100 of 5000, Stored update 2100 of 5000
-2 * Loglike = 2590.642406
Actual update 2150 of 5000, Stored update 2150 of 5000
-2 * Loglike = 2591.280187
Actual update 2200 of 5000, Stored update 2200 of 5000
-2 * Loglike = 2585.319839
Actual update 2250 of 5000, Stored update 2250 of 5000
-2 * Loglike = 2592.336145
Actual update 2300 of 5000, Stored update 2300 of 5000
-2 * Loglike = 2591.285905
Actual update 2350 of 5000, Stored update 2350 of 5000
-2 * Loglike = 2593.065948
Actual update 2400 of 5000, Stored update 2400 of 5000
-2 * Loglike = 2591.052767
Actual update 2450 of 5000, Stored update 2450 of 5000
-2 * Loglike = 2592.235381
Actual update 2500 of 5000, Stored update 2500 of 5000
-2 * Loglike = 2590.701800
Actual update 2550 of 5000, Stored update 2550 of 5000
-2 * Loglike = 2590.319799
Actual update 2600 of 5000, Stored update 2600 of 5000
-2 * Loglike = 2592.135881
Actual update 2650 of 5000, Stored update 2650 of 5000
-2 * Loglike = 2592.089608
Actual update 2700 of 5000, Stored update 2700 of 5000
-2 * Loglike = 2587.947954
Actual update 2750 of 5000, Stored update 2750 of 5000
-2 * Loglike = 2589.939323
Actual update 2800 of 5000, Stored update 2800 of 5000
-2 * Loglike = 2589.739719
Actual update 2850 of 5000, Stored update 2850 of 5000
-2 * Loglike = 2590.722767
Actual update 2900 of 5000, Stored update 2900 of 5000
-2 * Loglike = 2590.541238
Actual update 2950 of 5000, Stored update 2950 of 5000
-2 * Loglike = 2592.168696
Actual update 3000 of 5000, Stored update 3000 of 5000
-2 * Loglike = 2591.036347
Actual update 3050 of 5000, Stored update 3050 of 5000
-2 * Loglike = 2590.307136
Actual update 3100 of 5000, Stored update 3100 of 5000
-2 * Loglike = 2591.845229
Actual update 3150 of 5000, Stored update 3150 of 5000
-2 * Loglike = 2591.013616
Actual update 3200 of 5000, Stored update 3200 of 5000
-2 * Loglike = 2597.385478
Actual update 3250 of 5000, Stored update 3250 of 5000
-2 * Loglike = 2590.345531
Actual update 3300 of 5000, Stored update 3300 of 5000
-2 * Loglike = 2595.514683
Actual update 3350 of 5000, Stored update 3350 of 5000
-2 * Loglike = 2589.426875
Actual update 3400 of 5000, Stored update 3400 of 5000
-2 * Loglike = 2592.861416
Actual update 3450 of 5000, Stored update 3450 of 5000
-2 * Loglike = 2592.858537
Actual update 3500 of 5000, Stored update 3500 of 5000
-2 * Loglike = 2588.138266
Actual update 3550 of 5000, Stored update 3550 of 5000
-2 * Loglike = 2589.571574
Actual update 3600 of 5000, Stored update 3600 of 5000
-2 * Loglike = 2589.433225
Actual update 3650 of 5000, Stored update 3650 of 5000
-2 * Loglike = 2590.601317
Actual update 3700 of 5000, Stored update 3700 of 5000
-2 * Loglike = 2589.952170
Actual update 3750 of 5000, Stored update 3750 of 5000
-2 * Loglike = 2590.395170
Actual update 3800 of 5000, Stored update 3800 of 5000
-2 * Loglike = 2589.826861
Actual update 3850 of 5000, Stored update 3850 of 5000
-2 * Loglike = 2588.102185
Actual update 3900 of 5000, Stored update 3900 of 5000
-2 * Loglike = 2590.297300
Actual update 3950 of 5000, Stored update 3950 of 5000
-2 * Loglike = 2592.796900
Actual update 4000 of 5000, Stored update 4000 of 5000
-2 * Loglike = 2589.577441
Actual update 4050 of 5000, Stored update 4050 of 5000
-2 * Loglike = 2593.019622
Actual update 4100 of 5000, Stored update 4100 of 5000
-2 * Loglike = 2589.501620
Actual update 4150 of 5000, Stored update 4150 of 5000
-2 * Loglike = 2588.790930
Actual update 4200 of 5000, Stored update 4200 of 5000
-2 * Loglike = 2594.319619
Actual update 4250 of 5000, Stored update 4250 of 5000
-2 * Loglike = 2585.809142
Actual update 4300 of 5000, Stored update 4300 of 5000
-2 * Loglike = 2590.817792
Actual update 4350 of 5000, Stored update 4350 of 5000
-2 * Loglike = 2592.989553
Actual update 4400 of 5000, Stored update 4400 of 5000
-2 * Loglike = 2590.839469
Actual update 4450 of 5000, Stored update 4450 of 5000
-2 * Loglike = 2590.779146
Actual update 4500 of 5000, Stored update 4500 of 5000
-2 * Loglike = 2593.870488
Actual update 4550 of 5000, Stored update 4550 of 5000
-2 * Loglike = 2587.615504
Actual update 4600 of 5000, Stored update 4600 of 5000
-2 * Loglike = 2586.987066
Actual update 4650 of 5000, Stored update 4650 of 5000
-2 * Loglike = 2591.021368
Actual update 4700 of 5000, Stored update 4700 of 5000
-2 * Loglike = 2588.872930
Actual update 4750 of 5000, Stored update 4750 of 5000
-2 * Loglike = 2590.412870
Actual update 4800 of 5000, Stored update 4800 of 5000
-2 * Loglike = 2592.596384
Actual update 4850 of 5000, Stored update 4850 of 5000
-2 * Loglike = 2591.742665
Actual update 4900 of 5000, Stored update 4900 of 5000
-2 * Loglike = 2596.675991
Actual update 4950 of 5000, Stored update 4950 of 5000
-2 * Loglike = 2588.900644
Actual update 5000 of 5000, Stored update 5000 of 5000
-2 * Loglike = 2594.443487
PUPN c1003 c1004
AVER c1091 b99 b100

Count   =           5000
Average =          2590.9
S.D.    =          2.5187
S.E.M.  =        0.035619

ECHO 0
Execution completed


-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
MLwiN (version: 3.09)  multilevel model (Binomial) 
Estimation algorithm:  MCMC      Elapsed time : 7.24s 
Number of obs:  1934 (from total 1934)          Number of iter.: 5000  Chains: 1  Burn-in: 500 
Deviance statistic:  NA 
--------------------------------------------------------------------------------------------------- 
The model formula:
logit(use) ~ 1 + obsage
Level 1: l1id      
--------------------------------------------------------------------------------------------------- 
The fixed part estimates:  
               Coef.   Std. Err.      z    Pr(>|z|)       [95% Cred.   Interval]    ESS 
Intercept   -0.43839     0.04655  -9.42   4.625e-21  ***    -0.53559    -0.35113   1084 
obsage       0.00674     0.00591   1.14      0.2542         -0.00467     0.01836    843 
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1  
--------------------------------------------------------------------------------------------------- 
The random part estimates at the l1id level: 
                Coef.   Std. Err.   [95% Cred.   Interval]    ESS 
var_bcons_1   1.00000       1e-05      1.00000     1.00000   5000 
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- 
> 
> # 14.4 Measurement errors in more than one variable and
> #      misclassifications . . . . . . . . . . . . . . . . . . . . . . . .211
> 
> # Chapter learning outcomes . . . . . . . . . . . . . . . . . . . . . . .128
> 
> 
> 
> 
> 
> ############################################################################
> 
> proc.time()
   user  system elapsed 
   4.81    0.75   34.20