Я оценил SUR
модель с systemfit
(пакет R).С приблизительными результатами я пытаюсь получить logLik, AIC и BIC для каждой модели SUR.Для всего logLik я могу просто использовать функцию 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
)?