J'ai changé de PC et version de R ou rstudio et mon ancien script qui utilisait le package "games" avec la fonction Vuong ne peut plus marcher car la package n'existe pas pour cette version de R. Je compare deux régressions linéaires non-nested faites avec lm.
J'ai essayé une alternative. J'ai trouvé le package "nonnest2" avec la fonction testvuong. Mon problème c'est que j'ai a priori des résultats divergents. Que font-ils de différent ou comment puis je mettre les bonnes options pour retrouver la même chose ??? ça m'inquiète un peu.
Je ne vois pas trop comment donner du code reproductible mais voilà mon ancien code et le nouveau :
Code : Tout sélectionner
library(games) vuongtest(a, b)
et ça donnait ça comme résultat
##
## Vuong test for non-nested models
##
## Model 1 log-likelihood: -20709
## Model 2 log-likelihood: -20696
## Observations: 4382
## Test statistic: -1.2
##
## Neither model is significantly preferred (p = 0.24)
et maintenant :
Code : Tout sélectionner
library(nonnest2)
vuongtest(a, b)
avec ceci comme résultat :
Variance test
H0: Model 1 and Model 2 are indistinguishable
H1: Model 1 and Model 2 are distinguishable
w2 = 0.011, p = 0.0237
Non-nested likelihood ratio test
H0: Model fits are equal for the focal population
H1A: Model 1 fits better than Model 2
z = -1.789, p = 0.963
H1B: Model 2 fits better than Model 1
z = -1.789, p = 0.03678
Sinon, y a-t-il un remplaçant exact de "games" ?
Merci d'avance pour vos éclaircissements