IELogitNormalMR {RMark} | R Documentation |
Example of Immigration-Emigration LogitNormal Mark-Resight model
Description
Data and example illustrating Immigration-Emigration LogitNormal Mark-Resight model
Format
A data frame with 34 observations on the following variable.
- ch
a character vector
Examples
# This example is excluded from testing to reduce package check time
data(IELogitNormalMR)
IElogitNor.proc=process.data(IELogitNormalMR,model="IELogitNormalMR",
counts=list("Marked Superpopulation"=c(28, 29, 30, 30, 30, 33, 33, 33, 33, 34, 34, 34),
"Unmarked Seen"=c(264, 161, 152, 217, 217, 160, 195, 159, 166, 152, 175, 190),
"Marked Unidentified"=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)),
time.intervals=c(0,0,0,1,0,0,0,1,0,0,0))
IElogitNor.ddl=make.design.data(IElogitNor.proc)
mod1=mark(IElogitNor.proc,IElogitNor.ddl,
model.parameters=list(p=list(formula=~-1+session),
sigma=list(formula=~session),
alpha=list(formula=~-1+session:time),
Nstar=list(formula=~session),
Nbar=list(formula=~session)),delete=TRUE)
summary(mod1)
# You can use the initial value to get a better estimate.
mod2=mark(IElogitNor.proc,IElogitNor.ddl,
model.parameters=list(p=list(formula=~-1+session:time),
sigma=list(formula=~session),
alpha=list(formula=~-1+session:time),
Nstar=list(formula=~session),
Nbar=list(formula=~session)),
initial=mod1,delete=TRUE)
summary(mod2)
[Package RMark version 3.0.0 Index]