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

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

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

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

> ############################################################################
> #     MLwiN MCMC Manual
> #
> # 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.11/
To change this use: options(MLwiN_path="<path to MLwiN>")
> # MLwiN folder
> mlwin <- getOption("MLwiN_path")
> while (!file.access(mlwin, mode = 1) == 0) {
+   cat("Please specify the root MLwiN folder or the full path to the MLwiN executable:\n")
+   mlwin <- scan(what = character(0), sep = "\n")
+   mlwin <- gsub("\\", "/", mlwin, fixed = TRUE)
+ }
> options(MLwiN_path = mlwin)
> 
> ## 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.11)  multilevel model (Normal) 
Estimation algorithm:  IGLS        Elapsed time : 0.08s 
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.11)  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.11)  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.11)  multilevel model (Normal) 
Estimation algorithm:  MCMC      Elapsed time : 3.89s 
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.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 2.26s 
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.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 2.17s 
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.11)  multilevel model (Normal) 
        N min     mean max N_complete min_complete mean_complete max_complete
school 65   2 62.44615 198         65            2      62.44615          198
Estimation algorithm:  MCMC      Elapsed time : 7.53s 
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.11)  multilevel model (Binomial) 
Estimation algorithm:  MCMC      Elapsed time : 2.8s 
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.11)  multilevel model (Binomial) 
Estimation algorithm:  MCMC      Elapsed time : 2.92s 
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.11)  multilevel model (Binomial) 
Estimation algorithm:  MCMC      Elapsed time : 7.92s 
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 
   5.21    0.79   37.34