Как сохранить BI C, AI C и maxlog правдоподобие в переменных при использовании метода maxLogL - PullRequest
0 голосов
/ 29 марта 2020

Я подгоняю некоторое распределение методом Максимизации Вероятности. Метод, реализованный в R "maxLogL" - это удивительный инструмент, который прекрасно работает. Документация:

AI C и BI C печатаются с помощью функции итога, но я хочу сохранить эти значения в переменных. Вот код, который вы можете легко воспроизвести:

library(EstimationTools)

set.seed(10)
z <- rnorm(n = 1000, mean = 0.1, sd = 1)
fit1 <- maxlogL(x = z, dist = 'dnorm', start=c(0, 2), lower = 0, upper = 2)
a <-summary(fit1)

, который печатает следующее:

Optimization routine: nlminb 
Standard Error calculation: Hessian from optim 

       AIC      BIC
  2824.494 2820.494

     Estimate  Std. Error
mean  0.011375     0.0313
sd    0.991346     0.0222

Мой вопрос:

  1. Как сохранить ИИ C и BI C сохраняются в переменной?
  2. Сохраняется ли где-нибудь максимальное значение правдоподобия? Я мог бы использовать реверс-инжиниринг с BI C и AI C для его получения, но я бы хотел этого избежать.

Спасибо

PD: я особенно хочу использовать функцию maxLogL.

1 Ответ

0 голосов
/ 29 марта 2020

Более близкое решение, которое я достиг для вопроса 1:

results <- capture.output(summary(fit1))

, которое сохраняет вывод консоли на переменную. Однако это раздражающее решение, потому что тогда мне нужно обработать текст. Должно быть что-то простое и лучшее.

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