Cette étude a été menée conformément à la Déclaration d’Helsinki et approuvée par le Comité d’éthique de la recherche de l’Hôpital général du Hebei (numéro d’approbation IRB : LW-104). En raison de la nature rétrospective de l’étude et de l’utilisation d’informations anonymes des patients, le comité d’éthique de l’hôpital général du Hebei a levé l’exigence de consentement éclairé. Les outils de recherche utilisés dans ce protocole sont listés dans le tableau des matériaux.
1. Conception de l’étude
Cette enquête constituait une étude caso-témoin rétrospective unicentre, menée à l’hôpital général du Hebei. Les cas comprenaient des personnes ayant subi un premier AVC ischémique, et les témoins ont été recrutés parmi des patients d’examen médical contemporains ou des patients externes non cérébrovasculaires du même établissement. L’appariement individuel a été réalisé par âge (±3 ans) et sexe, avec un ratio cas/témoins de 1:2. Lorsque l’appariement individuel était irréalisable, le raccordement fréquentiel était utilisé comme stratégie complémentaire, les facteurs d’appariement étant intégrés dans les modèles statistiques ultérieurs.
Les sources de données comprenaient les dossiers médicaux électroniques, les systèmes d’archivage et de communication d’images pour la neuroimagerie, les systèmes de gestion des informations de laboratoire et les bases de données d’examens de santé. Pour garantir une taille d’échantillon adéquate et un nombre suffisant d’événements par variable (EPV) pour l’inférence statistique, une période de collecte des données du 1er janvier 2018 au 30 juin 2025 a été définie. Cette période a été calibrée pour produire au moins 300 cas et 600 témoins, répondant aux exigences statistiques prédéterminées et aux seuils de VE pour le modèle analytique primaire.
2. Participants à l’étude
Les cas ont été définis comme des individus subissant un premier AVC ischémique, avec des manifestations cliniques compatibles avec un déficit neurologique focal aigu de l’étiologie vasculaire, des preuves de neuroimagerie confirmatoires via tomodensitométrie (CT) ou imagerie par résonance magnétique (IRM) avec imagerie pondérée par diffusion (DWI) démontrant un infarctus cérébral aigu, et aucun antécédent antérieur d’AVC ischémique basé sur la revue du dossier médical et l’entretien par procuration du patient ou par procuration. Les sous-types d’AVC ischémique ont été classés selon les critères de l’essai Org 10172 dans le traitement de l’AVC aigu (TOAST) par des neurologues jugeants.
Les critères d’exclusion comprenaient un AVC hémorragique (hémorragie intracérébrale ou hémorragie sous-arachnoïdienne) ou une crise ischémique transitoire (AIT) sans infarctus à l’imagerie, une thrombose du sinus veineux cérébral, ainsi que de nombreuses données manquantes sur des variables critiques (>30 % de manque sur des mesures clés d’exposition ou de covariables qui ne pouvaient pas être imputées de manière fiable).
Deux neurologues certifiés ont jugé chaque cas potentiel de manière indépendante, les divergences étant résolues par consultation avec un troisième neurologue senior. Ce processus d’adjudication standardisé garantissait la cohérence diagnostique et minimisait les biais de mauvaise classification. La date d’indice pour les cas était définie comme la date d’apparition des symptômes d’AVC ou, lorsque le moment précis d’apparition n’était pas disponible, la date de présentation médicale initiale avec symptômes d’AVC.
Les témoins ont été prélevés auprès de personnes présentant le centre d’examen de santé de l’hôpital ou présentant des consultations externes non cérébrovasculaires durant la même période d’étude. Les critères d’éligibilité exigeaient l’absence de tout antécédents de maladie cérébrovasculaire, y compris un AVC ou un AIT, tel que documenté dans les dossiers médicaux et confirmé par un entretien structuré. Les témoins ont été individuellement appariés selon l’âge (±3 ans) et le sexe selon un ratio de 1:2. Lorsque l’appariement individuel n’était pas possible, l’appariement fréquentiel a été utilisé pour garantir une répartition globale d’âge et de sexe comparable entre les groupes. Dans de tels cas, les facteurs d’appariement étaient explicitement inclus comme covariables dans les modèles de régression logistique conditionnelle ou inconditionnelle pour contrôler la confusion résiduelle.
Des recherches antérieures ont documenté des problèmes omniprésents de taille d’échantillon insuffisante et d’EPV insuffisant dans des études épidémiologiques utilisant la régression logistique25. Pour répondre à ces limitations, la règle empirique de Peduzzi, recommandant un EPV minimum de 10, a été appliquée, avec des directives contemporaines suggérant un EPV ≥ 15 pour les modèles incorporant des termes d’interaction. Le modèle analytique principal a été spécifié a priori pour inclure environ 12 paramètres, comprenant des variables d’exposition catégorielles telles que le statut de tabagisme et la consommation d’alcool, des facteurs de confusion prédéfinis tels que l’âge avec splines, l’indice de masse corporelle (IMC), la pression artérielle systolique (SBP) avec splines, le diabète, la dyslipidémie, l’utilisation de médicaments, la fonction rénale et l’homocystéine, ainsi qu’un terme d’interaction principal entre hypertension et tabagisme actuel.
Les catégories par année de paquet et les catégories de statut de fumeur n’étaient pas incluses simultanément dans le modèle primaire pour éviter la linéarité ; au lieu de cela, ceux-ci ont été examinés dans des modèles séparés. L’inclusion simultanée du diagnostic d’hypertension et de la spline SBP a été justifiée car la variable binaire hypertension capture le statut de traitement, tandis que la spline SBP modélise la relation continue, potentiellement non linéaire, entre la pression artérielle mesurée et le risque d’AVC. L’évaluation du facteur d’inflation de la variance (VIF) a confirmé une colinéarité acceptable (VIF < 3,5 pour les deux variables). En appliquant le seuil conservateur de l’EPV ≥ 15, le nombre minimum de cas requis de 180 a été calculé.
La taille cible de l’échantillon a été fixée à ≥300 cas avec ≥600 contrôles appariés. Cette taille d’échantillon fournissait également une puissance statistique suffisante pour détecter les effets principaux et les interactions clés. En utilisant la méthode Hsieh26 pour les conceptions cas-témoins appariées, des calculs formels de puissance ont été réalisés pour l’exposition primaire au tabagisme, en supposant une prévalence d’exposition témoin de 30 % pour le tabagisme actuel, un rapport de chances (OR) anticipé de 1,6, α = 0,05 (bilatéral), un rapport de correspondance de 1:2, et une puissance cible de 90 %.
3. Variables et mesures
Le critère principal a été la première apparition d’un AVC ischémique (oui/non). La détermination des cas nécessitait une évaluation par un double neurologue basée sur la présentation clinique et la neuroimagerie confirmatoire. Le moment de l’AVC était défini comme la date d’apparition des symptômes ou, lorsque l’absence possible, la date de la première évaluation médicale documentant des déficits neurologiques aigus.
Les variables d’exposition comprenaient le tabagisme et la consommation d’alcool. Le comportement tabagiste était classé comme n’ayant jamais fumé, ancien fumeur ou fumeur actuel. Ne jamais fumer était défini comme une consommation à vie de moins de 100 cigarettes, un arrêt de fumée ancienne au moins 6 mois avant l’indice ou la date de référence, et le tabagisme actuel comme fumeur actif dans les 6 mois suivant l’indice ou la date de référence. L’intensité du tabagisme était quantifiée en années-paquets, calculées comme (moyenne de cigarettes par jour ÷ 20) × années de tabagisme. Des catégories ordinales (<10, 10–20, >20 années-paquet) ont été utilisées pour les analyses dose-réponse.
La consommation d’alcool a été standardisée en grammes d’éthanol pur par semaine (g/semaine), sur la base des hypothèses standards de teneur en éthanol pour la bière (5 %), le vin (12 %) et les spiritueux (40 %). Une boisson standard était définie comme environ 10 g d’éthanol. Les niveaux de consommation étaient classés comme nul (0 g/semaine), légers à modérés (1–100 g/semaine) et élevés (>100 g/semaine).
Toutes les données d’exposition ont été collectées par entretiens structurés complétés par une abstraction des dossiers médicaux. Dans les cas, les expositions reflétaient des schémas habituels l’année précédant le début de l’AVC ; Pour les témoins, les expositions correspondaient à la période précédant l’examen médical ou la visite en clinique. Lorsque cela était nécessaire, des entretiens téléphoniques de suivi avec les patients ou les mandataires étaient menés.
La Figure 1 présente un graphe orienté acyclique illustrant les relations hypothétées entre expositions, facteurs de confusion et résultats, qui ont influencé la sélection des covariables. L’ensemble de covariables pré-spécifiées comprenait l’âge, le sexe, l’IMC, la SBP, le diabète sucré, la dyslipidémie, l’utilisation de médicaments, le profil lipidique (LDL-C, HDL-C), la fonction rénale (eGFR), l’homocystéine, la fibrillation auriculaire et les antécédents familiaux. Les variables étaient codées de manière cohérente : 1 indique la présence ou l’élévation, et 0 indique l’absence ou les niveaux normaux. Le tableau 1 fournit des définitions détaillées et des spécifications de mesure pour les variables.

