Por favor, tenha paciência comigo se isso for bastante tênue e fique à vontade para fazer perguntas se eu tiver deixado alguma coisa de fora ...

Estou tentando fazer alguns cálculos de vento extremo de 50 anos com base no link a seguir

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

les parecem usar a distribuição gumbel, então usei a função gumbel no pacote "evir" para ajustar a distribuição aos dados e a função dgumbel no pacote "evd" como a função de plotage

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

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

Tentei plotar isso usando a função stat_function do ggplot2, assim (exceto agora eu coloquei valores fictícios para loc e scale.

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

Estou tendo o erro a seguir

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

Não tenho certeza se estou fazendo isso da maneira certa. Qualquer dica será muito apreciada