医学统计学实习课：STATA软件的应用(4)

1、 医学统计学实习课 STATA软件的应用(4) 邮箱账号：stata_密码：stata123456主要内容1.秩和检验；2.直线相关分析；3.等级相关分析；4.直线回归分析；5.复习秩和检验 秩变换genrank 新变量=原变量例 ex1.dta对数据 0 -2 -3 4 ? 5 47 0 进行编秩，其中?为缺失值。 x rx 1. -3 1 2. -2 2 3. 0 3.5 4. 0 3.5 5. 4 5 6. 5 6 7. 47 7 8. . . genrank rx=xegen rx=rank(x)list配对秩和检验signrank 变量1=变量2signrank 差值变量=0例 ex

2、2.dta某医院检验科试用新旧两种方法监测谷-草转氨酶，新方法的检测时间由20分钟缩短为10分钟。用两种方法测一份血清，结果见表，问两法所得结果有无差别？样品号旧法新法14060213214232122104808253825621224372302378951009236200103843signrank a=bsignrank a=bgen diff=a-bsignrank diff=0两样本比较秩和检验ranksum 观察值变量，by(分组变量)例 ex3.dta测得铅作业和非铅作业工人的血铅值(g/100mg)，问两组工人的血铅值有无差别？小样本定量资料不满足t检验要求时（方差不齐），

3、可以考虑采用秩和检验非铅作业组56791213151821铅作业组 17182025344344方差齐性检验方差不齐两样本秩和例 ex4.dta用复方猪胆胶囊治疗老年性慢性支气管炎患者403例，疗效见下表。问该药对此两型支气管炎疗效是否相同？疗效人数合计喘息型单纯型治愈236083显效8398181好转6551116无效111223expand franksum x,by(g)两组的定性资料秩和多样本秩和检验kwallis 观察值变量, by(分组变量)例 ex5.dta三种复方制剂治疗慢性胃炎的疗效比较疗效例数复方复方复方痊愈4256显效1861720有效753626无效504231expa

4、nd fkwallis x,by(g)多组的定性资料秩和相关与回归线性相关线性回归等级相关列联相关线性相关correlate 变量名，选择项注意：这里的变量可以是两个或多个，correlate命令给出的是变量间的两两的简单相关系数矩阵。选择项： means /*同时输出均数、标准差等统计量 covariance /*不输出相关系数矩阵，而输出协方差矩阵相关系数假设检验pwcorr 变量名，选择项obs 打印样本含量sig 打印假设检验之P值star(#) 如果相关系数的假设检验之P值小于#，则在相关系数旁打星号 例 ex6.dta以下资料选自Galton的一项研究，目的是探讨成年时的身高是否与

5、2岁时的身高（单位：英寸）有关。2岁时的身高3930323435363630成年身高7163636768687064分析步骤(1) 画散点图scatter y x(2) 计算相关系数corr y x(3) 相关系数的假设检验pwcorr y x,sig star(0.05) obspwcorr y x,sig star(0.05)(是否为0的假设检验)(区间估计)线性回归回归分析：regress 因变量 自变量回归预测：predict 新变量分析步骤（续 ex6.dta）(1) 画散点图scatter y x最小二乘法（SS剩最小）(2) 回归分析SS回归SS剩余SS总reg y x建立回归方

6、程后，可立即用该回归方程进行回代预测，并作回归线图。执行该命令后，stata将产生一个新变量yhat，并将根据回归方程估计的值写入该变量中。(3) 回归预测predict yhatline yhat x|scatter y x该命令作出的是yhat与x的回归图及y与x的散点图(4) 画回归直线等级相关spearman x y注意：该命令首先求出等级相关系数，然后进行假设检验。例 ex7.dta某研究所用野百合治疗白血病，并作抗白血病指数（简称抗白指数）及疗效的分析，结果见表，问抗白指数与临床疗效间有无关系？抗白指数X临床疗效Y223292103113300350450551701880900直线回归与t检验、F检验的关系方案血红蛋白增加量(g)A1.80.52.33.72.42.01.52.71.10.9(n=20)1.41.22.30.70.51.41.73.03.22.5B5.00.20.50.31.91.02.4-0.42.01.6(n=20)2.00.01.63.01.60.03.00.71.20.7ex8.dtat检验：ttest x,by(group)方差分析：oneway x group回归分析：reg x groupt=1.1981P=0.2383F=1.44P=0.2383F=1.44P=0.2383F=t2复习掌握输出结