[求助]如何用R软件求由delta法得出分位数的置信区间？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有一个问题，如何用R软件求由delta法得出分位数的置信区间？分位数的估计值x(p)我能算出了，用evir包中的qgev或qgpd得出。

>qgpd(0.95,xi,mu,sigma)

但是var( x(p) )用R软件算不出来，请求帮助。谢谢!

希望100金币奖励!

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Delta 置信区间 r软件 ELT del 软件 位数 置信区间 Delta

金币在哪里,呵呵!

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") 这能进行区间估计吗

那个嗨 发表于 2012-9-25 21:20
predict(fit,data.frame(x=10),intval="confident") 这能进行区间估计吗

有些小错误，其一，是interval 其二，是confidence(预测期望值），若预测实际值，用prediction

楼主的问题是一个函数在某一点的取值，这个不能算是统计推断吧。