回复 第1楼 的 zhangchaoting:可以从数学上先考虑下:
<br />
x=matrix(rbinom(338*21,1,.5),338)<br />
ind=rowSums(x)<br />
y=ind[ind>=10]<br />
tt='y'<br />
for(i in 1:9){<br />
assign(paste('y',i,sep=''),y-i)<br />
tt=paste(tt,paste('y',i,sep=''),sep='*')<br />
}<br />
sum(eval(parse(text=tt)))/388/gamma(11)<br />
不知道是不是你的意思
</p>