Моя работа заключается в составлении графика вероятности функции. Я новичок в этой теме, особенно в R. К счастью, мне дали пример кода (для другой функции), и я попытался изменить код и посмотреть, смогу ли я получить некоторые результаты.
Это таблица формул для оценки максимального правдоподобия
Проблема, над которой я работаю
9,52
А вот код, который у меня есть.
set.seed(4302019)
theta_true <- 2
n <- 30
U <- runif(n)
Y <- theta_true * U ^ (1 / 3)
y <- rbeta(n, U, Y)
theta <- seq(1.9, 2.1, length = 1000)
maxlike <- (theta_true + 1)^n + sum(y)^(theta_true)
plot(theta_true, maxlike, xlab = expression(theta_true),
ylab = expression(L(theta_true)), type = "l", lwd = 3)
theta_mome <- (2 * mean(y) - 1) / (1 - mean(y))
theta_mle <- -n / sum(y) - 1
abline(v = theta_mle, col = "red", lwd = 2)
abline(v = theta_mome, col = "blue", lwd = 2)
abline(v = theta_true, col = "black", lty = 2, lwd = 2)
Что я получаю, запустив код выше
Я должен получить несколько линий с разными цветами, но я не вижу никаких линий! Любая помощь будет оценена.