複数座位の線形和効果

k<-50
Ns<-10000
a<-array(0,c(Ns,2,k))
f<-runif(k)
for(i in 1:k){
	a[,,i]<-sample(0:1,Ns*2,replace=TRUE,prob=c(f[i],1-f[i]))
}

d<-a[,1,]+a[,2,]
d<-apply(d,1,cumsum)

counter<-10*k+1

for(i in 1:k){
	for(j in 1:10){
		
		lab=paste("No. markers = ",i,sep="")
		filename=paste("test",counter,".png",sep="")
		counter<-counter+1
		png(filename)
		#hist(d[i,],main=lab,xlab="",ylab="",axes=FALSE,breaks=0:(2*i))
		plot(density(d[i,]),axes=FALSE,main=lab,xlab="",ylab="")
		#par(new=TRUE)
		dev.off()
	}
}