remGRM: remGRM

View source: R/remGRM.R

remGRMR Documentation

remGRM

Description

remGRM fits the generalized removal model to data collected from repeated removal episodes from M sites over T primary periods with each primary consisting of J secondary periods. The model also facilitates the analysis of index (count) data collected in conjunction with the removal data to make joint inference on abundance.

Usage

remGRM(lamformula, phiformula, detformula, mdetformula, data, mixture = c("P", "NB"), K,
                  starts, method="BFGS", se=TRUE, ...)

Arguments

lamformula

formula for the latent abundance component.

detformula

formula for the removal detection component.

mdetformula

formula for the index detection component.

data

A eFrameGRM object containing the response (counts, index) and site-level covariates. see eFrameGRM for how to format the required data.

K

upper bound for superpopulation abundance

starts

Initial values for parameters

method

Optimisation method

se

flag to return the standard error (hessian).

model

for abundance, either Poisson 'P' or negative binomial 'NB'

Value

a efitGRM model object.

Examples

 rem<- san_nic_rem$rem
 ym<- san_nic_rem$ym
 emf <- eFrameGRM(rem, ym, numPrimary=1)
 mod <- remGRM(~1, ~1, ~1, ~1, data=emf)
 Nhat<- calcN(mod)


dslramsey/eradicate documentation built on March 16, 2024, 1:40 p.m.