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