$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
NHANES protokolü, Ulusal Sağlık İstatistikleri Merkezi (NCHS) Araştırma Etik İnceleme Kurulu tarafından onaylanmış ve tüm katılımcılardan yazılı bilgilendirilmiş onayı alınmıştır. Bu çalışma, kimliği silinmiş kamu kullanımına ait verilerin ikincil bir analiziydi; bu nedenle ek kurumsal etik onay gerekmiyordu. Tüm yazarlar nihai makaleyi okuyup onayladı.
1. Çalışma tasarımı ve veri kaynağı
Bu çalışma, ABD Hastalık Kontrol ve Önleme Merkezleri tarafından yönetilen ve NCHS Araştırma Etik İnceleme Kurulu tarafından denetlenen, ulusal çapraz ve ulusal çapraz anketler dizisi olan NHANES'in ikincil analizi olarak gerçekleştirilmiştir. Kamu kullanımına yönelik NHANES veri setleri tamamen kimliği çıkarıldı ve ikincil analiz için erişildi. 1999–2000, 2001–2002, 2003–2004 ve 2005–2006 döngülerinden veriler kullanılmıştır.
Her döngü için kamuya açık NHANES bileşen dosyaları indirildi; bunlar arasında (i) katılımcı tanımlayıcı (SEQN) ve anket tasarım değişkenlerini içeren demografik dosyalar, (ii) endometriozisin kendi raporunu içeren üreme sağlığı anket dosyaları ve (iii) bileşik indeks hesaplaması için gerekli laboratuvar dosyaları (C-reaktif protein, trigliseridler ve açlık plazma glikozu) yer aldı. Bu indeksler analiz edildiğinde karşılaştırıcı indeksler (örneğin, nötrofiller, lenfositler, trombositler) için gerekli muayene/antropometri dosyaları (örneğin vücut kitle indeksi) ve laboratuvar ölçümleri de elde edildi. Her 2 yıllık döngü içinde, bileşen dosyaları SEQN kullanılarak birleştirilir ve birleştirilen veri seti her SEQN başına bir kayıt olduğundan emin olmak için kontrol edilirdi. Daha sonra 1999–2006 analitik dosyasının birleşmesi için döngü düzeyinde veri setleri eklendi.
Analitik örneklem 20–54 yaş arası kadınlarla sınırlıydı. Endometriozis durumu eksikse, herhangi bir kompozit indeks bileşenleri (C-reaktif protein, trigliseridler veya açlık plazma glukozu) eksikse veya tam ayarlanmış model için gerekli temel kovaryatlar eksikse katılımcılar dışlandı. Karmaşık anket tasarım değişkenleri (katmanlar ve birincil örnekleme birimleri) ve oruç ölçümlerini içeren analizler için gereken oruç laboratuvarı alt örnek ağırlıkları korundu. Birden fazla NHANES döngüsü birleştirildiğinde, NHANES analitik rehberliğine göre 2 yıllık alt örneklem ağırlığı birleşik döngü sayısına bölünerek çok döngülü ağırlıklar oluşturuldu ve ortaya çıkan ağırlık, tabakalar ve PSU değişkenleri tüm analizlerde uygulandı. Katılımcı dahil etme ve dışlama adımları bir akış diyagramında belgelenmiştir (Şekil 1).
2. Endometriozisin tanımı
Endometriozis durumu, üreme sağlığı anketi maddesi kullanılarak tanımlandı: "Bir doktor veya başka bir sağlık profesyoneli size hiç endometriozis olduğunu söyledi mi?" "Evet" yanıtı veren katılımcılar endometriozis vakaları olarak, "Hayır" yanıtı verenler ise kontrol vakası olarak sınıflandırıldı. Bu tanım laparoskopik veya histolojik doğrulama yerine kendi raporuna dayandığı için, olası yanlış sınıflandırma çalışma sınırlaması olarak ele alındı.
3. Bileşik Indeksin (CTI) Tanımı
C-reaktif protein–trigliserid–glukoz bileşik indeksi, sistemik iltihabı ve metabolik bozukluğu birlikte yansıtacak şekilde operasyonel hale getirilmiştir. C-reaktif protein (mg/L), trigliseridler (mg/dL) ve açlık plazma glikozu (mg/dL) laboratuvar ölçümleri NHANES laboratuvar dosyalarından çıkarıldı. Trigliserid–glukoz indeksi, [trigliseridler × açlık plazma glukozu/2]'nin doğal logaritması olarak hesaplandı. CTI şu formülle hesaplandı: CTI = 0.412 × LN(CRP) + TyG. Daha yüksek CTI değerleri, düşük dereceli iltihap ve insülin direncinin birleşimde daha yüksek yükünügösterir 8.
Log dönüşümden önce herhangi bir C-reaktif protein değeri (örneğin, tespit sınırının altında veya üzerinde değerler) işlenilmesi gerekiyorsa, tüm döngülerde tutarlı şekilde tek bir önceden belirlenmiş kural uygulanmış ve tekrarlanabilirliği desteklemek için belgelenmiştir (örneğin, log dönüşümünden önce gözlemlenen en küçük pozitif ölçülebilir değerle pozitif olmayan değerlerin yerine koyulmak). Çeyrek kesme noktaları, tam analitik örneklemdeki ağırlıklı dağılımdan belirlenmiş ve kategorik analizler boyunca tutarlı şekilde uygulanmıştır; referans kategorisi olarak 1. Dörtre kullanılmıştır.
4. Ortak değişkenler
Kovariatlar, epidemiyolojik akıl yürütme ve önceki literatüre dayanarak karıştırıcı unsurları azaltmak için önceden belirlenmiştir. Demografik değişkenler arasında yaş, ırk/etnik köken, eğitim seviyesi ve medeni durum yer alıyordu. Yaşam tarzı değişkenleri arasında sigara geçmişi (≥ömür boyu 100 sigara vs. <100) ve alkol tüketimi (≥yılda 12 içkiye karşı <12) yer alıyordu. Komorbidite öyküsü, kendi bildirdiği hipertansiyon, diyabet, felç, koroner kalp hastalığı ve kanseri içeriyordu. Antropometrik ve laboratuvar değişkenleri arasında vücut kitle indeksi, hemoglobin, nötrofil sayısı, lenfosit sayısı ve trombosit sayısı yer alıyordu; Bu ölçümler ayrıca karşılaştırıcı inflamatuar indekslerin hesaplanmasını da destekledi (örneğin, nötrofil-lenfosit oranı, trombosit-lenfosit oranı, sistemik bağışıklık enflamasyon indeksi, sistemik iltihap yanıt indeksi). Üreme değişkenleri (örneğin, gravidity ve parite) seçilen döngülerde mevcut olduğunda dahil edilmiş ve NHANES dokümantasyonuna göre kodlanmıştır. Kategorik kovaryantlar, model girişinden önce gösterge değişkenlerine dönüştürüldü.
C-reaktif protein bileşik indeksin bir bileşeni olduğundan, aşırı ayarlama ve kollineerliği önlemek için çok değişkenli regresyon modellerinde bağımsız bir kovaryat olarak girilmemiştir. Bunun yerine, ayırtım analizlerinde C-reaktif protein ve trigliserid–glikoz indeksi karşılaştırıcı belirteçler olarak değerlendirildi.
5. İstatistiksel analiz
Tüm analizler, ulusal çapta temsil eden tahminler oluşturmak için NHANES'in karmaşık anket tasarımını hesaba kattı. Anket tasarımı, çok döngülü alt örneklem ağırlığı, tabakalar ve PSU değişkenlerinin analitik veri setine bağlanmasıyla belirlendi. Sürekli değişkenler standart sapmalarla ağırlıklı ortalamalar olarak özetlendi, kategorik değişkenler ise ağırlıklı sayılar ve yüzdelerle özetlendi. Başlangıç özellikleri, NHANES için uygun anket ağırlıklı prosedürler kullanılarak vakalar ve kontrol sistemleri arasında karşılaştırıldı ve başlangıç özellikleri Tablo 1'de özetlendi.
Bileşik indeks ile endometriozis arasındaki ilişkiler, anket ağırlıklı lojistik regresyon kullanılarak değerlendirildi. Uyumu göstermek için üç ardışık model uyguladı: ayarlanmamış bir model, yaş ve ırk/etnik kökene göre uyarlanmış bir model ve demografik faktörler, yaşam tarzı değişkenleri, komorbidite geçmişi, antropometrik/laboratuvar kovaryatları ve üreme geçmişi değişkenlerini içeren tam uyarlanmış bir model. Bileşik endeks hem sürekli (1 birim artış başlıkta) hem de kategorik olarak (dörtlükler, referans olarak 1. Dörtlük) analiz edildi ve regresyon tahminleri Tablo 2'de özetlendi. Dörtlüler arasında doğrusal eğilim, her çeyreğe ağırlıklı medyan değeri atanarak ve bu terimi sürekli modelleyerek test edildi.
Doğrusal olmayan doz-yanıt ilişkileri, önceden belirlenmiş düğüm yerleşimiyle anket ağırlıklı kısıtlı kübik splinlar kullanılarak değerlendirildi ve spline eğrileri Şekil 2'de çizildi. Eşik etkileri, segmentli ve tek eğimli spesifikasyonlar arasındaki model uyumu karşılaştırılarak anket ağırlıklı segmentli (parça parça) lojistik regresyon kullanılarak değerlendirildi ve dönüş noktasının her iki tarafındaki tahmini dönüş noktası ve eğim parametreleri Tablo 3'te bildirildi.
Alt grup analizleri, önceden belirlenmiş faktörlere (örneğin, yaş grubu, ırk/etnik köken, eğitim seviyesi, medeni durum ve seçilmiş yaşam tarzı faktörleri) göre etki değişikliğini incelemek için yapıldı. Etkileşim, sürekli bileşik endeks ile alt grup göstergeleri arasındaki çapraz ürün terimlerinin anket ağırlıklı çerçeve içinde dahil edilmesiyle test edilmiş ve alt grup ilişkileri Şekil 3'te özetlenmiştir.
Ayırtıcı performans, anket ağırlıklı lojistik modellerden türetilen model tahmin edilen olasılıklara dayanan alıcı işletim karakteristik analizleri kullanılarak değerlendirildi. Bileşik endeks, yaygın kullanılan inflamatuar indeksler ve bileşen belirteçleri için eğri altındaki alan tahminleri alındı ve AUC özetleri Ek Tablo 1'de sağlandı; genişletilmiş bir ROC karşılaştırması Ek Şekil 1'de sunulmuştur.
Eksik veriler, eksik endometriozis durumu, eksik kompozit indeks bileşenleri veya tam ayarlanmış model için gerekli temel kovaryantları eksik olan katılımcılar hariç tutulduktan sonra tam vaka analiziyle ele alındı. Bir dayanıklılık değerlendirmesi yapıldığında, önceden belirlenmiş bir ödüm modeli altında eksikliği olan kovaryatlar için birden fazla ödünç uygulandı ve tahminler tam vaka tahminleriyle karşılaştırıldı.
Analizler, R (sürüm 4.4.1) ve Materyal Tablosu'nda listelenen ek istatistiksel yazılımlar kullanılarak gerçekleştirildi. Anket çıkarımı, spline modellemesi, segmentli regresyon ve ROC tahmini için kullanılan anahtar paketler kaydedildi ve oturum bilgileri (işletim sistemi ve R oturum detayları) çoğaltmayı desteklemek için korundu.
6. İşlem bitiş noktası ve çıktılar
Uyumlu çok döngülü veri seti önceden belirlenmiş dahil etme/dışlama kriterleriyle oluşturulduktan (Şekil 1), bileşik indeks ve kovaryatlar belgelenmiş kurallara göre oluşturulduktan ve önceden belirlenmiş anket ağırlıklı regresyon, doğrusal olmayan/eşik değerlendirmesi, alt grup ve ayrımcılık analizleri aynı anket tasarım spesifikasyonu altında yürütüldükten sonra analitik iş akışı tamamlanmış sayıldı. Bu iş akışının birincil çıktıları temel özet (Tablo 1), ardışık ayarlama modelleri boyunca regresyon tahminleri (Tablo 2), eşik model parametreleri (Tablo 3), spline görselleştirme (Şekil 2), alt grup özet görselleştirmesi (Şekil 3) ve ayrımcılık özetleri (Ek Tablo 1 ve Ek Şekil 1) olarak düzenlenmiştir.
7. İçsel bağımsız doğrulama
Dahili bağımsız doğrulama, birleşik veri setinin NHANES döngülerine dayalı bir türetim kohortu ve örtüşmeyen bir doğrulama kohortu olarak bölünmesiyle gerçekleştirildi. 1999–2000 ve 2001–2002 döngülerinden katılımcılar türetim kohortuna atanmış, 2003–2004 ve 2005–2006 döngülerinden katılımcılar doğrulama kohortuna atanmıştır. Aynı dahil etme/dışlama kriteri, bileşik indeks hesaplaması, eşdeğer kodlama kuralları ve anket ağırlıklandırma stratejisi her kohort içinde bağımsız olarak uygulandı.
Türetim kohortu içinde, anket ağırlıklı lojistik regresyon modelleri tam ayarlanmış spesifikasyon kullanılarak yerleştirildi. Ana analizde kullanılan doğrusal olmayan ve eşik değerlendirme prosedürleri türetim kohortunda uygulandı ve ayrımcılık, model tahmin edilen olasılıklara dayalı ROC/AUC yöntemleriyle değerlendirildi. Aynı modelleme stratejisi, değişken tanımları, kodlama kuralları veya ağırlıklandırma spesifikasyonları değiştirilmeden doğrulama kohortunda tekrarlandı. Türetim-doğrulama tahminleri, ilişki tahminlerinin kohort-kohort karşılaştırması (Şekil 4) ve türetim-doğrulama (Şekil 5) olarak özetlenmiş, ilgili sayısal özetler ise Tablo 4'te sunulmuştur.