→ celestialgod:如果需要更進一步的程式諮詢 可能需要多一點的資訊 06/23 14:11
: N = 500
這些變數應該是你要儲存的資訊
長度應是你大迴圈的長度吧？
L0 = NULL; length(L0) = N
L1 = NULL; length(L1) = N
L2 = NULL; length(L2) = N
Alpha1 = NULL; length(alpha1) = N
Alpha2 = NULL; length(alpha2) = N
: Z1 = matrix(,ncol=2,nrow=length(t1))
: Z2 = matrix(,ncol=2,nrow=length(t2))
: for(i in 1:length(t1)){
: Z1[i,] = rmultinom(1,1,c(g1(t1[i]),(1-g1(t1[i]))))
: }
: for(i in 1:length(t2)){
: Z2[i,] = rmultinom(1,1,c(g2(t2[i]),(1-g2(t2[i]))))
: }
: }
這裡可以先做簡單的修改
Z1 = sapply(1:length(t1), function(k) rmultinom(1,1,c(g1(t1[k]),(1-g1(t1[k])))))
Z2 = sapply(1:length(t2), function(k) rmultinom(1,1,c(g2(t2[k]),(1-g2(t2[k]))))
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 36.239.245.94