Nothing
"simCOPv" <-
function(u, cop=NULL, para=NULL,
reflect=c("cop", "surv", "acute", "grave",
"1", "2", "3", "4"), ...) {
simCOPmicro(u, cop=cop, para=para, reflect=reflect, ...)
}
"simCOPmicro" <-
function(u, cop=NULL, para=NULL, seed=NULL,
reflect=c("cop", "surv", "acute", "grave",
"1", "2", "3", "4"), ...) {
reflect <- match.arg(reflect)
if(! is.null(seed)) set.seed(seed)
n <- length(u); t <- runif(n); v <- vector(mode="numeric", length=n)
v <- switch(reflect,
cop = sapply(1:n, function(i) { derCOPinv(cop=cop, u[i], t[i], para=para, ...) }),
surv = sapply(1:n, function(i) { 1 - derCOPinv(cop=cop, 1-u[i], 1-t[i], para=para, ...) }),
acute = sapply(1:n, function(i) { derCOPinv(cop=cop, 1-u[i], t[i], para=para, ...) }),
grave = sapply(1:n, function(i) { 1 - derCOPinv(cop=cop, u[i], 1-t[i], para=para, ...) }),
"1" = sapply(1:n, function(i) { derCOPinv(cop=cop, u[i], t[i], para=para, ...) }),
"2" = sapply(1:n, function(i) { 1 - derCOPinv(cop=cop, 1-u[i], 1-t[i], para=para, ...) }),
"3" = sapply(1:n, function(i) { derCOPinv(cop=cop, 1-u[i], t[i], para=para, ...) }),
"4" = sapply(1:n, function(i) { 1 - derCOPinv(cop=cop, u[i], 1-t[i], para=para, ...) }))
if(any(is.na(v))) warning("could not uniroot at least for one element in derCOPinv")
return(v)
}
# surv is a reflection on the horizontal AND vertical axes
# acute is a reflection on the horizontal axis
# grave is a reflection on the verical axis
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.