STATA软件应用(三)定量资料的统计分析

1、 医学统计学实习课 STATA软件的应用(3) 定量资料的统计分析 定量资料的统计分析样本均数与总体均数比较的t 检验配对设计t 检验 成组设计t 检验单因素方差分析两因素方差分析方差齐性检验正态性检验样本均数与总体均数比较的t检验o仅有数据之统计量时ttesti #obs #mean #sd #valo仅有原始数据时ttest 变量名= #val o注： #val /*总体均数 命令可以加一些if和in条件限制样本均数与总体均数比较的t检验o例4.4 命令：ttesti 25 5.1 0.88 4.6样本均数与总体均数比较的t检验 例：某医生随机抽查10名某病患者的血红蛋白，如下：问：该病患

2、者的平均Hb含量是否与正常人的平均Hb含量相同 (正常人的平均Hb含量为14.02(g/dl) )。data:ttest_1病例号病例号12345678910血红蛋白血红蛋白(x,g/dl)(x,g/dl)11.3151513.512.81011121312.3问题：问题： 014.02 1？ 正常人 患者样本含量：10 均 数: 12.59 g/dl标准差: 1.63 g/dl统计量与参数不同的两种可能统计量与参数不同的两种可能o其一：抽样误差 (偶然的、随机的、较小的)o其二：本质上的差别 (必然的、大于随机误差)样本均数与总体均数比较的t检验ttest x=14.02One-sample

3、 t test-Variable | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+- x | 10 12.59 .5162794 1.632619 11.42209 13.75791- mean = mean(x) t = -2.7698Ho: mean = 14.02 degrees of freedom = 9 Ha: mean 14.02 Pr(T |t|) = 0.0218 Pr(T t) = 0.9891配对设计t检验ottest 变量1=变量2o变量1与变量2必须成对输入，样本含量必须相等，如有缺省值，则用小数点表示，但与之对应

4、的记录在计算时被忽略配对设计t检验例：10例矽肺患者经某药治疗，其血红蛋白(g/dl)如下：data:ttest_2病例号:12345678910治疗前(x1):11.315.015.013.512.810.011.012.013.012.3治疗后(x2):14.013.814.013.513.512.014.711.413.812.0配对设计t检验. ttest x1=x2Paired t test-Variable | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+- x1 | 10 12.59 .5162794 1.632619 11.

5、42209 13.75791 x2 | 10 13.27 .3415813 1.080175 12.49729 14.04271-+- diff | 10 -.6799999 .5204272 1.645735 -1.857288 .4972881- mean(diff) = mean(x1 - x2) t = -1.3066 Ho: mean(diff) = 0 degrees of freedom = 9 Ha: mean(diff) 0 Pr(T |t|) = 0.2237 Pr(T t) = 0.8881配对设计t检验另一种思路： gen d=x0-x1 ttest d=0 结果一致成

6、组设计t检验ottest 变量1=变量2, unpaired unequalottest 变量,by（分组变量）unequalottesti #obs1 #mean1 #sd1 #obs2 #mean2 #sd2 ,unequal unpaired 表示非配对的，如不选就作配对检验 unequal 表示假设两组方差不齐，如不选表示假设两组方差达到齐性 成组设计t检验例：分别测得14例老年人煤饼病人及11例正常人的尿中17 酮类固醇排出量(mg/dl)如下，试比较两组的均数有无差别 病人病人： 2.905.415.484.604.035.104.974.244.362.722.37 2.097.

7、105.92健康人健康人：5.188.793.146.463.726.645.604.577.714.994.01data:ttest_3.dta ttest_4.dta成组设计t检验ttest x1=x2, unpairedttest x, by(g)q独立性q正态性q方差齐性 方差齐：t检验 方差不齐:t检验 均数的假设检验应用条件两组资料间的方差齐性检验osdtest 变量名1 = 变量名2osdtest 变量,by（分组变量）osdtesti #obs1 #mean1 #sd1 #obs2 #mean2 #sd2 sdtest x,by(g)Variance ratio test- G

