Mercurial > repos > greg > insect_phenology_model

comparison insect_phenology_model.R @ 30:9146801c5ba1 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 12:02:37 -0500
parents 4d38c1789ae0
children e56529d02c3e
comparison
equal deleted inserted replaced
29:4d38c1789ae0 30:9146801c5ba1
576 g1a.se <- apply(g1a.rep, 1, sd) / sqrt(opt$replications) 576 g1a.se <- apply(g1a.rep, 1, sd) / sqrt(opt$replications)
577 # SE for F2 adult 577 # SE for F2 adult
578 g2a.se <- apply(g2a.rep, 1, sd) / sqrt(opt$replications) 578 g2a.se <- apply(g2a.rep, 1, sd) / sqrt(opt$replications)
579 579
580 dev.new(width=20, height=20) 580 dev.new(width=20, height=20)
581 sink("/dev/null")
581 582
582 # Start PDF device driver to save charts to output. 583 # Start PDF device driver to save charts to output.
583 pdf(file=opt$output, height=20, width=20, bg="white") 584 pdf(file=opt$output, height=20, width=20, bg="white")
584 585
585 par(mar = c(5, 6, 4, 4), mfrow=c(3, 1)) 586 par(mar = c(5, 6, 4, 4), mfrow=c(3, 1))
647 # SE for eggs 648 # SE for eggs
648 lines (day.all, g2a + g2a.se, col = 4, lty = 2) 649 lines (day.all, g2a + g2a.se, col = 4, lty = 2)
649 lines (day.all, g2a - g2a.se, col = 4, lty = 2) 650 lines (day.all, g2a - g2a.se, col = 4, lty = 2)
650 } 651 }
651 652
653 sink()
652 # Turn off device driver to flush output. 654 # Turn off device driver to flush output.
653 dev.off() 655 dev.off()