金币在哪里,呵呵!

qci=function(x=x1,q=.5,alpha=0.05)
{n=length(x)
m=floor(n*q);s1=pbinom(0:m,n,q);s2=pbinom(m:(n-1),n,q,low=F)
ss=c(s1,s2);nn=length(ss)
if (nn/2==floor(nn/2)) ss1=ss[1:nn/2]
else ss1=ss[1:(floor(nn/2)+1)]
n1=length(ss1);a=NULL;for(i in 1:n1)for (j in i:n1)
a=rbind(a,c(ss1[i]+ss1[j],1-(ss1[i]+ss1[j],i,j))
a=a[(a[,1]<1),];id=order(a[,1]); a=a[id,];ii=nrow(a)
if (a[1,1]>alpha)out="alpha is too small,CI=full range"
else
  {for (i in 1:ii) if (a[i,1]>alpha)
  {out=c(a[i-1,],x[a[i-1,3]],x[n-a[i-1,4]+1]);break}
   if (out[3]!=out[4])
     {out=c(out[1:2], out[4],out[3],
       x[out[4]],x[n-out[3]+1])
        if  ((out[6]-out[5])>(out[6]-out[5]))out=out1}
out=out[c(2,5:6)]}
out}

3楼是吴喜之的分位数区间估计代码吧

有命令吗

predict(fit,data.frame(x=10),intval="confident") 这能进行区间估计吗