ls()
source("runms.R")
source("runms.R")
source("runms.R")
save(stat.voight, stat.3pops.sim, par.italy.sim, models, file="human.rda", compress=T)
q()
