Table XTMIXED Procedure in STATA with Output Systolic Blood Pressure, 2001. use "k:mydirectory,. xtmixed sbp nage20 nage30 nage40 nage50 nage70 nage80 nage90 winter male dept2 edu_bachelor median_household_income new_hyper asian hispanic black unknown nfac_id: npcp_id:,variance Performing EM optimization: Performing gradient-based optimization: Comment [JS1]: Age 60-69, female sex, and white race are omitted because they are the reference groups Comment [JS2]: This syntax specifies random effects for facility and primary care physician. By default, the second random effect is nested within the first in XTMIXED. Iteration 0: log restricted-likelihood = -1016708.8 Iteration 1: log restricted-likelihood = -1016708.8 Computing standard errors: Mixed-effects REML regression Number of obs = 232053 ----------------------------------------------------------- No. of Observations per Group Group Variable Groups Minimum Average Maximum ----------------+------------------------------------------ nfac_id 35 502 6630.1 18158 npcp_id 1049 26 221.2 793 ----------------------------------------------------------- Wald chi2(17) = 4083.54 Log restricted-likelihood = -1016708.8 Prob > chi2 = 0.0000 syst Coef. Std. Err. z P> z [95% Conf. Interval] -------------+---------------------------------------------------------------- nage20-5.467134.5624796-9.72 0.000-6.569574-4.364695 nage30-4.78875.2387905-20.05 0.000-5.25677-4.320729 nage40-4.008498.1449572-27.65 0.000-4.292608-3.724387 nage50-2.167071.1164383-18.61 0.000-2.395286-1.938856 nage70 1.216917.1139581 10.68 0.000.9935635 1.440271 nage80 1.423788.1475498 9.65 0.000 1.134595 1.71298 nage90 1.28458.3845358 3.34 0.001.5309031 2.038256 winter 1.249117.0881677 14.17 0.000 1.076312 1.421923 male -2.066664.086103-24.00 0.000-2.235423-1.897905 dept2-2.458218.1255729-19.58 0.000-2.704336-2.212099 edu_bachelor -.6203133.1330001-4.66 0.000 -.8809886 -.3596379 median_hou~e.5075523.1270388 4.00 0.000.2585608.7565439 new_hyper 1.494117.0822566 18.16 0.000 1.332897 1.655336 asian -1.110282.1369311-8.11 0.000-1.378662 -.8419017 hispanic.5464698.1438946 3.80 0.000.2644415.8284981 black 2.149833.142312 15.11 0.000 1.870907 2.42876 unknown 1.604231.1421661 11.28 0.000 1.32559 1.882871 _cons 143.5912.4358851 329.42 0.000 142.7369 144.4455
Random-effects Parameters Estimate Std. Err. [95% Conf. Interval] -----------------------------+------------------------------------------------ nfac_id: Identity var(_cons) 5.460574 1.46575 3.226675 9.241051 -----------------------------+------------------------------------------------ npcp_id: Identity var(_cons) 7.712068.4409736 6.894446 8.626652 -----------------------------+------------------------------------------------ var(residual) 371.3487 1.092643 369.2133 373.4964 LR test vs. linear regression: chi2(2) = 5353.29 Prob > chi2 = 0.0000 Comment [JS3]: Note that Confidence Intervals are asymmetric Note: LR test is conservative and provided only for reference.
GLIMMIX Procedure in STATA for Mammography Screening in Women 52-69 Yrs, 2001. use "c:\mamm01.dta", clear gllamm f_mamm age55 age65 age70 r_other r_black r_asian r_hisp edu_bachelor median_household_income, i(npcp_id nfac_id) fam(binom) link(logit) number of level 1 units = 236360 number of level 2 units = 1088 number of level 3 units = 35 Condition Number = 3646.3242 gllamm model log likelihood = -120450.06 f_mamm Coef. Std. Err. z P> z [95% Conf. Interval] -------------+---------------------------------------------------------------- age55 -.065825.0137513-4.79 0.000 -.0927771 -.038873 age65.0564901.0138102 4.09 0.000.0294227.0835576 age70.0515401.0145857 3.53 0.000.0229526.0801275 r_other -.8605286.0139808-61.55 0.000 -.8879305 -.8331267 r_black -.0319111.0217078-1.47 0.142 -.0744577.0106354 r_asian.0015076.0180698 0.08 0.934 -.0339086.0369238 r_hisp.0571813.0193694 2.95 0.003.0192179.0951447 edu_bachelor.5486208.0422151 13.00 0.000.4658808.6313609 median_hou~e 263.0839 30.26775 8.69 0.000 203.7602 322.4076 _cons.8255854.0285977 28.87 0.000.7695349.8816358 Variances of random effects ***level 2 (npcp_id) var(1):.11756377 (.00576969) ***level 3 (nfac_id) Comment [JS4]: Standard errors indicating uncertainty surrounding the variance estimates. var(1):.02052715 (.00433772)
Table FAC LEVEL:ADJUSTED SBP YR 2001 11:40 Wednesday, July 1, 2009 1 The UNIVARIATE Procedure Variable: R Moments N 35 Sum Weights 35 Mean 142.885429 Sum Observations 5000.99001 Std Deviation 2.58775442 Variance 6.69647296 Skewness 0.69461782 Kurtosis 0.80450613 Uncorrected SS 714796.283 Corrected SS 227.680081 Coeff Variation 1.8110695 Std Error Mean 0.43741033 Basic Statistical Measures Location Variability Mean 142.8854 Std Deviation 2.58775 Median 142.5621 Variance 6.69647 Mode. Range 12.21885 Interquartile Range 3.84372 Tests for Location: Mu0=0 Test -Statistic- -----p Value------ Student's t t 326.6622 Pr > t <.0001 Sign M 17.5 Pr >= M <.0001 Signed Rank S 315 Pr >= S <.0001 Tests for Normality Test --Statistic--- -----p Value------ Shapiro-Wilk W 0.966003 Pr < W 0.3432 Kolmogorov-Smirnov D 0.090118 Pr > D >0.1500 Cramer-von Mises W-Sq 0.049807 Pr > W-Sq >0.2500 Anderson-Darling A-Sq 0.352241 Pr > A-Sq >0.2500 Quantiles (Definition 5) Quantile Estimate 100% Max 150.396 99% 150.396 95% 147.573 90% 146.350
75% Q3 144.742 50% Median 142.562 25% Q1 140.899 10% 140.214 5% 138.808 1% 138.177 0% Min 138.177
FAC LEVEL:ADJUSTED SBP YR 2001 11:40 Wednesday, July 1, 2009 2 The UNIVARIATE Procedure Variable: R Extreme Observations ------Lowest----- -----Highest----- Value Obs Value Obs 138.177 7 145.326 19 138.808 3 146.350 35 140.166 32 146.824 13 140.214 20 147.573 1 140.379 6 150.396 26 Stem Leaf # Boxplot Normal Probability Plot 150 4 1 150.5+ * 149 ++ 148 ++++ 147 6 1 +*++ 146 38 2 **+ 145 23 2 +**+ 144 013789 6 +-----+ 144.5+ ****** 143 0344 4 **** 142 233467 6 *--+--* ***** 141 0024 4 ++** 140 2244579 7 +-----+ ** ****** 139 ++++ 138 28 2 138.5+ * +*+ ----+----+----+----+ +----+----+----+----+----+----+----+----+----+----+ -2-1 0 +1 +2
PROVIDER LEVEL:ADJUSTED SBP YR 2001 11:40 Wednesday, July 1, 2009 3 The UNIVARIATE Procedure Variable: R Moments N 1049 Sum Weights 1049 Mean 142.561263 Sum Observations 149546.765 Std Deviation 3.80450105 Variance 14.4742282 Skewness -0.1321348 Kurtosis 0.39300123 Uncorrected SS 21334744.8 Corrected SS 15168.9912 Coeff Variation 2.66867798 Std Error Mean 0.1174654 Basic Statistical Measures Location Variability Mean 142.5613 Std Deviation 3.80450 Median 142.5250 Variance 14.47423 Mode. Range 30.05818 Interquartile Range 4.95907 Tests for Location: Mu0=0 Test -Statistic- -----p Value------ Student's t t 1213.645 Pr > t <.0001 Sign M 524.5 Pr >= M <.0001 Signed Rank S 275362.5 Pr >= S <.0001 Tests for Normality Test --Statistic--- -----p Value------ Shapiro-Wilk W 0.996812 Pr < W 0.0328 Kolmogorov-Smirnov D 0.017759 Pr > D >0.1500 Cramer-von Mises W-Sq 0.048125 Pr > W-Sq >0.2500 Anderson-Darling A-Sq 0.359705 Pr > A-Sq >0.2500 Quantiles (Definition 5) Quantile Estimate 100% Max 154.969 99% 151.035 95% 148.775 90% 147.637
75% Q3 145.069 50% Median 142.525 25% Q1 140.110 10% 137.856 5% 136.561 1% 133.141 0% Min 124.911
PROVIDER LEVEL:ADJUSTED SBP YR 2001 11:40 Wednesday, July 1, 2009 4 The UNIVARIATE Procedure Variable: R Extreme Observations ------Lowest----- -----Highest----- Value Obs Value Obs 124.911 278 152.060 87 129.557 614 152.418 157 130.817 541 152.653 131 130.930 306 153.060 159 130.953 893 154.969 645 Histogram # Boxplot Normal Probability Plot 155+* 1 0 155+ *.* 5 0 *.**** 16 ****** 149+************* 64 149+ *******.********************** 109 ******.********************************* 165 +-----+ ****** 143+********************************************* 222 *--+--* 143+ ******.******************************************* 211 +-----+ *******.***************************** 144 ****** 137+*************** 71 137+ *******.***** 25 +****.** 10 0 **** 131+* 4 0 131+*.* 1 0 *. 125+* 1 * 125+* ----+----+----+----+----+----+----+----+----+ +----+----+----+----+----+----+----+----+----+----+ * may represent up to 5 counts -2-1 0 +1 +2