Figure 1 : Graphe acyclique orienté (DAG). Ce schéma présente les relations prédéfinies entre le tabagisme, la consommation d’alcool, les facteurs de risque cardiovasculaires traditionnels, les facteurs de confusion mesurés et l’AVC ischémique. Les flèches dirigées indiquent la structure analytique supposée utilisée pour identifier l’ensemble minimal d’ajustement suffisant avec l’algorithme de la dagitty. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
| Variable | Définition et mesure | Encodage | Notes |
| AVC ischémique | Toute première de l’histoire, imagerie confirmée | 1/0 | Résultat des primaires |
| Statut fumeur | Jamais/Ancien/Actuel | 0/1/2 | Aussi des années de meute record |
| Années de meute | (Cigarettes/jour ÷ 20) × ans | Continu | Modélisation RCS |
| Alcool (g/semaine) | Grammes d’éthanol purs/semaine | Continu | Catégories : 0, 1-100, >100 |
| TA systolique | mmHg | Continu | Modélisation RCS |
| Hypertension | TA ≥140/90 ou médicaments | 1/0 | Binaire |
| Diabète | Diagnostic ou médicaments | 1/0 | Binaire |
| Dyslipidémie | Diagnostic ou médicaments | 1/0 | Binaire |
| LDL-C | mmol/L | Continu | Plus haut = pire |
| HDL-C | mmol/L | Continu | Plus bas = pire |
| eGFR | CKD-EPI, mL/min/1,73 m² | Continu | Plus bas = pire |
| Homocystéine | μmol/L | Continu | Plus haut = pire |
| Utilisation des statines | Historique des médicaments | 1/0 | Binaire |
| Utilisation antihypertensive | Historique des médicaments | 1/0 | Binaire |
| HTN × Smoking | Terme produit | -- | Interaction primaire |
| La tension artérielle, la tension artérielle ; CKD-EPI, Collaboration sur l’épidémiologie des maladies rénales chroniques ; HDL-C, cholestérol lipoprotéiné à haute densité ; HTN, hypertension ; LDL-C, cholestérol lipoprotéiné de basse densité ; RCS, splines cubiques restreintes |
| (Uniforme 1=Présent/Haut, 0=Absent/Bas) |
Tableau 1 : Dictionnaire de codage variable (Uniforme 1=Présent/Haut, 0=Absent/Bas). Ce tableau définit le cadre de codage utilisé pour toutes les variables binaires, catégoriques et continues et documente les règles de mesure harmonisées appliquées lors de l’abstraction des données.
4. Sources de données, nettoyage et gestion des données manquantes
Pour minimiser les erreurs de transcription, deux chercheurs indépendants ont réalisé une abstraction parallèle des données pour tous les cas et un échantillon aléatoire de 10 % de témoins, les divergences étant résolues par un arbitrage tiers. L’assurance qualité des données comprenait des contrôles de plage pour des valeurs invraisemblables, des vérifications de cohérence logique et une résolution de doublons d’enregistrements. La fiabilité inter-évaluateurs a été évaluée à l’aide de coefficients de corrélation intraclasse pour les variables continues et des statistiques kappa de Cohen pour les variables catégorielles, toutes les valeurs dépassant 0,85.
Les données manquantes ont été traitées par imputation multiple par équations enchaînées (MICE). L’appariement prédictif des moyennes a été utilisé pour les variables continues, la régression logistique pour les variables binaires, et la régression logistique multinomiale ou ordinale pour les variables catégorielles selon le cas. Vingt ensembles de données imputés ont été générés pour améliorer la précision de l’estimation. Le modèle d’imputation comprenait toutes les variables d’analyse, les variables auxiliaires liées à l’absence et le résultat. Les variables dérivées, y compris les termes d’interaction, étaient imputées passivement pour maintenir la cohérence.
5. Analyse statistique
La bonté de l’ajustement du modèle et la performance prédictive ont été évaluées à l’aide de plusieurs métriques. La multicolinéarité a été évaluée à l’aide de facteurs d’inflation de variance (VIF), le VIF >10 indiquant une multicolinéarité problématique nécessitant une correction par réduction variable ou pénalisation de la crête. La discrimination a été quantifiée à l’aide de la surface sous la courbe (AUC). Il convient de noter que, puisqu’il s’agissait d’une étude cas-témoins appariée, les risques absolus ne pouvaient pas être estimés directement, et l’AUC a été dérivé d’un modèle de régression logistique inconditionnelle incluant l’âge, le sexe et toutes les covariables pré-spécifiées comme évaluation secondaire de la capacité discriminatoire plutôt qu’une revendication d’étalonnage au niveau de la population. Le score de Brier rapporté reflète la différence quadratique moyenne entre les probabilités prédites et le statut cas-témoins observé dans l’échantillon analytique et est présenté comme une mesure relative de performance du modèle plutôt qu’un indicateur calibré au niveau de la population, étant donné que le ratio cas/contrôle de 1:2 ne reflète pas la véritable prévalence de la maladie. L’analyse de la courbe de décision a été réalisée comme évaluation secondaire ; Les métriques numériques et les plages de seuils sont indiquées dans le texte.
Concernant les mesures d’interaction additive (RERI, AP, S), il est reconnu que ces indicateurs sont formellement définis en termes de risques ou de risques relatifs (RR). Dans le design actuel cas–témoin apparié, les rapports de chances (OR) estimés à partir de la régression logistique conditionnelle servent d’approximations aux RR sous l’hypothèse de maladie rare. D’après les données de la population source (3 847 cas cérébrovasculaires potentiels dépistés auprès d’une population clinique plus large avec 18 456 témoins éligibles), la prévalence brute des AVC dans ce cadre institutionnel était d’environ 17 % à 18 %. Bien que cela dépasse le seuil conventionnel pour l’hypothèse de maladie rare, des travaux méthodologiques récents ont montré que les mesures d’interaction additive dérivées des OR restent informatives et directionnellement cohérentes avec celles basées sur les RR même lorsque la prévalence de la maladie est modérée, bien que l’ampleur des RERI puisse être quelque peu surestimée17,18. Les OR conditionnels de l’analyse appariée ont été utilisés pour calculer RERI, AP et S, avec des intervalles de confiance démarrés (1 000 réplications) pour fournir une inférence valide. Ces résultats doivent être interprétés comme des mesures approximatives de la force de l’interaction additive plutôt que comme des attributions exactes de risque au niveau de la population.
Les six cas appariés en fréquence (1,9 % du total) ont été traités en incluant des variables d’appariement (âge et sexe) comme covariables dans le modèle de régression logistique conditionnelle. Une analyse de sensibilité excluant ces six cas et leurs témoins appariés a donné des résultats pratiquement identiques (données non montrées), confirmant que ce léger écart par rapport à l’appariement individuel n’a pas influencé les conclusions.
Des observations influentes et des valeurs aberrantes ont été identifiées à l’aide des statistiques delta–bêta et de la distance de Cook, et des analyses de sensibilité excluant les valeurs extrêmes ont été réalisées pour évaluer les résultats. Les analyses de sensibilité pré-spécifiées suivantes ont été réalisées : réanalyse en utilisant uniquement le fumeur actuel (contre jamais ou ancien combiné) comme exposition binaire ; l’évaluation des associations avec les sous-types d’AVC ischémique classés selon les critères TOAST lorsque les données le permettaient ; comparaison de modèles utilisant différents nombres et positions de nœuds à splines cubiques restreintes (RCS) pour des covariables continues ; analyse complète de cas restreinte aux individus disposant de données complètes sur toutes les variables du modèle ; suppression séquentielle des observations avec des valeurs de covariables extrêmes (au-delà des 1er et 99e percentiles) ; stratification par sexe et catégories d’âge (<60 et ≥60 ans), incluant une analyse distincte des jeunes patients victimes d’AVC (≤55 ans) ; et, si un déséquilibre résiduel de covariables persistait après correspondance (différence moyenne standardisée ≥0,10), l’inverse de la probabilité de pondération du traitement (IPTW) en utilisant les scores de propension pour les expositions clés comme analyses auxiliaires.
Concernant l’estimation du score de propension dans ce modèle cas-témoins, le modèle de propension a été adapté pour prédire la probabilité de tabagisme actuel (exposition), conditionnée aux covariables avant exposition, plutôt que de prédire le statut cas–témoins ; par conséquent, le ratio d’échantillonnage de 1:2 n’affecte pas la validité du score de propension. Toutes les analyses ont été réalisées avec la version R 4.3.x ou ultérieure (R Foundation for Statistical Computing, Vienne, Autriche). Les packages clés comprenaient survival (régression logistique conditionnelle), logistf (régression logistique pénalisée de Firth), rms (RCS et diagnostic de modèles), epiR (mesures d’interaction additives), mice (imputation multiple), boot (intervalles de confiance bootstrap) et ggplot2 (visualisation des données).