Найти (или рассчитать) логарифмическую вероятность, AIC и BIC для модели SUR (для каждого уравнения) с помощью systemfit
Я оценил SUR
модель с systemfit
(R пакет). С приблизительными результатами я пытаюсь получить logLik, AIC и BIC для каждой модели SUR. Для всего журнала я могу просто использовать logLik
функция. Однако, как я могу получить соответствие качества (logLik, AIC и BIC) для каждого уравнения модели SUR?
Вот пример оценки SUR с systemfit
,
data( "Kmenta" )
eqDemand <- consump ~ price + incomeeq
Supply <- consump ~ price + farmPrice + trend
system <- list( demand = eqDemand, supply = eqSupply )
fitsur <- systemfit( system, "SUR", data = Kmenta )
## residuals of all equations
logLik( fitsur )
Как я могу получить logLik
, AIC и BIC для каждого уравнения (eqDemand
,Supply
)?