do uphi=p to 0.9999 by 0.00001; utailp=probbnml(uphi,n,m);
if abs(utailp-alpha/2)<0.00001 then goto upper; end;
upper:put m n p utailp uphi; proc print; id m;
var n p ltailp utailp lphi uphi;
title 'Confidence Limits for Binomial Population'; run;
结果如下: 项目(1):
Confidence Limits for Binomial Population
M N P LTAILP UTAILP LPHI 35 60 0.58333 0.024993 0.025006 0.44883
0.70931
UPHI
项目(2):
Confidence Limits for Binomial Population
M N P LTAILP UTAILP LPHI
UPHI
17 60 0.28333 0.024993 0.025002 0.1745 0.41443
项目(3):
UPHI 0.20505
Confidence Limits for Binomial Population
M N P LTAILP UTAILP LPHI 6 60 0.1 0.024996 0.025008 0.03759
6.4 乳腺癌患者有着沉重的心理负担,主要表现为:焦虑、怀疑和否认、恐惧、依赖、自私、悲观失望等。经心理护理后,在很多方面都到改善,护理前和护理后的评分(y?s)见下表[29]:
表 现 心理护理前/评分 心理护理后/评分 样本含量(n) 健康 25.71±2.31 28.51±9.04 62 抑郁 10.54±4.92 4.61±1.63 62 恐怖 12.10±2.24 10.76±3.08 62 强迫 15.29±4.53 10.33±4.83 62 焦虑 14.41±2.47 9.26±3.35 62 怀疑 16.75±3.69 4.92±2.11 62 神经衰弱 11.73±2.46 8.84±3.83 62
计算上述各种表现平均数差的0.95置信区间。(注意方差不具齐性的情况。)
答:结果见下表:
F FUTAILP ALPHA LCLDMSEQ UCLDMSEQ LCLDMSUN UCLDMSUN
15.3149 0 0.05 0.45423 5.14577 健康
0.43601 5.16399
9.11077 1.3323E-15 0.05 4.62695 7.23305
抑郁
4.61850 7.24150
1.89063 .0070120 0.05 0.38253 2.29747
恐怖 0.38162 2.29838
1.13684 0.30903 0.05 3.29519 6.62481
强迫
3.29512 6.62488
1.83948 .0093243 0.05 4.10360 6.19640
焦虑
4.10268 6.19732
3.05835 .000011068 0.05 10.7613 12.8987
怀疑
10.7586 12.9014
2.42397 .00034969 0.05 1.74559 4.03441
神经衰弱 1.74360 4.03640 表现
6.5 紫杉烷类物质具有优良的抗癌作用,近年来已成功地开发出紫杉烷类抗癌新药紫杉醇和多烯紫杉醇。由此也引起人们对药源植物云南红豆杉的关注,测定了紫杉烷类物质在不同类型云南红豆杉中的含量。下面给出其中的两种物质的测定结果[30]:
种类 紫杉醇 三尖杉宁碱 y/% s/% 0.006 2 0.005 1 0.005 0 0.004 5 203 203 n 计算两种物质平均数差的0.95置信区间,并以H0:μ1-μ2=0,HA:μ1-μ2 ≠ 0的假设推断两者间的差异显著性。
答:结果见下表:
Confidence Limits on the Difference of Means
for Non-Primal Data
F FUTAILP ALPHA LCLDMSEQ UCLDMSEQ
LCLDMSUN UCLDMSUN
1.28444 0.037999
0.05 .00026156 .0021384 .00026152 .0021385
F的显著性概率P=0.037 999,P >0.025,方差具齐性。方差具齐性时的0.95置信区间为:0.000 261 56 ~ 0.002 138 4。在置信区间内不包括0,因此紫杉醇和三尖杉宁碱的含量差异显著。
6.6 流行病学调查表明,高同型半胱氨酸(Homocysteine,Hcy)是导致动脉粥样硬化性血管病的一个新的独立危险因素。测定了脑梗死组和对照组的Hcy,结果(y?s)如下表[31]:
Hcy/(μ mol · L-1) 21.10±5.58 15.83±3.40
组 别 n 95 90
脑梗死 对 照
计算两组平均数差的0.95置信区间,并解释所计算的结果。
答:结果如下:
Confidence Limits on the Difference of Means
for Non-Primal Data
F FUTAILP ALPHA LCLDMSEQ
UCLDMSEQ LCLDMSUN UCLDMSUN
2.69346 .0000020008 0.05 3.92068 6.61932
3.93589 6.60411
可以很明显看出,方差是不具齐性的。0.95置信区间为:3.935 89 ~ 6.604 11。置信区间内不包含0,因此脑梗死病人的同型半胱氨酸显著高于对照组。高同型半胱氨酸很可能是动脉粥样硬化性血管病的危险因素之一。
6.7 30名受试者同时采取两份静脉血,分别用传统的魏氏法和自动血沉仪测定血沉[32],结果为分别:
y魏氏法-血沉仪法=0.2667(mm/h),
s魏氏法-血沉仪法=2.9935(mm/h)。在α = 0.05水平上,通过置信区间
检验两种方法的差异显著性。
答:所用程序如下:
options nodate; data esr;
input n mean std ; alpha=0.05;
talpha=-tinv(alpha/2,n-1); lclm=mean-talpha*std/sqrt(n); uclm=mean+talpha*std/sqrt(n); cards;
30 0.2667 2.9935 ;
proc print;
id n;
var mean std alpha lclm uclm ; title1 'Confidence Limits for Mu'; title2 'Sigma Is Unknown'; run;
结果见下表:
Confidence Limits for Mu
Sigma Is Unknown
N MEAN STD ALPHA LCLM
UCLM
30 0.2667 2.9935 0.05 -0.85109
1.38449
在置信区间内包含0,因此传统魏氏法和自动血沉仪法测得的结果差异不显著。
6.8 生长激素缺乏症的患儿,在用生长激素治疗前和治疗6个月后的身高和体重数据如下表[33]:
项 目 身高/cm 体重/kg
治疗前(y?s) 108±12 20.9±2.2
治疗后(y?s) 114±13 24.2±4.3
样本含量n 20 20
先用t检验,推断治疗前和治疗后的平均身高和平均体重在α = 0.05水平上的差异显著性,再用治疗前和治疗后的平均数差数的0.95置信区间验证。你认为这是一种很好的实验设计吗?怎样做检验的效果可能会更好?
答:1. 先做成组数据t检验: (1)身高:
T-Test for Non-Primal Data
F FUTAILP T DF
TUTAILP
1.17361 0.36536 1.51668 38.0000
0.068812
1.17361 0.36536 1.51668 37.7591 0.068838
(2)体重:
TUTAILP
T-Test for Non-Primal Data
F FUTAILP T DF 3.82025 .0026673 3.05542
38.0000 .0020482
3.82025 .0026673 3.05542
28.3091 .0024304
2. 计算置信区间: