gistfile1.r

simulationOneTrial <- function(start.pop, knumGenerations, trial.index) {  
  df.allele <- NULL
  df.gen <- NULL
  
  #Keep track of the individuals in each generation
  df.gen <- rbind(df.gen, start.pop)
  
  x <- start.pop
  for ( gen in 1:knumGenerations) {
    #print(unlist(nex))
    #print(paste("Gen:", gen))
    numA <- calcAllelesInGeneration(x)
    df.gen <- rbind(df.gen, x)
    df.allele <- rbind(df.allele, c(gen, trial.index, numA))
    nex <- getNextGen(x)
    x <- sample(nex) #shuffle the population for breeding  
  }