# \donttest{
# 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)
# }
Run the code above in your browser using DataLab