Я пытаюсь подогнать усеченное нормальное распределение к данным, используя fitdistrplus::fitdistr
и указав верхнюю и нижнюю границы. Однако при сравнении параметров MLE-подгонки с параметрами MLE-подгонки без ограничений они кажутся одинаковыми.
library(fitdistrplus)
library(MASS)
dt <- rnorm(100, 1, 0.5)
cat("truncated:", fitdistr(dt, "normal", lower = 0, upper = 1.5, method = "mle")$estimate,
"original:", fitdist(dt, "norm", method = "mle")$estimate, sep = "\n")
truncated:
1.034495
0.4112629
original:
1.034495
0.4112629
Я не гений статистики, но я почти уверен, что параметры должны быть другими, потому что усечение распределения изменит как среднее, так и sd (потому что распределение масштабируется). Это правильно?
Спасибо за совет
Ваше здоровье,
Саймон
lower=
иupper=
не создают усеченный дистрибутив. Они просто передаются вoptim()
, где они контролируют, где искать оптимальное значение. Это не меняет дистрибутив. - person MrFlick   schedule 09.02.2016optim
выдавал ошибки, предположительно потому, что я не указал нижнюю границу стандартного отклонения, и он начал тестировать отрицательные значения. Я вернусь к этому и посмотрю, решит ли указание границ эту проблему. - person e-cology   schedule 10.02.2016