Сейчас я работаю над подбором множественных распределений вероятностей с помощью пакета «fitdistrplus», и я знаю, что в этом пакете он предложил нам подобрать некоторые распространенные распределения вероятностей, такие как нормальное распределение, распределение Коши ......, и вы также можете напишите свой собственный дистрибутив, например:
dgumbel <- function(x, a, b) 1/b*exp((a-x)/b)*exp(-exp((a-x)/b))
pgumbel <- function(q, a, b) exp(-exp((a-q)/b))
qgumbel <- function(p, a, b) a-b*log(-log(p))
fitgumbel <- fitdist(serving, "gumbel", start=list(a=10, b=10))
summary(fitgumbel)
plot(fitgumbel)
Приведенный выше код скопирован из примера, предлагаемого пакетом «fitdistrplus». И теперь я хочу соответствовать распределению levy, и этот пакет не предоставляет связанных функций, поэтому мне нужно сначала написать его. Затем я обнаружил, что есть еще один пакет под названием «rmutil», который позволяет нам создавать дистрибутив levy, например:
dlevy(y, m=0, s=1, log=FALSE)
plevy(q, m=0, s=1)
qlevy(p, m=0, s=1)
rlevy(n, m=0, s=1)
. Поэтому я скопировал исходный код и попытался приспособить дистрибутив levy, и он дал мне предупреждение сообщение: «Функция dlevy должна возвращать вектор со значениями NaN, когда входные данные имеют противоречивые значения и не вызывает ошибку»
Итак, я хочу спросить, как построить распределение levy, которое может запускать функцию fitdist ()? Спасибо.