Table. XTMIXED Procedure in STATA with Output Systolic Blood Pressure, use "k:mydirectory,

Transcription

1 Table XTMIXED Procedure in STATA with Output Systolic Blood Pressure, 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 = Iteration 1: log restricted-likelihood = Computing standard errors: Mixed-effects REML regression Number of obs = No. of Observations per Group Group Variable Groups Minimum Average Maximum nfac_id npcp_id Wald chi2(17) = Log restricted-likelihood = Prob > chi2 = syst Coef. Std. Err. z P> z [95% Conf. Interval] nage nage nage nage nage nage nage winter male dept edu_bachelor median_hou~e new_hyper asian hispanic black unknown _cons

2 Random-effects Parameters Estimate Std. Err. [95% Conf. Interval] nfac_id: Identity var(_cons) npcp_id: Identity var(_cons) var(residual) LR test vs. linear regression: chi2(2) = Prob > chi2 = Comment [JS3]: Note that Confidence Intervals are asymmetric Note: LR test is conservative and provided only for reference.

3 GLIMMIX Procedure in STATA for Mammography Screening in Women Yrs, 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 = number of level 2 units = 1088 number of level 3 units = 35 Condition Number = gllamm model log likelihood = f_mamm Coef. Std. Err. z P> z [95% Conf. Interval] age age age r_other r_black r_asian r_hisp edu_bachelor median_hou~e _cons Variances of random effects ***level 2 (npcp_id) var(1): ( ) ***level 3 (nfac_id) Comment [JS4]: Standard errors indicating uncertainty surrounding the variance estimates. var(1): ( )

4 Table FAC LEVEL:ADJUSTED SBP YR :40 Wednesday, July 1, The UNIVARIATE Procedure Variable: R Moments N 35 Sum Weights 35 Mean Sum Observations Std Deviation Variance Skewness Kurtosis Uncorrected SS Corrected SS Coeff Variation Std Error Mean Basic Statistical Measures Location Variability Mean Std Deviation Median Variance Mode. Range Interquartile Range Tests for Location: Mu0=0 Test -Statistic p Value Student's t t 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 Pr < W Kolmogorov-Smirnov D Pr > D > Cramer-von Mises W-Sq Pr > W-Sq > Anderson-Darling A-Sq Pr > A-Sq > Quantiles (Definition 5) Quantile Estimate 100% Max % % %

5 75% Q % Median % Q % % % % Min

6 FAC LEVEL:ADJUSTED SBP YR :40 Wednesday, July 1, The UNIVARIATE Procedure Variable: R Extreme Observations Lowest Highest----- Value Obs Value Obs Stem Leaf # Boxplot Normal Probability Plot * * ** ** ****** **** *--+--* ***** ** ** ****** * +*

7 PROVIDER LEVEL:ADJUSTED SBP YR :40 Wednesday, July 1, The UNIVARIATE Procedure Variable: R Moments N 1049 Sum Weights 1049 Mean Sum Observations Std Deviation Variance Skewness Kurtosis Uncorrected SS Corrected SS Coeff Variation Std Error Mean Basic Statistical Measures Location Variability Mean Std Deviation Median Variance Mode. Range Interquartile Range Tests for Location: Mu0=0 Test -Statistic p Value Student's t t Pr > t <.0001 Sign M Pr >= M <.0001 Signed Rank S Pr >= S <.0001 Tests for Normality Test --Statistic p Value Shapiro-Wilk W Pr < W Kolmogorov-Smirnov D Pr > D > Cramer-von Mises W-Sq Pr > W-Sq > Anderson-Darling A-Sq Pr > A-Sq > Quantiles (Definition 5) Quantile Estimate 100% Max % % %

8 75% Q % Median % Q % % % % Min

9 PROVIDER LEVEL:ADJUSTED SBP YR :40 Wednesday, July 1, The UNIVARIATE Procedure Variable: R Extreme Observations Lowest Highest----- Value Obs Value Obs Histogram # Boxplot Normal Probability Plot 155+* *.* 5 0 *.**** 16 ****** 149+************* *******.********************** 109 ******.********************************* ****** 143+********************************************* 222 *--+--* 143+ ******.******************************************* *******.***************************** 144 ****** 137+*************** *******.***** 25 +****.** 10 0 **** 131+* *.* 1 0 *. 125+* 1 * 125+* * may represent up to 5 counts