複数座位の線形和効果
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() } }