Ich bin relativ neu zu R und versuche, eine Gumbel Copula mit zu simulieren und einer dieser Ränder ist Gumbel pdf und der andere ist exponentiell.Fehler in eval (expr, envir, enclos): konnte Funktion "qGumbel" nicht finden
durch diesen Code
G3 <- gumbelCopula(1.37, dim=2)
gMvd2 <- mvdc(G3, c("Gumbel","exp"), param = list(list(shape=10.2988298881251, scale=1.02463492397923), list(rate=4)))
set.seed(11)
# n <- if(Xtras) 10000 else 200 # sample size (realistic vs short for example)
x <- rMvdc(100, gMvd2)
head(x)
Der Ausgang war:
Error in eval(expr, envir, enclos): could not find function "qGumbel"
Traceback:
1. rMvdc(100, gMvd2)
2. eval(qdf.expr, list(x = u[, i]))
3. eval(expr, envir, enclos)
Sie erwähnen auch Bemerkung
Warnmeldung in mvdCheckM (Ränder, "p"):
"margins correct? Currently, have no function(s) named: pGumbel"Warning message in mvdCheckM(margins, "d"):
"margins correct? Currently, have no function(s) named: dGumbel"
Ich denke, der Profi In den Parameterwerten der Randverteilungen liegt ein Fehler, denn wenn ich tippe
G3 <- gumbelCopula(1.37, dim=2)
gMvd2 <- mvdc(G3, c("exp","exp"), param = list(list(rate=2), list(rate=4)))
set.seed(11)
# n <- if(Xtras) 10000 else 200 # sample size (realistic vs short for example)
x <- rMvdc(100, gMvd2)
head(x)
Kein Problem !!
Versuchen Klein G. http://artax.karlin.mff.cuni.cz/r-help/library/ActuDistns/html/qgumbel.html – Elin
versuchte es auch nicht, daher denke ich, dass ich die Funktion selbst programmieren muss. – rsc05