Estoy probando varios modelos de panel basados ​​en el mismo conjunto de datos, que es grande con 265.904 observaciones.

Digamos que estoy comparando los modelos (1), (2) y (3) de la Tabla 4.

Table_4_1 <- plm(DV ~ X1   +    X3 + X4 + X5 + X6, index="firm_id", model="within", data=Survivors_upto_2019) 
Table_4_2 <- plm(DV ~      X2 + X3 + X4 + X5 + X6, index="firm_id", model="within", data=Survivors_upto_2019)  
Table_4_3 <- plm(DV ~      X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9, index="firm_id", model="within", data=Survivors_upto_2019)

summary(Table_4_1)  
summary(Table_4_2)  
summary(Table_4_3)

Table_4_1 y Table_4_2 no están anidados, pero Table_4_2 y Table_4_3 están anidados. Por lo tanto, utilicé los siguientes comandos,

jtest(Table_4_1,Table_4_2)  
pFtest(Table_4_2, Table_4_3)  

Para la prueba J, obtuve los siguientes resultados.

Estimate Std. Error t value  Pr(>|t|)    
M1 + fitted(M2)  0.95029  0.0034222 277.683 < 2.2e-16 ***  
M2 + fitted(M1)  0.35485  0.0075335  47.103 < 2.2e-16 ***

Para pFtest, recibí un error.

Warning message:  
In pf(stat, df1, df2, lower.tail = FALSE) : NaNs produced

Mis preguntas son:

1. ¿Cómo interpreto los resultados de la estadística J?

2. Para una comparación anidada de paneles, ¿cómo debo solucionar pFtest?