Niniejsze retrospektywne badanie kohortowe w jednym ośrodku zostało zaprojektowane i opisane zgodnie z wytycznymi STROBE Statement14. W badaniu wykorzystano wcześniej zgromadzoną dokumentację medyczną pacjentów ambulatoryjnych, zdjęcia wewnątrzustne oraz dane obrazowania cyfrowego, bez konieczności przeprowadzania dodatkowych badań lub zabiegów. Protokół badania został zatwierdzony przez Komisję Etyczną Szpitala Stomatologicznego Sinopharm Dongfeng przy Uniwersytecie Medycznym w Hubei (numer zatwierdzenia LW-2026-01). Ponieważ wszystkie dane zostały zanonimizowane przed analizą, komisja etyczna zwolniła badaczy z obowiązku uzyskania świadomej zgody.
Źródło przypadków i populacja badawcza
Badanie zostało przeprowadzone w Katedrze Stomatologii Zachowawczej w Szpitalu Stomatologicznym Sinopharm Dongfeng, Uniwersytetu Medycznego w Hubei. Okres badania retrospektywnego objął kolejnych pacjentów, którzy otrzymali analizowaną terapię w okresie od 1 stycznia 2022 r. do 31 grudnia 2023 r. Dane z obserwacji zbierano do 28 lutego 2025 r., a baza danych została zamknięta 31 marca 2025 r. Dane pozyskano z elektronicznego systemu dokumentacji medycznej, systemu rejestracji zabiegów, systemu kodowania rozliczeń, cyfrowego archiwum obrazowego oraz systemu fotografii wewnątrzustnej. Datą indeksową było wskazanie daty zakończenia kwalifikującego się mikrowinwazyjnego leczenia. Jednostką analizy była analizowana zmiana w obrębie powierzchni okluzjalnej zęba, a nie pacjent.
Przydział do grup zabiegowych nie był randomizowany i nie został ustalony w celach badawczych. Zamiast tego, wybór metody leczenia został dokonany przez dentystę prowadzącego w ramach rutynowej praktyki klinicznej i nie był systematycznie rejestrowany jako niezależna zmienna decyzyjna. Na wybór leczenia mogły wpłynąć cechy zmiany, status wyrzynania zęba, możliwość kontroli wilgoci, preferencje operatora oraz inne rozważania kliniczne. W konsekwencji różne kwalifikujące się zęby u tego samego pacjenta mogły zostać objęte różnymi strategiami leczenia, a każdy ząb został zaklasyfikowany zgodnie z faktycznie zastosowanym zabiegiem. Nierówna liczebność grup odzwierciedlała względną częstotliwość stosowania obu strategii leczenia w rutynowej praktyce klinicznej, a nie zaplanowany przydział w stosunku 1:1. Analizy statystyczne uwzględniały grupowanie na poziomie pacjenta.
Populacja badana składała się z trwałych pierwszych lub drugich zębów trzonowych u pacjentów w wieku 6–18 lat. Zęby zakwalifikowane do badania spełniały wszystkie poniższe kryteria włączenia: podczas wizyty indeksowej, po oczyszczeniu zęba, izolacji i ciągłym osuszaniu powietrzem przez 5 s, docelowa powierzchnia okluzyjna została zdiagnozowana jako kod 1 lub 2 według Międzynarodowego Systemu Wykrywania i Oceny Próchnicy (ICDAS); wyjściowe cyfrowe zdjęcie rentgenowskie typu bitewing nie wykazało prześwietleń w obrębie zębiny na powierzchni okluzyjnej; oraz podczas tej samej wizyty wykonano albo infiltrację żywiczną połączoną z uszczelnianiem bruzd i dołków płynnym kompozytem żywicznym, albo samo uszczelnianie bruzd i dołków płynnym kompozytem żywicznym.
W przypadku zębów, w których nie wystąpiło zdarzenie końcowe w ciągu 12 miesięcy, wymagana była wizyta kontrolna w tym samym szpitalu między 10. a 14. miesiącem po leczeniu indeksowym, przy dostępności możliwych do interpretacji dokumentacji klinicznej, zdjęć wewnątrzustnych oraz cyfrowych radiogramów kątowych (bitewing). W przypadku zębów, które poddano restauracji operacyjnej z powodu postępu próchnicy na tej samej powierzchni okluzyjnej w ciągu 12 miesięcy, dokumentacja leczenia potwierdzająca restaurację służyła jako dowód wystąpienia zdarzenia i nie było wymaganych dodatkowych obrazów punktu końcowego.
Kryteria wykluczenia były następujące: wcześniejsza infiltracja żywicą, uszczelnianie bruzd i dołków, restauracja, korona protetyczna lub leczenie kanałowe obejmujące badaną powierzchnię okluzyjną; miejscowe uszkodzenie szkliwa, wyraźna próchnica, widoczny cień zębiny lub radiologiczna prześwietlność zębiny w badaniu wyjściowym; współistniejąca hipoplazja szkliwa, umiarkowane lub ciężkie wady strukturalne związane z fluorozą lub urazowe uszkodzenia korony; elementy ortodontyczne pokrywające badaną powierzchnię okluzyjną; brak możliwości potwierdzenia pierwotnej ekspozycji lub nieinterpretowalne dane z obserwacji. Jeśli w okresie badania zidentyfikowano wiele kwalifikujących się zapisów leczenia dla tego samego zęba, za leczenie indeksowe uznano wyłącznie najwcześniejszą kwalifikującą się procedurę.
Kryteria diagnostyki i klasyfikacji zmian
Wstępna diagnoza zmian opierała się na zestandaryzowanym badaniu klinicznym, fotografiach wewnątrzustnych oraz cyfrowych zdjęciach rentgenowskich typu bitewing wykonanych w dniu wizyty indeksowej. Cyfrowe zdjęcia rentgenowskie typu bitewing w badaniu wyjściowym i kontrolnym wykonano przy użyciu tego samego modelu wewnątrzustnego sensora CMOS (complementary metal-oxide-semiconductor) zgodnie ze zestandaryzowanym protokołem klinicznym (70 kVp, 7 mA oraz czas ekspozycji 0,10–0,14 s, w zależności od rozmiaru zęba i łuku). W celu standaryzacji pozycjonowania sensora i wyrównania wiązki promieni rentgenowskich zastosowano system pozycjonujący do zdjęć typu bitewing, składający się z uchwytu sensora, bloku przygryzowego, ramienia wyrównawczego i pierścienia celowniczego. Sensor ustawiono równolegle do koron zębów trzonowych objętych badaniem, a centralną wiązkę promieni skierowano przez punkty styku przyzębne. Podczas wizyt kontrolnych, o ile było to technicznie możliwe, stosowano ten sam typ sensora, system pozycjonujący, protokół ekspozycji oraz warunki wyświetlania obrazów w oprogramowaniu. Obrazy wykazujące znaczne nakładanie się styków przyzębnych, zniekształcenia geometryczne lub niepełną wizualizację powierzchni okluzyjnej badanych zębów uznano za nieinterpretowalne. Przez cały okres trwania badania przeprowadzano rutynowe codzienne kontrole systemu oraz comiesięczne procedury kontroli jakości.
Przed badaniem klinicznym powierzchnię zęba oczyszczono pastą profilaktyczną bez zawartości fluorków. Po izolacji ząb był suszony strumieniem powietrza przez 5 s, a ocena zmiany została przeprowadzona w standardowym oświetleniu gabinetowym przez lekarzy stomatologów przeszkolonych w Międzynarodowym Systemie Wykrywania i Oceny Próchnicy (ICDAS). W razie potrzeby użyto sondy z kulistym końcem wyłącznie w celu usunięcia zanieczyszczeń i potwierdzenia ciągłości powierzchni; nie wykonywano sondowania ostrym końcem. Kod ICDAS 1 zdefiniowano jako pierwszą widoczną zmianę w szkliwie, która stała się widoczna dopiero po osuszeniu powietrzem, w tym zmianę o barwie kredowo-białej lub przebarwioną. Kod ICDAS 2 zdefiniowano jako wyraźną zmianę koloru lub przezroczystości szkliwa, widoczną w stanie zwilżonym i stającą się bardziej oczywistą po osuszeniu powietrzem15. Oba kody ICDAS wymagały nienaruszonej powierzchni szkliwa, bez miejscowego ubytku szkliwa lub wyraźnej kawitacji.
Ponieważ zdjęcia rentgenowskie typu bitewing nie zapewniają takiej samej precyzji anatomicznej w przypadku zmian okkluzyjnych, jak w przypadku zmian proksymalnych, kategorie radiologiczne E0, E1, E2, D1, D2 i D3 zostały wykorzystane wyłącznie jako uprzednio określone, specyficzne dla badania porządkowe deskryptory obrazu radiologicznego, a nie jako bezpośrednie miary histologicznej głębokości zmiany. E0 oznaczało brak widocznej prześwietlenia w obrębie powierzchni okluzyjnej; E1 i E2 oznaczały prześwietlenie ograniczone odpowiednio do zewnętrznej i wewnętrznej połowy szkliwa; natomiast D1, D2 i D3 oznaczały prześwietlenie rozszerzające się odpowiednio do zewnętrznej, środkowej i wewnętrznej jednej trzeciej zębiny16. Kwalifikujące się zmiany zostały sklasyfikowane klinicznie jako kod ICDAS 1 lub 2 i nie wykazywały radiologicznych oznak zajęcia zębiny, co odpowiadało specyficznym dla badania kategoriom radiologicznym E0–E2.
W przypadku wystąpienia rozbieżności pomiędzy dokumentacją kliniczną, zdjęciami wewnątrzjuszczowymi a wynikami radiologicznymi, dwóch oceniających niezależnie analizowało dany przypadek. Jeśli po niezależnej analizie nadal występowała niezgodność, przypadek był przekazywany do rozstrzygnięcia przez trzeciego, starszego lekarza stomatologa, zgodnie z wcześniej zdefiniowaną procedurą rozstrzygania opisaną poniżej. Przypadki, które po rozstrzygnięciu nadal pozostawały nieinterpretowalne, zostały wykluczone z badania.
Protokoły ekspozycji i leczenia
Analizowaną ekspozycją była strategia leczenia mikroinwazyjnego zastosowana podczas wizyty indeksowej. Wszystkie zabiegi zostały wykonane zgodnie ze zstandaryzowanym protokołem oddziałowym przez trzech stomatologów z ponad 5-letnim doświadczeniem w stomatologii zachowawczej i regeneracyjnej: operatora 1 (X.Q.), operatora 2 (Z.Z.) oraz operatora 3 (J.Z.). Przed rozpoczęciem badania wszyscy operatorzy przeszli zstandaryzowane szkolenie z zakresu procedur klinicznych i protokołów leczenia stosowanych w niniejszym badaniu.
Przed rozpoczęciem leczenia powierzchnię zęba oczyszczono pastą profilaktyczną bez dodatku fluoru. Izolację uzyskano za pomocą koferdamu lub wałków bawełnianych w połączeniu z odsysaniem wysokoprzepływowym. Izolację koferdamem preferowano w przypadkach, gdy stopień wyrzutu i pozycja zęba pozwalały na stabilne osadzenie klamry oraz pełną kontrolę pola zabiegowego. Wałków bawełnianych w połączeniu z odsysaniem wysokoprzepływowym używano wtedy, gdy częściowy wyrzut, pozycja zęba lub niewystarczająca stabilność klamry uniemożliwiały niezawodne założenie koferdamu. Metoda izolacji była rejestrowana dla każdego leczonego zęba i uwzględniona jako zdefiniowana wcześniej współzmienna, ponieważ kontrola wilgotności mogła wpływać zarówno na wybór metody leczenia, jak i na właściwości materiałów. Rozkład metod izolacji pomiędzy grupami badawczymi przedstawiono w Tabeli 1. Wszystkie procedury światłoutwardzania przeprowadzono przy użyciu lampy LED o natężeniu wyjściowym co najmniej 1000 mW/cm2.
W grupie z leczeniem łączonym do zmiany w jamkach i szczelinach powierzchni okluzyjnej zastosowano żel do wstępnego trawienia zawierający 15% kwas solny przez 120 s, a następnie opłukano wodą przez 30 s i osuszono strumieniem powietrza. Następnie zmianę odwodniono za pomocą 99% etanolu przez 30 s. Kolejno nałożono niskolepki, utwardzany światłem infiltrant żywiczny na 3 min, a następnie utwardzano go światłem przez 40 s. Aplikację infiltranta powtórzono przez 1 min, po czym zastosowano drugi 40-sekundowy cykl utwardzania światłem. Po zakończeniu infiltracji żywicą wejście do szczeliny oraz przyległy szkliwo trawiono 37% kwasem fosforowym przez 20 s, dokładnie opłukano i osuszono powietrzem. Następnie nałożono uniwersalny system wiążący i utwardzano go światłem przez 20 s, po czym system szczelin uszczelniono płynnym kompozytem żywicznym. Aby ułatwić penetrację materiału do szczelin i usunąć wszelkie pęcherzyki powietrza przed utwardzaniem światłem przez 40 s, zastosowano sondę stomatologiczną. Na koniec oceniono adaptację brzeżną oraz kontakty okluzyjne i w razie potrzeby dokonano korekty.
W grupie z samym uszczelnianiem nie przeprowadzono infiltracji żywicą. Zamiast tego bruzdy oraz przyległe szkliwo poddano bezpośredniemu trawieniu 37% kwasem fosforowym przez 20 s, a następnie opłukano i wysuszono strumieniem powietrza. Zastosowano ten sam uniwersalny adhesyw, który utwardzono światłem przez 20 s, a do wykonania uszczelniania bruzd i dołków użyto tego samego płynnego kompozytu żywicznego. Końcowe utwardzanie światłem, ocena adaptacji brzeżnej oraz korekta okluzji zostały przeprowadzone zgodnie z procedurami zastosowanymi w grupie z leczeniem łączonym.
Po leczeniu pacjenci w obu grupach otrzymali identyczne instrukcje dotyczące higieny jamy ustnej i zostali pouczeni o konieczności szczotkowania zębów dwa razy dziennie pastą z zawartością 1450 ppm fluorku. Przestrzeganie tych instrukcji higienicznych nie było systematycznie monitorowane ze względu na retrospektywny charakter badania. Dodatkowe uszczelnianie wykonane w ciągu 30 dni po leczeniu indeksowym w celu skorygowania pęcherzyków brzeżnych lub miejscowego niedopełnienia uznano za część pierwotnej procedury leczniczej i nie sklasyfikowano go jako zdarzenia końcowego.
Dane wyjściowe i potencjalne czynniki zakłócające
Potencjalne czynniki zakłócające zostały określone przed ekstrakcją danych. Zmienne na poziomie pacjenta obejmowały wiek, płeć, całkowity wskaźnik dmft+DMFT oraz wskaźnik widocznej płytki nazębnej. Całkowity wskaźnik dmft+DMFT uzyskano na podstawie pełnego zapisu zębów próchniczych, usuniętych i plombowanych w dniu wizyty indeksowej, a wskaźnik widocznej płytki nazębnej uzyskano podczas badania pełnej jamy ustnej przeprowadzonego podczas tej samej wizyty.
Zmienne na poziomie zęba obejmowały łuk zębowy (szczękowy lub żuchwowy), rodzaj zęba trzonowego (pierwszy lub drugi trwały ząb trzonowy), status wyrzynania, wyjściową klasyfikację ICDAS oraz wyjściową kategorię obrazu radiologicznego specyficzną dla badania. Zmienne na poziomie operatora obejmowały tożsamość operatora oraz metodę izolacji.
Pełne wyrżnięcie zdefiniowano jako stan, w którym cała powierzchnia okluzyjna znajduje się powyżej brzegu dziąsła, co pozwala na całkowitą izolację. Wyrżnięcie częściowe zdefiniowano jako utrzymywanie się operculum dziąsłowego przykrywającego część powierzchni okluzyjnej. Wskaźnik widocznej płytki obliczono jako proporcję zbadanych powierzchni zębów, na których podczas wizyty indeksowej stwierdzono obecność widocznej płytki. Wszystkie zmienne towarzyszące określono w czasie leczenia indeksowego; do zdefiniowania zmiennych wyjściowych nie wykorzystano informacji z okresu obserwacji.
Obserwacja i ocena wyników
Rutynowe kontrole ambulatoryjne były przeprowadzane przez lekarzy prowadzących mniej więcej 6 i 12 miesięcy po leczeniu w ramach standardowej opieki klinicznej. Okres obserwacji w niniejszym badaniu wynosił 12 miesięcy po leczeniu indeksowym. W przypadku zębów objętych badaniem, w których nie stwierdzono progresji zmiany w okresie obserwacji, wizytę kontrolną najbliższą 12. miesiącowi w określonym przedziale czasowym od 10 do 14 miesięcy uznano za wizytę końcową. W przypadku zębów objętych badaniem, w których wykonano wypełnienie operacyjne z powodu progresji ubytków próchnicowych na tej samej powierzchni okluzyjnej w ciągu 12 miesięcy, pierwszą wizytę dokumentującą progresję zmiany uznano za punkt czasowy ustalenia wyniku. Zęby objęte badaniem, dla których nie odnotowano wyniku w ciągu 12 miesięcy i które nie miały kwalifikującej się wizyty kontrolnej w 10.–14. miesiącu, zaklasyfikowano jako utracone z obserwacji. Wizyta w 6. miesiącu służyła przede wszystkim rutynowemu przeglądowi klinicznemu. Wyniki tej wizyty wpłynęły na dalsze postępowanie kliniczne i zostały wykorzystane do ustalenia wyników, jeśli zidentyfikowano progresję zmiany; w przeciwnym razie wynik główny określano na podstawie wizyty kontrolnej najbliższej 12. miesiącowi.
Głównym punktem końcowym była progresja zmiany w ciągu 12 miesięcy. Progresję zmiany zdefiniowano jako wystąpienie któregokolwiek z następujących zdarzeń w okresie obserwacji: (1) powstanie wyraźnej kawerni lub miejscowego uszkodzenia szkliwa na brzegu uszczelniacza lub na odpowiadającej mu powierzchni okluzyjnej pod uszczelniaczem, przy potwierdzeniu w dokumentacji klinicznej, że zmiana ma charakter próchnicowy; (2) pojawienie się nowego prześwietlenia rozciągającego się poza szkliwo do zębiny na cyfrowym zdjęciu RTG typu bitewing powierzchni okluzyjnej indeksowej, sklasyfikowanego na potrzeby badania jako D1, D2 lub D3; lub (3) wykonanie wypełnienia operacyjnego na powierzchni okluzyjnej z powodu progresji zmiany próchnicowej.
Częściowa lub całkowita utrata uszczelniacza w przypadku braku któregokolwiek z powyższych znalezisk nie była uznawana za progresję zmiany. Podobnie, wypełnienia wykonane z powodu urazu, interferencji okluzyjnej lub złamania materiału niezwiązanego z próchnicą nie zostały zaklasyfikowane jako zdarzenia końcowe.
Ekstrakcja danych i kontrola jakości
Ekstrakcja danych została przeprowadzona z wykorzystaniem uprzednio określonego słownika danych oraz ustandaryzowanych formularzy ekstrakcji. Dwóch badaczy (X.Q. i Z.Z.) niezależnie ekstrahowało dane z elektronicznego systemu dokumentacji medycznej, systemu zapisów leczenia, systemu kodów rozliczeniowych, cyfrowego archiwum obrazów oraz systemu fotografii wewnątrzustnej. Baza danych badania została utworzona przy zastosowaniu podwójnego wprowadzania danych, a wszelkie rozbieżności rozstrzygano poprzez weryfikację oryginalnych zapisów źródłowych.
Przypisanie do grupy terapeutycznej musiało zostać potwierdzone zarówno na podstawie dokumentacji leczenia, jak i odpowiadających im kodów rozliczeniowych. W przypadku stwierdzenia rozbieżności między tymi źródłami analizowano dokumentację procedur oraz zdjęcia śródoperacyjne. Wykluczono przypadki, w których po analizie nie udało się potwierdzić przypisania do grupy terapeutycznej.
Przed oceną wszystkie pliki obrazowe zostały zanonimizowane, ponownie oznaczone przy użyciu numerów identyfikacyjnych badania oraz oddzielone od dokumentacji leczenia i informacji rozliczeniowych. Wstępna interpretacja obrazów została przeprowadzona bez dostępu do identyfikatorów pacjentów lub udokumentowanego przydziału do grup leczenia. W przypadkach, gdy do rozstrzygnięcia wyników wymagana była analiza dokumentacji klinicznej, dostęp do tych rekordów uzyskano dopiero po zakończeniu wstępnej interpretacji obrazów.
Dwaj równomiernie przeszkoleni stomatolodzy (X.Q. i Z.Z.) niezależnie ocenili wyjściową klasyfikację ICDAS, kategorię obrazu radiologicznego specyficzną dla badania oraz wyniki obserwacji. Obaj oceniający pełnili również rolę lekarzy prowadzących w kohorcie klinicznej. W związku z tym, mimo że informacje o leczeniu zostały udokumentowane, nie można było zagwarantować pełnego maskowania przydziału do grup terapeutycznych podczas wstępnej analizy obrazów.
Przed formalną interpretacją obrazów przeprowadzono kalibrację oceniających, wykorzystując 50 przypadków spoza badania. Formalna ocena rozpoczęła się dopiero po osiągnięciu wartości ważonych współczynników κ na poziomie 0,80 zarówno dla klasyfikacji ICDAS, jak i dla kategorii radiologicznego wyglądu specyficznej dla badania. W zbiorze kalibracyjnym końcowe ważone współczynniki κ wynosiły 0,84 dla klasyfikacji ICDAS oraz 0,88 dla kategorii radiologicznego wyglądu specyficznej dla badania.
W przypadku wystąpienia rozbieżności podczas formalnej oceny, dwóch oceniających przeprowadzało najpierw wspólną ponowną ocenę. Jeśli nie udało się osiągnąć konsensusu, przypadek był rozstrzygany przez trzeciego, starszego lekarza dentystę (J.L.), który nie brał udziału w początkowej ocenie parzystej. Rozstrzygnięcie w ramach trzeciej recenzji było wymagane dla 58 z 1 117 ocen wyjściowych (5,19%) oraz 46 z 989 ocen wyników w obserwacji (4,65%).
Z analizy wykluczono przypadki z brakującymi danymi dotyczącymi głównej ekspozycji lub głównego punktu końcowego. Brakujące dane ograniczały się do całkowitej oceny dmft+DMFT oraz indeksu widocznej płytki nazębnej i zostały uzupełnione za pomocą wielokrotnej imputacji metodą łańcuchowych równań z zastosowaniem 20 imputacji. Model imputacji obejmował grupę terapeutyczną, progresję zmiany po 12 miesiącach, wiek, płeć, całkowitą ocenę dmft+DMFT, indeks widocznej płytki nazębnej, łuk zębowy, typ zęba trzonowego, status wyrzynania, klasyfikację ICDAS wyjściową, wyjściową kategorię wyglądu radiologicznego specyficzną dla badania, operatora oraz metodę izolacji.
Analiza statystyczna
Nie przeprowadzono a priori obliczeń wielkości próby. Zamiast tego do analizy włączono wszystkie kolejne kwalifikujące się przypadki zidentyfikowane w okresie badania. Jako jednostkę analizy przyjęto zęba objętego procesem chorobowym, a wszystkie modele statystyczne uwzględniały grupowanie na poziomie pacjenta w celu skorygowania korelacji między wieloma objętymi procesem zębami u tej samej osoby.
W pierwszej kolejności oceniono rozkład zmiennych ciągłych. Zmienne o rozkładzie normalnym przedstawiono jako średnią ± odchylenie standardowe, natomiast zmienne o rozkładzie innym niż normalny przedstawiono jako medianę i rozstęp międzykwartylowy. Zmienne kategoryczne przedstawiono w formie liczności i procentów. Równowagę współzmiennych wyjściowych oceniono za pomocą standaryzowanych różnic średnich (SMDs), przyjmując, że bezwzględna wartość SMD <0,10 wskazuje na odpowiednią równowagę. Testowanie hipotez wyjściowych nie było stosowane do określania włączenia współzmiennych do modeli regresji.
Zdefiniowana a priori hipoteza zerowa zakładała, że 12-miesięczne ryzyko progresji zmiany nie będzie różnić się między dwiema grupami terapeutycznymi, co odpowiada względnemu ryzyku (RR) wynoszącemu 1,00. W analizie podstawowej w pierwszej kolejności oszacowano skumulowaną zapadalność oraz surowe RR progresji zmiany po 12 miesiącach w obu grupach, wraz z 95% przedziałami ufności (CIs). Analizy niepoprawione przeprowadzono za pomocą zmodyfikowanej regresji Poissona z zastosowaniem uogólnionych równań estymujących (GEEs), przyjmując grupę terapeutyczną jako jedyną zmienną niezależną. Analiza skorygowana wykorzystywała ten sam model i dodatkowo uwzględniała wiek, płeć, całkowity wskaźnik dmft+DMFT, wskaźnik widocznej płytki nazębnej, typ zęba trzonowego, łuk zębowy, status wyrzędu, klasyfikację ICDAS w punkcie wyjścia, kategorię radiologicznego wyglądu specyficznego dla badania w punkcie wyjścia, operatora oraz metodę izolacji. Zastosowano robustą estymację wariancji z wymienną strukturą korelacji roboczej i podano skorygowane wartości RR wraz z odpowiednimi 95% CIs. Zmienne towarzyszące wybrano a priori na podstawie istotności klinicznej i sekwencji czasowej, a nie wartości p z analiz jednowymiarowych.
W przypadku głównej ekspozycji oraz głównego punktu końcowego nie zastosowano imputacji. Brakujące dane dotyczące współzmiennych uzupełniono za pomocą wielokrotnej imputacji metodą łańcuchowych równań, generując 20 imputowanych zbiorów danych, a szacunki połączono zgodnie z regułami Rubina. Przeprowadzono cztery zaplanowane analizy wrażliwości w celu oceny odporności głównych wyników. Po pierwsze, analizę główną powtórzono, stosując analizę przypadków pełnych. Po drugie, analizę główną powtórzono, stosując odwrócone prawdopodobieństwo ważenia leczenia (IPTW) w oparciu o wskaźniki skłonności. Model wskaźnika skłonności obejmował te same współzmienne wyjściowe co główna analiza skorygowana, a równowagę współzmiennych po ważeniu ponownie oceniono za pomocą standaryzowanych różnic średnich. Po trzecie, analizę powtórzono po ograniczeniu populacji badania wyłącznie do zębów z ubytkami o kodzie ICDAS 2 w punkcie wyjścia. Po czwarte, analizę powtórzono po wykluczeniu częściowo wyrżniętych zębów trzonowych.
Wszystkie testy statystyczne były obustronne, a wartość p < 0,05 uznano za statystycznie istotną. Analizy statystyczne przeprowadzono przy użyciu oprogramowania R, wersja 4.4.2.