第十章直线回归和相关Stata实现.doc

第十章 直线回归和相关的Stata实现本章使用的Stata命令为：直线回归regress 应变量 自变量残差计算在执行regress命令后，执行下列命令：predict 变量名,residualPearson相关pwcorr 变量1 变量n, sigSpearman相关spearman 变量1 变量2 例101 为了研究血清胆固醇含量与舒张压之间是否存在依存关系，2006年在郑州某大学随机抽取10名成年男性，测得他们的血清胆固醇(mg/dL)含量和舒张压（mmHg）如表10-1，请作统计分析。表10-1 10名成人的血清胆固醇(mg/dL)含量和舒张压（mmHg）指标编 号12345678910胆固醇307259341237254416267320374316舒张压807590707510570858878：总体回归方程不成立（）：总体回归方程成立（）=0.05Stata数据：xy13078022597533419042377052547564161057267708320859374881031678Stata命令为：reg y x结果： Source | SS df MS Number of obs = 10-+- F( 1, 8) = 65.11 Model | 928.343132 1 928.343132 Prob F = 0.0000 Residual | 114.056868 8 14.2571086 R-squared = 0.8906-+- Adj R-squared = 0.8769 Total | 1042.4 9 115.822222 Root MSE = 3.7759- y | Coef. Std. Err. t P|t| 95% Conf. Interval-+- x | .177864 .0220419 8.07 0.000 .1270352 .2286927 _cons | 26.62225 6.916995 3.85 0.005 10.67163 42.57287-t=8.07 ,P值0.001(Stata输出值0.000),构建直线回归方程 将和代入式（10-2），可知， 本例中，的统计学意义为：血清胆固醇含量每增加1mg/dL，总体中舒张压平均增加0.178mmHg。总体均数的区间估计 给定时，的总体均数的点估计，例101中，当自变量取值为307 mg/dL时。个体值的容许区间估计 给定值时，估计总体中个体值的波动范围，以例10-1中第一个样本点的数据（307,80）为例。Stata命令：predict yhatpredict stdp, stdppredict stdf, stdfgenerate clm1 = yhat - invttail(7,0.025)*stdpgenerate clm2 = yhat + invttail(7,0.025)*stdpgenerate clp1 = yhat - invttail(7,0.025)*stdfgenerate clp2 = yhat + invttail(7,0.025)*stdftwoway (lfitci y x, level(95) (scatter y x) (line clm1 clm2 clp1 clp2 x, pstyle(p2 p2 p3 p3) sort) 结果：xyyhatstdpstdfclm1clm2clp1clp23078081.231.193.9678.484.0571.8690.592597572.691.634.1168.8476.5362.9782.413419087.271.394.028490.5577.7696.782377068.781.994.2764.0873.4858.6978.872547571.81.74.1467.7775.836281.59416105100.612.644.6194.37106.8689.72111.512677074.111.514.0770.5477.6964.4983.733208583.541.223.9780.6686.4274.1692.923748893.141.864.2188.7497.5583.19103.13167882.831.23.9679.9885.6773.4692.2例102 某发热门诊医生根据患者就诊顺序随机抽取12名2040岁发热患者，试探讨体温与脉搏之间的伴随关系，数据见表10-4。表10-4 12名发热患者的体温（）与脉搏（次/分）指标编 号123456789101112体温 41.037.540.738.438.840.038.039.039.539.840.839.6脉搏 13585134108120124106117125130138135Stata数据xy141135237.585340.7134438.4108538.8120640124738106839117939.51251039.81301140.81381239.6135H0: 0 H1: 0Stata命令：pwcorr y x,sig结果： |