R/mcCovMat.Bbdpars.R

In dbd: Discretised Beta Distribution

mcCovMat.Bbdpars <- function(object,nsim=100,seed=NULL,maxit=1000) {
    if(is.null(seed)) seed <- sample(1:1e5,1)
    set.seed(seed)
    par0   <- unlist(object[1:2])
    simdat <- simulate(object,nsim=nsim)
    fitz   <- lapply(simdat,mleBb,size=attr(object,"size"),
                     par0=par0,maxit=maxit,covmat=FALSE)
    mcCovMatEngine(fitz,par0,seed)
}