Извлечение остаточных степеней свободы из объекта pgmm
У меня проблемы с извлечением количества степеней свободы из pgmm
объект. Команда df.resudual()
раньше работал просто отлично, а сейчас нет. Вместо того, чтобы возвращать число, он возвращает df.residual()
функционировать сам.
Вот минимальный пример использования plm
документация:
data("EmplUK", package="plm")
emp.gmm <- pgmm(log(emp)~lag(log(emp), 1:2)+lag(log(wage), 0:1)+log(capital)+
lag(log(output), 0:1)|lag(log(emp), 2:99),
data = EmplUK, effect = "twoways", model = "twosteps")
df.residual(emp.gmm)
> function (object, ...)
> UseMethod("df.residual")
> <bytecode: 0x39a9370>
> <environment: namespace:stats>
Я попытался заглянуть внутрь команды pgmm, чтобы выяснить, что происходит, и я думаю, что проблема в этой строке:
result <- list(coefficients = coefficients, residuals = residuals,
vcov = vcov, fitted.values = fitted.values, df.residual = df.residual,
model = yX, W = W, A1 = A1, A2 = A2, call = cl, args = args)
Так как df.residual
ранее не было определено, тогда запись списка df.residual = df.residual просто возвращает глобальную запись для (универсальной) функции.
Это ошибка в пакете? Или здесь что-то еще происходит? Я был в состоянии извлечь DoFs из объектов pgmm, прежде чем просто отлично.
PS. Это было закрыто на Cross Validated. Первоначально здесь: https://stats.stackexchange.com/questions/339239/extracting-residual-degrees-of-freedom-from-pgmm-object