第十七章 多因素回归分析的Stata 实现 本章使用的Stata 命令: 多因素回归 regress depvar [indepvars] 逐步回归 stepwise [, options ] : command Logistic 回归 logistic depvar indepvars [weight] 生存时间数据设定 stset timevar [weight] [, failure(failvar[==numlist])] Cox回归 stcox [varlist] 例17-4 某研究者为了研究某种避孕药对人体血糖的影响,分别在正在使用这种避孕药的人群、6 个月前曾经使用过这种避孕药的人群、从未使用过避孕药的人群中各随机抽取 6 人。考虑到血糖可能与年龄有关,所以该研究者不仅测定了这18 位对象的血糖,而且也记录了这18 位对象的年龄,具体资料见表 17-4。请根据研究问题作统计分析。 表 17-4 三种避孕药使用情况下的年龄( X ,岁)与血糖水平(Y ,mg%) 现服药者 曾服药者 从未服药者 X Y X Y X Y 20 120 24 126 28 135 21 122 26 130 32 137 23 124 27 132 34 138 23 126 29 131 35 137 24 125 29 134 35 139 24 127 30 136 37 144 本研究的问题是比较三种用药情况下的血糖平均水平是否不同,因此首先考虑以下总体均数的情况。 解:Stata 数据如下: x y g1 g2 20 120 0 1 21 122 0 1 23 124 0 1 23 126 0 1 24 125 0 1 24 127 0 1 24 126 1 0 26 130 1 0 27 132 1 0 29 131 1 0 29 134 1 0 30 136 1 0 28 135 0 0 32 137 0 0 34 138 0 0 35 137 0 0 35 139 0 0 37 144 0 0 Stata 命令如下: reg y x g1 g2 结果: Source | SS df MS Number of obs = 18 -------------+------------------------------ F( 3, 14) = 83.88 Model | 717.684557 3 239.228186 Prob > F = 0.0000 Residual | 39.9265537 14 2.85189669 R-squared = 0.9473 -------------+------------------------------ Adj R-squared = 0.9360 Total | 757.611111 17 44.5653595 Root MSE = 1.6888 ------------------------------------------------------------------------------ y | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- x | 1.067797 .1795127 5.95 0.000 .6827801 1.4...