0

我想将伽玛和帕累托分布拟合到我的数据中。这是示例数据。在这一点上，我只想拟合分布，而不是它们是否是适合的分布

  x <- c(3449.48886699726, 2890.01223895794, 2868.97132193186, 2559.69741012966, 
         2999.1270689335, 3455.34137281101, 3772.37075245719, 3744.05575693446, 
         3555.42346782386, 3453.59115457264, 3488.64443911085, 3949.28561615232, 
         3794.83245108087, 4242.77755251986, 3894.6542925462, 3497.71878257211, 
         3158.22673257603, 2733.45331438388, 2410.35327717621, 2372.67597041366, 
         3538.97187741642, 2435.33759646514, 3532.46932136698, 3649.01008976503, 
         2574.36498196242, 3583.20376922858, 3755.22256460135, 3417.08141838327, 
         3578.08810809661)

  library(MASS)
  library(fitdistrplus)      

  gammaMLE <- fitdist(x, 'gamma', method = 'mle')
  
  dpareto <- function(x, a, l){
    
    fx <- ((a) * (l ^ a)) / ((l + x) ^ (a + 1))
    }
  
  paretoMLE <- fitdist(x, 'pareto', method = 'mle')
  
  Error in checkparamlist(arg_startfix$start.arg, arg_startfix$fix.arg,  : 
                            'start' must specify names which are arguments to 'distr'

我无法理解这个错误，想知道这里是否有人有将帕累托分布拟合到数据的经验。

4

0 回答 0