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

Я пытаюсь подогнать усеченное нормальное распределение к данным, используя 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 (потому что распределение масштабируется). Это правильно?

Спасибо за совет

Ваше здоровье,

Саймон


person e-cology    schedule 09.02.2016    source источник
comment
lower= и upper= не создают усеченный дистрибутив. Они просто передаются в optim(), где они контролируют, где искать оптимальное значение. Это не меняет дистрибутив.   -  person MrFlick    schedule 09.02.2016
comment
Спасибо, это очень полезно! Я пытался подогнать усеченные распределения размера огня в соответствии с этим примером: ссылка Но optim выдавал ошибки, предположительно потому, что я не указал нижнюю границу стандартного отклонения, и он начал тестировать отрицательные значения. Я вернусь к этому и посмотрю, решит ли указание границ эту проблему.   -  person e-cology    schedule 10.02.2016