On Anticipation Effect in Stepped Wedge Cluster Randomized Trials
Wang, Chen, Courtright et al.
In stepped wedge cluster randomized trials (SW-CRTs), the intervention is rolled out to clusters over multiple periods. A standard approach for analyzing SW-CRTs utilizes the linear mixed model, where the treatment effect is only present after the treatment adoption, under the assumption of no anticipation. This assumption, however, may not always hold in practice because stakeholders, providers, or individuals who are aware of the treatment adoption timing (especially when blinding is challenging or infeasible) can inadvertently change their behaviors in anticipation of the forthcoming intervention. We provide an analytical framework to address the anticipation effect in SW-CRTs and study its impact. We derive expectations of the estimators based on a collection of linear mixed models and demonstrate that when the anticipation effect is ignored, these estimators give biased estimates of the treatment effect. We also provide updated sample size formulas that explicitly account for anticipation effects, exposure-time heterogeneity, or both in SW-CRTs and illustrate their impact on study power. Through simulation studies and empirical analyses, we compare the treatment effect estimators with and without adjusting for anticipation, and provide some practical considerations.
academic
Sur l'Effet d'Anticipation dans les Essais Randomisés en Grappe à Coin Progressif
Dans les essais randomisés en grappe à coin progressif (SW-CRT), l'intervention est progressivement mise en œuvre dans plusieurs grappes au cours de plusieurs périodes. Les méthodes d'analyse standard utilisent des modèles linéaires mixtes, en supposant que l'effet du traitement n'apparaît que après la mise en œuvre de l'intervention et qu'aucun effet d'anticipation n'existe. Cependant, cette hypothèse peut ne pas tenir en pratique, car les parties prenantes, prestataires ou individus conscients du calendrier d'intervention (particulièrement lorsque l'insu est difficile ou impossible) peuvent involontairement modifier leur comportement en anticipant l'intervention imminente. Cet article fournit un cadre d'analyse pour traiter les effets d'anticipation dans les SW-CRT et étudier leurs impacts. Les auteurs dérivent l'espérance des estimateurs basés sur des modèles linéaires mixtes et prouvent que ces estimateurs produisent des estimations biaisées de l'effet du traitement lorsque l'effet d'anticipation est ignoré. L'article fournit également des formules de calcul de la taille d'échantillon qui considèrent explicitement les effets d'anticipation, l'hétérogénéité du temps d'exposition ou les deux, et illustre leurs impacts sur la puissance de l'étude.
L'essai randomisé en grappe à coin progressif (SW-CRT) est un type particulier de conception d'essai où toutes les grappes commencent initialement dans une condition de contrôle, puis reçoivent progressivement l'intervention selon un ordre aléatoire au cours de plusieurs périodes. L'analyse traditionnelle suppose que l'effet du traitement n'apparaît que après la mise en œuvre de l'intervention, sans effet d'anticipation (anticipation effect).
Prévalence en pratique: Dans de nombreux SW-CRT, l'insu présente des défis ou est impossible, et les participants peuvent connaître à l'avance le calendrier d'intervention
Modification du comportement: Les membres des grappes conscients de l'intervention imminente peuvent modifier leur comportement, produisant un effet d'anticipation
Biais d'estimation: Ignorer l'effet d'anticipation peut conduire à des biais graves dans l'estimation de l'effet du traitement
Hypothèse d'absence d'anticipation: Pratiquement tous les travaux antérieurs supposent l'absence d'effet d'anticipation
Erreur de spécification du modèle: La littérature existante se concentre principalement sur l'hétérogénéité du temps d'exposition, sans étudier systématiquement l'effet d'anticipation
Calcul de la taille d'échantillon: Les formules de taille d'échantillon existantes ne considèrent pas l'impact de l'effet d'anticipation
Cadre théorique: Fournit pour la première fois un cadre d'analyse systématique pour l'effet d'anticipation dans les SW-CRT, définissant quatre modèles linéaires mixtes (HH, HH-ANT, ETI, ETI-ANT)
Formules de biais: Dérive des formules de biais en forme fermée pour 16 combinaisons de scénarios de modèle vrai et modèle de travail, quantifiant l'impact d'ignorer l'effet d'anticipation et/ou l'hétérogénéité du temps d'exposition
Formules de variance: Fournit des expressions de variance pour l'estimateur de l'effet du traitement considérant l'effet d'anticipation, utilisées pour le calcul de la taille d'échantillon
Méthode de taille d'échantillon: Développe des formules de taille d'échantillon considérant explicitement l'effet d'anticipation, l'hétérogénéité du temps d'exposition ou les deux
Preuve empirique: Valide les résultats théoriques par des études de simulation et une réanalyse de l'essai REDAPS, fournissant des recommandations pratiques
Outil logiciel: Fournit une fonction R (Trt.Ant.Power) pour l'analyse de puissance considérant l'effet d'anticipation
3. Modèle ETI (Exposure Time Indicator) - Considérant l'hétérogénéité du temps d'exposition, sans anticipation:
Yijk=μ+αj+θ(sij)Zij+bi+ϵijk
4. Modèle ETI-ANT - Considérant simultanément l'effet d'anticipation et l'hétérogénéité du temps d'exposition:
Yijk=μ+αj+γAij+θ(sij)Zij+bi+ϵijk
Où:
μ: moyenne globale
αj: paramètre de tendance à long terme pour la période j
θ ou θ(s): effet du traitement
γ: effet d'anticipation
bi∼N(0,τ2): effet aléatoire de grappe
ϵijk∼N(0,σ2): résidu
Définition de l'indicateur d'anticipation:
Aij={1,0,j∗−l≤j<j∗sinon
où j∗ est la période d'adoption du traitement et l est l'ordre de l'effet d'anticipation (l'article considère principalement l=1).
où les constantes de conception dépendent de l'allocation du traitement et de l'indicateur d'anticipation.
Découvertes importantes:
V(θ^HH)<V(θ^HH−ANT)
L'inflation de variance augmente avec l'ICC et diminue avec le nombre de périodes
Pour J=9 périodes, l'inflation de variance peut atteindre 40-50%
4. Calcul de la Taille d'Échantillon et de la Puissance
Basé sur les formules de variance, fournit la formule de puissance pour le test de Wald bilatéral:
Puissance=Φ{V(θ^)∣θ∗∣−z1−α/2}
Cependant, il faut noter que lorsque le modèle de travail est mal spécifié, θ^ converge vers une valeur biaisée plutôt que vers le vrai effet du traitement, ce qui entraîne une inflation de l'erreur de type I.
Test d'hétérogénéité du temps d'exposition du modèle ETI: p=0.15
Test d'hétérogénéité du temps d'exposition du modèle ETI-ANT: p=0.36
Conclusion: Pas de preuve forte soutenant l'hétérogénéité du temps d'exposition, se concentrer sur HH et HH-ANT
Estimations de l'effet du traitement:
Modèle HH: θ^=−1.23 (IC 95%: -1.52 à -0.93)
Modèle HH-ANT: θ^=−1.43 (IC 95%: -1.76 à -1.10)
Différence: L'estimation HH-ANT de l'effet absolu est plus grande (réduction du temps de consultation de 1.43 jours vs 1.23 jours)
Estimations de l'effet d'anticipation:
Modèle HH-ANT: γ^=−0.59 (IC 95%: -1.01 à -0.17), statistiquement significatif
Interprétation: Le temps de consultation diminue de 0.59 jours avant l'intervention, indiquant que les cliniciens ont déjà ajusté leur comportement avant la mise en œuvre formelle
Signification pratique:
Ignorer l'effet d'anticipation négatif entraîne une sous-estimation de l'effet du traitement d'environ 14% (1.23 vs 1.43)
La préparation éducative a produit des bénéfices précoces
L'effet s'est renforcé après la mise en œuvre formelle
Wang et al. (2024): Prouvent que le modèle linéaire mixte est robuste à l'erreur de spécification des covariables, des effets aléatoires et de la structure d'erreur, mais exige une spécification stricte de la structure d'effet de traitement
Chen & Li (2025): Proposent l'estimation des effets de variation du temps calendaire et de la taille de grappe informative
Systématicité: Énumère pour la première fois 16 combinaisons de scénarios de modèles, fournissant un cadre théorique complet
Rigueur: Dérive des formules de biais et de variance en forme fermée, avec des dérivations mathématiques détaillées (50 pages d'appendice)
Innovativité: Découvre que les scénarios de non-correspondance peuvent être pires que la sous-spécification, remettant en question la sagesse conventionnelle
2. Valeur Méthodologique:
Opérabilité: Fournit un package R et une fonction de calcul de taille d'échantillon
Praticabilité: Basé sur le modèle linéaire mixte standard, facile à mettre en œuvre
Généralité: Le cadre peut s'étendre à d'autres conceptions (SW-CRT incomplètes)
3. Support Empirique:
Simulation suffisante: 4 études de simulation couvrant les scénarios clés, 2000 répétitions chacune
Données réelles: L'analyse de l'essai REDAPS fournit des preuves empiriques convaincantes
Visualisation: 20+ figures présentent clairement les modèles de biais et les fonctions de poids
4. Qualité de Rédaction:
Structure claire, logique rigoureuse de la motivation à la théorie à l'application
Équilibre bien entre détails techniques et explications intuitives
Le Tableau 2 fournit une excellente feuille de route de recherche
Structure de corrélation: Considère uniquement la structure échangeable, peut ne pas être réaliste en pratique
Résultats continus: Pas de solution en forme fermée pour les résultats binaires (comme la mortalité), limitant l'applicabilité
Conception standard: Les résultats principaux sont basés sur les SW-CRT standard, les allocations inégales et les conceptions irrégulières nécessitent une recherche supplémentaire
2. Insuffisances Empiriques:
Cas unique: Une seule application de données réelles, nécessite plus de preuves empiriques
Détection d'anticipation: Ne fournit pas de test statistique pour l'existence de l'effet d'anticipation
Analyse de sensibilité: N'explore pas suffisamment la sensibilité à l'ICC et à la taille de l'effet d'anticipation mal spécifiée
3. Orientations Pratiques:
Règles de décision: Manque de critères clairs pour savoir quand considérer l'effet d'anticipation au stade de la conception
Analyse coûts-bénéfices: Ne discute pas du compromis entre l'inflation de variance et la réduction du biais
Stratégies d'insu: Moins de recommandations sur comment éviter l'effet d'anticipation en pratique
4. Lacunes Théoriques:
Cadre causal: Ne formalise pas l'effet d'anticipation dans le cadre des résultats potentiels
Dépendance temporelle: Ne considère pas l'effet d'anticipation variant dans le temps
Interactions: N'explore pas les interactions entre l'effet d'anticipation et les caractéristiques individuelles
Hussey MA, Hughes JP (2007). Design and analysis of stepped wedge cluster randomized trials. Contemporary Clinical Trials, 28(2):182-191.
Article fondateur des SW-CRT
Kenny A, et al. (2022). Analysis of stepped wedge cluster randomized trials in the presence of a time-varying treatment effect. Statistics in Medicine, 41(22):4311-4339.
Travail précurseur important sur l'hétérogénéité du temps d'exposition
Maleyeff L, et al. (2022). Assessing exposure-time treatment effect heterogeneity in stepped-wedge cluster randomized trials. Biometrics, 79(3):2551-2564.
Modèle ETI à effets aléatoires
Courtright KR, et al. (2024). Default palliative care consultation for seriously ill hospitalized patients: a pragmatic cluster randomized trial. JAMA, 331(3):224-232.
Résultats principaux de l'essai REDAPS
Hemming K, Taljaard M (2020). Reflection on modern methods: when is a stepped-wedge cluster randomized trial a good study design choice? International Journal of Epidemiology, 49(3):1043-1052.
Guide de choix de conception pour SW-CRT
Évaluation Globale: Cet article est une contribution méthodologique de haute qualité qui comble une lacune importante dans l'analyse des SW-CRT. La théorie est rigoureuse, les preuves empiriques sont suffisantes et la valeur pratique est élevée. Les contributions principales résident dans la quantification systématique de l'impact de l'effet d'anticipation et la fourniture de solutions opérationnelles. Les limitations résident principalement dans la restrictivité des hypothèses du modèle et l'étendue des preuves empiriques. Il est recommandé d'appliquer prudemment en pratique, en combinant les connaissances du domaine pour évaluer la probabilité d'effet d'anticipation et en effectuant une analyse de sensibilité.