Mercurial > repos > greg > bmsb

changeset 11:99a386ac1f5b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Tue, 16 Aug 2016 09:52:16 -0400
parents 303a662c7881
children 7b0c47f9436a
files bmsb.R
diffstat 1 files changed, 5 insertions(+), 11 deletions(-) [+]
line wrap: on
line diff
--- a/bmsb.R	Fri Aug 12 14:45:12 2016 -0400
+++ b/bmsb.R	Tue Aug 16 09:52:16 2016 -0400
@@ -1,13 +1,7 @@
 #!/usr/bin/env Rscript
 
-
-suppressPackageStartupMessages(require("hash"))
-suppressPackageStartupMessages(require("KernSmooth"))
-suppressPackageStartupMessages(require("methods"))
-suppressPackageStartupMessages(library("optparse"))
-suppressPackageStartupMessages(library("rjson"))
-suppressPackageStartupMessages(require("utils"))
-
+#suppressPackageStartupMessages(library("optparse"))
+#suppressPackageStartupMessages(library("rjson"))
 
 option_list <- list(
     make_option(c("-s", "--save_log"), action="store_true", default=FALSE, help="Save R logs"),
@@ -150,7 +144,7 @@
 
 
 # model initialization
-# TODO: setwd(“/home/lunarmouse/Nielsen's project/”)
+# TODO: add tool params for the following options.
 # start with 1000 individuals
 n <- 1000
 # Generation, Stage, DD, T, Diapause
@@ -436,9 +430,9 @@
         N.death[day] <- n.death
         N.adult[day] <- n.adult
     }   
-print(n.sim)
+#print(n.sim)
 }
 
 proc.time() - ptm
 dd.cum <- cumsum(dd.day)
-save(dd.day, dd.cum, S0, S1, S2, S3, S4, S5, N.newborn, N.death, N.adult, tot.pop, gen0.pop, gen1.pop, gen2.pop, g0.adult, g1.adult, g2.adult, file="asheville2014sim.Rdat")
+save(dd.day, dd.cum, S0, S1, S2, S3, S4, S5, N.newborn, N.death, N.adult, tot.pop, gen0.pop, gen1.pop, gen2.pop, g0.adult, g1.adult, g2.adult, file=opt$output)