Mercurial > repos > greg > insect_phenology_model

comparison insect_phenology_model.R @ 27:24949e72f7ec draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Thu, 09 Nov 2017 11:35:44 -0500
parents 3c06cab3db2c
children 3ef94dce7454
comparison
equal deleted inserted replaced
26:3c06cab3db2c 27:24949e72f7ec
213 cat("Adult mortality rate: ", opt$adult_mort, "\n") 213 cat("Adult mortality rate: ", opt$adult_mort, "\n")
214 cat("Min clutch size: ", opt$min_clutch_size, ", ") 214 cat("Min clutch size: ", opt$min_clutch_size, ", ")
215 cat("Max clutch size: ", opt$max_clutch_size, "\n") 215 cat("Max clutch size: ", opt$max_clutch_size, "\n")
216 cat("(egg->young nymph): ", opt$young_nymph_accum, ", ") 216 cat("(egg->young nymph): ", opt$young_nymph_accum, ", ")
217 cat("(young nymph->old nymph): ", opt$old_nymph_accum, "\n") 217 cat("(young nymph->old nymph): ", opt$old_nymph_accum, "\n")
218 cat("(old nymph->adult): ", opt$adult_accum 218 cat("(old nymph->adult): ", opt$adult_accum)
219 219
220 # Initialize matrix for results from all replications 220 # Initialize matrix for results from all replications.
221 S0.rep <- S1.rep <- S2.rep <- S3.rep <- S4.rep <- S5.rep <- matrix(rep(0, opt$num_days * opt$replications), ncol = opt$replications) 221 S0.rep <- S1.rep <- S2.rep <- S3.rep <- S4.rep <- S5.rep <- matrix(rep(0, opt$num_days * opt$replications), ncol = opt$replications)
222 newborn.rep <- death.rep <- adult.rep <- pop.rep <- g0.rep <- g1.rep <- g2.rep <- g0a.rep <- g1a.rep <- g2a.rep <- matrix(rep(0, opt$num_days * opt$replications), ncol=opt$replications) 222 newborn.rep <- death.rep <- adult.rep <- pop.rep <- g0.rep <- g1.rep <- g2.rep <- g0a.rep <- g1a.rep <- g2a.rep <- matrix(rep(0, opt$num_days * opt$replications), ncol=opt$replications)
223 223
224 # loop through replications 224 # loop through replications
225 for (N.rep in 1:opt$replications) { 225 for (N.rep in 1:opt$replications) {