library(RMark)
data(mallard)
mallard$habitat <- as.factor(mallard$habitat)
mallard.pr <- process.data(mallard,
nocc=90,
model="Nest",
groups=("habitat"))
ddl=make.design.data(mallard.pr)
df=data.frame(precip=1,group="Native",time=1:89)
df=rbind(df,data.frame(precip=2,group="Wetland",time=1:89))
df=rbind(df,data.frame(precip=3,group="Planted",time=1:89))
df=rbind(df,data.frame(precip=3,group="Roadside",time=1:89))
merge_design.covariates(ddl$S,df,bygroup=TRUE,bytime=TRUE)
ddl$S=merge_design.covariates(ddl$S,df,bygroup=TRUE,bytime=TRUE)
mark(mallard.pr,ddl,model.parameters=list(S=list(formula=~precip)))
Users browsing this forum: No registered users and 17 guests