R version 4.4.1 (2024-06-14 ucrt) -- "Race for Your Life" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ############################################################################ > # MLwiN MCMC Manual > # > # 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.11/ 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.11) multilevel model (Normal) Estimation algorithm: MCMC Elapsed time : 13.58s 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.11) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: IGLS Elapsed time : 0.08s 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.11) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.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 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.11) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.25s 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.11) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.15s 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.11) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.21s 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.11) multilevel model (Normal) N min mean max N_complete min_complete mean_complete max_complete school 65 2 62.44615 198 65 2 62.44615 198 Estimation algorithm: MCMC Elapsed time : 2.16s 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 15.53 2.51 47.04