18/06/30 13:18:23.33 j2CU1Lw0.net
>>341
このシミュレーションは少し難渋。
重複順列を返す関数があれば楽なのだが仕方ないので手作り。
a=3 ; b=2
(idx=combn(1:(a+b),a))
n=ncol(idx)
mat=matrix(0,nrow=n,ncol=a+b)
for(i in 1:n) mat[i,idx[,i]]=1
mat
tally <- function (x){
n=length (x)
t1=NULL
t0=NULL
for (i in 1:n){
if (x[i]){t1[i]=1;t0[i]=0}
else {t1[i]=0;t0[i]=1}
}
T1=cumsum (t1)
T0=cumsum (t0)
any(T1==T0)
}
mean(apply (mat,1,tally))