##single-season occupancy model example modified from ?occu
if (FALSE) {
require(unmarked)
##single season
data(frogs)
pferUMF <- unmarkedFrameOccu(pfer.bin)
## add some fake covariates for illustration
siteCovs(pferUMF) <- data.frame(sitevar1 = rnorm(numSites(pferUMF)),
sitevar2 = rnorm(numSites(pferUMF)))
## observation covariates are in site-major, observation-minor order
obsCovs(pferUMF) <- data.frame(obsvar1 = rnorm(numSites(pferUMF) *
obsNum(pferUMF)))
##run model set
fm1 <- occu(~ obsvar1 ~ sitevar1, pferUMF)
##extract log-likelihood
extractLL(fm1)
detach(package:unmarked)
}
Run the code above in your browser using DataLab