《医学统计学实习课：STATA软件的应用(4)》由会员分享，可在线阅读，更多相关《医学统计学实习课：STATA软件的应用(4)（37页珍藏版）》请在金锄头文库上搜索。
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=xgenrank rx=xegen rx=rank(x)egen rx=rank(x)listlist配
2、对秩和检验signrank 变量1=变量2signrank 差值变量=0例 ex2.dta某医院检验科试用新旧两种方法监测谷-草转氨酶，新方法的检测时间由20分钟缩短为10分钟。用两种方法测一份血清，结果见表，问两法所得结果有无差别？样品号旧法新法14060213214232122104808253825621224372302378951009236200103843signrank a=bsignrank a=bsignrank a=bsignrank a=bgen diff=a-bgen diff=a-bsignrank diff=0signrank diff=0两样本比较秩和检验ran
3、ksum 观察值变量，by(分组变量)例 ex3.dta测得铅作业和非铅作业工人的血铅值(g/100mg)，问两组工人的血铅值有无差别？小样本定量资料不满足t检验要求时（方差不齐），可以考虑采用秩和检验非铅作业组56791213151821铅作业组 17182025344344方差齐性检验方差不齐两样本秩和例 ex4.dta用复方猪胆胶囊治疗老年性慢性支气管炎患者403例，疗效见下表。问该药对此两型支气管炎疗效是否相同？疗效人数合计喘息型单纯型治愈236083显效8398181好转6551116无效111223expand fexpand franksum x,by(g)ranksum x,b
4、y(g)两组的定性资料秩和多样本秩和检验kwallis 观察值变量, by(分组变量)例 ex5.dta三种复方制剂治疗慢性胃炎的疗效比较疗效例数复方复方复方痊愈4256显效1861720有效753626无效504231expand fexpand fkwallis x,by(g)kwallis x,by(g)多组的定性资料秩和相关与回归线性相关线性回归等级相关列联相关线性相关correlate 变量名，选择项注意：这里的变量可以是两个或多个，correlate命令给出的是变量间的两两的简单相关系数矩阵。选择项：n means /*同时输出均数、标准差等统计量n covariance /*不输
5、出相关系数矩阵，而输出协方差矩阵相关系数假设检验pwcorr 变量名，选择项nobs 打印样本含量nsig 打印假设检验之P值nstar(#) 如果相关系数的假设检验之P值小于#，则在相关系数旁打星号 例 ex6.dta以下资料选自Galton的一项研究，目的是探讨成年时的身高是否与2岁时的身高（单位：英寸）有关。2岁时的身高3930 323435363630成年身高7163 636768687064分析步骤(1) 画散点图scatter y xscatter y x(2) 计算相关系数corrcorr y x y x(3) 相关系数的假设检验pwcorr y x,sig star(0.05)
6、 obspwcorr y x,sig star(0.05) obspwcorr y x,sig star(0.05)pwcorr y x,sig star(0.05)(是否为0的假设检验)(区间估计)线性回归回归分析：regress 因变量 自变量回归预测：predict 新变量分析步骤（续 ex6.dta）(1) 画散点图scatter y xscatter y x最小二乘法（SS剩最小）(2) 回归分析SS回归回归SS剩余剩余SS总总regreg y x y xo建立回归方程后，可立即用该回归方程进行回代预测，并作回归线图。n执行该命令后，stata将产生一个新变量yhat，并将根据回归方
7、程估计的值写入该变量中。(3) 回归预测predict yhatpredict 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
8、.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复习掌握输出结果中各项的意思重点数据描述中sum、ci、mean等命令定量资料ttest、anova、oneway等命令定性资料tab等命令等级资料ranksum、kwallis等命令相关回归corr、pwcorr、reg等命令其他
