18/06/30 13:30:18.13 j2CU1Lw0.net
>>348
for loop無しで書ける。
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
tally <- function (x) any(cumsum(x)==cumsum(1-x))
mean(apply (mat,1,tally))
18/06/30 13:30:18.13 j2CU1Lw0.net
>>348
for loop無しで書ける。
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
tally <- function (x) any(cumsum(x)==cumsum(1-x))
mean(apply (mat,1,tally))