if (FALSE) {
library(Distance)
library(dsm)
# load the Gulf of Mexico dolphin data (see ?mexdolphins)
data(mexdolphins)
# fit a detection function and look at the summary
hr.model <- ds(distdata, truncation=6000,
key = "hr", adjustment = NULL)
summary(hr.model)
# fit a simple smooth of x and y
mod1 <- dsm(count~s(x, y), hr.model, segdata, obsdata)
summary(mod1)
# calculate the variance by 500 moving block bootstraps
mod1.movblk <- dsm_var_movblk(mod1, preddata, n.boot = 500,
block.size = 3, samp.unit.name = "Transect.Label",
off.set = preddata$area,
bar = TRUE, bs.file = "mexico-bs.csv", ds.uncertainty = TRUE)
}
Run the code above in your browser using DataLab