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

1、 医学统计学实习课 STATA软件的应用(2) 数值变量资料的统计分析 样本均数与总体均数比较的t检验配对设计t检验 成组设计t检验及方差齐性检验单因素方差分析及方差齐性检验 两因素的方差分析样本均数与总体均数比较的t检验例 ex2-1.dta样本：某医生随机抽查10名某病患者的血红蛋白。问题：该病患者的平均Hb含量是否与正常人的平均Hb含量相同 (正常人的平均Hb含量为14.02(g/dl) )。病例号12345678910血红蛋白11.315.015.013.512.810.011.012.013.012.3问题： 014.02 1？ 正常人 患者？样本含量：10 均 数: 12.59 g

2、/dl标准差: 1.63 g/dl统计量与参数不同的两种可能其一：抽样误差 (偶然的、随机的、较小的)其二：本质上的差别 (必然的、大于随机误差)样本：某医生随机抽查10名某病患者的血红蛋白，求得其均数为12.59 (g/dl),标准差为1.632619 (g/dl) 。问题：该病患者的平均Hb含量是否与正常人的平均Hb含量相同 (正常人的平均Hb含量为14.02 (g/dl)。样本均数与总体均数比较的t检验(1) ttest 变量名= #val(总体均数) ttest x=14.02(2) ttesti #obs #mean #sd #val ttesti 10 12.59 1.632619

3、 14.02 样本均数与总体均数比较的t检验配对设计检验 例 ex2-2.dta10例矽肺患者经某药治疗，其血红蛋白(g/dl) ：病例号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配对设计检验 ttest 变量1=变量2成组设计例 ex2-3.dta分别测得14例老年人煤饼病人及11例正常人的尿中17 酮类固醇排出量(mg/dl)如下，试比较两组的均数有无差别 病人：2.905.415.484.604.035.104.974.2

4、44.362.722.372.097.105.92健康人：5.188.793.146.463.726.645.604.577.714.994.01成组设计 t 检验 ttest 变量1=变量2, unpaired unequalttest x, by(group) unequalttesti #obs1 #mean1 #sd1 #obs2 #mean2 #sd2 ,unequal unpaired 表示非配对的，如不选就作配对检验 unequal 表示假设两组方差不齐，如不选表示假设两组方差达到齐性 成组设计 t 检验成组设计 t 检验stack x2 x1,into(x)rename _st

5、ack group应用条件(成组) 独立性 正态性 方差齐性 （1）方差齐： t检验 （2）方差不齐: t检验方差齐性检验sdtest 变量名1 = 变量名2sdtesti #obs1 #mean1 #sd1 #obs2 #mean2 #sd2 方差齐性检验练习（样本与总体）测得25例某病女性患者血红蛋白(Hb),其均数为150(g/l),标准差为16.5(g/l),而正常成年女性的Hb均数为132(g/l),问该地区女性患者的Hb含量是否与正常成年女性含量不同？练习（配对设计）某医生研究脑缺氧对脑组织中生化指标的影响，将乳猪按出生体重配成7对，一组为对照组，一组为脑缺氧模型组，结果见表第（2

6、）、（3）栏。试比较两组猪脑组织钙泵的含量有无差别。练习（成组设计）某医生研究转铁蛋白对病毒性肝炎诊断的临床意义，测得12名正常人和15名病毒性肝炎患者血清转铁蛋白含量（ug/dl）,问患者和正常人转铁蛋白含量是否有差别？正常人：265.4 271.5 284.6291.3254.8275.9281.7268.6264.4273.2270.8260.5患者：235.9218.8215.4233.8251.8230.9224.7240.7228.3256.9231.1260.7253.0224.4221.7单因素方差分析例 ex2-4.dta 三组小鼠的FDP酶活力 对照组水层RNA组酚层RNA

7、组2.79 3.83 5.412.69 3.15 3.473.11 4.70 4.923.47 3.97 4.071.77 2.03 2.182.44 2.87 3.132.83 3.65 3.772.52 5.09 4.26单因素方差分析oneway 因变量 分组变量，选择项 tabulate /* 打印各组的基本统计量表 scheffe /* Scheffe法进行两两比较 （bonferroni, sidak）单因素方差分析(1) oneway x g(2) oneway x g , t(3) oneway x g , t sch单因素方差分析两两比较（ Scheffe ）配伍组a1a2a3a4b10.800.360.170.28b20.740.500.420.36b30.310.200.380.25b40.480.180.440.22b50.760.260.280.13两因素方差分析例 ex2-5.dta 四种抗癌药物抑瘤效果 两因素方差分析anova 因变量 分组变量1 分组变量2anova x a b 两因素方差分析样本均数与总体均数比较的 t 检验 ttest x= 14.02 配对设计 t