Mercurial > repos > greg > insect_phenology_model

comparison insect_phenology_model.R @ 4:4114785abbe7 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Tue, 08 Aug 2017 12:38:07 -0400
parents f20a67192cfa
children 80d38a26b2e8
comparison
equal deleted inserted replaced
3:574a052becbe 4:4114785abbe7
591 pdf(file=opt$output, height=20, width=20, bg="white") 591 pdf(file=opt$output, height=20, width=20, bg="white")
592 592
593 par(mar = c(5, 6, 4, 4), mfrow=c(3, 1)) 593 par(mar = c(5, 6, 4, 4), mfrow=c(3, 1))
594 594
595 # Subfigure 2: population size by life stage 595 # Subfigure 2: population size by life stage
596 plot(day.all, sa, main = "Total Population Size by Life Stage", type = "l", ylim = c(0, max(se + se.se, sn + sn.se, sa + sa.se)), axes = F, lwd = 2, xlab = "", ylab = "Number", cex = 2, cex.lab = 2, cex.axis = 2, cex.main = 2) 596 plot(day.all, sa, main = "BSMB Total Population Size by Life Stage", type = "l", ylim = c(0, max(se + se.se, sn + sn.se, sa + sa.se)), axes = F, lwd = 2, xlab = "", ylab = "Number", cex = 2, cex.lab = 2, cex.axis = 2, cex.main = 2)
597 # Young and old nymphs 597 # Young and old nymphs
598 lines(day.all, sn, lwd = 2, lty = 1, col = 2) 598 lines(day.all, sn, lwd = 2, lty = 1, col = 2)
599 # Eggs 599 # Eggs
600 lines(day.all, se, lwd = 2, lty = 1, col = 4) 600 lines(day.all, se, lwd = 2, lty = 1, col = 4)
601 axis(1, at = c(1:12) * 30 - 15, cex.axis = 2, labels = c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec")) 601 axis(1, at = c(1:12) * 30 - 15, cex.axis = 2, labels = c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec"))
614 lines (day.all, se + se.se, col = 4, lty = 2) 614 lines (day.all, se + se.se, col = 4, lty = 2)
615 lines (day.all, se - se.se, col = 4, lty = 2) 615 lines (day.all, se - se.se, col = 4, lty = 2)
616 } 616 }
617 617
618 # Subfigure 3: population size by generation 618 # Subfigure 3: population size by generation
619 plot(day.all, g0, main = "Total Population Size by Generation", type = "l", ylim = c(0, max(g2)), axes = F, lwd = 2, xlab = "", ylab = "Number", cex = 2, cex.lab = 2, cex.axis = 2, cex.main = 2) 619 plot(day.all, g0, main = "BSMB Total Population Size by Generation", type = "l", ylim = c(0, max(g2)), axes = F, lwd = 2, xlab = "", ylab = "Number", cex = 2, cex.lab = 2, cex.axis = 2, cex.main = 2)
620 lines(day.all, g1, lwd = 2, lty = 1, col = 2) 620 lines(day.all, g1, lwd = 2, lty = 1, col = 2)
621 lines(day.all, g2, lwd = 2, lty = 1, col = 4) 621 lines(day.all, g2, lwd = 2, lty = 1, col = 4)
622 axis(1, at = c(1:12) * 30 - 15, cex.axis = 2, labels = c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec")) 622 axis(1, at = c(1:12) * 30 - 15, cex.axis = 2, labels = c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec"))
623 axis(2, cex.axis = 2) 623 axis(2, cex.axis = 2)
624 leg.text <- c("P", "F1", "F2") 624 leg.text <- c("P", "F1", "F2")
635 lines (day.all, g2 + g2.se, col = 4, lty = 2) 635 lines (day.all, g2 + g2.se, col = 4, lty = 2)
636 lines (day.all, g2 - g2.se, col = 4, lty = 2) 636 lines (day.all, g2 - g2.se, col = 4, lty = 2)
637 } 637 }
638 638
639 # Subfigure 4: adult population size by generation 639 # Subfigure 4: adult population size by generation
640 plot(day.all, g0a, ylim = c(0, max(g2a) + 100), main = "Adult Population Size by Generation", type = "l", axes = F, lwd = 2, xlab = "Year", ylab = "Number", cex = 2, cex.lab = 2, cex.axis = 2, cex.main = 2) 640 plot(day.all, g0a, ylim = c(0, max(g2a) + 100), main = "BSMB Adult Population Size by Generation", type = "l", axes = F, lwd = 2, xlab = "Year", ylab = "Number", cex = 2, cex.lab = 2, cex.axis = 2, cex.main = 2)
641 lines(day.all, g1a, lwd = 2, lty = 1, col = 2) 641 lines(day.all, g1a, lwd = 2, lty = 1, col = 2)
642 lines(day.all, g2a, lwd = 2, lty = 1, col = 4) 642 lines(day.all, g2a, lwd = 2, lty = 1, col = 4)
643 axis(1, at = c(1:12) * 30 - 15, cex.axis = 2, labels = c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec")) 643 axis(1, at = c(1:12) * 30 - 15, cex.axis = 2, labels = c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec"))
644 axis(2, cex.axis = 2) 644 axis(2, cex.axis = 2)
645 leg.text <- c("P", "F1", "F2") 645 leg.text <- c("P", "F1", "F2")