生物统计学课后习题答案(杜荣骞第三版)

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

联系客服:779662525#qq.com(#替换为@)