$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Protokół NHANES został zatwierdzony przez Radę ds. Etyki Etyki Badań Narodowego Centrum Statystyki Zdrowia (NCHS), a wszyscy uczestnicy uzyskali pisemną świadomą zgodę. Praca ta była analizą drugorzędną danych z użytku publicznego bez identyfikacji; dlatego nie było wymagane dodatkowe instytucjonalne zatwierdzenie etyczne. Wszyscy autorzy przeczytali i zatwierdzili ostateczny manuskrypt.
1. Projekt badania i źródło danych
Badanie to zostało przeprowadzone jako analiza wtórna NHANES, serii przekrojowych, reprezentatywnych na poziomie krajowym ankiet prowadzonych przez amerykańskie Centra Kontroli i Prewencji Chorób i nadzorowanych przez Radę ds. Etyki Badań NCHS. Publiczne zbiory danych NHANES zostały w pełni zdeidentyfikowane i udostępnione do analizy wtórnej. Wykorzystano dane z cykli 1999–2000, 2001–2002, 2003–2004 oraz 2005–2006.
Dla każdego cyklu pobrano publiczne pliki komponentów NHANES, w tym (i) pliki demograficzne zawierające identyfikator uczestnika (SEQN) i zmienne projektowe badania, (ii) pliki kwestionariuszy zdrowia reprodukcyjnego zawierające samoraport o endometriozie oraz (iii) pliki laboratoryjne wymagane do obliczenia indeksu złożonego (białko C-reaktywne, trójglicerydy i glukoza na czczo). Dodatkowo podczas analizy tych indeksów uzyskano pliki badań/antropometryczne (np. wskaźnik masy ciała) oraz pomiary laboratoryjne wymagane dla indeksów porównawczych (np. neutrofile, limfocyty, płytki krwi). W każdym dwuletnim cyklu pliki komponentów były łączone za pomocą SEQN, a połączony zbiór danych sprawdzano, aby zapewnić jeden rekord na SEQN. Następnie dołączono zbiory danych na poziomie cyklicznym, aby stworzyć połączony plik analityczny z lat 1999–2006.
Próba analityczna była ograniczona do kobiet w wieku 20–54 lat. Uczestniczki były wykluczone, jeśli brakowało statusu endometriozy, brakowało jakichkolwiek składników indeksu złożonego (białko C-reaktywne, trójglicerydy lub glukoza na czczo) lub jeśli w przypadku strategii pełnej sprawy brakowało istotnych kowarianów wymaganych dla w pełni skorygowanego modelu. Zachowano złożone zmienne projektowe ankiety (warstwy i jednostki pierwotnego próbkowania), a także wagi podpróbek laboratoryjnych na czczo wymagane do analiz uwzględniających miary na czczo. Gdy łączono wiele cykli NHANES, według wytycznych analitycznych NHANES utworzono wagi wielocyklowe, dzieląc wagę podpróby 2 lat przez liczbę cykli łączonych, a we wszystkich analizach zastosowano wartości wagi, warstw i PSU. Kroki dotyczące włączenia i wykluczenia uczestników zostały udokumentowane na diagramie przepływowym (Rysunek 1).
2. Definicja endometriozy
Status endometriozy definiowano za pomocą pytania z kwestionariusza zdrowia reprodukcyjnego: "Czy lekarz lub inny pracownik zdrowia kiedykolwiek powiedział Ci, że masz endometriozę?" Uczestniczki, które odpowiedziały "Tak", były klasyfikowane jako osoby z endometriozą, a te, które odpowiedziały "Nie", jako grupy kontrolne. Ponieważ definicja ta opierała się na samoocenach, a nie na potwierdzeniu laparoskopowym czy histologicznym, potencjalna błędna klasyfikacja została potraktowana jako ograniczenie badania.
3. Definicja indeksu złożonego (CTI)
Indeks kompozytowy białko-trójgliceryd-glukoza C-reaktywny został operacjonalizowany, aby wspólnie odzwierciedlać stany zapalne i zaburzenia metaboliczne. Laboratoryjne pomiary białka C-reaktywnego (mg/L), trójglicerydów (mg/dL) oraz glukozy na czczo w osocze (mg/dL) zostały wyekstrahowane z dokumentacji laboratoryjnej NHANES. Indeks trójglicerydów–glukozy został obliczony jako naturalny logarytm [trójglicerydów × glukozy na czczo w osoczu/2]. CTI obliczono według następującego wzoru: CTI = 0,412 × ln(CRP) + TyG. Wyższe wartości CTI wskazują na wyższy łączny ciężar niskiego stopnia zapalenia i insulinooporności8.
Jeśli jakiekolwiek wartości białek C-reaktywnych wymagały obsługi przed transformacją logarytmu (np. wartości na lub poniżej limitu wykrywania), stosowano jedną z góry określoną regułę konsekwentnie we wszystkich cyklach i dokumentowano ją w celu wspierania replikowalności (na przykład zastępowanie wartości niedodatnich najmniejszą dodatnią mierzalną wartością obserwowaną przed transformacją logarytmu). Punkty cięcia kwartyla były wyznaczane na podstawie rozkładu ważonego w pełnej próbie analitycznej i stosowane konsekwentnie we wszystkich analizach kategorycznych, przy czym Kwartyl 1 był używany jako kategoria odniesienia.
4. Kowaryjne
Kowarianity zostały wcześniej określone w celu ograniczenia zakłóceń na podstawie rozumowania epidemiologicznego i wcześniejszej literatury. Zmienne demograficzne obejmowały wiek, rasę/pochodzenie etniczne, poziom wykształcenia oraz stan cywilny. Zmienne stylu życia obejmowały historię palenia (≥100 papierosów w życiu vs. <100) oraz spożycie alkoholu (≥12 drinków/rok vs. <12). Historia współistniejących zachorowań obejmowała samozgłaszane nadciśnienie, cukrzycę, udar, chorobę wieńcową serca oraz nowotwory. Zmienne antropometryczne i laboratoryjne obejmowały wskaźnik masy ciała, hemoglobinę, liczbę neutrofyli, limfocyty oraz liczbę płytek krwi; Te miary wspierały również obliczanie wskaźników zapalnych porównywalnych, gdzie to możliwe (np. stosunek neutrofilów do limfocytów, stosunek płytek krwi do limfocytów, ogólnoustrojowy wskaźnik zapalenia immunologicznego, wskaźnik ogólnoustrojowej odpowiedzi zapalnej). Zmienne rozrodcze (np. gravidity i parity) zostały uwzględnione, gdy były dostępne w wybranych cyklach i kodowane zgodnie z dokumentacją NHANES. Kategoryczne kowarianty były przekształcane w zmienne wskaźnikowe przed wejściem do modelu.
Ponieważ białko C-reaktywne było składnikiem indeksu złożonego, nie było ono wpisywane jako niezależny kowariant w modelach regresji wielowymiarowej, aby uniknąć nadmiernej adaptacji i kolinearyczności. Zamiast tego w analizach dyskryminacyjnych oceniono białko C-reaktywne oraz indeks trójgliceryd-glukozy jako markery porównawcze.
5. Analiza statystyczna
Wszystkie analizy uwzględniały złożony projekt badań NHANES, aby generować reprezentatywne szacunki na poziomie krajowym. Projekt badania został określony poprzez powiązanie wagi podpróby, strat i zmiennych PSU z analitycznym zbiorem danych. Zmienne ciągłe podsumowano jako średnie ważone z odchyleniami standardowymi, a zmienne kategoryczne jako ważone liczby i procenty. Charakterystyki bazowe porównano między przypadkami a kontrolami za pomocą procedur ważonych ankietą odpowiednich dla NHANES, a charakterystyki bazowe podsumowano w Tabeli 1.
Powiązania między indeksem złożonym a endometriozą oceniono za pomocą regresji logistycznej ważonej przez ankietę. Zastosowano trzy modele sekwencyjne, aby wykazać dostosowanie się: model nieskorygowany, model dostosowany pod względem wieku i rasy/pochodzenia etnicznego oraz w pełni skorygowany model uwzględniający czynniki demograficzne, zmienne stylu życia, historię współwystępów, antropometryczne/laboratoryjne współzmienne oraz zmienne historii rozrodu. Indeks złożony był analizowany zarówno ciągłie (na wzrost o 1 jednostkę), jak i kategorycznie (kwartyle, z punktem odniesienia kwartyl 1), a szacunki regresji podsumowano w Tabeli 2. Liniowy trend w kwartylach testowano poprzez przypisanie każdemu kwartylowi jego ważonej wartości mediany i ciągłe modelowanie tego wyrazu.
Nieliniowe relacje dawka–odpowiedź oceniano za pomocą ograniczonych sześciennych spline'ów ważonych ankietą z wcześniej określonym położeniem węzłów, a krzywe spline przedstawiono na Rysunku 2. Efekty progowe oceniono za pomocą segmentowanej regresji logistycznej ważonej przez ankietę, porównując dopasowanie modelu między specyfikacjami segmentowanymi a pojedynczymi, a szacowane parametry punktu zwrotnego i nachylenia po obu stronach punktu zwrotnego zostały przedstawione w Tabeli 3.
Analizy podgrupowe zostały przeprowadzone w celu zbadania modyfikacji efektów według wcześniej określonych czynników (np. grupy wiekowej, rasy/pochodzenia etnicznego, poziomu wykształcenia, stanu cywilnego oraz wybranych czynników stylu życia). Interakcja została przetestowana poprzez uwzględnienie składników międzyproduktowych między ciągłym indeksem złożonym a wskaźnikami podgrup w ramach ważonych przez badania, a powiązania podgrup podsumowano na Rysunku 3.
Dyskryminacyjne wyniki oceniano za pomocą analiz charakterystyk operacyjnych odbiorców opartych na prawdopodobieństwach przewidzianych przez model, pochodzących z modeli logistycznych ważonych przez ankiety. Obszar pod krzywą oszacowane zostały dla indeksu złożonego, powszechnie stosowanych wskaźników zapalnych oraz markerów składowych, a podsumowania AUC przedstawiono w Tabeli Uzupełniającej 1; rozszerzone porównanie ROC zostało przedstawione w Rysunku Uzupełniającym 1.
Brakujące dane zostały obsłużone za pomocą analizy kompletnego przypadku po wykluczeniu uczestników z brakującym statusem endometriozy, brakującymi komponentami indeksu złożonego lub nieistotnymi istotnymi kowarianami wymaganymi dla w pełni dostosowanym modelem. Podczas oceny odporności stosowano wielokrotne imputacje dla kowarianów z brakiem w uprecyzowanym modelu imputacji, a imputowane oszacowania porównywano z pełnymi przypadkami.
Analizy przeprowadzano przy użyciu R (wersja 4.4.1) oraz dodatkowego oprogramowania statystycznego wymienionego w Tabeli Materiałów. Rejestrowano kluczowe pakiety używane do wnioskowania ankietowego, modelowania spline, regresji segmentowanej oraz szacowania ROC, a informacje o sesji (szczegóły systemu operacyjnego i sesji R) zostały zachowane w celu wsparcia replikacji.
6. Punkt końcowy procedury i wyjścia
Analityczny przepływ pracy uznano za kompletny po stworzeniu zharmonizowanego wielocyklowego zbioru danych z wcześniej określonymi kryteriami włączenia/wykluczenia (Rysunek 1), wygenerowano indeks złożony i kowarianity zgodnie z udokumentowanymi regułami, a regresja ważona ankietą, ocena nieliniowości/progu, podgrupy oraz analiza dyskryminacji została wykonana według tej samej specyfikacji projektu badania. Główne wyniki tego przepływu pracy zostały zorganizowane jako podsumowanie bazowe (Tabela 1), szacunki regresji dla modeli korekt sekwencyjnych (Tabela 2), parametry modelu progowego (Tabela 3), wizualizację spline (Rysunek 2), wizualizację podsumowań podgrup (Rysunek 3) oraz podsumowania dyskryminacji (Tabela Uzupełniająca 1 i Rysunek Uzupełniający 1).
7. Wewnętrzna niezależna walidacja
Wewnętrzna niezależna walidacja przeprowadzono poprzez podział połączonego zbioru danych na kohortę derywacyjną oraz kohortę walidacyjną nienakładającą się na podstawie cykli NHANES. Uczestnicy z cykli 1999–2000 i 2001–2002 zostali przypisani do kohorty derywacyjnej, a uczestnicy z cykli 2003–2004 oraz 2005–2006 do kohorty walidacyjnej. Te same kryteria włączenia/wykluczenia, obliczenia indeksu złożonego, reguły kodowania kowarianckiego oraz strategia ważenia ankiet były stosowane niezależnie w każdej kohorcie.
W kohorcie wyprowadzającej dostosowano modele regresji logistycznej ważone według ankiet, korzystając z w pełni dostosowanej specyfikacji. Procedury oceny nieliniowości i progowej stosowane w głównej analizie zastosowano w kohorcie wyprowadzania, a rozróżnienie oceniano metodami ROC/AUC opartymi na prawdopodobieństwach przewidzianych przez model. Ta sama strategia modelowania była następnie powtarzana w kohorcie walidacyjnej bez zmiany definicji zmiennych, reguł kodowania czy specyfikacji wagowania. Szacunki wyprowadzania względem walidacji podsumowano jako porównanie szacunków asocjacji skojarzeń kohorta do kohorty (rysunek 4) oraz jako krzywe ROC pochodzenia względem walidacji (rysunek 5), z odpowiadającymi im podsumowaniami liczbowymi w Tabeli 4.