Comment puis-je regrouper les valeurs p amorcées dans plusieurs ensembles de données imputées?

Je suis préoccupé par le problème que j'aimerais amorcer la valeur de p pour une estimation de partir de données multipliées imputées (MI), mais qu'il n'est pas clair pour moi comment combiner les valeurs de p entre les ensembles d'IM.$\theta$

Pour les ensembles de données MI, l'approche standard pour arriver à la variance totale des estimations utilise les règles de Rubin. Voir ici pour un examen de la mise en commun des ensembles de données MI. La racine carrée de la variance totale sert d'estimation d'erreur standard de . Cependant, pour certains estimateurs, la variance totale n'a pas de forme fermée connue ou la distribution d'échantillonnage n'est pas normale. La statistique θ / s e ( θ $\theta$${\theta}/{se(\theta)}$ peut alors ne pas être distribuée t, pas même asymptotiquement.

Par conséquent, dans le cas des données complètes, une option alternative consiste à amorcer la statistique pour trouver la variance, une valeur de p et un intervalle de confiance, même si la distribution d'échantillonnage n'est pas normale et sa forme fermée inconnue. Dans le cas MI, il y a alors deux options:

• Mettre en commun la variance amorcée entre les ensembles de données MI
• Mettre en commun la valeur de p ou les limites de confiance entre les ensembles de données MI

La première option utiliserait alors à nouveau les règles de Rubin. Cependant, je pense que cela pose problème si a une distribution d'échantillonnage non normale. Dans cette situation (ou plus généralement, dans toutes les situations), la valeur de p amorcée peut être utilisée directement. Cependant, dans le cas de l'IM, cela conduirait à plusieurs valeurs de p ou intervalles de confiance, qui doivent être regroupés entre les ensembles de données de l'IM.$\theta$

Ma question est donc la suivante: comment regrouper plusieurs valeurs de p (ou intervalles de confiance) amorcées sur plusieurs ensembles de données imputées?

Je serais heureux de recevoir toute suggestion sur la façon de procéder, merci.

Je pense que les deux options donnent la bonne réponse. En général, je préférerais la méthode 1 car elle préserve la distribution entière.

$k$$m$$m$$k \times m$ échantillons qui incluent la variation entre imputations. Traitez ensuite cela comme un échantillon bootstrap conventionnel pour obtenir des intervalles de confiance. Utilisez le bootstrap bayésien pour les petits échantillons. Je ne connais aucun travail de simulation qui étudie cette procédure, et c'est en fait un problème ouvert à étudier.

Pour la méthode 2, utilisez la procédure Licht-Rubin. Voir Comment obtenir des valeurs de p regroupées sur des tests effectués dans plusieurs jeux de données imputés?

Ce n'est pas une littérature que je connais, mais une façon d'aborder cela pourrait être d'ignorer le fait qu'il s'agit de valeurs p bootstrapées et de regarder la littérature sur la combinaison de valeurs p à travers des ensembles de données imputées multipliées.

Dans ce cas, Li, Meng, Raghunathan et Rubin (1991) s'appliquent. La procédure est basée sur des statistiques de chacun des ensembles de données imputées, pondérées à l'aide d'une mesure de la perte d'informations due à l'imputation. Ils rencontrent des problèmes liés à la distribution conjointe des statistiques entre les imputations et émettent des hypothèses simplificatrices.

Meng (1994) présente un intérêt connexe .

Mise à jour

Une procédure pour combiner des valeurs de p à travers des ensembles de données imputées à plusieurs reprises est décrite dans la dissertation de Christine Licht, Ch. 4 . L'idée, qu'elle attribue à Don Rubin, est essentiellement de transformer les valeurs de p à distribuer normalement, qui peuvent ensuite être combinées entre les ensembles de données MI en utilisant les règles standard pour la combinaison de statistiques z.