R version 3.6.3 (2020-02-29) -- "Holding the Windsock" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ############################################################################ > # MLwiN MCMC Manual > # > # 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-2017 Zhengzheng Zhang, Christopher M. J. Charlton, Richard M. A. Parker, William J. Browne and George Leckie Support provided by the Economic and Social Research Council (ESRC) (Grants RES-149-25-1084, RES-576-25-0032 and ES/K007246/1) To cite R2MLwiN in publications use: Zhengzheng Zhang, Richard M. A. Parker, Christopher M. J. Charlton, George Leckie, William J. Browne (2016). R2MLwiN: A Package to Run MLwiN from within R. Journal of Statistical Software, 72(10), 1-43. doi:10.18637/jss.v072.i10 A BibTeX entry for LaTeX users is @Article{, title = {{R2MLwiN}: A Package to Run {MLwiN} from within {R}}, author = {Zhengzheng Zhang and Richard M. A. Parker and Christopher M. J. Charlton and George Leckie and William J. Browne}, journal = {Journal of Statistical Software}, year = {2016}, volume = {72}, number = {10}, pages = {1--43}, doi = {10.18637/jss.v072.i10}, } The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.05/ To change this use: options(MLwiN_path="") > # 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.05) 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.05) multilevel model (Normal) Estimation algorithm: IGLS Elapsed time : 0.05s 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.05) multilevel model (Normal) Estimation algorithm: IGLS Elapsed time : 0.05s 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 for 50 iterations out of 500 -2 * Loglike = 9973.123915 Burning in for 100 iterations out of 500 -2 * Loglike = 9773.645338 Burning in for 150 iterations out of 500 -2 * Loglike = 9836.680390 Burning in for 200 iterations out of 500 -2 * Loglike = 9902.678258 Burning in for 250 iterations out of 500 -2 * Loglike = 9882.388699 Burning in for 300 iterations out of 500 -2 * Loglike = 9810.770999 Burning in for 350 iterations out of 500 -2 * Loglike = 9867.724807 Burning in for 400 iterations out of 500 -2 * Loglike = 9818.774370 Burning in for 450 iterations out of 500 -2 * Loglike = 9818.068350 Burning in for 500 iterations out of 500 -2 * Loglike = 9799.382025 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 = 9781.244691 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9752.228870 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9814.378974 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9752.635050 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9785.182159 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9850.599135 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9833.927048 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9810.716117 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9854.570985 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9870.894503 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9857.733413 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9798.816474 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9791.906287 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9870.674753 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9826.987486 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9784.336287 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9796.037524 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9803.527817 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9825.219009 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9882.614832 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9762.567739 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9851.999867 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9816.027586 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9819.890513 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9794.445332 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9832.736911 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9823.976816 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9803.015078 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9889.513035 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9811.139520 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9772.377675 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9864.801490 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9810.611717 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9850.534842 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9808.494448 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9879.278502 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9838.799561 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9873.424801 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9827.016801 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9855.329727 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9814.788326 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9823.390745 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9847.456335 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9865.245926 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9816.840221 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9833.022073 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9793.182547 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9916.821790 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9812.998741 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9882.033478 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9881.986822 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9898.947755 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9904.201909 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9791.212674 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9898.290017 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9776.170938 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9823.732559 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9832.787440 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9831.074783 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9885.489347 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9830.228442 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9914.820803 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9839.744360 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9929.118328 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9898.284711 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9829.582693 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9898.523978 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9764.458796 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9877.961912 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9800.133825 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9820.204186 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9807.172311 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9836.474018 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9892.254142 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9864.372643 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9858.586871 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9897.579628 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9834.419321 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9876.760126 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9857.945675 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9908.496879 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9857.243708 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9919.866119 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9776.363477 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9804.700535 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9856.272757 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9811.202091 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9857.958731 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9802.562444 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9822.545879 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9805.894682 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9881.076264 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9906.104797 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9770.269127 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9881.345269 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9797.872894 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9842.065541 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9801.744945 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9893.289343 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9820.367283 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9835.7 S.D. = 41.602 S.E.M. = 0.58835 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.05) multilevel model (Normal) Estimation algorithm: MCMC Elapsed time : 4.65s 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.00078 0.01345 -0.06 0.9539 -0.02667 0.02567 4113 obslrt 0.58239 0.01485 39.22 0 *** 0.55423 0.61137 2826 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.66097 0.01615 0.62976 0.69297 3397 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > # 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 for 50 iterations out of 500 -2 * Loglike = 9145.871584 Burning in for 100 iterations out of 500 -2 * Loglike = 9132.292720 Burning in for 150 iterations out of 500 -2 * Loglike = 9116.964735 Burning in for 200 iterations out of 500 -2 * Loglike = 9109.819138 Burning in for 250 iterations out of 500 -2 * Loglike = 9142.444254 Burning in for 300 iterations out of 500 -2 * Loglike = 9126.443022 Burning in for 350 iterations out of 500 -2 * Loglike = 9158.382599 Burning in for 400 iterations out of 500 -2 * Loglike = 9130.050027 Burning in for 450 iterations out of 500 -2 * Loglike = 9132.164613 Burning in for 500 iterations out of 500 -2 * Loglike = 9095.782185 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 = 9141.251505 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9092.112338 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9114.453660 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9130.042851 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9112.630195 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9127.189612 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9121.760133 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9095.484118 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9119.696741 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9097.278119 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9091.835441 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9147.067486 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9115.848152 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9136.364432 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9122.999236 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9125.626814 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9124.942111 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9123.715072 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9139.162535 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9093.087478 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9103.780217 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9125.714840 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9117.856745 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9104.224610 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9136.726333 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9115.553162 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9135.082659 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9106.260473 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9113.023248 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9128.720286 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9128.243265 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9112.524919 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9116.163849 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9100.141526 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9112.566069 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9129.956343 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9131.034399 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9121.937589 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9094.252149 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9117.545274 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9099.416940 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9132.162100 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9135.159139 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9110.982424 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9112.731006 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9129.335122 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9135.832789 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9115.598454 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9119.806787 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9106.961083 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9094.496029 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9114.022952 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9122.506561 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9108.681836 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9110.341205 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9131.645275 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9143.193063 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9096.461754 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9142.098601 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9090.543006 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9123.431925 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9100.574500 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9135.650108 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9131.254028 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9130.803281 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9137.320115 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9151.055699 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9113.381797 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9136.395805 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9119.182591 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9099.498807 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9094.669026 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9124.954760 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9113.843599 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9122.879346 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9138.609563 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9106.548693 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9104.254128 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9117.331443 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9094.789350 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9135.693654 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9101.392385 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9114.650394 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9113.173542 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9126.374967 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9142.771045 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9133.755790 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9146.895681 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9130.279064 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9126.639139 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9112.546091 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9130.121874 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9102.731900 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9101.131722 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9081.592796 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9107.770592 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9129.044185 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9118.474018 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9118.583370 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9132.471942 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9122.7 S.D. = 16.894 S.E.M. = 0.23892 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.05) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 1.96s 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 9122.672 9031.177 91.496 9214.168 --------------------------------------------------------------------------------------------------- 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.01325 0.03984 -0.33 0.7395 -0.08901 0.07336 225 standlrt 0.55687 0.02033 27.39 4.031e-165 *** 0.51599 0.59630 563 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.09706 0.02004 0.06544 0.14260 3136 cov_Intercept_standlrt 0.01955 0.00737 0.00647 0.03595 1644 var_standlrt 0.01549 0.00484 0.00804 0.02681 930 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. [95% Cred. Interval] ESS var_Intercept 0.55432 0.01247 0.53027 0.57958 4575 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > (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 for 50 iterations out of 500 -2 * Loglike = 9563.551994 Burning in for 100 iterations out of 500 -2 * Loglike = 9551.869706 Burning in for 150 iterations out of 500 -2 * Loglike = 9540.142431 Burning in for 200 iterations out of 500 -2 * Loglike = 9529.532910 Burning in for 250 iterations out of 500 -2 * Loglike = 9565.554519 Burning in for 300 iterations out of 500 -2 * Loglike = 9547.181186 Burning in for 350 iterations out of 500 -2 * Loglike = 9581.360863 Burning in for 400 iterations out of 500 -2 * Loglike = 9556.335165 Burning in for 450 iterations out of 500 -2 * Loglike = 9553.004192 Burning in for 500 iterations out of 500 -2 * Loglike = 9519.383411 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 = 9567.138560 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9510.322185 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9540.128711 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9556.863836 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9533.140188 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9548.930489 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9540.267426 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9511.213504 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9539.121927 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9511.377877 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9513.618935 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9565.434894 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9536.710660 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9553.272149 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9539.491449 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9541.214855 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9546.394392 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9541.206324 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9558.995171 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9510.058732 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9522.336055 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9552.663122 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9543.151269 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9528.914709 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9549.666090 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9533.860812 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9552.501978 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9529.248909 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9539.545067 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9552.254348 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9556.977622 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9540.656971 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9533.443215 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9527.842694 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9536.620302 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9548.861469 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9556.116918 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9540.985751 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9519.432812 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9527.122047 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9521.620350 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9553.034826 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9557.490654 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9529.135733 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9532.092805 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9544.196905 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9557.878075 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9533.469696 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9535.903317 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9531.140925 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9517.062521 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9538.034397 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9543.757858 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9525.729624 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9530.017209 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9549.205874 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9565.466221 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9516.686220 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9568.974273 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9510.431241 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9547.673224 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9519.123479 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9550.481660 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9549.280389 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9547.539287 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9559.387211 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9565.804720 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9529.451395 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9565.709388 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9537.295629 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9521.661811 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9519.091064 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9551.781795 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9531.282340 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9546.119340 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9562.439026 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9525.432763 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9522.752703 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9538.576919 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9515.999066 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9556.598016 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9522.649596 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9532.958609 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9533.789928 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9544.539244 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9562.923948 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9552.283234 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9560.913974 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9551.833434 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9543.183140 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9534.137845 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9549.803868 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9524.314743 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9519.928389 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9505.597018 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9526.670073 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9552.302909 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9544.637780 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9545.238917 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9547.621105 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9543.1 S.D. = 17.079 S.E.M. = 0.24154 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.05) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 1.97s 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.105 9453.964 89.141 9632.246 --------------------------------------------------------------------------------------------------- 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.01443 0.04222 -0.34 0.7325 -0.09475 0.07697 222 obslrt 0.44706 0.01815 24.63 5.829e-134 *** 0.41024 0.48176 678 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.11011 0.02258 0.07411 0.16072 3187 cov_Intercept_obslrt 0.01686 0.00687 0.00462 0.03195 1458 var_obslrt 0.01124 0.00390 0.00532 0.02044 743 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. [95% Cred. Interval] ESS var_Intercept 0.61481 0.01383 0.58817 0.64291 4569 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > (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 for 50 iterations out of 500 -2 * Loglike = 9304.594936 Burning in for 100 iterations out of 500 -2 * Loglike = 9142.963882 Burning in for 150 iterations out of 500 -2 * Loglike = 9223.639991 Burning in for 200 iterations out of 500 -2 * Loglike = 9077.813808 Burning in for 250 iterations out of 500 -2 * Loglike = 9126.651636 Burning in for 300 iterations out of 500 -2 * Loglike = 9206.756947 Burning in for 350 iterations out of 500 -2 * Loglike = 9140.663043 Burning in for 400 iterations out of 500 -2 * Loglike = 9060.584955 Burning in for 450 iterations out of 500 -2 * Loglike = 9185.922804 Burning in for 500 iterations out of 500 -2 * Loglike = 9200.587822 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.646844 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9163.661841 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9161.941051 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9109.203855 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9112.079916 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9157.513451 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9127.602336 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9071.730703 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9179.602694 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9099.570732 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9157.243076 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9161.637534 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9110.876319 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9144.423149 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9135.893895 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9124.894502 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9129.942191 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9120.786838 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9093.433572 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9103.421697 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9105.698191 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9283.230257 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9162.351491 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9186.635968 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9145.152566 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9163.925509 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9177.123521 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9112.319794 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9109.167638 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9132.284819 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9146.217163 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9018.938462 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9167.716585 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9146.712835 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9165.320042 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9191.209884 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9039.447638 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9104.808624 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9135.413319 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9215.787723 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9101.081884 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9143.765793 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9171.958242 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9131.212001 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9098.579377 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9070.268990 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9019.692183 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9155.669370 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9115.650391 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9134.096415 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9055.297621 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9161.189830 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9143.521445 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9191.333757 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9147.813086 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9102.893528 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9063.826064 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9180.618616 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9148.251794 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9169.093911 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9083.307666 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9181.544045 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9138.599441 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9131.672831 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9243.618561 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9220.741287 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9149.695346 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9117.109333 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9148.993357 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9151.487883 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9136.152414 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9151.518971 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9154.099566 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9237.968724 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9116.878881 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9220.455956 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9110.498814 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9244.481229 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9153.048657 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9110.848027 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9168.632670 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9196.249179 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9117.137655 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9153.306973 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9029.657950 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9137.165515 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9223.198738 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9113.596257 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9129.847243 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9211.235763 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9157.244759 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9160.408190 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9143.180550 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9107.878221 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9212.615510 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9109.064464 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9179.576662 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9173.206628 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9102.269183 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9158.071413 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9142.4 S.D. = 49.264 S.E.M. = 0.69669 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.05) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 7s 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.01761 0.04418 -0.40 0.6902 -0.10123 0.07178 179 obslrt 0.53636 0.02463 21.77 4.07e-105 *** 0.48947 0.58733 388 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.10981 0.02321 0.07246 0.16314 3227 cov_Intercept_obslrt 0.02276 0.00915 0.00695 0.04337 1747 var_obslrt 0.02149 0.00677 0.01072 0.03680 868 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. [95% Cred. Interval] ESS var_Intercept 0.55711 0.01437 0.52966 0.58476 3093 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > # 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 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.05) multilevel model (Binomial) Estimation algorithm: MCMC Elapsed time : 2.03s 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 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.05) multilevel model (Binomial) Estimation algorithm: MCMC Elapsed time : 2.05s 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 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.05) multilevel model (Binomial) Estimation algorithm: MCMC Elapsed time : 7.06s 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.32 1.03 34.61