作者granzi (乌木)
看板Statistics
标题Re: [R] 想问如何用R将模拟後资料分组
时间Tue Apr 8 00:08:56 2008
※ 引述《freefeel (自由感觉)》之铭言:
: 不知道在这里问合不合适..
: 我想问的是
: 我最近用R模拟一组资料,想将所得的结果分组
: 但是如果所分到该组的样本数若不够的话.再继续跑回圈
: 不过想归想... 实在是不太知道程式该怎麽写
: 举例:
: 假设 X ~ Bin(n=50, p =1/2)
: 用rbinom随机先产生一笔x资料
: 将所得结果分组
: 若x范围介於 [0~10],则分到第一组;否则舍弃(x落在范围外),直到得到10笔样本资料
: 接着同样的步骤继续抽样
: 只是x范围改变
: 若x范围介於 [0~25]分到第二组;否则舍弃,直到得到10笔样本资料
: 若x范围介於 [0~50]分到第三组;否则舍弃,直到得到10笔样本资料。
: 我想得到三组不同范围的 x ,各10笔资料... 可是不太知道程式该如何写 ><
: 不知道板上有人能帮帮小弟吗
: 谢谢 ^^
x <- rbino(50, p=0.5)
y <- array(NA, c(10, 3))
y[,1] <- (x[((x <= 10) & (x >= 0)])[1:10]
y[,2] <- (x[((x <= 25) & ( x >= 0)])[1:10]
y[,3] <- (x[((x <= 50) & ( x >= 0)])[1:10]
即为所求
原PO请确你的分组范围是否正确,这样很有可能得到三组一模一样的list
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 220.135.68.240
1F:推 freefeel:thx :p 07/01 15:19