|
|
|
||||
|
|
|||||
use http://www.ats.ucla.edu/stat/examples/imm/imm23, clear
describe
Contains data from imm23.dta
obs: 519
vars: 18 19 Mar 2002 10:39
size: 39,444 (99.8% of memory free)
-------------------------------------------------------------------------------
storage display value
variable name type format label variable label
-------------------------------------------------------------------------------
schid float %9.0g School ID
stuid float %9.0g Student ID
ses float %9.0g Socioecnonomic Status
meanses float %9.0g Mean SES for the school
homework float %9.0g Time spent on math homework
each week
white float %9.0g Race: 1=white, 0=non-white
parented float %9.0g Parents highest education level
public float %9.0g Public school: 1=public,
0=non-public
ratio float %9.0g Student-Teacher ratio
percmin float %9.0g Percent minority in school
math float %9.0g Math score
sex float %9.0g Sex: 1=male, 2=female
race float %9.0g race of student, 1=asian,
2=Hispanic, 3=Black, 4=White,
5=Native American
sctype float %9.0g Type of school, 1=public,
2=catholic, 3=Private other
religious, 4=Private non-r
cstr float %9.0g
scsize float %9.0g
urban float %9.0g
region float %9.0g
-------------------------------------------------------------------------------
Sorted by: schid stuid
Page 64, 4.2.2, The Null Model, Model 0.
xtmixed math || schid:, variance
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1899.3393
Iteration 1: log restricted-likelihood = -1899.3393
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(0) = .
Log restricted-likelihood = -1899.3393 Prob > chi2 = .
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_cons | 50.75885 1.151181 44.09 0.000 48.50258 53.01513
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Identity |
var(_cons) | 26.12419 8.987209 13.31109 51.27103
-----------------------------+------------------------------------------------
var(Residual) | 81.24387 5.153894 71.74516 92.00015
------------------------------------------------------------------------------
LR test vs. linear regression: chibar2(01) = 134.06 Prob >= chibar2 = 0.0000
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1899.339 3 3804.679 3817.434
------------------------------------------------------------------------------
Page 65, 4.2.3 'Homework and 'MathAchievement', Model 1.
xtmixed math homework || schid:, variance
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1864.6598
Iteration 1: log restricted-likelihood = -1864.6598
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(1) = 74.95
Log restricted-likelihood = -1864.6598 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 2.399867 .2771976 8.66 0.000 1.856569 2.943164
_cons | 46.35575 1.162776 39.87 0.000 44.07676 48.63475
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Identity |
var(_cons) | 21.34157 7.574593 10.64428 42.78943
-----------------------------+------------------------------------------------
var(Residual) | 71.28374 4.530462 62.93498 80.74002
------------------------------------------------------------------------------
LR test vs. linear regression: chibar2(01) = 96.43 Prob >= chibar2 = 0.0000
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1864.66 4 3737.32 3754.327
------------------------------------------------------------------------------
Pages 66 and 67, 4.2.4 Random slope for 'Homework', model 2.
xtmixed math homework || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1817.781
Iteration 1: log restricted-likelihood = -1817.781
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(1) = 4.55
Log restricted-likelihood = -1817.781 Prob > chi2 = 0.0329
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.980164 .9284474 2.13 0.033 .1604405 3.799888
_cons | 46.32561 1.758932 26.34 0.000 42.87816 49.77305
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 17.72647 6.260262 8.871821 35.41863
var(_cons) | 62.42438 21.38146 31.90087 122.1535
cov(homework,_cons) | -27.59383 10.56622 -48.30324 -6.884415
-----------------------------+------------------------------------------------
var(Residual) | 53.29463 3.465963 46.91659 60.53973
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 190.19 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1817.781 6 3647.562 3673.074
------------------------------------------------------------------------------
Page 69, 4.2.5 Adding 'ParentEducation', model 3.
xtmixed math homework parented || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1800.0217
Iteration 1: log restricted-likelihood = -1800.0217
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(2) = 44.61
Log restricted-likelihood = -1800.0217 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.881702 .830247 2.27 0.023 .2544478 3.508956
parented | 1.841414 .2959197 6.22 0.000 1.261422 2.421406
_cons | 40.85465 1.790225 22.82 0.000 37.34588 44.36343
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 13.8786 5.039229 6.812086 28.27557
var(_cons) | 47.86951 17.01638 23.84947 96.08136
cov(homework,_cons) | -21.95137 8.502263 -38.6155 -5.287237
-----------------------------+------------------------------------------------
var(Residual) | 50.77825 3.311781 44.68502 57.70234
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 119.94 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1800.022 7 3614.043 3643.807
------------------------------------------------------------------------------
Page 70 and 71, 4.2.6 Traditional regression model.
xtmixed math homework parented , variance
Mixed-effects REML regression Number of obs = 519
Wald chi2(2) = 264.49
Log restricted-likelihood = -1859.9932 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 2.335361 .2683967 8.70 0.000 1.809313 2.861409
parented | 3.003987 .2765059 10.86 0.000 2.462046 3.545929
_cons | 37.23915 .9962789 37.38 0.000 35.28648 39.19182
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
var(Residual) | 76.11628 4.725077 67.39652 85.96421
------------------------------------------------------------------------------
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1859.993 4 3727.986 3744.994
------------------------------------------------------------------------------
/* using the regress command */
regress math homework parented
Source | SS df MS Number of obs = 519
-------------+------------------------------ F( 2, 516) = 132.25
Model | 20132.0446 2 10066.0223 Prob > F = 0.0000
Residual | 39276.0016 516 76.1162822 R-squared = 0.3389
-------------+------------------------------ Adj R-squared = 0.3363
Total | 59408.0462 518 114.687348 Root MSE = 8.7245
------------------------------------------------------------------------------
math | Coef. Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 2.335361 .2683967 8.70 0.000 1.808077 2.862646
parented | 3.003987 .2765059 10.86 0.000 2.460771 3.547203
_cons | 37.23915 .9962789 37.38 0.000 35.28189 39.19641
------------------------------------------------------------------------------
Page 73/74, 4.3.2 A model with 'SchoolSize' (Model 2).
xtmixed math homework scsize || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1817.1153
Iteration 1: log restricted-likelihood = -1817.1153
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(2) = 5.02
Log restricted-likelihood = -1817.1153 Prob > chi2 = 0.0812
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.981179 .9298449 2.13 0.033 .1587166 3.803642
scsize | .4259322 .6409595 0.66 0.506 -.8303253 1.68219
_cons | 44.97198 2.710876 16.59 0.000 39.65876 50.2852
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 17.78479 6.281797 8.900087 35.53883
var(_cons) | 66.42771 23.33067 33.37269 132.2231
cov(homework,_cons) | -28.75128 11.07166 -50.45133 -7.051229
-----------------------------+------------------------------------------------
var(Residual) | 53.3058 3.467273 46.92541 60.55372
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 189.33 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1817.115 7 3648.231 3677.994
------------------------------------------------------------------------------
Page 74/75, 4.3.3 Changing 'SchoolSize' to 'Public' (Model 3).
xtmixed math homework public || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1814.194
Iteration 1: log restricted-likelihood = -1814.194
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(2) = 8.92
Log restricted-likelihood = -1814.194 Prob > chi2 = 0.0116
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.97527 .9184314 2.15 0.031 .1751778 3.775363
public | -4.060999 1.977904 -2.05 0.040 -7.93762 -.1843776
_cons | 49.05035 2.184881 22.45 0.000 44.76807 53.33264
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 17.31097 6.113316 8.664078 34.58762
var(_cons) | 60.17005 20.89702 30.46202 118.8508
cov(homework,_cons) | -27.45916 10.38128 -47.8061 -7.112219
-----------------------------+------------------------------------------------
var(Residual) | 53.34333 3.471559 46.95527 60.60046
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 157.64 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1814.194 7 3642.388 3672.151
------------------------------------------------------------------------------
Page 77, 4.3.4 Adding a cross level interaction with 'Public', (Model 4).
generate homewrkXpublic = homework*public
xtmixed math homework public homewrkXpublic || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1812.5753
Iteration 1: log restricted-likelihood = -1812.5753
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(3) = 8.72
Log restricted-likelihood = -1812.5753 Prob > chi2 = 0.0333
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 2.292791 1.591376 1.44 0.150 -.8262482 5.411831
public | -3.261944 3.714642 -0.88 0.380 -10.54251 4.018621
homewrkXpu~c | -.4956655 1.972783 -0.25 0.802 -4.362249 3.370918
_cons | 48.52887 3.016064 16.09 0.000 42.6175 54.44025
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 18.25189 6.551994 9.031228 36.88663
var(_cons) | 62.55576 22.05624 31.3434 124.85
cov(homework,_cons) | -28.95327 11.08199 -50.67357 -7.232962
-----------------------------+------------------------------------------------
var(Residual) | 53.33265 3.470174 46.94705 60.58679
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 159.63 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1812.575 8 3641.151 3675.166
------------------------------------------------------------------------------
Page 80, 4.3.5 Model 4 will full NELS-88 data (we don't have these data, so this is omitted).
Page 80/82, 4.3.6 Deleting 'HomePublic' and adding 'White' (Model 5).
xtmixed math homework public white || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1807.7266
Iteration 1: log restricted-likelihood = -1807.7266
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(3) = 21.73
Log restricted-likelihood = -1807.7266 Prob > chi2 = 0.0001
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.899176 .9050575 2.10 0.036 .125296 3.673056
public | -3.882385 1.798686 -2.16 0.031 -7.407744 -.3570262
white | 3.309409 .9698517 3.41 0.001 1.408534 5.210283
_cons | 46.62829 2.192207 21.27 0.000 42.33164 50.92493
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 16.77652 5.931869 8.389431 33.54834
var(_cons) | 56.03024 19.63033 28.19694 111.3379
cov(homework,_cons) | -26.81495 10.00251 -46.41951 -7.210387
-----------------------------+------------------------------------------------
var(Residual) | 52.72817 3.436984 46.40436 59.91377
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 133.03 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1807.727 8 3631.453 3665.468
------------------------------------------------------------------------------
Page 82/83, 4.3.7 Adding a random part for 'White' (Model 6).
xtmixed math homework public white || schid: homework white, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1805.255
Iteration 1: log restricted-likelihood = -1805.2505
Iteration 2: log restricted-likelihood = -1805.2505
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(3) = 15.98
Log restricted-likelihood = -1805.2505 Prob > chi2 = 0.0011
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.940045 .9004501 2.15 0.031 .1751951 3.704894
public | -4.929207 1.650688 -2.99 0.003 -8.164497 -1.693917
white | 2.598035 1.553804 1.67 0.095 -.4473641 5.643434
_cons | 48.22768 2.340021 20.61 0.000 43.64132 52.81403
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 16.63131 5.874699 8.322529 33.23514
var(white) | 26.50044 22.54574 5.001239 140.4199
var(_cons) | 69.16463 31.24256 28.53549 167.642
cov(homework,white) | 2.923541 7.944215 -12.64683 18.49392
cov(homework,_cons) | -28.54352 12.26504 -52.58255 -4.504491
cov(white,_cons) | -22.08364 20.96947 -63.18305 19.01577
-----------------------------+------------------------------------------------
var(Residual) | 51.15698 3.384706 44.9352 58.24023
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(6) = 137.98 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1805.251 11 3632.501 3679.272
------------------------------------------------------------------------------
Page 85, 4.3.8 Making the coefficient of 'White' fixed and adding 'MeanSES' (Model 7).
xtmixed math homework public white meanses || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1803.2599
Iteration 1: log restricted-likelihood = -1803.2599
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(4) = 31.50
Log restricted-likelihood = -1803.2599 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.924179 .8964219 2.15 0.032 .1672248 3.681134
public | .1604258 2.270945 0.07 0.944 -4.290545 4.611396
white | 3.096913 .9630887 3.22 0.001 1.209294 4.984532
meanses | 4.977274 1.950878 2.55 0.011 1.153624 8.800924
_cons | 44.61584 2.235094 19.96 0.000 40.23513 48.99654
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 16.44364 5.821773 8.21552 32.91249
var(_cons) | 54.13689 19.09921 27.11402 108.0918
cov(homework,_cons) | -27.0358 9.902771 -46.44487 -7.626721
-----------------------------+------------------------------------------------
var(Residual) | 52.79435 3.444721 46.45669 59.99659
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 109.59 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1803.26 9 3624.52 3662.787
------------------------------------------------------------------------------
Page 86, 4.3.9 Deleting the school characteristic 'Public' (Model 8).
xtmixed math homework white meanses || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1804.9907
Iteration 1: log restricted-likelihood = -1804.9907
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(3) = 32.52
Log restricted-likelihood = -1804.9907 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.925085 .8951596 2.15 0.032 .1706048 3.679566
white | 3.114922 .9570099 3.25 0.001 1.239217 4.990627
meanses | 4.892483 1.340592 3.65 0.000 2.264971 7.519995
_cons | 44.7022 1.787226 25.01 0.000 41.1993 48.2051
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 16.39899 5.804124 8.195013 32.8159
var(_cons) | 53.57389 18.91041 26.8224 107.0061
cov(homework,_cons) | -26.99907 9.864397 -46.33293 -7.665203
-----------------------------+------------------------------------------------
var(Residual) | 52.79042 3.444296 46.45353 59.99175
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 108.55 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1804.991 8 3625.981 3659.997
------------------------------------------------------------------------------
Page 87/88, 4.3.10 Adding an interaction between 'HomeWork' and 'MeanSES' (Model 9).
generate homewrkXmeanses = homework*meanses
xtmixed math homework white meanses homewrkXmeanses || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1803.5271
Iteration 1: log restricted-likelihood = -1803.5271
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(4) = 32.31
Log restricted-likelihood = -1803.5271 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.974678 .9291318 2.13 0.034 .1536126 3.795742
white | 3.115026 .9571186 3.25 0.001 1.239108 4.990944
meanses | 3.977613 3.022684 1.32 0.188 -1.946739 9.901965
homewrkXme~s | .5530935 1.64573 0.34 0.737 -2.672479 3.778666
_cons | 44.61137 1.836096 24.30 0.000 41.01268 48.21005
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 17.23849 6.210693 8.508005 34.92775
var(_cons) | 55.88547 20.01922 27.6938 112.7756
cov(homework,_cons) | -28.38963 10.52383 -49.01596 -7.763291
-----------------------------+------------------------------------------------
var(Residual) | 52.78088 3.443055 46.4462 59.97952
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 111.28 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1803.527 9 3625.054 3663.321
------------------------------------------------------------------------------
Page 88/89, 4.3.11 Adding another student-level variable (Model 10).
xtmixed math homework white meanses ses || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1796.3777
Iteration 1: log restricted-likelihood = -1796.3777
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(4) = 50.25
Log restricted-likelihood = -1796.3777 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.825613 .8504767 2.15 0.032 .1587094 3.492517
white | 2.170787 .9732755 2.23 0.026 .2632017 4.078372
meanses | 2.945142 1.428078 2.06 0.039 .1461609 5.744123
ses | 2.205733 .5356316 4.12 0.000 1.155914 3.255551
_cons | 45.67481 1.752034 26.07 0.000 42.24089 49.10874
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 14.64975 5.247581 7.259843 29.56194
var(_cons) | 49.89752 17.73365 24.86341 100.1376
cov(homework,_cons) | -24.37951 9.046941 -42.11119 -6.647832
-----------------------------+------------------------------------------------
var(Residual) | 51.29601 3.350503 45.13211 58.30175
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 103.65 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1796.378 9 3610.755 3649.023
------------------------------------------------------------------------------
Page 88/89, 4.3.12: Analyses with NELS-88 (we don't have these data, so these analyses are omitted).
Page 91, 4.4.1 'SES' as a student-level explanatory variable (Model 1).
xtmixed math ses || schid:, variance
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1873.0887
Iteration 1: log restricted-likelihood = -1873.0887
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(1) = 58.53
Log restricted-likelihood = -1873.0887 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
ses | 4.332311 .5662715 7.65 0.000 3.222439 5.442182
_cons | 51.20085 .8507325 60.18 0.000 49.53345 52.86826
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Identity |
var(_cons) | 12.6317 5.015357 5.800912 27.506
-----------------------------+------------------------------------------------
var(Residual) | 75.32781 4.78723 66.50585 85.32001
------------------------------------------------------------------------------
LR test vs. linear regression: chibar2(01) = 43.58 Prob >= chibar2 = 0.0000
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1873.089 4 3754.177 3771.185
------------------------------------------------------------------------------
Page 92, 4.4.2 Adding a random slope (Model 2). Model never converged.
set maxiter 100
xtmixed math ses || schid: ses, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1873.3707
Iteration 1: log restricted-likelihood = -1873.0375
[output omitted]
Iteration 99: log restricted-likelihood = -1872.9972 (not concave)
Iteration 100: log restricted-likelihood = -1872.9972 (not concave)
convergence not achieved
Computing standard errors:
standard error calculation failed
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(1) = 58.02
Log restricted-likelihood = -1872.9972 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
ses | 4.32386 .5676435 7.62 0.000 3.211299 5.436421
_cons | 51.24839 .8562642 59.85 0.000 49.57015 52.92664
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(ses) | .0962659 . . .
var(_cons) | 12.89715 . . .
cov(ses,_cons) | -1.114251 . . .
-----------------------------+------------------------------------------------
var(Residual) | 75.19665 . . .
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 43.77 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1872.997 2 3749.994 3758.498
------------------------------------------------------------------------------
Page 93, 4.4.3 Adding 'PercentMinorities' (Model 3), SAS Program.
xtmixed math ses percmin || schid:, variance
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1870.8054
Iteration 1: log restricted-likelihood = -1870.8054
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(2) = 66.23
Log restricted-likelihood = -1870.8054 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
ses | 4.298805 .5618229 7.65 0.000 3.197652 5.399958
percmin | -.8093957 .3647249 -2.22 0.026 -1.524243 -.094548
_cons | 53.12666 1.178066 45.10 0.000 50.81769 55.43562
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Identity |
var(_cons) | 10.69835 4.337275 4.833108 23.68137
-----------------------------+------------------------------------------------
var(Residual) | 75.17249 4.769303 66.38266 85.1262
------------------------------------------------------------------------------
LR test vs. linear regression: chibar2(01) = 40.65 Prob >= chibar2 = 0.0000
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1870.805 5 3751.611 3772.87
------------------------------------------------------------------------------
Page 95, 4.4.4 Adding 'MeanSES' (Model 4).
xtmixed math ses percmin meanses || schid:, variance
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1867.8075
Iteration 1: log restricted-likelihood = -1867.8075
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(3) = 73.80
Log restricted-likelihood = -1867.8075 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
ses | 3.88476 .6097731 6.37 0.000 2.689627 5.079894
percmin | -.6921886 .3456913 -2.00 0.045 -1.369731 -.014646
meanses | 2.805075 1.479369 1.90 0.058 -.0944343 5.704584
_cons | 53.0969 1.102666 48.15 0.000 50.93572 55.25809
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Identity |
var(_cons) | 8.846737 3.936787 3.698277 21.16249
-----------------------------+------------------------------------------------
var(Residual) | 75.24051 4.777618 66.4358 85.2121
------------------------------------------------------------------------------
LR test vs. linear regression: chibar2(01) = 25.48 Prob >= chibar2 = 0.0000
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1867.808 6 3747.615 3773.126
------------------------------------------------------------------------------
Page 95, 4.4.5 Analyses with NELS-88, models 2 and 3 (we do not have these data, so these analyses are omitted).
Page 99, 4.5.1 Analysis with class size and a cross level interaction (Model 1).
xtmixed math homework ratio || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1818.3188
Iteration 1: log restricted-likelihood = -1818.3188
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(2) = 4.74
Log restricted-likelihood = -1818.3188 Prob > chi2 = 0.0934
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 1.979276 .9288626 2.13 0.033 .1587383 3.799813
ratio | -.0944215 .2123894 -0.44 0.657 -.510697 .3218541
_cons | 47.9616 4.083297 11.75 0.000 39.95848 55.96471
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 17.74195 6.268408 8.876944 35.46005
var(_cons) | 63.48001 21.7818 32.40153 124.368
cov(homework,_cons) | -27.70352 10.64132 -48.56012 -6.846919
-----------------------------+------------------------------------------------
var(Residual) | 53.30429 3.467218 46.92401 60.5521
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 183.07 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1818.319 7 3650.638 3680.401
------------------------------------------------------------------------------
Page 100, 4.5.2 Interaction between 'Ratio' and 'HomeWork' (Model 2), SAS Program.
generate homewrkXratio = homework*ratio
xtmixed math homework homewrkXratio || schid: homework, variance covar(un)
Performing EM optimization:
Performing gradient-based optimization:
Iteration 0: log restricted-likelihood = -1818.9455
Iteration 1: log restricted-likelihood = -1818.9455
Computing standard errors:
Mixed-effects REML regression Number of obs = 519
Group variable: schid Number of groups = 23
Obs per group: min = 5
avg = 22.6
max = 67
Wald chi2(2) = 4.66
Log restricted-likelihood = -1818.9455 Prob > chi2 = 0.0974
------------------------------------------------------------------------------
math | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
homework | 2.884085 2.156817 1.34 0.181 -1.3432 7.111369
homewrkXra~o | -.0525855 .1122818 -0.47 0.640 -.2726539 .1674828
_cons | 46.33064 1.760591 26.32 0.000 42.87995 49.78134
------------------------------------------------------------------------------
------------------------------------------------------------------------------
Random-effects Parameters | Estimate Std. Err. [95% Conf. Interval]
-----------------------------+------------------------------------------------
schid: Unstructured |
var(homework) | 18.09523 6.390325 9.056545 36.15477
var(_cons) | 62.55502 21.42087 31.97298 122.3887
cov(homework,_cons) | -27.80043 10.65508 -48.68401 -6.916853
-----------------------------+------------------------------------------------
var(Residual) | 53.29771 3.46642 46.91886 60.5438
------------------------------------------------------------------------------
LR test vs. linear regression: chi2(3) = 180.57 Prob > chi2 = 0.0000
Note: LR test is conservative and provided only for reference
estat ic
------------------------------------------------------------------------------
Model | Obs ll(null) ll(model) df AIC BIC
-------------+----------------------------------------------------------------
| 519 . -1818.945 7 3651.891 3681.654
------------------------------------------------------------------------------
Page 100, 4.5.3 Reporting the modeling session with NELS-88 (we do not have these data, so these analyses are omitted).
UCLA Researchers are invited to our Statistical Consulting Services
We recommend others to our list of Other Resources for Statistical Computing Help
These pages are Copyrighted (c) by UCLA Academic Technology Services