Hola de nuevo, aquí va algo de código de información de los datos: estoy
usando ismev, evd y una fuente (Rstufflee.r) que podeis encontrar aquí:
http://www.mas.ncl.ac.uk/~ndw/Rstufflee.r
Se supone que gev.ret te devuelve el dato para el periodo que le indiques….
Gracias!
> q_arenas_mean
[1] 509.00 153.00 76.00 307.00 435.00 252.00 325.00 313.00 303.00 204.00
216.00 288.00 314.00 261.00 413.00
[16] 376.00 333.00 84.65 28.52 81.24 41.79 302.59 84.17 46.16 19.22
32.22
> fit <- gev.fit(q_arenas_mean)
$conv
[1] 0
$nllh
[1] 165.185
$mle
[1] 174.5922200 139.8051747 -0.3049459
$se
[1] 31.2325845 23.6186751 0.1666812
> gev.ret(q_arenas_mean,100)
$conv
[1] 0
$nllh
[1] 165.185
$mle
[1] 174.5922200 139.8051747 -0.3049459
$se
[1] 31.2325845 23.6186751 0.1666812
Error en return(ret.level, ret.se) :
multi-argument returns are not permitted
> gev.ret(q_arenas_mean)
$conv
[1] 0
$nllh
[1] 165.185
$mle
[1] 174.5922200 139.8051747 -0.3049459
$se
[1] 31.2325845 23.6186751 0.1666812
Error en 1/period : period'' esta perdido
De: Carlos Ortega [mailto:cof@qualityexcellence.es]
Enviado el: jueves, 07 de junio de 2012 0:24
Para: Virginia Ruiz
CC: r-help-es@r-project.org
Asunto: Re: [R-es] estimar cuantiles de un ajuste GEV
Hola Virginia,
Si añades en tu pregunta algo más de contexto, un conjunto de datos y algo
de código que podamos reproducir aumentan las posibilidades de que podamos
ayudarte de una forma más efectiva.
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 6 de junio de 2012 19:14, Virginia Ruiz <v.ruiz@igme.es> escribió:
Hola,
Quisiera saber como se puede extraer un cuantil (un caudal para un
determinado periodo de retorno) una vez que ya he ajustado mi serie de
caudales máximos con una función de distribución (eg: GEV). He encontrado
una función (gev.ret) dentro del paquete extremes, pero me da error:
Error en return(ret.level, ret.se) :
multi-argument returns are not permitted
lo que me gustaría es conocer el caudal asociado a T=100 años, y si es
posible con un intervalo de confianza.
Gracias!
Saludos,
Virginia
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]