export(IDEA)
export(WP)
export(seqB)
+importFrom(stats,pexp)
+importFrom(stats,pgamma)
+importFrom(stats,qexp)
+importFrom(stats,qgamma)
#' 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)) {
#' \code{J0} (the location for the maximum likelihood estimator \code{Rhat}). If \code{JJ} and \code{J0}
#' are not the same, this means that the maximum likelihood estimator is not unique.
#'
+#' @importFrom stats pgamma qgamma
+#'
#' @keywords internal
WP_unknown <- function(NT, B=100, shape.max=10, scale.max=10, tol=0.999) {
shape <- seq(0, shape.max, length.out=B+1)