set.seed(87654)
x <- rnorm(1000, mean=0, sd=20)
x <- x[x >= 0]
x <- units::set_units(x, "mi")
dfunc <- dfuncEstim(x~1, likelihood="halfnorm")
plot(dfunc, nbins = 40, col="lightgrey", border=NA, vertLines=FALSE)
lines(dfunc, col="grey", lwd=15)
lines(dfunc, col="black", lwd=5, lty = 2)
# Multiple lines
data(sparrowDetectionData)
data(sparrowSiteData)
dfuncObs <- dfuncEstim(formula = dist ~ observer
, likelihood = "halfnorm"
, detectionData = sparrowDetectionData
, siteData = sparrowSiteData)
plot(dfuncObs
, vertLines = FALSE
, lty = 0
, col = c("grey","lightgrey")
, border=NA
, main="Detection by observer"
, legend = FALSE)
y <- lines(dfuncObs
, newdata = data.frame(observer = levels(sparrowSiteData$observer))
, col = palette.colors(length(levels(sparrowSiteData$observer)))
, lty = 1
, lwd = 4)
head(y) # values returned, same as predict method
Run the code above in your browser using DataLab