统计学5相关与回归

医学统计学实习课STATA软件的应用(5)--相关与回归相关与回归分析线性相关分析线性回归分析等级相关分析列联相关分析

21线性相关分析例1以下资料选自Galton的一项研究。探讨成年时的身高是否与2岁时的身高存在相关。(ex_1.dta)单位：英寸

32岁身高3930323435363630成年身高71636367686870641线性相关分析考察变量间线性相关趋势绘制两个变量间的散点图scatter变量1变量2scatteryx

42岁身高与成年身高的散点图

51线性相关分析计算变量间的线性相关系数correlate变量列表correlateyx

61线性相关分析线性相关系数的假设检验pwcorr变量列表,[选项]pwcorryx,sigstar(0.05)sig /*显示假设检验P值*/star(0.05) /*星号标记P≤0.05者*/

71线性相关分析correlate变量名[，选择项]注意：这里的变量可以是两个或多个，corrlate命令给出的是变量间的两两的简单相关系数矩阵。选择项：

means /*同时输出均数、标准差等统计量

covariance /*不输出相关系数矩阵，而输出协方差

矩阵pwcorr变量名,[选择项]sig打印假设检验之P值star(#)如果相关系数的假设检验之P值小于#，

则在相关系数旁打星号

2线性回归分析例2请按例1提供的数据。能否利用2岁时的身高预测成年后的身高。已经判断两个变量间存在线性相关建立模型建立2岁身高与成年身高间的回归模型模型预测将2岁身高纳入回归模型预测成年身高

92线性回归分析建立模型regress因变量自变量regressyx

102线性回归分析

112线性回归分析模型预测predict新变量,[选项]predictyhat,xb /*线性预测值*/predictresidual,residuals /*残差*/predictse,stdp /*预测值标准误*/predictsd,stdf /*观察值标准差*/

122线性回归分析估计预测值的95%可信区间genCI_L=yhat-invttail(6,0.025)*segenCI_U=yhat+invttail(6,0.025)*se估计观察值的95%容许区间genR_L=yhat-invttail(6,0.025)*sdgenR_U=yhat+invttail(6,0.025)*sd

132线性回归分析绘制回归分析五线谱sortxscatteryx||lineyhatx||lineCI_Lx||lineCI_Ux||lineR_Lx||lineR_Ux

14回归分析五线谱

153等级相关分析例3某研究所用野百合治疗白血病，并作抗白指数及疗效的分析，结果见表，问抗白指数与临床疗效间有无关系？(ex_3.dta)

163等级相关分析spearmanxy

174列联相关（定性资料的相关）tabvar1var2[fw=频数变量][,V]4列联相关例4某县防疫站观察三种药物驱钩虫的疗效，在服药后7天得粪检钩虫卵阴转率（%）如下，问三药疗效是否不同？(ex_4.dta)3种不同药物与是否转阴之间有无相关？tabdrugeffect[fw=f],V

19药物阴转例数未阴转例数合计阴转率(%)复方敌百虫片2893775.7纯敌百虫片18203847.4灭虫灵10243429.4合计565310951.4直线回归与t检验、F检验的关系(ex_5.dta)方案血红蛋白增加量(g)A(n=20)1.80.52.33.72.42.01.52.71.10.91.41.22.30.70.51.41.73.03.22.5B(n=20)5.00.20.50.31.91.02.4-0.42.01.62.00.01.63.01.60.03.00.71.20.7用A和B组研究t检验和直线回归的关系ttestx,by(group)regxgroup

22练习(lx.dta)44名女医学生的肺活量(L)和身高(cm)资料，试作如下分析:（1）以肺活量作为应变量Y，身高作为自变量X，做回归分析；（2）求相当于各身高值的肺活量总体均数的95%可信区间及个体值的95%容许区间，并绘制回归“五线谱”。

scatteryxregyxpredictyhatpredictseyhat,stdpgenl1=yhat-invttail(42,0.025)*seyhatgenl2=yhat+invttail(42,0.025)*seyhatpredictsey,stdfgenl3=yhat-invttail(42,0.025)*seygenl4=yhat+invtt