Archive/NIMBioS.code/fastsimcoal/fsc.scenarioCheck.R

In christianparobek/skeleSim: Genetic Simulation Engine

fsc.scenarioCheck <- function(params) {
  # check that sample times and growth rates are of length number of populations, and
  #   that historical events matrix converges 
  results <- sapply(params@scenarios, function(sc) {
    
    c(
      samptime.eq.npops = length(sc@simulator.params@sample.times) == sc@num.pops,
      growth.eq.npops = length(sc@simulator.params@growth.rate) == sc@num.pops,
      hist.ev.good = fsc.histEvConverges(
                            sc@simulator.params@hist.ev,
                            sc@pop.size,
                            sc@simulator.params@growth.rate,
                            length(sc@migration))
    )
  })
  return(results)
}