Максимальная вероятность равномерного распределения - PullRequest
0 голосов
/ 17 апреля 2020

Я решал книгу, в которой нашел один пример, в котором мне нужно смоделировать 100 выборок из плотности, когда тета = 2, f (x) = 3 * x * x на [0, 1].

Это то, что я пробовал:
RandomVariable U = X ^ 3 = F (X) => X = U ^ 1/3 и U непрерывно равномерно на [0, 1].

library(bbmle)
x <- runif(100, 0, 1)
LL <- function(min, max) {
  R <- suppressWarnings(dunif(x, min, max))
  R <- R ^ (1 / 3)
  -sum(log(R))
}
mle2(LL, start = list(min = 0, max = 1), methos = "L-BFGS-B")

Но, похоже, это не работает, я получаю сообщение об ошибке:
Ошибка в optim (par = c (min = 0, max = 1), fn = function (p): не конечное разностное значение [2].
Я также попробовал mle пакета stats4 вместо mle2, но я получаю ту же ошибку.
Я также проверил это:
Как найти MLE равномерного распределения?
Пожалуйста, помогите мне в этом.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...