+source("computeLL.R")
+source("WP_known.R")
+
#' WP method background function WP_unknown
#'
#' This is a background/internal function called by \code{WP}. It computes the maximum likelihood estimator
p <- diff(pgamma(0:range.max, shape=shape[i], scale=scale[j]))
p <- p / sum(p)
mle <- WP_known(NT, p)
- resLL[i,j] <- computeLL(p, NT, mle$R)
- resR0[i,j] <- mle$R
+ resLL[i,j] <- computeLL(p, NT, mle)
+ resR0[i,j] <- mle
}
}