Подгонка усеченного логнормального распределения в R - PullRequest
0 голосов
/ 26 марта 2019

Ссылаясь на тему:

Подгонка логнормального распределения к усеченным данным в R

Я пытаюсь оценить параметры усеченного логнормального распределения:

library(fitdistrplus)
library(truncdist)


D <- rlnorm(1000,meanlog = -0.75, sdlog = 1.5)
# Censor data #
min <- 0.10
max <- 20
Dt <- D[D > min]
Dt <- Dt[Dt <= max]

#fitt <- fitdist(Dt, "lognormal", lower = min, upper = max)

dtruncated_log_normal <- function(x, meanlog, sdlog) 
  dtrunc(x, "lnorm", a=.10, b=20, meanlog=meanlog, sdlog=sdlog)
ptruncated_log_normal <- function(x, meanlog, sdlog) 
  ptrunc(x, "lnorm", a=.10, b=20, meanlog=meanlog, sdlog=sdlog)

fitdist(Dt, "truncated_log_normal", start = c(meanlog=0, sdlog=1))

Но я получаю следующую ошибку:

Error in manageparam(start.arg = start, fix.arg = fix.arg, obs = data,  : 
  Wrong type of argument for start

Где проблема? Речь идет о начальных значениях. Но я не могу решить это.

Буду очень рад любой помощи.

Большое спасибо.

...