8、roup | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+- 1 | 14 4.377857 .3875 1.449892 3.540714 5.215 2 | 11 5.528182 .5232431 1.735401 4.362324 6.69404-+-combined | 25 4.884 .3306453 1.653227 4.201582 5.566418- ratio = sd(1) / sd(2) f = 0.6980 Ho: ratio = 1 degrees of freedom = 13, 10 Ha: ratio 1

9、 Pr(F f) = 0.2675 2*Pr(F f) = 0.7325单因素方差分析ooneway 因变量 分组变量，选择项noanova /* 不打印方差分析表missing /* 将缺省值作为单独的一组tabulate /* 打印各组的基本统计量表 简写：tscheffe /* Scheffe法 简写：schBonferroni /* Bonferroni法 简写：bonsidak /* Sidak法 简写：si 单因素方差分析o三组小鼠的三组小鼠的FDP酶活力酶活力 对照组水层RNA组酚层RNA组2.79 3.83 5.412.69 3.15 3.473.11 4.70 4.923.4

10、7 3.97 4.071.77 2.03 2.182.44 2.87 3.132.83 3.65 3.772.52 5.09 4.26单因素方差分析ooneway x groupooneway x group , tooneway x group ,t schdata:oneway.dtao. oneway x g,t scho | Summary of xo g | Mean Std. Dev. Freq.o-+-o 1 | 2.7025 .50013569 8o 2 | 3.66125 .98508069 8o 3 | 3.9012501 1.0164425 8o-+-o Total |

11、3.4216667 .98273207 24 Analysis of Variance Source SS df MS F Prob F-Between groups 6.43680897 2 3.21840449 4.28 0.0275 Within groups 15.7757246 21 .751224983- Total 22.2125336 23 .965762331Bartletts test for equal variances: chi2(2) =3.4559 Probchi2=0.178 Comparison of x by g (Scheffe)Row Mean-|Col

12、 Mean | 1 2-+- 2 | .95875 | 0.111 | 3 | 1.19875 .24 | 0.038 0.859两因素方差分析o四种抗癌药物抑瘤效果配伍组a1a2a3a4b10.800.360.170.28b20.740.500.420.36b30.310.200.380.25b40.480.180.440.22b50.760.260.280.13qanova 因变量 分组变量1 分组变量2 两因素的方差分析qegen a=seq(), b(5)qegen b=seq(), f(1) t(5)qanova x a bdata：anova.dta两因素的方差分析. anova

13、x group block Number of obs = 20 R-squared = 0.7058 Root MSE = .134818 Adj R-squared = 0.5341 Source | Partial SS df MS F Prob F -+- Model | .523170004 7 .074738572 4.11 0.0157 | group | .410839998 3 .136946666 7.53 0.0043 block | .112330006 4 .028082501 1.55 0.2514 | Residual | .218110001 12 .01817

14、5833 -+- Total | .741280005 19 .039014737 正态性检验o命令：sktest 变量 该命令要求资料的样本含量至少为8 Skewness/Kurtosis tests for Normality - joint -Variable | Pr(Skewness) Pr(Kurtosis) adj chi2(2) Probchi2-+- db | 0.304 0.304 2.53 0.2829小结o样本均数与总体均数比较的t检验 ttest x= 14.02 o配对设计检验 ttest x1=x2o成组设计检验及方差齐性检验 sdtest x1=x2 ttest x1=x2, unpaired unequal ttest x,by（group）unequalo单因素方差分析 oneway x g,t scho两因素方差分析 anova x a b 练习例4.5（P55）:配对t检验例4.7（P58）：成组t检验例5.1（ P68）：单因素方差分析思考题 某医师为评价两种药物对血沉的影响，用药前后结果如下，某医师为评价两种药物对血沉的影响，用药前后结果如下，问两种药物对血沉影响有无差别。问两种药物对血沉影响有无差别。 A药组 B药组前 后前后40 3639