mledist для распределения Парето

Пытаюсь найти дистрибутив. с помощью функции mledist. Но каждый раз, когда я его использую, я нахожу сообщение об ошибке:

[1] "Error in optim(par = vstart, fn = fnobj, fix.arg = fix.arg, obs = data,  : \n  non-finite finite-difference value [2]\n"
attr(,"class")
[1] "try-error"
attr(,"condition")
<simpleError in optim(par = vstart, fn = fnobj, fix.arg = fix.arg, obs = data,     ddistnam = ddistname, hessian = TRUE, method = meth, lower = lower,     upper = upper, ...): non-finite finite-difference value [2]>

Дело в том, что я знаю, что у "Алаэ" распределение по Парето, но не могу понять, где моя ошибка. Пожалуйста, помогите :) Это мой код

library(fitdistrplus)
library(CASdatasets)
data(lossalae)
pertes<-lossalae
Alae <- pertes$ALAE   
Alae_par <- mledist(Alae, 'pareto',start=c(shape=5,scale=2))

TNX


r distribution pareto-chart
person NANA    schedule 17.02.2015    source источник
comment
Я получаю другую ошибку: Ошибка в mledist(Alae, pareto, start = c(shape = 5, scale = 2)) : функция dpareto должна быть определена   -  person IRTFM    schedule 17.02.2015

Ответы (1)


arrow_upward
0
arrow_downward

Учитывая, что источником CASdatasets является канадский автор страховых услуг, я предполагаю, что функция dpareto должна исходить от pkg:actuar:

> library(actuar)

Attaching package: ‘actuar’

The following object is masked from ‘package:grDevices’:

    cm

> library(fitdistrplus)
> library(CASdatasets)
> data(lossalae)
> pertes<-lossalae
> Alae <- pertes$ALAE   
> Alae_par <- mledist(Alae, 'pareto',start=c(shape=5,scale=2))
Warning messages:
1: In dpareto(c(3806L, 5658L, 321L, 305L, 758L, 8768L, 1805L, 78L,  :
  NaNs produced
snipped 6 other warnings

> Alae_par
$estimate
       shape        scale 
    2.221351 15106.971534 

$convergence
[1] 0

$loglik
[1] -15413.45

$hessian
             shape         scale
shape 303.98838362 -3.083551e-02
scale  -0.03083551  3.637979e-06

$optim.function
[1] "optim"

$fix.arg
NULL
person IRTFM    schedule 17.02.2015