ÃÑ 327ÆäÀÌÁö
240ÆäÀÌÁö º»¹®½ÃÀÛ
¥±. À§ÇèºÐ¼®ÀÇ Åë°èÀû ÀÌ·Ð
install.packages('mc2d')
library(mc2d)
i
j
ij
v1
A
#º¯¼öpa¸¦±¸ÇϱâÀ§Çѹݺ¹¹®
for(nin1:7){
for(kin1:100){
a
A[k,ij[2,n]]
#pa°ªÁßÇϳªÀÇ°ªÃßÃâ
pa
Nii
B
# Nii¸¦±¸ÇϱâÀ§Çѹݺ¹¹®
for(nin1:100){Ni
if(Ni==0)
Nti=0
elseNti=rhyper(1,Ni,50-Ni,2)
if(Nti==0)
Np
elseNp=rbinom(1,Nti,0.85)
Ps=rpert(1,0.2,0.4,0.5)
Pd=rpert(1,0.1,0.3,0.4)
if(Ni==0)
Nin=0
elseNin=rbinom(1,Ni,Pd)
if(Nin+Np==0)
Nif
elseNif=0
B[i,]
#È÷½ºÅä±×·¥±×¸®±â
hist(B$Nii)
±×¸² 3-33 100°³ ´ß´Ù¸®·Î ±¸¼ºµÈ 1°³ ¼öÀÔÆÐŶ¿¡¼ °¨¿°µÈ ´ß ´Ù¸® ¦ÀÇ ¼ö
(É®??)¿¡ ´ëÇÑ È÷½ºÅä±×·¥
225
install.packages('mc2d')
library(mc2d)
i
j
ij
v1
A
#º¯¼öpa¸¦±¸ÇϱâÀ§Çѹݺ¹¹®
for(nin1:7){
for(kin1:100){
a
A[k,ij[2,n]]
#pa°ªÁßÇϳªÀÇ°ªÃßÃâ
pa
Nii
B
# Nii¸¦±¸ÇϱâÀ§Çѹݺ¹¹®
for(nin1:100){Ni
if(Ni==0)
Nti=0
elseNti=rhyper(1,Ni,50-Ni,2)
if(Nti==0)
Np
elseNp=rbinom(1,Nti,0.85)
Ps=rpert(1,0.2,0.4,0.5)
Pd=rpert(1,0.1,0.3,0.4)
if(Ni==0)
Nin=0
elseNin=rbinom(1,Ni,Pd)
if(Nin+Np==0)
Nif
elseNif=0
B[i,]
#È÷½ºÅä±×·¥±×¸®±â
hist(B$Nii)
±×¸² 3-33 100°³ ´ß´Ù¸®·Î ±¸¼ºµÈ 1°³ ¼öÀÔÆÐŶ¿¡¼ °¨¿°µÈ ´ß ´Ù¸® ¦ÀÇ ¼ö
(É®??)¿¡ ´ëÇÑ È÷½ºÅä±×·¥
225
240ÆäÀÌÁö º»¹®³¡
¸Þ´º
ÇöÀç Æ÷Ä¿½ºÀÇ ¾Æ·¡³»¿ëµéÀº µ¿ÀÏÇÑ ÄÁÅÙÃ÷¸¦ °¡Áö°í ÆäÀÌÁö³Ñ±è È¿°ú¹× ½Ã°¢Àû È¿°ú¸¦ Á¦°øÇÏ´Â ÆäÀÌÁöÀ̹ǷΠ½ºÅ©¸°¸®´õ »ç¿ëÀÚ´Â ¿©±â±îÁö¸¸ ³¶µ¶ÇϽðí À§ÀÇ ÆäÀÌÁöÀ̵¿ ¸µÅ©¸¦ »ç¿ëÇÏ¿© ´ÙÀ½ÆäÀÌÁö·Î À̵¿ÇϽñ⠹ٶø´Ï´Ù.