H0:?1??2HA:?1??2
首先检验数据的正态性,利用正态概率图发现分布有些正偏,经对数变换后,得到了较好的正态性。程序和结果如下:
options linesize=76 nodate; data wool;
input treat y @@; fat=log(y); cards;
1 0.19 1 0.18 1 0.21 1 0.30 1 0.66 1 0.42 1 0.08 1 0.12 1 0.30 1 0.27 2 0.15 2 0.13 2 0.07 2 0.24 2 0.19 2 0.08 2 0.20 ;
proc ttest; class treat; var fat;
title 'T-Test for Pooled Data'; run;
T-Test for Pooled Data
TTEST PROCEDURE
Variable: FAT
TREAT N Mean Std Dev Std Error
---------------------------------------------------------------------------- 1 10 -1.45824645 0.60395399 0.19098702
2 7 -1.97423071 0.46804661 0.17690499
Variances T DF Prob>|T| ---------------------------------------
Unequal 1.9820 14.8 0.0664 Equal 1.8913 15.0 0.0781
For H0: Variances are equal, F' = 1.67 DF = (9,6) Prob>F' = 0.5506
方差齐性检验的结果,方差具齐性。在方差具齐性时得到的t=1.891 3,其单侧显著性概率P=0.039 1,小于0.05,拒绝H0。因此,羊毛在处理前后的含脂率有显著性差异。
5.33 为了检验正常成人男、女性血液红细胞平均数之间的差异,随机抽取60对健康青年夫妻进行检测,得到以下结果:
夫/(104·
mm-3)
妻/(104·
mm-3)
d /(104· mm-3)(夫
-妻) d?42 2sd?6000
y1?465 s12?2040 y2?423 2s2?1800
n1= 60 n2= 60 nd= 60
使用适当的检验方法,检验其差异显著性。
答:夫妻之间是不能配成对子的,因此只能按成组数据处理。假定数据是从正态总体中获得的。 方差齐性检验的统计假设为:
H0:?1??2HA:?1??2
根据题意,本题之平均数差的显著性检验是双侧检验,统计假设为:
H0:?1??2HA:?1??2
结果如下:
T-Test for Non-Primal Data
F FUTAILP T DF
TUTAILP
1.13332 0.31615 5.24957 118.000 .00000034122
1.13332 0.31615 5.24957 117.541 .00000034296
方差是具齐性的。t=5.249 57,t的显著性概率P=0.000 000 341 22,远远小于0.005,因此拒绝H0。结论是男女两性血红细胞平均含量差异极显著。
5.34 随机选择11名25岁健康男性,抽取静脉血并制成血滤液。随机抽取其中1人作为对照,制成血滤液后马上测定血糖浓度,结果为102(10-2mg/mL)。另外10份滤液放置2小时后再测定,结果为:83、84、90、90、88、94、88、80、81、87(10-2mg/mL)。问酵解作用对血糖浓度是否有影响?
答:对照组无重复,无法计算s2,故不能进行统计推断。
第六章 参数估计
6.1 以每天每千克体重52 表[9]:
对照组 5-羟色胺处理组
mol 5-羟色胺处理家兔 14天后,对血液中血清素含量的影响如下
y/(g · L-1) 4.20 8.49
g · L-1) 0.35 0.37
s/(n 12 9
建立对照组和5-羟色胺处理组平均数差的0.95置信限。
答:程序如下:
options nodate; data common;
alpha=0.05;
input n1 m1 s1 n2 m2 s2; dfa=n1-1; dfb=n2-1; vara=s1**2; varb=s2**2;
if vara>varb then F=vara/varb; else F=varb/vara;
if vara>varb then Futailp=1-probf(F,dfa,dfb); else Futailp=1-probf(F,dfb,dfa); df=n1+n2-2; t=tinv(1-alpha/2,df); d=abs(m1-m2);
lcldmseq=d-t*sqrt(((dfa*vara+dfb*varb)/(dfa+dfb))*(1/n1+1/n2)); ucldmseq=d+t*sqrt(((dfa*vara+dfb*varb)/(dfa+dfb))*(1/n1+1/n2));
k=vara/n1/(vara/n1+varb/n2); df0=1/(k**2/dfa+(1-K)**2/dfb); t0=tinv(1-alpha/2,df0); lcldmsun=d-t0*sqrt(vara/n1+varb/n2); ucldmsun=d+t0*sqrt(vara/n1+varb/n2); cards;
12 4.20 0.35 9 8.49 0.37 ;
proc print; id f;
var Futailp alpha lcldmseq ucldmseq lcldmsun ucldmsun; title1 'Confidence Limits on the Difference of Means'; title2 'for Non-Primal Data'; run;
结果见下表:
Confidence Limits on the Difference of Means
for Non-Primal Data
F FUTAILP ALPHA LCLDMSEQ UCLDMSEQ
LCLDMSUN UCLDMSUN
1.11755 0.42066 0.05 3.95907 4.62093 3.95336
4.62664
首先,方差是具齐性的。在方差具齐性的情况下,平均数差的0.95置信下限为3.959 07,置信上限为4.620 93。0.95置信区间为3.959 07 ~ 4.620 93。
6.2 不同年龄的雄岩羊角角基端距如下表[27]:
年龄/a 4.5 5.5 y/cm 28.92 31.81 s/cm 2.17 2.44 13 11 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.26433 0.34528 0.05 0.93873 4.84127 0.90910
4.87090
因为方差具齐性,所以平均数差的0.95置信区间为:0.938 73 ~ 4.841 27。置信区间内不包括0,因此两者间的差异是显著的。
6.3 了解我国风险识别、风险评价和风险缓解的现状,对于应对突发事件有重要作用。以下是关于应对突发公共卫生事件能力调查(共调查了60个单位)的部分数据[28]:
项 目
识别了当地可能发生的突发公共卫生事件 对所识别的突发公共卫生事件进行了风险评价 根据风险评价结果确定了当地突发公共卫生事件的工作重点
分别计算上述三个项目的0.95置信区间。
答:程序如下:
options nodate; data clbi;
n=60; m=35; p=m/n; alpha=0.05; do lphi=0.0001 to p by 0.00001; ltailp=1-probbnml(lphi,n,m-1);
if abs(ltailp-alpha/2)<0.00001 then goto lower; end;
lower:put m n p ltailp lphi;
单位数 35 17
6
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置信区间。(注意方差不具齐性的情况。)