R version 4.3.2 (2023-10-31 ucrt) -- "Eye Holes" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ############################################################################ > # MLwiN MCMC Manual > # > # 6 Random Slopes Regression Models . . . . . . . . . . . . . . . . . . 71 > # > # Browne, W.J. (2009) MCMC Estimation in MLwiN, v2.13. Centre for > # Multilevel Modelling, University of Bristol. > ############################################################################ > # R script to replicate all analyses using R2MLwiN > # > # Zhang, Z., Charlton, C., Parker, R, Leckie, G., and Browne, W.J. > # Centre for Multilevel Modelling, 2012 > # http://www.bristol.ac.uk/cmm/software/R2MLwiN/ > ############################################################################ > > library(R2MLwiN) R2MLwiN: A package to run models implemented in MLwiN from R Copyright 2013-2024 Zhengzheng Zhang, Christopher M. J. Charlton, Richard M. A. Parker, William J. Browne and George Leckie Support provided by the Economic and Social Research Council (ESRC) (Grants RES-149-25-1084, RES-576-25-0032 and ES/K007246/1) To cite R2MLwiN in publications use: Zhengzheng Zhang, Richard M. A. Parker, Christopher M. J. Charlton, George Leckie, William J. Browne (2016). R2MLwiN: A Package to Run MLwiN from within R. Journal of Statistical Software, 72(10), 1-43. doi:10.18637/jss.v072.i10 A BibTeX entry for LaTeX users is @Article{, title = {{R2MLwiN}: A Package to Run {MLwiN} from within {R}}, author = {Zhengzheng Zhang and Richard M. A. Parker and Christopher M. J. Charlton and George Leckie and William J. Browne}, journal = {Journal of Statistical Software}, year = {2016}, volume = {72}, number = {10}, pages = {1--43}, doi = {10.18637/jss.v072.i10}, } The MLwiN_path option is currently set to C:/Program Files/MLwiN v3.09/ To change this use: options(MLwiN_path="") > # 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") > > ## Choose MCMC algoritm for estimation (IGLS will be used to obtain starting values for MCMC) > (mymodel <- runMLwiN(normexam ~ 1 + standlrt + school + school:standlrt + (1 | student), estoptions = list(EstM = 1), + data = tutorial)) MLwiN is running, please wait...... /nogui option ignored ECHO 0 Echoing is ON BATC 1 Batch mode is ON MAXI 2 STAR iteration 0 iteration 1 Convergence achieved TOLE 2 MAXI 20 BATC 1 Batch mode is ON NEXT iteration 2 Convergence achieved ECHO 0 Echoing is ON MCMC 0 500 1 5.8 50 10 1 1 1 1 1 1 BURNING IN... Burning in for 50 iterations out of 500 -2 * Loglike = 9093.418531 Burning in for 100 iterations out of 500 -2 * Loglike = 9092.900586 Burning in for 150 iterations out of 500 -2 * Loglike = 9118.380097 Burning in for 200 iterations out of 500 -2 * Loglike = 9118.685771 Burning in for 250 iterations out of 500 -2 * Loglike = 9105.378461 Burning in for 300 iterations out of 500 -2 * Loglike = 9105.934099 Burning in for 350 iterations out of 500 -2 * Loglike = 9134.500789 Burning in for 400 iterations out of 500 -2 * Loglike = 9138.704459 Burning in for 450 iterations out of 500 -2 * Loglike = 9159.735979 Burning in for 500 iterations out of 500 -2 * Loglike = 9129.797786 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 = 9091.857742 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9139.149573 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9114.754009 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9110.714071 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9126.173028 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9116.028105 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9113.275752 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9124.253577 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9117.518915 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9131.493118 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9105.840549 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9123.911967 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9103.578211 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9116.991012 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9117.854430 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9106.765475 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9144.235795 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9099.419419 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9127.688079 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9158.443020 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9094.865749 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9116.082855 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9117.468759 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9127.819119 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9123.786680 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9097.689286 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9128.028155 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9128.326314 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9135.997273 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9140.903782 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9118.373033 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9104.554069 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9123.293877 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9107.257709 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9116.612875 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9109.452977 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9133.951513 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9124.577192 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9116.301573 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9105.496177 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9102.734252 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9158.727444 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9115.842466 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9130.662720 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9115.131284 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9116.122359 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9118.500282 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9112.228719 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9100.720416 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9120.453749 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9107.097140 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9096.640975 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9111.706882 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9112.883092 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9084.434828 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9110.722011 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9128.702644 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9114.949033 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9122.983417 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9124.000837 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9113.898441 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9096.052350 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9096.346731 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9108.993970 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9084.693941 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9095.251089 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9120.634017 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9134.141289 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9098.777982 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9126.291347 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9115.912032 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9109.901836 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9100.639733 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9114.264070 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9109.996807 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9121.732075 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9118.397513 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9129.093108 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9125.064921 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9135.931942 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9135.830370 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9118.841312 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9129.628546 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9158.583421 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9128.597314 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9141.338552 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9105.516208 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9118.578704 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9143.639765 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9121.966083 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9111.260866 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9084.455123 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9107.838627 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9143.512938 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9135.113191 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9105.161130 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9106.430748 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9113.083779 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9090.183744 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9089.552195 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9118.2 S.D. = 16.250 S.E.M. = 0.22980 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.09) multilevel model (Normal) Estimation algorithm: MCMC Elapsed time : 13.39s 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 9118.244 8987.053 131.191 9249.435 --------------------------------------------------------------------------------------------------- The model formula: normexam ~ 1 + standlrt + school + school:standlrt + (1 | student) Level 1: student --------------------------------------------------------------------------------------------------- The fixed part estimates: Coef. Std. Err. z Pr(>|z|) [95% Cred. Interval] ESS Intercept 0.38271 0.08830 4.33 1.463e-05 *** 0.20691 0.55172 5000 standlrt 0.70866 0.08359 8.48 2.286e-17 *** 0.54932 0.87178 4405 school2 0.10153 0.13727 0.74 0.4595 -0.15955 0.37598 5000 school3 0.17442 0.14755 1.18 0.2372 -0.11027 0.46765 5218 school4 -0.37771 0.12128 -3.11 0.001843 ** -0.61238 -0.14086 5000 school5 -0.12408 0.15923 -0.78 0.4358 -0.42466 0.19109 4765 school6 0.22047 0.13557 1.63 0.1039 -0.05312 0.49925 5000 school7 0.01391 0.11747 0.12 0.9057 -0.21128 0.24506 5000 school8 -0.40751 0.11446 -3.56 0.0003702 *** -0.63371 -0.18199 5000 school9 -0.61942 0.16529 -3.75 0.0001786 *** -0.93911 -0.29622 4799 school10 -0.70629 0.13734 -5.14 2.706e-07 *** -0.97614 -0.43423 5000 school11 -0.11689 0.14193 -0.82 0.4102 -0.39734 0.15586 5000 school12 -0.45226 0.13969 -3.24 0.001206 ** -0.72349 -0.17598 5000 school13 -0.53425 0.12955 -4.12 3.725e-05 *** -0.78283 -0.27958 4802 school14 -0.56553 0.10432 -5.42 5.915e-08 *** -0.76867 -0.35989 5000 school15 -0.62306 0.12054 -5.17 2.353e-07 *** -0.85031 -0.38458 5000 school16 -0.76907 0.12358 -6.22 4.862e-10 *** -1.01683 -0.53227 5000 school17 -0.56996 0.11110 -5.13 2.898e-07 *** -0.78873 -0.35469 5000 school18 -0.43912 0.11322 -3.88 0.0001051 *** -0.65609 -0.21582 5000 school19 -0.48033 0.14524 -3.31 0.0009422 *** -0.76378 -0.19227 4839 school20 -0.11167 0.15746 -0.71 0.4782 -0.42450 0.20309 5000 school21 -0.11323 0.12796 -0.88 0.3762 -0.36221 0.13685 4912 school22 -0.84839 0.11804 -7.19 6.624e-13 *** -1.07823 -0.62132 4532 school23 -1.02078 0.16900 -6.04 1.541e-09 *** -1.35310 -0.69083 5518 school24 -0.20243 0.15984 -1.27 0.2054 -0.51407 0.10435 5000 school25 -0.65219 0.13593 -4.80 1.603e-06 *** -0.91858 -0.38347 5000 school26 -0.41732 0.13773 -3.03 0.002445 ** -0.69391 -0.15562 5000 school27 -0.35167 0.16693 -2.11 0.03514 * -0.67248 -0.02402 5278 school28 -1.15058 0.13629 -8.44 3.111e-17 *** -1.42227 -0.89001 5000 school29 -0.17831 0.12647 -1.41 0.1586 -0.42835 0.06649 5000 school30 -0.26454 0.14924 -1.77 0.07629 . -0.55834 0.02926 5000 school31 -0.42110 0.16026 -2.63 0.0086 ** -0.72622 -0.10387 4681 school32 -0.32350 0.15693 -2.06 0.03927 * -0.63081 -0.01003 5000 school33 -0.34249 0.12076 -2.84 0.004568 ** -0.58280 -0.11190 5000 school34 -0.49188 0.17873 -2.75 0.005923 ** -0.84768 -0.13590 4633 school35 -0.25388 0.15104 -1.68 0.09279 . -0.55939 0.04208 5000 school36 -0.59038 0.12710 -4.64 3.402e-06 *** -0.83935 -0.34724 5000 school37 -0.86927 0.21226 -4.10 4.215e-05 *** -1.28587 -0.45619 5000 school38 -0.53762 0.13743 -3.91 9.155e-05 *** -0.80212 -0.27754 5000 school39 -0.25183 0.14214 -1.77 0.07644 . -0.53325 0.02748 5000 school40 -0.63313 0.12551 -5.04 4.546e-07 *** -0.87256 -0.38803 5000 school41 -0.16816 0.13081 -1.29 0.1986 -0.42451 0.09057 5281 school42 -0.30032 0.13435 -2.24 0.02539 * -0.56023 -0.03373 5000 school43 -0.55221 0.13998 -3.94 7.987e-05 *** -0.82695 -0.27764 5302 school44 -0.72795 0.16690 -4.36 1.291e-05 *** -1.05648 -0.40477 5212 school45 -0.50330 0.13843 -3.64 0.0002771 *** -0.77450 -0.23208 4225 school46 -0.77132 0.12093 -6.38 1.793e-10 *** -1.00380 -0.53146 5000 school47 -0.41330 0.12130 -3.41 0.0006564 *** -0.64475 -0.17346 5000 school48 -3.68837 5.32817 -0.69 0.4888 -14.11963 6.84850 5218 school49 -0.33544 0.11323 -2.96 0.003051 ** -0.55636 -0.11335 5000 school50 -0.70863 0.12291 -5.77 8.139e-09 *** -0.94798 -0.46708 4861 school51 -0.52511 0.13712 -3.83 0.0001284 *** -0.78197 -0.25762 5000 school52 0.00464 0.13169 0.04 0.9719 -0.25335 0.27030 5820 school53 0.21353 0.13141 1.62 0.1042 -0.04643 0.46998 5224 school54 -1.04554 0.38320 -2.73 0.006364 ** -1.78573 -0.28862 5000 school55 0.16860 0.14045 1.20 0.23 -0.09525 0.45207 5000 school56 -0.34410 0.14914 -2.31 0.02104 * -0.63614 -0.03996 5000 school57 -0.34234 0.13084 -2.62 0.008882 ** -0.59461 -0.08727 5697 school58 -0.17603 0.15290 -1.15 0.2496 -0.46937 0.13334 5000 school59 -1.23429 0.15498 -7.96 1.661e-15 *** -1.53503 -0.92901 5000 school60 -0.13018 0.12290 -1.06 0.2895 -0.36911 0.10737 5000 school61 -0.42106 0.12852 -3.28 0.001052 ** -0.66691 -0.17119 5000 school62 -0.43537 0.12516 -3.48 0.0005044 *** -0.67947 -0.19233 5000 school63 0.30099 0.16364 1.84 0.06587 . -0.02332 0.61095 4434 school64 -0.34527 0.14024 -2.46 0.01382 * -0.61800 -0.07258 4974 school65 -0.55862 0.12088 -4.62 3.814e-06 *** -0.79556 -0.31851 5634 standlrt:school2 0.05226 0.12267 0.43 0.6701 -0.19517 0.28807 4735 standlrt:school3 -0.12948 0.14564 -0.89 0.374 -0.41738 0.15885 4282 standlrt:school4 0.05016 0.11754 0.43 0.6696 -0.18315 0.27786 5000 standlrt:school5 -0.02736 0.19718 -0.14 0.8897 -0.41440 0.37434 5000 standlrt:school6 -0.17468 0.12761 -1.37 0.171 -0.41653 0.07711 5000 standlrt:school7 -0.46607 0.12446 -3.74 0.0001805 *** -0.71707 -0.22575 5000 standlrt:school8 -0.14259 0.10633 -1.34 0.1799 -0.34966 0.06581 4387 standlrt:school9 -0.30727 0.13924 -2.21 0.02733 * -0.57804 -0.03744 5000 standlrt:school10 -0.41294 0.16498 -2.50 0.01231 * -0.73154 -0.09508 5000 standlrt:school11 -0.24883 0.12726 -1.96 0.05054 . -0.50192 -0.00755 5000 standlrt:school12 -0.23097 0.15459 -1.49 0.1352 -0.53862 0.06827 5287 standlrt:school13 -0.08266 0.12292 -0.67 0.5013 -0.32303 0.15919 5000 standlrt:school14 -0.10264 0.10097 -1.02 0.3094 -0.30301 0.09606 4768 standlrt:school15 0.02731 0.12869 0.21 0.8319 -0.22596 0.27944 5000 standlrt:school16 -0.30086 0.12884 -2.34 0.01953 * -0.55951 -0.04604 5030 standlrt:school17 -0.21106 0.10317 -2.05 0.04079 * -0.41235 -0.01114 4475 standlrt:school18 -0.34935 0.13606 -2.57 0.01024 * -0.62009 -0.08228 5000 standlrt:school19 0.08373 0.17520 0.48 0.6327 -0.26263 0.43062 5000 standlrt:school20 -0.18960 0.14532 -1.30 0.192 -0.47697 0.09308 4814 standlrt:school21 -0.15535 0.12924 -1.20 0.2293 -0.40602 0.09962 5000 standlrt:school22 -0.18216 0.11219 -1.62 0.1045 -0.39863 0.03815 4771 standlrt:school23 -0.31964 0.15845 -2.02 0.04367 * -0.62885 -0.00364 5000 standlrt:school24 -0.31331 0.15838 -1.98 0.04791 * -0.62482 -0.00934 5000 standlrt:school25 -0.17697 0.11786 -1.50 0.1332 -0.41141 0.05372 5000 standlrt:school26 -0.16377 0.12511 -1.31 0.1905 -0.40904 0.07793 4598 standlrt:school27 -0.13750 0.14882 -0.92 0.3555 -0.43241 0.14757 5000 standlrt:school28 -0.42468 0.14009 -3.03 0.002434 ** -0.69977 -0.15780 5000 standlrt:school29 -0.30777 0.12403 -2.48 0.01308 * -0.54425 -0.06700 5000 standlrt:school30 0.09813 0.13098 0.75 0.4537 -0.15524 0.35442 5000 standlrt:school31 -0.30210 0.17715 -1.71 0.08813 . -0.65689 0.04046 4711 standlrt:school32 -0.04495 0.13342 -0.34 0.7362 -0.31170 0.21727 5630 standlrt:school33 -0.20371 0.13476 -1.51 0.1306 -0.46684 0.06263 5000 standlrt:school34 0.02214 0.13675 0.16 0.8714 -0.24143 0.29129 4926 standlrt:school35 -0.31009 0.17941 -1.73 0.08392 . -0.66894 0.03555 5000 standlrt:school36 -0.27160 0.11192 -2.43 0.01523 * -0.49529 -0.05391 4534 standlrt:school37 -0.47284 0.16961 -2.79 0.005306 ** -0.80965 -0.14025 5000 standlrt:school38 -0.08749 0.12810 -0.68 0.4946 -0.33842 0.16468 5000 standlrt:school39 -0.25462 0.12840 -1.98 0.04736 * -0.49816 -0.00278 5000 standlrt:school40 0.01589 0.11842 0.13 0.8933 -0.22020 0.24658 5215 standlrt:school41 -0.22346 0.12989 -1.72 0.08537 . -0.48211 0.02702 5000 standlrt:school42 -0.29495 0.14374 -2.05 0.04017 * -0.57601 -0.01537 5000 standlrt:school43 0.02279 0.14718 0.15 0.8769 -0.26333 0.30472 5000 standlrt:school44 -0.33957 0.15752 -2.16 0.0311 * -0.64449 -0.03109 5273 standlrt:school45 -0.14094 0.13705 -1.03 0.3038 -0.41090 0.12930 4856 standlrt:school46 -0.24479 0.12353 -1.98 0.04752 * -0.48822 -0.00520 5156 standlrt:school47 -0.04506 0.12764 -0.35 0.7241 -0.29423 0.20342 5000 standlrt:school48 -7.68741 12.79794 -0.60 0.5481 -32.71448 17.90470 5222 standlrt:school49 -0.22227 0.11456 -1.94 0.05235 . -0.44621 -0.00070 4646 standlrt:school50 -0.03298 0.12346 -0.27 0.7894 -0.27853 0.21116 5000 standlrt:school51 -0.36964 0.11934 -3.10 0.001952 ** -0.60903 -0.13759 5000 standlrt:school52 0.03102 0.12579 0.25 0.8052 -0.21896 0.27835 4772 standlrt:school53 0.36912 0.12898 2.86 0.004212 ** 0.12021 0.62302 3643 standlrt:school54 -0.65808 0.46565 -1.41 0.1576 -1.59369 0.24131 4934 standlrt:school55 -0.08303 0.13494 -0.62 0.5383 -0.34431 0.18281 5000 standlrt:school56 0.16623 0.15945 1.04 0.2972 -0.14419 0.47501 5000 standlrt:school57 -0.09614 0.12963 -0.74 0.4583 -0.34438 0.16009 5000 standlrt:school58 -0.35361 0.16785 -2.11 0.03514 * -0.68121 -0.02570 5000 standlrt:school59 -0.34581 0.14469 -2.39 0.01684 * -0.63027 -0.06631 5000 standlrt:school60 -0.07452 0.12492 -0.60 0.5508 -0.32392 0.17191 5000 standlrt:school61 -0.07192 0.12526 -0.57 0.5659 -0.31168 0.17209 5000 standlrt:school62 -0.16317 0.13989 -1.17 0.2435 -0.43722 0.11427 5000 standlrt:school63 -0.37337 0.19961 -1.87 0.06142 . -0.75748 0.02473 5000 standlrt:school64 -0.00309 0.13196 -0.02 0.9813 -0.26321 0.25776 5000 standlrt:school65 -0.14124 0.11689 -1.21 0.2269 -0.37096 0.08993 4708 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.55376 0.01250 0.52986 0.57890 5000 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > ## Define the model Choose IGLS algoritm for estimation Fit the model > (mymodel0a <- runMLwiN(normexam ~ 1 + standlrt + (1 + standlrt | school) + (1 | student), data = tutorial)) /nogui option ignored ECHO 0 Echoing is ON BATC 1 Batch mode is ON MAXI 2 STAR iteration 0 iteration 1 Convergence not achieved TOLE 2 MAXI 20 NEXT iteration 2 iteration 3 Convergence achieved ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.09) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: IGLS Elapsed time : 0.1s Number of obs: 4059 (from total 4059) The model converged after 4 iterations. Log likelihood: -4658.4 Deviance statistic: 9316.9 --------------------------------------------------------------------------------------------------- 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% Conf. Interval] Intercept -0.01151 0.03978 -0.29 0.7724 -0.08948 0.06647 standlrt 0.55673 0.01994 27.92 1.344e-171 *** 0.51765 0.59581 Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 --------------------------------------------------------------------------------------------------- The random part estimates at the school level: Coef. Std. Err. var_Intercept 0.09044 0.01792 cov_Intercept_standlrt 0.01804 0.00672 var_standlrt 0.01454 0.00441 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. var_Intercept 0.55366 0.01248 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > ## Choose MCMC algoritm for estimation (IGLS will be used to obtain starting values for MCMC) > (mymodel0 <- runMLwiN(normexam ~ 1 + standlrt + (1 + standlrt | school) + (1 | student), estoptions = list(EstM = 1), + data = tutorial)) MLwiN is running, please wait...... /nogui option ignored ECHO 0 Echoing is ON BATC 1 Batch mode is ON MAXI 2 STAR iteration 0 iteration 1 Convergence not achieved TOLE 2 MAXI 20 BATC 1 Batch mode is ON NEXT iteration 2 iteration 3 Convergence achieved ECHO 0 Echoing is ON MCMC 0 500 1 5.8 50 10 G30[1] G30[2] 1 1 1 1 1 1 BURNING IN... Burning in for 50 iterations out of 500 -2 * Loglike = 9136.615416 Burning in for 100 iterations out of 500 -2 * Loglike = 9129.472590 Burning in for 150 iterations out of 500 -2 * Loglike = 9137.752388 Burning in for 200 iterations out of 500 -2 * Loglike = 9149.641457 Burning in for 250 iterations out of 500 -2 * Loglike = 9113.737768 Burning in for 300 iterations out of 500 -2 * Loglike = 9127.854908 Burning in for 350 iterations out of 500 -2 * Loglike = 9099.670195 Burning in for 400 iterations out of 500 -2 * Loglike = 9107.133017 Burning in for 450 iterations out of 500 -2 * Loglike = 9108.686421 Burning in for 500 iterations out of 500 -2 * Loglike = 9142.796685 ERAS G30 LINK 0 G30 ERAS c1090 c1091 MCMC 1 5000 1 c1090 c1091 c1003 c1004 1 1 Actual update 50 of 5000, Stored update 50 of 5000 -2 * Loglike = 9114.789855 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9128.936360 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9120.626663 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9140.221260 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9133.366931 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9112.633267 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9110.267312 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9134.271816 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9146.903828 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9118.578020 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9132.630258 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9105.115711 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9126.654041 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9138.820724 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9119.239813 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9145.879653 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9132.483550 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9094.657080 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9118.502105 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9094.902922 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9135.645832 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9128.926707 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9104.509733 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9108.989181 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9111.759420 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9103.538840 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9120.142978 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9117.533977 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9114.717113 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9118.929584 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9146.845851 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9126.654978 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9127.854549 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9114.945398 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9156.130168 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9134.732656 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9142.867734 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9137.471434 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9127.702153 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9133.939751 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9107.640871 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9137.519015 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9147.472132 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9104.767457 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9118.396071 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9143.494727 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9130.468641 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9137.527730 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9106.747459 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9154.348069 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9142.340919 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9100.831693 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9131.653716 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9151.916582 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9165.855130 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9125.508254 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9146.519009 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9122.524343 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9128.808395 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9125.880481 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9130.006670 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9116.446060 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9105.444348 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9141.749060 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9122.819747 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9140.592905 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9119.521430 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9137.983211 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9128.855413 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9130.638160 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9119.619021 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9124.980114 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9135.552922 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9127.477138 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9114.227379 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9110.953194 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9118.839082 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9118.763419 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9117.590737 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9106.053640 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9112.274103 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9115.311806 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9122.276612 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9120.226668 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9140.403968 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9153.535294 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9117.253885 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9116.415887 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9136.007855 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9119.650550 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9108.406619 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9104.027872 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9104.517145 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9102.506754 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9125.791725 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9121.949129 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9123.201133 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9116.574595 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9120.080554 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9120.628553 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9123.2 S.D. = 16.792 S.E.M. = 0.23747 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.09) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.22s 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 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > # 6.1 Prediction intervals for a random slopes regression model . . . . . 75 > > ## Save level 2 residual chains > (mymodel <- runMLwiN(normexam ~ 1 + standlrt + (1 + standlrt | school) + (1 | student), estoptions = list(EstM = 1, + mcmcMeth = list(iterations = 5001), resi.store.levs = 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 = 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 5001 1 c1090 c1091 c1003 c1004 1 1 Actual update 50 of 5001, Stored update 50 of 5001 -2 * Loglike = 9114.789855 Actual update 100 of 5001, Stored update 100 of 5001 -2 * Loglike = 9128.936360 Actual update 150 of 5001, Stored update 150 of 5001 -2 * Loglike = 9120.626663 Actual update 200 of 5001, Stored update 200 of 5001 -2 * Loglike = 9140.221260 Actual update 250 of 5001, Stored update 250 of 5001 -2 * Loglike = 9133.366931 Actual update 300 of 5001, Stored update 300 of 5001 -2 * Loglike = 9112.633267 Actual update 350 of 5001, Stored update 350 of 5001 -2 * Loglike = 9110.267312 Actual update 400 of 5001, Stored update 400 of 5001 -2 * Loglike = 9134.271816 Actual update 450 of 5001, Stored update 450 of 5001 -2 * Loglike = 9146.903828 Actual update 500 of 5001, Stored update 500 of 5001 -2 * Loglike = 9118.578020 Actual update 550 of 5001, Stored update 550 of 5001 -2 * Loglike = 9132.630258 Actual update 600 of 5001, Stored update 600 of 5001 -2 * Loglike = 9105.115711 Actual update 650 of 5001, Stored update 650 of 5001 -2 * Loglike = 9126.654041 Actual update 700 of 5001, Stored update 700 of 5001 -2 * Loglike = 9138.820724 Actual update 750 of 5001, Stored update 750 of 5001 -2 * Loglike = 9119.239813 Actual update 800 of 5001, Stored update 800 of 5001 -2 * Loglike = 9145.879653 Actual update 850 of 5001, Stored update 850 of 5001 -2 * Loglike = 9132.483550 Actual update 900 of 5001, Stored update 900 of 5001 -2 * Loglike = 9094.657080 Actual update 950 of 5001, Stored update 950 of 5001 -2 * Loglike = 9118.502105 Actual update 1000 of 5001, Stored update 1000 of 5001 -2 * Loglike = 9094.902922 Actual update 1050 of 5001, Stored update 1050 of 5001 -2 * Loglike = 9135.645832 Actual update 1100 of 5001, Stored update 1100 of 5001 -2 * Loglike = 9128.926707 Actual update 1150 of 5001, Stored update 1150 of 5001 -2 * Loglike = 9104.509733 Actual update 1200 of 5001, Stored update 1200 of 5001 -2 * Loglike = 9108.989181 Actual update 1250 of 5001, Stored update 1250 of 5001 -2 * Loglike = 9111.759420 Actual update 1300 of 5001, Stored update 1300 of 5001 -2 * Loglike = 9103.538840 Actual update 1350 of 5001, Stored update 1350 of 5001 -2 * Loglike = 9120.142978 Actual update 1400 of 5001, Stored update 1400 of 5001 -2 * Loglike = 9117.533977 Actual update 1450 of 5001, Stored update 1450 of 5001 -2 * Loglike = 9114.717113 Actual update 1500 of 5001, Stored update 1500 of 5001 -2 * Loglike = 9118.929584 Actual update 1550 of 5001, Stored update 1550 of 5001 -2 * Loglike = 9146.845851 Actual update 1600 of 5001, Stored update 1600 of 5001 -2 * Loglike = 9126.654978 Actual update 1650 of 5001, Stored update 1650 of 5001 -2 * Loglike = 9127.854549 Actual update 1700 of 5001, Stored update 1700 of 5001 -2 * Loglike = 9114.945398 Actual update 1750 of 5001, Stored update 1750 of 5001 -2 * Loglike = 9156.130168 Actual update 1800 of 5001, Stored update 1800 of 5001 -2 * Loglike = 9134.732656 Actual update 1850 of 5001, Stored update 1850 of 5001 -2 * Loglike = 9142.867734 Actual update 1900 of 5001, Stored update 1900 of 5001 -2 * Loglike = 9137.471434 Actual update 1950 of 5001, Stored update 1950 of 5001 -2 * Loglike = 9127.702153 Actual update 2000 of 5001, Stored update 2000 of 5001 -2 * Loglike = 9133.939751 Actual update 2050 of 5001, Stored update 2050 of 5001 -2 * Loglike = 9107.640871 Actual update 2100 of 5001, Stored update 2100 of 5001 -2 * Loglike = 9137.519015 Actual update 2150 of 5001, Stored update 2150 of 5001 -2 * Loglike = 9147.472132 Actual update 2200 of 5001, Stored update 2200 of 5001 -2 * Loglike = 9104.767457 Actual update 2250 of 5001, Stored update 2250 of 5001 -2 * Loglike = 9118.396071 Actual update 2300 of 5001, Stored update 2300 of 5001 -2 * Loglike = 9143.494727 Actual update 2350 of 5001, Stored update 2350 of 5001 -2 * Loglike = 9130.468641 Actual update 2400 of 5001, Stored update 2400 of 5001 -2 * Loglike = 9137.527730 Actual update 2450 of 5001, Stored update 2450 of 5001 -2 * Loglike = 9106.747459 Actual update 2500 of 5001, Stored update 2500 of 5001 -2 * Loglike = 9154.348069 Actual update 2550 of 5001, Stored update 2550 of 5001 -2 * Loglike = 9142.340919 Actual update 2600 of 5001, Stored update 2600 of 5001 -2 * Loglike = 9100.831693 Actual update 2650 of 5001, Stored update 2650 of 5001 -2 * Loglike = 9131.653716 Actual update 2700 of 5001, Stored update 2700 of 5001 -2 * Loglike = 9151.916582 Actual update 2750 of 5001, Stored update 2750 of 5001 -2 * Loglike = 9165.855130 Actual update 2800 of 5001, Stored update 2800 of 5001 -2 * Loglike = 9125.508254 Actual update 2850 of 5001, Stored update 2850 of 5001 -2 * Loglike = 9146.519009 Actual update 2900 of 5001, Stored update 2900 of 5001 -2 * Loglike = 9122.524343 Actual update 2950 of 5001, Stored update 2950 of 5001 -2 * Loglike = 9128.808395 Actual update 3000 of 5001, Stored update 3000 of 5001 -2 * Loglike = 9125.880481 Actual update 3050 of 5001, Stored update 3050 of 5001 -2 * Loglike = 9130.006670 Actual update 3100 of 5001, Stored update 3100 of 5001 -2 * Loglike = 9116.446060 Actual update 3150 of 5001, Stored update 3150 of 5001 -2 * Loglike = 9105.444348 Actual update 3200 of 5001, Stored update 3200 of 5001 -2 * Loglike = 9141.749060 Actual update 3250 of 5001, Stored update 3250 of 5001 -2 * Loglike = 9122.819747 Actual update 3300 of 5001, Stored update 3300 of 5001 -2 * Loglike = 9140.592905 Actual update 3350 of 5001, Stored update 3350 of 5001 -2 * Loglike = 9119.521430 Actual update 3400 of 5001, Stored update 3400 of 5001 -2 * Loglike = 9137.983211 Actual update 3450 of 5001, Stored update 3450 of 5001 -2 * Loglike = 9128.855413 Actual update 3500 of 5001, Stored update 3500 of 5001 -2 * Loglike = 9130.638160 Actual update 3550 of 5001, Stored update 3550 of 5001 -2 * Loglike = 9119.619021 Actual update 3600 of 5001, Stored update 3600 of 5001 -2 * Loglike = 9124.980114 Actual update 3650 of 5001, Stored update 3650 of 5001 -2 * Loglike = 9135.552922 Actual update 3700 of 5001, Stored update 3700 of 5001 -2 * Loglike = 9127.477138 Actual update 3750 of 5001, Stored update 3750 of 5001 -2 * Loglike = 9114.227379 Actual update 3800 of 5001, Stored update 3800 of 5001 -2 * Loglike = 9110.953194 Actual update 3850 of 5001, Stored update 3850 of 5001 -2 * Loglike = 9118.839082 Actual update 3900 of 5001, Stored update 3900 of 5001 -2 * Loglike = 9118.763419 Actual update 3950 of 5001, Stored update 3950 of 5001 -2 * Loglike = 9117.590737 Actual update 4000 of 5001, Stored update 4000 of 5001 -2 * Loglike = 9106.053640 Actual update 4050 of 5001, Stored update 4050 of 5001 -2 * Loglike = 9112.274103 Actual update 4100 of 5001, Stored update 4100 of 5001 -2 * Loglike = 9115.311806 Actual update 4150 of 5001, Stored update 4150 of 5001 -2 * Loglike = 9122.276612 Actual update 4200 of 5001, Stored update 4200 of 5001 -2 * Loglike = 9120.226668 Actual update 4250 of 5001, Stored update 4250 of 5001 -2 * Loglike = 9140.403968 Actual update 4300 of 5001, Stored update 4300 of 5001 -2 * Loglike = 9153.535294 Actual update 4350 of 5001, Stored update 4350 of 5001 -2 * Loglike = 9117.253885 Actual update 4400 of 5001, Stored update 4400 of 5001 -2 * Loglike = 9116.415887 Actual update 4450 of 5001, Stored update 4450 of 5001 -2 * Loglike = 9136.007855 Actual update 4500 of 5001, Stored update 4500 of 5001 -2 * Loglike = 9119.650550 Actual update 4550 of 5001, Stored update 4550 of 5001 -2 * Loglike = 9108.406619 Actual update 4600 of 5001, Stored update 4600 of 5001 -2 * Loglike = 9104.027872 Actual update 4650 of 5001, Stored update 4650 of 5001 -2 * Loglike = 9104.517145 Actual update 4700 of 5001, Stored update 4700 of 5001 -2 * Loglike = 9102.506754 Actual update 4750 of 5001, Stored update 4750 of 5001 -2 * Loglike = 9125.791725 Actual update 4800 of 5001, Stored update 4800 of 5001 -2 * Loglike = 9121.949129 Actual update 4850 of 5001, Stored update 4850 of 5001 -2 * Loglike = 9123.201133 Actual update 4900 of 5001, Stored update 4900 of 5001 -2 * Loglike = 9116.574595 Actual update 4950 of 5001, Stored update 4950 of 5001 -2 * Loglike = 9120.080554 Actual update 5000 of 5001, Stored update 5000 of 5001 -2 * Loglike = 9120.628553 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5001 Average = 9123.2 S.D. = 16.791 S.E.M. = 0.23743 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.09) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.37s Number of obs: 4059 (from total 4059) Number of iter.: 5001 Chains: 1 Burn-in: 500 Bayesian Deviance Information Criterion (DIC) Dbar D(thetabar) pD DIC 9123.209 9031.714 91.495 9214.704 --------------------------------------------------------------------------------------------------- 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.01256 0.04088 -0.31 0.7586 -0.09054 0.06864 210 standlrt 0.55645 0.02084 26.70 5.159e-157 *** 0.51501 0.59701 597 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 847 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. [95% Cred. Interval] ESS var_Intercept 0.55440 0.01242 0.53053 0.57904 4750 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > predLines(mymodel, xname = "standlrt", lev = 2, selected = NULL, probs = c(0.025, 0.975), legend.space = "right", + legend.ncol = 2) > dev.new() dev.new(): using pdf(file="Rplots1.pdf") > predLines(mymodel, xname = "standlrt", lev = 2, selected = c(30, 44, 53, 59), probs = c(0.025, 0.975)) > > # 6.2 Alternative priors for variance matrices . . . . . . . . . . . . . .78 > > # 6.3 WinBUGS priors (Prior 2) . . . . . . . . . . . . . . . . . . . . . .78 > > ## Change the starting values for Level 2 variance matrix to .1 on diagonal 0 otherwise. > RP.b <- c(0.1, 0, 0.1, 0.554) > names(RP.b) <- c("RP2_var_Intercept", "RP2_cov_Intercept_standlrt", "RP2_var_standlrt", "RP1_var_Intercept") > > ## Fit the model > (mymodel1 <- runMLwiN(normexam ~ 1 + standlrt + (1 + standlrt | school) + (1 | student), estoptions = list(EstM = 1, + startval = list(RP.b = RP.b)), data = tutorial)) MLwiN is running, please wait...... /nogui option ignored ECHO 0 Echoing is ON BATC 1 Batch mode is ON MAXI 2 STAR iteration 0 iteration 1 Convergence not achieved JOIN 0.1 0 0.1 0.554 '_RP_b' 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 = 9100.597867 Burning in for 100 iterations out of 500 -2 * Loglike = 9114.859458 Burning in for 150 iterations out of 500 -2 * Loglike = 9131.060396 Burning in for 200 iterations out of 500 -2 * Loglike = 9139.251421 Burning in for 250 iterations out of 500 -2 * Loglike = 9103.324707 Burning in for 300 iterations out of 500 -2 * Loglike = 9111.514369 Burning in for 350 iterations out of 500 -2 * Loglike = 9095.770644 Burning in for 400 iterations out of 500 -2 * Loglike = 9098.416226 Burning in for 450 iterations out of 500 -2 * Loglike = 9104.599695 Burning in for 500 iterations out of 500 -2 * Loglike = 9136.912799 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 = 9108.783409 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9117.677290 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9110.662119 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9118.702238 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9129.997741 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9105.572549 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9101.071629 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9123.680481 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9138.690457 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9112.370849 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9119.786772 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9095.188890 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9113.284402 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9124.930549 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9116.363410 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9121.332262 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9124.130784 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9094.416030 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9107.261022 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9088.937646 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9117.450175 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9119.190132 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9098.201233 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9096.094239 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9106.264247 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9101.696788 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9113.184011 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9106.601132 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9112.889403 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9114.845772 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9144.617915 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9109.667293 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9121.148947 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9107.170897 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9141.810247 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9126.947055 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9131.819346 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9129.540429 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9118.420232 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9117.427898 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9101.300438 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9138.064959 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9133.628313 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9100.786852 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9116.828274 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9127.704594 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9122.857630 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9124.178782 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9101.451288 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9146.468731 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9128.020840 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9097.882656 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9108.402604 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9137.233359 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9158.027340 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9111.587321 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9131.869404 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9123.846017 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9126.986281 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9123.588155 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9119.549677 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9110.668430 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9100.801380 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9120.272246 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9118.804743 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9131.912345 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9108.860416 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9134.331512 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9118.443349 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9126.297678 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9113.868091 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9119.473195 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9136.180979 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9109.025668 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9106.066952 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9106.553866 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9107.832145 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9104.237202 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9101.244536 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9099.538602 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9103.409078 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9109.325336 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9119.525742 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9092.642044 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9131.730154 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9141.444076 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9117.231273 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9105.794832 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9120.185690 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9111.448137 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9095.347747 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9097.124331 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9098.910699 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9100.766044 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9116.543186 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9110.396105 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9114.922001 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9111.816747 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9116.523506 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9110.674661 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9115.1 S.D. = 15.607 S.E.M. = 0.22071 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.09) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.27s 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 9115.072 9016.075 98.997 9214.068 --------------------------------------------------------------------------------------------------- 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.01347 0.04086 -0.33 0.7416 -0.09129 0.06713 209 standlrt 0.55467 0.02375 23.36 1.223e-120 *** 0.50713 0.60061 548 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.09750 0.02008 0.06496 0.14273 2584 cov_Intercept_standlrt 0.01791 0.00818 0.00349 0.03540 2161 var_standlrt 0.02269 0.00581 0.01336 0.03613 1749 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. [95% Cred. Interval] ESS var_Intercept 0.55329 0.01237 0.52920 0.57779 4785 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > # 6.4 Uniform prior . . . . . . . . . . . . . . . . . . . . . . . . . . . 79 > > ## Diffuse priors (Uniform priors) > (mymodel2 <- runMLwiN(normexam ~ 1 + standlrt + (1 + standlrt | school) + (1 | student), estoptions = list(EstM = 1, + mcmcMeth = list(priorcode = 0)), data = tutorial)) MLwiN is running, please wait...... /nogui option ignored ECHO 0 Echoing is ON BATC 1 Batch mode is ON MAXI 2 STAR iteration 0 iteration 1 Convergence not achieved TOLE 2 MAXI 20 BATC 1 Batch mode is ON NEXT iteration 2 iteration 3 Convergence achieved ECHO 0 Echoing is ON MCMC 0 500 1 5.8 50 10 G30[1] G30[2] 1 1 1 1 0 1 BURNING IN... Burning in for 50 iterations out of 500 -2 * Loglike = 9135.903908 Burning in for 100 iterations out of 500 -2 * Loglike = 9104.652611 Burning in for 150 iterations out of 500 -2 * Loglike = 9125.882163 Burning in for 200 iterations out of 500 -2 * Loglike = 9123.153698 Burning in for 250 iterations out of 500 -2 * Loglike = 9104.706169 Burning in for 300 iterations out of 500 -2 * Loglike = 9122.541224 Burning in for 350 iterations out of 500 -2 * Loglike = 9097.078493 Burning in for 400 iterations out of 500 -2 * Loglike = 9103.644340 Burning in for 450 iterations out of 500 -2 * Loglike = 9114.741152 Burning in for 500 iterations out of 500 -2 * Loglike = 9139.496241 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 = 9115.081465 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9125.307254 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9116.740999 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9134.760861 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9121.194725 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9108.974087 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9105.654024 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9131.687439 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9108.355460 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9132.555414 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9115.943999 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9152.513070 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9141.354966 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9133.258306 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9116.810713 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9108.564120 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9121.420451 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9138.742381 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9114.736940 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9091.555016 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9131.363142 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9125.332568 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9097.844330 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9138.484394 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9082.156207 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9101.290608 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9124.239409 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9114.619063 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9114.747777 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9116.337528 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9144.766528 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9122.181723 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9124.973066 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9111.972089 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9151.472475 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9131.683764 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9138.257016 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9133.554004 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9123.409766 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9128.297585 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9104.296050 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9137.095718 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9143.896737 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9102.070659 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9117.396992 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9138.499239 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9127.133528 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9113.067866 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9103.903447 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9149.144987 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9138.087506 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9114.282994 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9126.694810 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9147.028098 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9161.879847 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9121.180491 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9142.555087 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9122.311287 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9127.227149 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9124.217925 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9125.753068 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9117.732995 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9102.999132 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9142.671662 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9122.920747 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9105.394085 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9109.937096 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9133.266238 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9128.653530 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9136.158971 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9105.982545 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9120.144608 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9114.575031 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9133.197650 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9130.344926 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9092.512702 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9114.278857 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9144.742337 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9120.743215 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9142.298781 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9118.469790 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9109.795673 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9088.073600 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9161.811465 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9137.419303 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9115.843319 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9115.648891 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9150.094864 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9100.919860 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9096.297411 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9131.796980 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9085.465876 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9140.397553 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9123.370402 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9100.746109 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9121.987877 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9135.811667 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9101.529840 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9089.833247 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9153.221559 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9120.0 S.D. = 16.246 S.E.M. = 0.22975 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.09) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.19s 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 9120.047 9025.797 94.250 9214.297 --------------------------------------------------------------------------------------------------- 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.01421 0.04113 -0.35 0.7297 -0.09113 0.06585 206 standlrt 0.55491 0.02162 25.67 2.801e-145 *** 0.51271 0.59711 641 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.10334 0.02219 0.06820 0.15409 3208 cov_Intercept_standlrt 0.02055 0.00834 0.00598 0.03916 1818 var_standlrt 0.01771 0.00561 0.00908 0.03075 949 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. [95% Cred. Interval] ESS var_Intercept 0.55424 0.01243 0.53049 0.57877 4456 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > # 6.5 Informative prior . . . . . . . . . . . . . . . . . . . . . . . . . 80 > > ## Informative normal prior for Sigma_u > (mymodel3 <- runMLwiN(normexam ~ 1 + standlrt + (1 + standlrt | school) + (1 | student), estoptions = list(EstM = 1, + mcmcMeth = list(priorParam = list(rp2 = list(estimate = matrix(c(0.09, 0.018, 0.018, 0.015), 2, 2), size = 65)))), + data = tutorial)) MLwiN is running, please wait...... /nogui option ignored ECHO 0 Echoing is ON BATC 1 Batch mode is ON MAXI 2 STAR iteration 0 iteration 1 Convergence not achieved TOLE 2 MAXI 20 BATC 1 Batch mode is ON NEXT iteration 2 iteration 3 Convergence achieved ECHO 0 Echoing is ON MCMC 0 500 1 5.8 50 10 G30[1] G30[2] c1092 1 1 1 1 1 1 BURNING IN... Burning in for 50 iterations out of 500 -2 * Loglike = 9122.363827 Burning in for 100 iterations out of 500 -2 * Loglike = 9130.332318 Burning in for 150 iterations out of 500 -2 * Loglike = 9131.643496 Burning in for 200 iterations out of 500 -2 * Loglike = 9126.868873 Burning in for 250 iterations out of 500 -2 * Loglike = 9104.697630 Burning in for 300 iterations out of 500 -2 * Loglike = 9115.331556 Burning in for 350 iterations out of 500 -2 * Loglike = 9119.458734 Burning in for 400 iterations out of 500 -2 * Loglike = 9100.655210 Burning in for 450 iterations out of 500 -2 * Loglike = 9121.330592 Burning in for 500 iterations out of 500 -2 * Loglike = 9107.375644 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 = 9128.993505 Actual update 100 of 5000, Stored update 100 of 5000 -2 * Loglike = 9119.505287 Actual update 150 of 5000, Stored update 150 of 5000 -2 * Loglike = 9148.452617 Actual update 200 of 5000, Stored update 200 of 5000 -2 * Loglike = 9132.150660 Actual update 250 of 5000, Stored update 250 of 5000 -2 * Loglike = 9154.761068 Actual update 300 of 5000, Stored update 300 of 5000 -2 * Loglike = 9103.272896 Actual update 350 of 5000, Stored update 350 of 5000 -2 * Loglike = 9118.421083 Actual update 400 of 5000, Stored update 400 of 5000 -2 * Loglike = 9110.961018 Actual update 450 of 5000, Stored update 450 of 5000 -2 * Loglike = 9143.142419 Actual update 500 of 5000, Stored update 500 of 5000 -2 * Loglike = 9142.754873 Actual update 550 of 5000, Stored update 550 of 5000 -2 * Loglike = 9118.419338 Actual update 600 of 5000, Stored update 600 of 5000 -2 * Loglike = 9136.143754 Actual update 650 of 5000, Stored update 650 of 5000 -2 * Loglike = 9144.192037 Actual update 700 of 5000, Stored update 700 of 5000 -2 * Loglike = 9137.121837 Actual update 750 of 5000, Stored update 750 of 5000 -2 * Loglike = 9120.656566 Actual update 800 of 5000, Stored update 800 of 5000 -2 * Loglike = 9138.366900 Actual update 850 of 5000, Stored update 850 of 5000 -2 * Loglike = 9126.876356 Actual update 900 of 5000, Stored update 900 of 5000 -2 * Loglike = 9117.054350 Actual update 950 of 5000, Stored update 950 of 5000 -2 * Loglike = 9124.748664 Actual update 1000 of 5000, Stored update 1000 of 5000 -2 * Loglike = 9139.747982 Actual update 1050 of 5000, Stored update 1050 of 5000 -2 * Loglike = 9126.389336 Actual update 1100 of 5000, Stored update 1100 of 5000 -2 * Loglike = 9107.797044 Actual update 1150 of 5000, Stored update 1150 of 5000 -2 * Loglike = 9095.840081 Actual update 1200 of 5000, Stored update 1200 of 5000 -2 * Loglike = 9141.822991 Actual update 1250 of 5000, Stored update 1250 of 5000 -2 * Loglike = 9121.696205 Actual update 1300 of 5000, Stored update 1300 of 5000 -2 * Loglike = 9147.667123 Actual update 1350 of 5000, Stored update 1350 of 5000 -2 * Loglike = 9108.237056 Actual update 1400 of 5000, Stored update 1400 of 5000 -2 * Loglike = 9113.271547 Actual update 1450 of 5000, Stored update 1450 of 5000 -2 * Loglike = 9117.493572 Actual update 1500 of 5000, Stored update 1500 of 5000 -2 * Loglike = 9122.095028 Actual update 1550 of 5000, Stored update 1550 of 5000 -2 * Loglike = 9107.090794 Actual update 1600 of 5000, Stored update 1600 of 5000 -2 * Loglike = 9122.623136 Actual update 1650 of 5000, Stored update 1650 of 5000 -2 * Loglike = 9114.323257 Actual update 1700 of 5000, Stored update 1700 of 5000 -2 * Loglike = 9121.046034 Actual update 1750 of 5000, Stored update 1750 of 5000 -2 * Loglike = 9114.603646 Actual update 1800 of 5000, Stored update 1800 of 5000 -2 * Loglike = 9087.809069 Actual update 1850 of 5000, Stored update 1850 of 5000 -2 * Loglike = 9121.917764 Actual update 1900 of 5000, Stored update 1900 of 5000 -2 * Loglike = 9110.638821 Actual update 1950 of 5000, Stored update 1950 of 5000 -2 * Loglike = 9116.615444 Actual update 2000 of 5000, Stored update 2000 of 5000 -2 * Loglike = 9118.590341 Actual update 2050 of 5000, Stored update 2050 of 5000 -2 * Loglike = 9121.614019 Actual update 2100 of 5000, Stored update 2100 of 5000 -2 * Loglike = 9146.183013 Actual update 2150 of 5000, Stored update 2150 of 5000 -2 * Loglike = 9110.755170 Actual update 2200 of 5000, Stored update 2200 of 5000 -2 * Loglike = 9119.857268 Actual update 2250 of 5000, Stored update 2250 of 5000 -2 * Loglike = 9099.723205 Actual update 2300 of 5000, Stored update 2300 of 5000 -2 * Loglike = 9105.841641 Actual update 2350 of 5000, Stored update 2350 of 5000 -2 * Loglike = 9118.414326 Actual update 2400 of 5000, Stored update 2400 of 5000 -2 * Loglike = 9098.804582 Actual update 2450 of 5000, Stored update 2450 of 5000 -2 * Loglike = 9111.643500 Actual update 2500 of 5000, Stored update 2500 of 5000 -2 * Loglike = 9117.582528 Actual update 2550 of 5000, Stored update 2550 of 5000 -2 * Loglike = 9123.440345 Actual update 2600 of 5000, Stored update 2600 of 5000 -2 * Loglike = 9096.578302 Actual update 2650 of 5000, Stored update 2650 of 5000 -2 * Loglike = 9094.830638 Actual update 2700 of 5000, Stored update 2700 of 5000 -2 * Loglike = 9116.414242 Actual update 2750 of 5000, Stored update 2750 of 5000 -2 * Loglike = 9134.420171 Actual update 2800 of 5000, Stored update 2800 of 5000 -2 * Loglike = 9097.513740 Actual update 2850 of 5000, Stored update 2850 of 5000 -2 * Loglike = 9128.829932 Actual update 2900 of 5000, Stored update 2900 of 5000 -2 * Loglike = 9117.332961 Actual update 2950 of 5000, Stored update 2950 of 5000 -2 * Loglike = 9125.885188 Actual update 3000 of 5000, Stored update 3000 of 5000 -2 * Loglike = 9140.476065 Actual update 3050 of 5000, Stored update 3050 of 5000 -2 * Loglike = 9120.044965 Actual update 3100 of 5000, Stored update 3100 of 5000 -2 * Loglike = 9128.496487 Actual update 3150 of 5000, Stored update 3150 of 5000 -2 * Loglike = 9092.561518 Actual update 3200 of 5000, Stored update 3200 of 5000 -2 * Loglike = 9111.864079 Actual update 3250 of 5000, Stored update 3250 of 5000 -2 * Loglike = 9091.474662 Actual update 3300 of 5000, Stored update 3300 of 5000 -2 * Loglike = 9141.544971 Actual update 3350 of 5000, Stored update 3350 of 5000 -2 * Loglike = 9127.760418 Actual update 3400 of 5000, Stored update 3400 of 5000 -2 * Loglike = 9128.264474 Actual update 3450 of 5000, Stored update 3450 of 5000 -2 * Loglike = 9100.419185 Actual update 3500 of 5000, Stored update 3500 of 5000 -2 * Loglike = 9105.780915 Actual update 3550 of 5000, Stored update 3550 of 5000 -2 * Loglike = 9153.061392 Actual update 3600 of 5000, Stored update 3600 of 5000 -2 * Loglike = 9120.717553 Actual update 3650 of 5000, Stored update 3650 of 5000 -2 * Loglike = 9149.681723 Actual update 3700 of 5000, Stored update 3700 of 5000 -2 * Loglike = 9147.668659 Actual update 3750 of 5000, Stored update 3750 of 5000 -2 * Loglike = 9135.198110 Actual update 3800 of 5000, Stored update 3800 of 5000 -2 * Loglike = 9133.315718 Actual update 3850 of 5000, Stored update 3850 of 5000 -2 * Loglike = 9132.389113 Actual update 3900 of 5000, Stored update 3900 of 5000 -2 * Loglike = 9143.346829 Actual update 3950 of 5000, Stored update 3950 of 5000 -2 * Loglike = 9111.110959 Actual update 4000 of 5000, Stored update 4000 of 5000 -2 * Loglike = 9169.782815 Actual update 4050 of 5000, Stored update 4050 of 5000 -2 * Loglike = 9108.794660 Actual update 4100 of 5000, Stored update 4100 of 5000 -2 * Loglike = 9141.687174 Actual update 4150 of 5000, Stored update 4150 of 5000 -2 * Loglike = 9105.278517 Actual update 4200 of 5000, Stored update 4200 of 5000 -2 * Loglike = 9119.949246 Actual update 4250 of 5000, Stored update 4250 of 5000 -2 * Loglike = 9125.489570 Actual update 4300 of 5000, Stored update 4300 of 5000 -2 * Loglike = 9117.859523 Actual update 4350 of 5000, Stored update 4350 of 5000 -2 * Loglike = 9110.987858 Actual update 4400 of 5000, Stored update 4400 of 5000 -2 * Loglike = 9116.800327 Actual update 4450 of 5000, Stored update 4450 of 5000 -2 * Loglike = 9133.244208 Actual update 4500 of 5000, Stored update 4500 of 5000 -2 * Loglike = 9131.761474 Actual update 4550 of 5000, Stored update 4550 of 5000 -2 * Loglike = 9104.211476 Actual update 4600 of 5000, Stored update 4600 of 5000 -2 * Loglike = 9110.829910 Actual update 4650 of 5000, Stored update 4650 of 5000 -2 * Loglike = 9116.094565 Actual update 4700 of 5000, Stored update 4700 of 5000 -2 * Loglike = 9117.727257 Actual update 4750 of 5000, Stored update 4750 of 5000 -2 * Loglike = 9133.727343 Actual update 4800 of 5000, Stored update 4800 of 5000 -2 * Loglike = 9130.106899 Actual update 4850 of 5000, Stored update 4850 of 5000 -2 * Loglike = 9111.917504 Actual update 4900 of 5000, Stored update 4900 of 5000 -2 * Loglike = 9118.645420 Actual update 4950 of 5000, Stored update 4950 of 5000 -2 * Loglike = 9124.067795 Actual update 5000 of 5000, Stored update 5000 of 5000 -2 * Loglike = 9120.404518 PUPN c1003 c1004 AVER c1091 b99 b100 Count = 5000 Average = 9121.3 S.D. = 15.747 S.E.M. = 0.22270 ECHO 0 Execution completed -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- MLwiN (version: 3.09) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.09s 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 9121.322 9030.486 90.837 9212.159 --------------------------------------------------------------------------------------------------- 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.01132 0.03953 -0.29 0.7746 -0.08603 0.06779 208 standlrt 0.55655 0.02053 27.11 6.997e-162 *** 0.51621 0.59658 600 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.09060 0.01195 0.06979 0.11647 3955 cov_Intercept_standlrt 0.01811 0.00413 0.01075 0.02696 3116 var_standlrt 0.01487 0.00233 0.01103 0.01996 2462 --------------------------------------------------------------------------------------------------- The random part estimates at the student level: Coef. Std. Err. [95% Cred. Interval] ESS var_Intercept 0.55410 0.01233 0.53063 0.57856 4559 -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- > > # 6.6 Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81 > > if (!require(texreg)) { + warning("texreg package required to use screenreg() function") + } else { + screenreg(list(mymodel0a, mymodel0, mymodel1, mymodel2, mymodel3), custom.model.names=c("IGLS", "default", "prior 2", "uniform", "prior 4"), groups = list("Fixed Part" = 1:2, "Level-2" = 3:5, "Level-1" = 6:6), + stars = numeric(0), include.nobs=FALSE, include.loglik=FALSE, include.deviance=FALSE, include.dbar=FALSE, include.dthetabar=FALSE, include.pd=FALSE, include.dic=FALSE) + } Loading required package: texreg Version: 1.39.3 Date: 2023-11-09 Author: Philip Leifeld (University of Essex) Consider submitting praise using the praise or praise_interactive functions. Please cite the JSS article in your publications -- see citation("texreg"). ================================================================================================== IGLS default prior 2 uniform prior 4 -------------------------------------------------------------------------------------------------- Fixed Part FP_Intercept -0.01 -0.01 -0.01 -0.01 -0.01 (0.04) [-0.09; 0.07] [-0.09; 0.07] [-0.09; 0.07] [-0.09; 0.07] FP_standlrt 0.56 0.56 * 0.55 * 0.55 * 0.56 * (0.02) [ 0.52; 0.60] [ 0.51; 0.60] [ 0.51; 0.60] [ 0.52; 0.60] Level-2 RP2_var_Intercept 0.09 0.10 * 0.10 * 0.10 * 0.09 * (0.02) [ 0.06; 0.14] [ 0.06; 0.14] [ 0.07; 0.15] [ 0.07; 0.12] RP2_cov_Intercept_standlrt 0.02 0.02 * 0.02 * 0.02 * 0.02 * (0.01) [ 0.01; 0.04] [ 0.00; 0.04] [ 0.01; 0.04] [ 0.01; 0.03] RP2_var_standlrt 0.01 0.02 * 0.02 * 0.02 * 0.01 * (0.00) [ 0.01; 0.03] [ 0.01; 0.04] [ 0.01; 0.03] [ 0.01; 0.02] Level-1 RP1_var_Intercept 0.55 0.55 * 0.55 * 0.55 * 0.55 * (0.01) [ 0.53; 0.58] [ 0.53; 0.58] [ 0.53; 0.58] [ 0.53; 0.58] ================================================================================================== * Null hypothesis value outside the confidence interval. > > # Chapter learning outcomes . . . . . . . . . . . . . . . . . . . . . . . 81 > > > > > > ############################################################################ > > proc.time() user system elapsed 16.68 2.45 46.51