Valeur p globale et valeurs p par paire?

11

J'ai ajusté un modèle linéaire général dont la probabilité logarithmique est .L u

y=β0+β1X1+β2X2+β3X3,
Lu

Maintenant, je souhaite tester si les coefficients sont les mêmes.

  • Tout d'abord, test global : la probabilité logarithmique du modèle réduit est . Par test de rapport de vraisemblance, le modèle complet est significativement meilleur que le modèle réduit avec .L r p = 0,02y=β0+β1(X1+X2+X3)Lrp=0,02
  • Ensuite, ? Le modèle réduit est . Le résultat est que n'est PAS différent de avec . y = β 0 + β 1( x 1 + x 2 ) + β 2 x 3 β 1 β 2 p = 0,15β1=β2y=β0+β1(X1+X2)+β2X3β1β2p=0,15
  • De même, ? Ils sont différents avec . p = 0,007β1=β3p=0,007
  • Enfin, ? Ils ne sont PAS différents avec . p = 0,12β2=β3p=0,12

C'est assez déroutant pour moi, car je m'attends à ce que le global soit inférieur à , car évidemment est un critère beaucoup plus strict que (qui génère ).0,007 β 1 = β 2 = β 3 β 1 = β 3 p = 0,007p0,007β1=β2=β3β1=β3p=0,007

C'est-à-dire, puisque je suis déjà " confiant" que ne tient pas, je devrais être "plus confiant" que ne tient pas. Donc mon devrait descendre.0,007β1=β3β1=β2=β3p

Suis-je les tester à tort? Sinon, où ai-je tort dans le raisonnement ci-dessus?

Sibbs Gambling
la source
Je suppose que x1, x2 et x3 sont différents niveaux d'un facteur similaire, codé fictivement. Ensuite, je pense que de tels résultats surprenants pourraient provenir du nombre différent de répétitions indépendantes (= unités expérimentales) à chaque niveau.
Rodolphe
La période de grâce de la prime touche à sa fin, n'hésitez pas à critiquer ou à demander des précisions si besoin.
brumar

Réponses:

7

C'est-à-dire, puisque je suis déjà "0.007 confiant" que ne tient pas, je devrais être "plus confiant" que β 1 = β 2 = β 3 ne tient pas. Donc mon p devrait descendre β1=β3β1=β2=β3

Réponse courte: Votre probabilité devrait diminuer. Mais ici, les valeurs de p ne mesurent pas la probabilité, mais si la libération de certaines contraintes fournit une amélioration significative de la probabilité. C'est pourquoi il n'est pas nécessairement plus facile de rejeter que de rejeter β 1 = β 3 car il faut montrer de bien meilleures améliorations de vraisemblance dans le modèle le plus contraint pour prouver que la libération de 2 degrés de liberté pour atteindre le modèle complet en "valait la peine".β1=β2=β3β1=β3

Élaboration: dessinons un graphique des améliorations de vraisemblance. graphique de vraisemblance
La seule contrainte pour éviter une contradiction est que les améliorations de vraisemblance doivent être égales à la somme des améliorations de vraisemblance du chemin indirect. Voilà comment j'ai trouvé la valeur de p à l'étape 1 du chemin indirect: Par améliorations de vraisemblance, je veux dire le rapport de vraisemblance logarithmique représenté par leΔChi au carré, c'est pourquoi ils sont additionnés dans le graphique. Avec ce schéma, on peut éliminer la contradiction apparente car une grande partie de l'amélioration de la probabilité du chemin direct provient de la libération d'un seul degré de liberté (β1=β3). Je suggérerais deux facteurs qui peuvent contribuer à ce modèle.

L3L1=L3L2×L2L1
Δβ1=β3
  • a un grand intervalle de confiance dans le modèle completβ2
  • est proche de la moyenne de β 3 et β 1 dans le modèle completβ2β3β1

β3=β1=β2β3=β1β2

β3+β12=β2

brumar
la source