作者celestialgod (攸蓝)
看板Statistics
标题Re: [问题] 统计问题-如何重抽
时间Sun Dec 29 03:27:05 2013
用R吧!!
简单的想法，瞧瞧可不可以用
pop_build = function(n.v){
Reduce(c, sapply(1:length(n.v), function(i) rep(LETTERS[i], n.v[i])))
}
n.v = c(6,3,7,4,10)
populations = pop_build(n.v)
samplesize = 15
rep_N = 10000
samples = replicate(rep_N, sample(populations, samplesize, rep = FALSE))
samples = t(sapply(apply(samples, 2, table), function(v){
out = rep(0, length(n.v))
out[match(names(v), LETTERS[1:length(n.v)])] = v
out
}))
table(apply(samples,1, paste, collapse = ","))
apply(samples, 2, mean)
apply(samples, 2, sd)
※ 引述《buzzmonkey (笑笑羊)》之铭言：
: 如果是跟统计软体有关请重发文章
: 如果跟论文有关也烦请您重发文章
: 文章类别是为了帮助大家搜寻资料与解答，造成不便之处请见谅
: 想请教一下群聚重抽，有什麽样的软体可以和要如何做
: 例如我有ABCDE五个物种，右边分别为所观察的到只数加起来共30只
: A 6
: B 3
: C 7
: D 4
: E 10
: 现在要进行重抽看如果是只剩15只的情况，那分别为哪几个物种各为几只?
: 谢谢大家
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 218.164.75.46
1F：推 anovachen:感谢分享!! 这是比较有效率的写法。 12/29 16:02
2F：→ buzzmonkey:感谢分享>< 12/29 17:12
3F：推 buzzmonkey:感谢分享>< 12/29 17:19
※ 编辑: celestialgod 来自: 140.116.152.221 (12/30 14:33)