# HG changeset patch
# User greg
# Date 1512139749 18000
# Node ID 89492de287a6309007b8c236b90bfa5a2bc447a6
# Parent  06acadd388242ebd16845be2ca325f6f803a4441
Uploaded

diff -r 06acadd38824 -r 89492de287a6 insect_phenology_model.R
--- 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