Трехпараметрическая функция плотности вероятности Вейбулла класса mle2

Question

Я подгоняю свои данные (данные о ветре) к различным распределениям (распределения Вейбулла и гамма с двумя и тремя параметрами). Я использовал описанный здесь метод (mle2, используя метод «L-BFGS-B») для оценки параметров трехпараметрического распределения Вейбулла: Ошибка оценки параметра распределения Вейбулла Я использовал fitdist (из библиотеки fitdistrplus) для оценки параметров распределения двух Вейбулла и гаммы. Теперь я хочу сравнить соответствие двухпараметрического распределения Вейбулла, трехпараметрического Вейбулла и гаммы с моими данными на одном графике. Я смог сравнить соответствие двухпараметрического Вейбулла и гаммы, но не смог включить трехпараметрический Вейбулл (это другой класс). Можно ли построить двухпараметрический и трехпараметрический Вейбулла на одном графике. Буду очень признателен за любые предложения. Спасибо.

оценка параметров трехпараметрического и двухпараметрического Вейбулла:

    library(bbmle)
    library(FAdist)

    y = my.data
    eps <- 1e-8
    weibull3 <- mle2(y ~ dweibull3(shape=shape, scale = scale, thres = thres), start=list(shape=1, scale=1, thres=-5),lower=c(shape=0,scale=0,thres=-Inf), upper=c(shape=Inf,scale=Inf,thres=min(y)-eps), method="L-BFGS-B", data=data.frame(y))

    weibull3
    call:

     mle2(minuslogl = y ~ dweibull3(shape = shape, scale = scale, thres = thres), start = list(shape = 1, scale = 1, thres = -5), method = "L-BFGS-B", data = data.frame(y), lower = c(shape = 0, scale = 0, thres = -Inf), upper = c(shape = Inf, scale = Inf, thres = min(y) - eps))

     coefficients:
      shape      scale      thre
    2.4828383  5.5248459 -0.2530033


    Log-likelihood: -1358214 


    library (fitdistrplus)
    weibull2 = fitdist(my.data, "weibull")
    weibull2

    Fitting of the distribution ' weibull ' by maximum likelihood 
       Parameters:
          estimate  Std. Error
    shape 2.324290 0.002344834
    scale 5.238753 0.002980665