第十七篇多因素回归分析的Stata实现.docVIP

第十七章 多因素回归分析的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。请根据研究问题作统计分析。表17- 三种避孕药使用情况下的年龄（，岁）与血糖水平(，mg％) 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 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.452813 g1 | -.4265537 1.452834 -0.29 0.773 -3.542572 2.689465 g2 | -2.587571 2.202234 -1.17 0.260 -7.310893 2.135752 _cons | 102.5621 6.053067 16.94 0.000 89.57961 115.5447 --------------------------------------------------------------------