[Pomp-commits] r1022 - pkg/pompExamples/inst/examples

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Wed Dec 17 20:12:15 CET 2014


Author: kingaa
Date: 2014-12-17 20:12:15 +0100 (Wed, 17 Dec 2014)
New Revision: 1022

Modified:
   pkg/pompExamples/inst/examples/parus.R
Log:
- include better documentation of Parus data

Modified: pkg/pompExamples/inst/examples/parus.R
===================================================================
--- pkg/pompExamples/inst/examples/parus.R	2014-12-17 19:12:10 UTC (rev 1021)
+++ pkg/pompExamples/inst/examples/parus.R	2014-12-17 19:12:15 UTC (rev 1022)
@@ -1,6 +1,31 @@
+#' Annual *Parus major* counts in Wytham Wood, Oxfordshire, England.
+#' 
+#' @param proc  the process model (Gompertz or Ricker are currently supported)
+#' @param meas  the measurement model
+#' (lognormal, negative binomial, and Poisson are currently supported)
+#' @references McCleery, R. & Perrins, C. (1991)
+#' Effects of predation on the numbers of Great Tits, Parus major.
+#' In: Bird Population Studies, 
+#' edited by Perrins, C.M., Lebreton, J.-D. & Hirons, G.J.M. 
+#' Oxford. Univ. Press. pp. 129--147.
+#' @author Aaron A. King \email{kingaa@@umich.edu}
+
 require(pomp)
 
-dat <- 'year,pop
+dat <- '
+## Parus major (Great Tit) census (all individuals)
+## Wytham Wood, Oxfordshire
+## Global Population Dynamics Database dataset #10163.
+## (NERC Centre for Population Biology, Imperial College (2010) 
+## The Global Population Dynamics Database Version 2. 
+## http://www.sw.ic.ac.uk/cpb/cpb/gpdd.html").
+## Original source: 
+## McCleery, R. & Perrins, C. (1991)
+## Effects of predation on the numbers of Great Tits, Parus major.
+## In: Bird Population Studies, 
+## edited by Perrins, C.M., Lebreton, J.-D. & Hirons, G.J.M. 
+## Oxford. Univ. Press. pp. 129--147.
+year,pop
 1960,148
 1961,258
 1962,185
@@ -30,7 +55,7 @@
 1986,211
 '
 
-dat <- read.csv(text=dat)
+dat <- read.csv(text=dat,comment.char="#")
 
 parus.example <- function (proc = c("Gompertz","Ricker"),
                            meas = c("lognormal","Poisson","negbin")) {



More information about the pomp-commits mailing list