Wpływ przepisów antynikotynowych na śmiertelność we wczesnym okresie życia i niską masę urodzeniową w Anglii

Podstawowe pytanie badawcze Czy po wprowadzeniu 1 lipca 2007 r. w Anglii zakazu palenia w miejscach publicznych nastąpiła zmiana liczby dzieci urodzonych z niską masą urodzeniową lub umierających w okresie okołoporodowym lub niemowlęcym?

Projekt badania Retrospektywne badanie kohortowe (z wykorzystaniem prospektywnych, rutynowo zbieranych danych dotyczących opieki zdrowotnej)

Badana populacja Wszystkie pojedyncze urodzenia w Anglii w okresie od 1 stycznia 1995 r. do 31 grudnia 2011 r.

Interwencja Badana interwencja to zakaz palenia w zamkniętych miejscach publicznych i miejscach pracy wprowadzony w Anglii z dnia na dzień 1 lipca 2007 roku.

Kryteria włączenia i wyłączenia Uwzględnimy wszystkie zarejestrowane pojedyncze narodziny w Anglii, które miały miejsce między 1 stycznia 1995 r. a 31 grudnia 2011 r. Jest to maksymalny okres wokół wprowadzenia zakazu, dla którego wymagane dane dotyczące urodzeń są dostępne w źródle danych. Dane zostały pierwotnie pobrane za okres od 1 stycznia 1993 r. do 31 grudnia 2011 r. Jednak kod pocztowy nie został zarejestrowany w latach 1993-1994, co doprowadziło do brakujących wartości dla kwintylu wskaźnika wielokrotnej deprywacji (IMD), regionu i poziomu urbanizacji w tym okresie. Ponieważ w pierwotnych analizach zmienne te uznano za kluczowe potencjalne czynniki zakłócające, podjęto decyzję o ograniczeniu modelowania do okresu od 1 stycznia 1995 r. do 31 grudnia 2011 r.

Kodowanie Międzynarodowej Klasyfikacji Chorób (ICD) zmieniło się z wersji 9 na 10 w styczniu 2001 r., co doprowadziło do znacznego spadku zarejestrowanych przypadków SIDS. Dlatego analizy SIDS są ograniczone do okresu 2001-2011.

Niemowlęta z anomaliami chromosomowymi zostaną wykluczone.

Wynik

Główne wyniki to:

• Niska masa urodzeniowa (waga urodzeniowa < 2500 gramów)
• Poród martwego płodu (śmierć wewnątrzmaciczna od 24+0 tygodnia ciąży)
• Śmiertelność noworodków (zgon w ciągu pierwszych 28 dni życia)
• Zespół nagłej śmierci niemowląt (SIDS; zgon w ciągu pierwszego roku życia z wzmianką w akcie zgonu o kodzie ICD-10-U R95 lub R99 bez innej specyfikacji)

Aby ocenić, czy przepisy antynikotynowe miały selektywny wpływ na niektóre podgrupy wyników, zidentyfikowaliśmy ponadto szereg drugorzędnych wyników:

• Bardzo niska masa urodzeniowa (waga urodzeniowa < 1500 gramów)
• Wczesna śmiertelność noworodków (zgon w pierwszym tygodniu życia)
• Późna śmiertelność noworodków (zgon między 7 a 28 dniem życia)
• Śmiertelność poporodowa (zgon między 28 dniem a 1 rokiem życia)
• Śmiertelność niemowląt (zgon w pierwszym roku życia)

Źródła danych Dane pozyskiwane są za pośrednictwem Office for National Statistics (ONS). Uwzględniono wszystkie zarejestrowane martwe i żywe urodzenia, które miały miejsce w Anglii między 1 stycznia 1995 r. a 31 grudnia 2011 r. Są one powiązane z aktami zgonu dla wszystkich zgonów, które miały miejsce przed pierwszymi urodzinami.

Ekstrakcja i przetwarzanie danych Indywidualne dane dotyczące okresu okołoporodowego i śmiertelności są łączone przez ONS w jedną bazę danych zawierającą następujące zmienne towarzyszące na poziomie indywidualnym: miesiąc urodzenia, rok urodzenia, miesiąc śmierci, rok śmierci, wiek w chwili śmierci, płeć, masa urodzeniowa, wiek matki, stan cywilny matki, liczba porodów, IMD, region, stopień urbanizacji.

Następujące współzmienne zostały podzielone na kategorie ze względu na zarządzanie informacjami:

