2017-03-11 2 views
0

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 !!

+0

Versuchen Klein G. http://artax.karlin.mff.cuni.cz/r-help/library/ActuDistns/html/qgumbel.html – Elin

+0

versuchte es auch nicht, daher denke ich, dass ich die Funktion selbst programmieren muss. – rsc05

Antwort

0

Als Elin im Kommentar erwähnt ich habe es durch das Kapital G zu klein g verändert und hinzugefügt meine eigene Funktion

qgumbel <- function(p,shape,scale) shape-scale *log(-log(p)) 

und jetzt der Code funktioniert

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) 
# parameter values of the marginal distributions 
Verwandte Themen