Por favor, tenga paciencia conmigo si esto es bastante tenue, y siéntase libre de hacer preguntas si he dejado algo fuera .

Estoy intentando hacer unos 50 años de cálculos de vientos extremos basados en el siguiente enlace

http: //www.wasp.dk/Products/weng/ExtremeWinds.ht

Parecen usar la distribución gumbel, así que he usado la función gumbel en el paquete "evir" para ajustar la distribución a los datos, y la función dgumbel en el paquete "evd" como la función de trazado.

package("evd")
package("evir")

speeds2 <- data.frame(speed=sample(10:50,1000,rep=TRUE))
gumbel(speeds2$speed)

uego, traté de trazar esto usando la función stat_función de ggplot2, así (excepto por ahora, he puesto valores ficticios para loc y scale.

library(ggplot2)
ggplot(data=speeds2, aes(x=speed)) + 
  stat_function(fun=dgumbel, args=list(loc=1, scale=0.5))

Obtuve el siguiente error

Error in dgev(x, loc = loc, scale = scale, shape = 0, log = log) : 
  unused argument(s) (loc = loc, scale = scale, shape = 0, log = log)

No estoy seguro si estoy haciendo esto de la manera correcta. Cualquier sugerencia sería muy apreciada.