Nothing
# These functions are
# Copyright (C) 1998-2023 T.W. Yee, University of Auckland.
# All rights reserved.
simulate.vlm <- function (object, nsim = 1, seed = NULL, ...) {
if (!exists(".Random.seed", envir = .GlobalEnv, inherits = FALSE))
runif(1)
if (is.null(seed)) {
RNGstate <- get(".Random.seed", envir = .GlobalEnv)
} else {
R.seed <- get(".Random.seed", envir = .GlobalEnv)
set.seed(seed)
RNGstate <- structure(seed, kind = as.list(RNGkind()))
on.exit(assign(".Random.seed", R.seed, envir = .GlobalEnv))
}
ftd <- fitted(object)
nm <- names(ftd)
n <- length(ftd)
ntot <- n * nsim
Fam <- if (inherits(object, "vlm")) {
object@family
} else {
stop("cannot get at the 'family' slot")
}
val <-
if (length(Fam@simslot) > 0) {
Fam@simslot(object, nsim)
} else {
stop(gettextf("family '%s' not implemented", Fam), domain = NA)
}
if (!is.list(val)) {
dim(val) <- c(n, nsim)
val <- as.data.frame(val)
} else {
class(val) <- "data.frame"
}
names(val) <- paste("sim", seq_len(nsim), sep = "_")
if (!is.null(nm))
row.names(val) <- nm
attr(val, "seed") <- RNGstate
val
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.