уйста, потерпите меня, если это довольно незначительно, и не стесняйтесь задавать вопросы, если я что-то упустил ...

Я пытаюсь сделать примерно 50-летний расчет экстремального ветра на основе следующей ссылки

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

Кажется, что они используют распределение gumbel, поэтому я использовал функцию gumbel в пакете «evir», чтобы подогнать распределение к данным, и функцию dgumbel в пакете «evd» в качестве функции построения графиков.

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

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

Затем я попытался построить это с помощью stat_function в ggplot2, вот так (за исключением того, что сейчас я ввел фиктивные значения для loc и scale.

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

Я получаю следующую ошибку:

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

Я не уверен, правильно ли я поступаю. Любые указатели будут высоко оценены.