TaskShare 发表于 2012-10-21 11:25 
我不是很明白你的意思。如你所说,已知lambda是个常数=2,那么一个常数如何会有什么置信区间呢?所以,我想 ...
谢谢你了,我最早也是理解错了,是用另外一个值去估算lambda。在问你一下,我在算t confidence interval你能帮我看一下我这个code哪里写错了吗
%let lambda = 2;
data exp;
do i=1 to 20;
E=ranexp(2345)/λ
output;
end;
run;
%macro analyze(data=,out=);
proc means noprint data=&data vardef=n;
output out=&out(drop=_freq_ _type_)
mean=mean_i mean_E;
var i E;
%bystmt;
run;
data &out;
set &out;
lambda1 =1/mean_E;
drop mean_i;
run;
%mend;
title2 "BC confidence interval"
%boot(data=exp,stat=lambda1,samples=200,random=123,alpha=0.10);
前半部分应该没错,后半部分
%macro analyze(data=,out=);
proc means noprint data=&data vardef=n;
output out=&out(drop=_freq_ _type_)
mean= mean_E; kurtosis=kurt_E;
var i E;
%bystmt;
run;
data &out;
set &out;
lambda1=1/mean_E
kurtl=1/kurt_E
studl=lambda1*sqrt(kurtl+2);
RUN;
%mend;
title2 'T Confidence Interval';
%bootci(t,stat=lambda1,student=studl,alpha=.10)