С небольшой модификацией вашего кода (добавление, geom) он отлично работает для меня.
уйста, потерпите меня, если это довольно незначительно, и не стесняйтесь задавать вопросы, если я что-то упустил ...
Я пытаюсь сделать примерно 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)
Я не уверен, правильно ли я поступаю. Любые указатели будут высоко оценены.