推 anovachen:感謝分享!! 這是比較有效率的寫法。 12/29 16:02
→ buzzmonkey:感謝分享>< 12/29 17:12
推 buzzmonkey:感謝分享>< 12/29 17:19
※ 編輯: celestialgod 來自: 140.116.152.221 (12/30 14:33)
用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