R code for moose popn model

sessionInfo() 

## R version 3.5.2 (2018-12-20)
## Platform: x86_64-w64-mingw32/x64 (64-bit)
## Running under: Windows 10 x64 (build 17763)
## 
## Matrix products: default
## 
## locale:
## [1] LC_COLLATE=English_United States.1252 
## [2] LC_CTYPE=English_United States.1252   
## [3] LC_MONETARY=English_United States.1252
## [4] LC_NUMERIC=C                          
## [5] LC_TIME=English_United States.1252    
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## loaded via a namespace (and not attached):
##  [1] compiler_3.5.2  magrittr_1.5    tools_3.5.2     htmltools_0.3.6
##  [5] yaml_2.2.0      Rcpp_1.0.0      stringi_1.3.1   rmarkdown_1.12 
##  [9] knitr_1.22      stringr_1.4.0   xfun_0.6        digest_0.6.18  
## [13] evaluate_0.13

library(population)

## Loading required package: parallel

## Loading required package: abind

n0<-4350 #2014 moose population estimate from DelGiudice (2014)
s0<-0.85 #median adult moose survival from DNR collaring study (Carstensen et al. 2016)
b1<-1.12 #median litter size at birth from DNR collaring study (Severud et al. 2016)
years<-30
runs<-1000
results <- project(
  years = years,
  runs = runs,
  initial_population = c(n0),
  survival = cbind(c(s0), 0.03), #stochasticity from standard deviation of adult survival from DNR collaring study (Carstensen et al. 2016)
  litter = cbind(c(b1), 0.32) #stochasticity from standard deviation of litter size at birth from DNR collaring study (Severud et al. 2016)
)

## 
## |**************************************************|

# Plot projection
plot_projection(results, "mean")