作者phil5566 (5566)
看板R_Language
标题[问题] 函数boot.ci的问题?
时间Wed Dec 21 18:03:27 2016
[问题类型]:
程式谘询: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
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.74.231
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1482314610.A.408.html
1F：→ andrew43: 直接把boot.obj$t拿出来做，例如 12/21 21:24
2F：→ andrew43: apply(boot.out$t, 2, quantile, c(0.025, 0.975)) 12/21 21:24
3F：→ phil5566: 谢谢a大,试过您提供的方法,有成功,又学到一招,感谢 12/23 19:33
4F：→ phil5566: 不知道是不是boot.ci这个函数一次只能算第一列的C.I,还 12/23 19:36
5F：→ phil5566: 是要作什麽调整之後,第二列2~第5列的C.I也可以计算 12/23 19:38
6F：→ andrew43: 我也没试出来才这样建议你的。 12/23 19:39
7F：推 andrew43: 反正最麻烦的部分boot()已经帮你做完了，要采用那种分 12/23 19:42
8F：推 andrew43: 数都只是简单的运算，自己写在apply里即可。 12/23 19:42
9F：→ phil5566: 了解了,感谢 12/23 19:51