Я хотел бы выполнить тест общего фактора, проводимый с использованием
тест отношения правдоподобия на модели пространственной ошибки и модели пространственного Дурбина
(оба в формате панели).
Я не нашел команду в R, которая помогла бы мне
провести тест. Поэтому я пытаюсь выполнить тест вручную, используя пакет splm и данные, доступные в пакетах splm и plm.
Вот мой код:
library(splm)
library(plm)
data(Produc, package="plm")
data(usaww)
Produc <- pdata.frame(Produc, index = c("state", "year"))
usaww<- mat2listw(usaww, style="W")
Produc$slagUnemp <- slag(Produc$unemp, listw = usaww)
durbin <- spml(gsp~unemp + slagUnemp,
data=Produc, listw=usaww, effect = "twoways", model="within", lag=TRUE, spatial.error = "none", quiet = FALSE)
spem <- spml(gsp ~ unemp,
data=Produc, listw=usaww, effect = "twoways", model="within", lag=FALSE, spatial.error = "b", quiet = FALSE)
Правильно ли брать последнее значение функции из вывода консоли для вычисления отношения правдоподобия?
Значение, в данном примере, для вычисления 2 * (функция (не вложенная) - функция (вложенная)) = 2 * (10261.79 - 10255.74) = 12,1?
Если это правильно, как я могу вычислить значения значимости?
Если это не правильно, есть ли лучший способ вычислить это?