第六章 参数估计

6.1 以每天每千克体重52 mol 5－羟色胺处理家兔 14天后，对血液中血清素含量的影响如下表[9]：

对照组

5－羟色胺处理组

/（ g · L） s/（ g · L）

4.20

8.49

0.35 0.37

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

/cm

s/cm n