STATA软件应用(三)定量资料的统计分析解析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
mean(diff) = mean(x1 - x2)
t = -1.3066
Ho: mean(diff) = 0
degrees of freedom =
9
Ha: mean(diff) < 0 Pr(T < t) = 0.1119
Ha: mean(diff) != 0 Pr(|T| > |t|) = 0.2237
[,unequal]
➢ unpaired 表示非配对的,如不选就作配对t检验 ➢ unequal 表示假设两组方差不齐,如不选表示假设两
组方差达到齐性
成组设计t检验
例:分别测得14例老年人煤饼病人及11例正常人的 尿中17 酮类固醇排出量(mg/dl)如下,试比较 两组的均数有无差别
病人: 2.90 5.41 5.48 4.60 4.03 5.10 4.97 4.24 4.36 2.72 2.37 2.09 7.10 5.92
仅有原始数据时
ttest 变量名= #val
注: #val /*总体均数 命令可以加一些if和in条件限制
样本均数与总体均数比较的t检验
例4.4 命令:ttesti 25 5.1 0.88 4.6
样本均数与总体均数比较的t检验
例:某医生随机抽查10名某病患者的血红蛋白,如下:
病例号
1 23 4
---------+--------------------------------------------------------------------
x1 | 10 12.59 .5162794 1.632619 11.42209 13.75791
x2 | 10 13.27 .3415813 1.080175 12.49729 14.04271
Ha: mean != 14.02 Pr(|T| > |t|) = 0.0218
Ha: mean > 14.02 Pr(T > t) = 0.9891
配对设计t检验
ttest 变量1=变量2
变量1与变量2必须成对输入,样本含量必须相 等,如有缺省值,则用小数点表示,但与之对 应的记录在计算时被忽略
样本含量:10 均 数: 12.59 g/dl 标准差: 1.63 g/dl
Leabharlann Baidu
统计量与参数不同的两种可能
其一:抽样误差 (偶然的、随机的、较小的)
其二:本质上的差别 (必然的、大于随机误差)
样本均数与总体均数比较的t检验
ttest x=14.02
One-sample t test
------------------------------------------------------------------------------
data:ttest_2
配对设计t检验
. ttest x1=x2
Paired t test
------------------------------------------------------------------------------
Variable | Obs
Mean Std. Err. Std. Dev. [95% Conf. Interval]
------------------------------------------------------------------------------
mean = mean(x)
t = -2.7698
Ho: mean = 14.02
degrees of freedom = 9
Ha: mean < 14.02 Pr(T < t) = 0.0109
医学统计学实习课
STATA软件的应用(3) ——定量资料的统计分析
定量资料的统计分析
➢ 样本均数与总体均数比较的t 检验 ➢ 配对设计t 检验 ➢ 成组设计t 检验 ➢ 单因素方差分析 ➢ 两因素方差分析 ➢ 方差齐性检验 ➢ 正态性检验
样本均数与总体均数比较的t检验
仅有数据之统计量时
ttesti #obs #mean #sd #val
---------+--------------------------------------------------------------------
diff | 10 -.6799999 .5204272 1.645735 -1.857288 .4972881
------------------------------------------------------------------------------
Ha: mean(diff) > 0 Pr(T > t) = 0.8881
配对设计t检验
另一种思路: gen d=x0-x1 ttest d=0 结果一致
成组设计t检验
ttest 变量1=变量2, unpaired [unequal] ttest 变量,by(分组变量)[unequal] ttesti #obs1 #mean1 #sd1 #obs2 #mean2 #sd2
Variable | Obs
Mean Std. Err. Std. Dev. [95% Conf. Interval]
---------+--------------------------------------------------------------------
x | 10 12.59 .5162794 1.632619 11.42209 13.75791
配对设计t检验
例:10例矽肺患者经某药治疗,其血红蛋白 (g/dl)如下:
病例号:
1
2
3
4
5
6
7
8
9 10
治疗前(x1): 11.3 15.0 15.0 13.5 12.8 10.0 11.0 12.0 13.0 12.3 治疗后(x2): 14.0 13.8 14.0 13.5 13.5 12.0 14.7 11.4 13.8 12.0
5
6 7 8 9 10
血红蛋白 11.3 15 15 13.5 12.8 10 11 12 13 12.3 (x,g/dl)
问:该病患者的平均Hb含量是否与正常人的平均Hb含量 相同 (正常人的平均Hb含量为14.02(g/dl) )。
data:ttest_1
问题:
正常人
0=14.02
患者
1=?