-source("WP_known.R")
-source("WP_unknown.R")
-
#' WP method
#'
#' This function implements an R0 estimation due to White and Pagano (Statistics in Medicine, 2008).
#' res4 <- WP(NT=NT[1:5], mu=5/7) # serial distribution has mean of five days
#' res4$Rhat
#'
+#' @importFrom stats pexp qexp
+#'
#' @export
WP <- function(NT, mu=NA, search=list(B=100, shape.max=10, scale.max=10), tol=0.999) {
if (is.na(mu)) {