• Wiek w chwili śmierci: wczesny noworodek, późny noworodek, postneonatalny
• Waga urodzeniowa: <1000 gramów, 1000-1499 gramów, 1500-2499 gramów, 2500-3999 gramów, ≥4000 gramów
• Wiek matki: <20 lat, 20-24 lata, 25-29 lat, 30-34 lata, 35-39 lat, >40 lat
• Parzystość: 0, 1, 2, ≥3
• IMD: kwintyle
• Region: 10 regionów
• Poziom urbanizacji: miejski, wiejski

Wielkość próby Obliczenie wielkości próby dla analiz zorientowanych czasowo jest skomplikowane ze względu na złożoność modeli. W bieżącym badaniu wykorzystamy dane krajowe, które – zgodnie z naszą najlepszą wiedzą – będą największą oceną wpływu przepisów antynikotynowych na zdrowie okołoporodowe, zarówno pod względem wielkości populacji, jak i przedziału czasowego. Ponieważ używamy maksymalnego przedziału czasowego i dostępnej populacji, obliczenie wielkości próby można w pewnym sensie uznać za zbędne.

Znamy tylko jedno opublikowane badanie dotyczące ustawodawstwa dotyczącego zakazu palenia i śmiertelności we wczesnym okresie życia (odnośnik 1). Ze względu na problemy projektowe nie jest możliwe wykorzystanie danych z tego badania do porównania z bieżącym badaniem.

W wielu badaniach oceniano wcześniej wpływ przepisów antynikotynowych na niską masę urodzeniową. Proponowane przez nas podejście jest najlepiej porównywalne z podejściem przeprowadzonym wcześniej w Szkocji (referencja 2). Korzystając ze szkockich danych dotyczących 757 795 porodów, które miały miejsce w latach 1996-2009, wykazano natychmiastowy spadek liczby noworodków o niskiej masie urodzeniowej o -9,9% (95% przedział ufności -14,2; -5,2). Biorąc pod uwagę dłuższy okres badania (1993-2011) i znacznie większą populację (n>10 milionów), można oczekiwać, że nasze badanie będzie miało wystarczającą moc, aby wykryć podobny spadek liczby dzieci z niską masą urodzeniową w Anglii, jeśli występuje.

Analiza statystyczna

Opisane zostaną istotne cechy populacji. Przeprowadzona zostanie analiza regresji logistycznej w celu zbadania związku między wprowadzeniem przepisów antynikotynowych a nagłymi („krokowymi”) i/lub stopniowymi („skok”) zmianami (odpowiednio) prawdopodobieństwa wystąpienia każdego wyniku. Analizy zostaną dostosowane do masy urodzeniowej, płci, wieku matki, stanu cywilnego matki, parzystości (tylko analizy wtórne, patrz poniżej), kwintylu IMD, regionu i poziomu urbanizacji. W stosownych przypadkach zostaną uwzględnione wzorce sezonowe i nieliniowość podstawowego trendu czasowego. Ostateczny wybór modelu będzie oparty na kryteriach informacyjnych Akaike'a i bayesowskich (AIC i BIC). Mianownik analiz będzie różny w zależności od różnych wyników:

• martwe urodzenia: wszystkie urodzenia w zbiorze danych
• niska masa urodzeniowa, bardzo niska masa urodzeniowa, śmiertelność noworodków, wczesna śmiertelność noworodków, śmiertelność niemowląt: wszystkie żywo urodzone w zbiorze danych
• późna śmiertelność noworodków: wszystkie żywe urodzenia w zbiorze danych, które przeżyły wczesny okres noworodkowy
• śmiertelność poporodowa: wszystkie żywe urodzenia w zbiorze danych, które przeżyły okres noworodkowy

Podstawowe analizy zostaną przeprowadzone na przypadkach z pełnymi danymi dotyczącymi wszystkich zmiennych towarzyszących. Parzystość jest jedyną zmienną, w której brak danych wynosi >10% (około 40-50%), ponieważ występuje tylko u zamężnych kobiet. Ponieważ nie oczekuje się, że parytet będzie kluczowym czynnikiem zakłócającym, podstawowe analizy przeprowadzimy bez uwzględniania parytetu w modelach, aby zmaksymalizować wielkość populacji.

Analizy wrażliwości Aby ocenić możliwe zakłócenia parzystości, zostaną przeprowadzone analizy wrażliwości uwzględniające parytet w modelu. W drugim zestawie analiz wrażliwości zostanie przeprowadzona imputacja w celu zbadania odporności wyników na brakujące dane. Aby zminimalizować problemy związane z wieloma testami, analizy wrażliwości zostaną przeprowadzone tylko dla głównych wyników.

Wszystkie analizy zostaną przeprowadzone przy użyciu Stata 12.0.