问题描述
我想从plm
对象中提取对数似然。
当将基本logLik
软件包中的stats
函数与felm
软件包中的lfe
或{{1 }}程序包,但不包含任何feols
对象,其中产生的错误消息是:
fixest
我检查了一下,发现plm
不是Error in UseMethod("logLik") :
no applicable method for 'logLik' applied to an object of class "c('plm','panelmodel')"
包中定义的类(请参见plm
)。
我在概念上是否缺少某些东西(意识到这是ML估算值)?例如,为什么lfe和fixest使其起作用而plm不能起作用?有解决方法吗?
谢谢!
stats
解决方法
我想我明白了:
object
是一个plm对象,就像上面的zz
。
logLik.plm <- function(object){
-plm:::nobs.panelmodel(summary(object)) * log(2 * var(object$residuals) * pi)/2 - plm:::deviance.panelmodel(summary(object))/(2 * var(object$residuals))
}