averageFixationTime <- c() for (N in 2:50) { print(N) time.fixation <- c() for (i in 1:5000) { population <- .5 generation <- 0 while (!(population %in% c(0,1))) { children <- rbinom(1,N,population) population <- children/N generation <- generation+1 } time.fixation <- c(time.fixation,generation) } averageFixationTime <- c(averageFixationTime, mean(time.fixation)) } #plot(fA,ylim=c(0,N),xlab='generation',ylab='frequeny of A',type='l') #hist(time.fixation,nclass=50)