Mercurial > repos > greg > insect_phenology_model

changeset 94:89492de287a6 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Fri, 01 Dec 2017 09:49:09 -0500
parents 06acadd38824
children 08571981adb3
files insect_phenology_model.R
diffstat 1 files changed, 7 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
--- a/insect_phenology_model.R	Fri Dec 01 09:43:56 2017 -0500
+++ b/insect_phenology_model.R	Fri Dec 01 09:49:09 2017 -0500
@@ -574,7 +574,7 @@
 
 # Subfigure 1: population size by life stage
 title <- paste(opt$insect, ": Total pop. by life stage :", opt$location, ": Lat:", latitude, ":", start_date, "-", end_date, sep=" ")
-plot(day.all, mean_value_adult, main=title, type="l", ylim=c(0, max(mean_value_egg + mean_value_egg.se, mean_value_nymph + mean_value_nymph.se, mean_value_adult + mean_value_adult.se)), axes=F, lwd=2, xlab="", ylab="", cex=3, cex.lab=3, cex.axis=3, cex.main=3)
+plot(day.all, mean_value_adult, main=title, type="l", ylim=c(0, max(mean_value_egg + mean_value_egg.std_error, mean_value_nymph + mean_value_nymph.std_error, mean_value_adult + mean_value_adult.std_error)), axes=F, lwd=2, xlab="", ylab="", cex=3, cex.lab=3, cex.axis=3, cex.main=3)
 # Young and old nymphs.
 lines(day.all, mean_value_nymph, lwd=2, lty=1, col=2)
 # Eggs
@@ -586,14 +586,14 @@
 if (opt$se_plot == 1) {
     # Add Standard error lines to plot
     # Standard error for adults
-    lines (day.all, mean_value_adult+mean_value_adult.se, lty=2)
-    lines (day.all, mean_value_adult-mean_value_adult.se, lty=2) 
+    lines (day.all, mean_value_adult+mean_value_adult.std_error, lty=2)
+    lines (day.all, mean_value_adult-mean_value_adult.std_error, lty=2) 
     # Standard error for nymphs
-    lines (day.all, mean_value_nymph+mean_value_nymph.se, col=2, lty=2)
-    lines (day.all, mean_value_nymph-mean_value_nymph.se, col=2, lty=2)
+    lines (day.all, mean_value_nymph+mean_value_nymph.std_error, col=2, lty=2)
+    lines (day.all, mean_value_nymph-mean_value_nymph.std_error, col=2, lty=2)
     # Standard error for eggs
-    lines (day.all, mean_value_egg+mean_value_egg.se, col=4, lty=2)
-    lines (day.all, mean_value_egg-mean_value_egg.se, col=4, lty=2)
+    lines (day.all, mean_value_egg+mean_value_egg.std_error, col=4, lty=2)
+    lines (day.all, mean_value_egg-mean_value_egg.std_error, col=4, lty=2)
 }
 
 # Subfigure 2: population size by generation