1 |
n |
|
g |
|
h |
|
p |
|
cmat |
|
SEED |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (n, g = 0, h = 0, p = 2, cmat = diag(rep(1, p)), SEED = FALSE)
{
x <- rmulnorm(n, p, cmat, SEED = SEED)
for (j in 1:p) {
if (g > 0) {
x[, j] <- (exp(g * x[, j]) - 1) * exp(h * x[, j]^2/2)/g
}
if (g == 0)
x[, j] <- x[, j] * exp(h * x[, j]^2/2)
}
x
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.