# R/summary.sightest.R In SightabilityModel: Wildlife Sightability Modeling

#### Documented in summary.sightest

#' Summarize sightability estimator
#'
#' Calculates confidence interval (based on asymptotic [normal or log-normal
#' assumption])
#'
#'
#' @aliases summary.sightest summary.sightest_ratio
#' @param object Sightability object, output from call to Sight.Est function.
#' @param ... arguments to be passed to or from other methods
#' @return \item{Nhat or Ratiohat}{Sightability population estimate}
#' \item{lcl}{Lower confidence limit} \item{ucl}{Upper confidence limit}
#' @author John Fieberg and Carl James Schwarz
#' @keywords summary
#' @export
summary.sightest <-
function(object,...){
# total number of animals seen
Tot.seen <- sum(object$odat$total)

# z statistic
z <- qnorm(1-object$alpha/2) # Normal assumption confidence interval if(object$CI.method == "normal"){
temp <- rep(object$est[1], 3)+c(0, -z, z)*sqrt(object$est[2])
}
else{
#CI under lognormal assumption (see Wong p. 65-67)
#Estimated number of animals that were not seen (assumed to be lognormally distributed)
tau.m.T <- object$est[1]-Tot.seen cv2 <- object$est[2]/(tau.m.T)^2
cfact <- exp(z*sqrt(log(1+cv2)))
temp <- rep(tau.m.T, 3)*c(1,(1/cfact)*sqrt(1+cv2), cfact*sqrt(1+cv2))+ rep(Tot.seen, 3)
}
names(temp) <- NULL
temp <- format(round(temp, 0), big.mark=",")
cat("\n")
temp2 <- paste("tau.hat = ", temp[1], ";  ",100*(1-object\$alpha), "% CI = (", temp[2], ", ", temp[3], ")")
print(format(temp2[1], big.mark=","), quote = FALSE)
out.summary <- list(tau.hat = temp[1], lcl = temp[2], ucl = temp[3])
}


## Try the SightabilityModel package in your browser

Any scripts or data that you put into this service are public.

SightabilityModel documentation built on Aug. 20, 2023, 1:08 a.m.