[问题] 函数boot.ci的问题?

楼主: phil5566 (5566)   2016-12-21 18:03:27
[问题类型]:
程式咨询:boot这个package里的boot.ci这个function如何一次
产生多个参数的confidence intervals
[软件熟悉度]:
新手(没写过程式)
[问题叙述]:
由于程式码太长,因此只附上关键部分的问题,
我先用boot()这个function把data T bootstrap 100次
并且把每bootstrap一次的结果执行
自定函数"theta.boot"这个运算
再用boot.ci()这个function对boot.obj(boot()这个function输出的结果)
运算出"basic", "norm", "perc"三种类型的confidence intervals
程式码如下:
boot.obj <- boot(T, statistic = theta.boot, R = 100)
boot.ci(boot.out = boot.obj, type = c("basic", "norm", "perc"))
结果如下:
Bootstrap Statistics :
original bias std. error
t1* -0.0972177 -0.03674216 0.1914715
t2* 1.3326862 0.02823161 0.7443144
t3* 4.8126388 0.11100370 0.3723238
t4* 1.0862814 -0.04447950 0.1474000
t5* 4.7145885 0.04016033 0.2725067
Intervals :
Level Normal Basic Percentile
95% (-0.4358, 0.3148 ) (-0.3816, 0.1833 ) (-0.3777, 0.1871 )
问题点就是:
算confidence intervals怎么只输出t1那一列的confidence intervals
t2,t3,t4,t5的confidence intervals到底要怎么产生?
求高手解救,谢谢!
[程式范例]:
http://codepad.org/aKtpe7FB
[环境叙述]:
win7,
R x64 3.3.1,
"boot"package
[关键字]:
boot,boot.ci
作者: andrew43 (讨厌有好心推文后删文者)   2016-12-21 21:24:00
直接把boot.obj$t拿出来做,例如apply(boot.out$t, 2, quantile, c(0.025, 0.975))
楼主: phil5566 (5566)   2016-12-23 19:33:00
谢谢a大,试过您提供的方法,有成功,又学到一招,感谢不知道是不是boot.ci这个函数一次只能算第一列的C.I,还是要作什么调整之后,第二列2~第5列的C.I也可以计算
作者: andrew43 (讨厌有好心推文后删文者)   2016-12-23 19:39:00
我也没试出来才这样建议你的。反正最麻烦的部分boot()已经帮你做完了,要采用那种分数都只是简单的运算,自己写在apply里即可。
楼主: phil5566 (5566)   2016-12-23 19:51:00
了解了,感谢

Links booklink

Contact Us: admin [ a t ] ucptt.com