作者pcshdan (我的人生就像冒险一样)
看板Statistics
标题[问题] 关於[R]的问题
时间Sun Jul 20 14:51:07 2008
我是R的新手
所以内容叙述不清楚请见谅
我想请问关於bootstrap resampling的问题
今天我有一组data 里面有1000笔资料
资料的尺度是intreval
请问我要如何抽出新的100组data (抽後放回)
让每一组data里面同样有1000笔资料的呢
又一组里面如果只要有500笔资料又要怎麽写呢
我看R的PDF档查了很久都没查到所以才来问问题
我有看到指令boot()和boot.ci()
但里面似乎都要写一个函数
我不懂到底是要如何使用?
试了很久都不行
因为我的最终目的是要形成一个新的矩阵
也就是那100组data我最後会写一个函数使每一组data都产生一个系数
然後做此系数的信赖区间
请问这样是不是使用loop会比较好
我是第一次学统计软体
问题叙述有点长
麻烦大家了
谢谢^^
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 219.91.78.147
※ 编辑: pcshdan 来自: 219.91.78.147 (07/20 14:55)
※ 编辑: pcshdan 来自: 219.91.78.147 (07/20 14:58)
1F：推 mrk:囧，有这个函数喔，我直接用runif取。Orz 07/20 18:02
2F：推 tew:用google搜寻 bootstrap R 前面几个网页就有范例 07/20 18:49
3F：推 tew:我常这样搜寻 学会自己要的东西 07/20 18:50
4F：推 waltervic:我也在学R= = 07/20 23:01
5F：→ clickhere:bootstrap要resampling 1000笔, 500是不对的. 07/21 02:23
6F：→ clickhere:把你写好的产生系数的函数丢给boot(),boot.ci()就行了. 07/21 02:24
7F：→ clickhere:sample的问题, boot()会处理. 07/21 02:25
8F：→ clickhere:用sample(..., replace=T)可以直接rexampling. 07/21 02:26
9F：→ pcshdan:感恩喔~谢谢各位的回答^^但我把产生系数的函数丢进boot() 07/21 18:08
10F：→ pcshdan:它一直说找不到引数"m" 不知道是哪步出问题@@ 07/21 18:09
11F：→ pcshdan:似乎在写boot()里面的函数跟平常的写法不一样 07/21 23:13