do.salamander=function()
{
data(salamander)
occ.p.dot=mark(salamander,model="Occupancy")
occ.p.time=mark(salamander,model="Occupancy",
model.parameters=list(p=list(formula=~time)))
occ.p.mixture=mark(salamander,model="OccupHet",
model.parameters=list(p=list(formula=~mixture)))
return(collect.models())
}
salamander.results=do.salamander()
print(salamander.results)
Run the code above in your browser using DataLab