Populacja badana i przydział do grup
Do analizy końcowej włączono łącznie 159 kobiet w połogu. Spośród nich 95 miało prawidłowy czas trwania drugiego okresu porodu, natomiast 64 doświadczyło przedłużonego drugiego okresu porodu (Tabela 1). W porównaniu z grupą prawidłową, w grupie z przedłużonym czasem trwania odnotowano znacznie wyższą częstość występowania cierpienia płodu oraz niskich wyników w skali Apgara (w obu przypadkach p < 0,05). Nie zaobserwowano istotnych różnic w śmiertelności okołoporodowej między grupami (p > 0,05).
| Zmienna | Wszystkie
(N = 159) | Prawidłowy drugi okres porodu
(N = 95) | Przedłużony drugi okres porodu
(N = 64) | wartość p |
| Czas trwania drugiego etapu (min) | 109 [40.5; 172] | 88.4 [40.5; 125] | 139 [100; 172] | <0.0001 |
| Hipotonia płodu | | | | 0.0347 |
| Nieobecny | 124 (78.0%) | 80 (84.2%) | 44 (68.8%) | |
| Przedstawić | 35 (22.0%) | 15 (15.8%) | 20 (31.2%) | |
| Niska punktacja w skali Apgara | | | | 0.0001 |
| Nie | 79 (49.7%) | 60 (63.2%) | 19 (29.7%) | |
| Tak | 80 (50.3%) | 35 (36.8%) | 45 (70.3%) | |
| Śmiertelność okołoporodowa | | | | 0.5653 |
| Nieobecny | 156 (98.1%) | 94 (98.9%) | 62 (96.9%) | |
| Prezentuj | 3 (1.89%) | 1 (1.05%) | 2 (3.12%) | |
Tabela 1: Wpływ przedłużonego drugiego okresu porodu na kliniczne wyniki noworodkowe. Zmienne ciągłe przedstawiono jako medianę [rozstęp międzykwartylowy (IQR)], a zmienne kategoryczne jako liczbę (%). Grupa z normalnym drugim okresem porodu obejmowała uczestniczki, u których czas trwania drugiego okresu wynosił ≤120 min, natomiast grupa z przedłużonym drugim okresem porodu obejmowała uczestniczki, u których czas trwania drugiego okresu wynosił >120 min. Niski wynik w skali Apgar (Appearance, Pulse, Grimace, Activity, and Respiration) zdefiniowano jako wynik Apgar w 1. minucie <7. Wartości P reprezentują porównania między grupą z normalnym a grupą z przedłużonym drugim okresem porodu.
Porównanie charakterystyki demograficznej i klinicznej obu grup przedstawiono w Tabeli 2. Rodzące kobiety z przedłużonym drugim okresem porodu charakteryzowały się wyższą medianą wieku matki oraz większym przyrostem BMI w czasie ciąży niż kobiety z prawidłowym drugim okresem porodu. Częstość występowania cukrzycy ciążowej, bezwładności macicy oraz zaburzeń stanu psychicznego była również istotnie wyższa w grupie z przedłużonym okresem porodu (wszystkie p < 0,05). Nie zaobserwowano istotnych różnic między grupami w zakresie nadciśnienia ciążowego ani wieku ciążowego w momencie porodu (oba p > 0,05).
| Zmienna | Wszystkie
(N = 159) | Prawidłowy drugi okres porodu
(N = 95) | Przedłużony drugi okres porodu
(N = 64) | Wartość p |
| Wiek (lat) | 28.8 [19.4; 4.9] | 26.3 [19.4; 3.2] | 32.4 [2.3; 4.9] | <0.01 |
| Wzrost BMI w czasie ciąży (kg/m²) | 1.4 ± 3.62 | 10.1 ± 3.3 | 13.3 ± 3.17 | <0.001 |
| Cukrzyca ciążowa | | | | 0.002 |
| Brak | 141 (8.7%) | 92 (96.8%) | 49 (76.6%) | |
| Występuje | 18 (1.3%) | 3 (3.16%) | 15 (23.4%) | |
| Nadciśnienie tętnicze w ciąży | | | | 0.8723 |
| Brak | 124 (78.0%) | 75 (78.9%) | 49 (76.6%) | |
| Występuje | 35 (22.0%) | 20 (21.1%) | 15 (23.4%) | |
| Masa urodzeniowa noworodka (kg) | 2.61 [1.68; 4.60] | 2.46 [1.68; 2.85] | 3.39 [1.71; 4.60] | <0.001 |
| Stan psychiczny | | | | 0.017 |
| Zaburzony | 139 (87.4%) | 90 (94.7%) | 49 (76.6%) | |
| Zdrowy | 20 (12.6%) | 5 (5.26%) | 15 (23.4%) | |
| Inercja macicy | | | | <0.01 |
| Brak | 123 (7.4%) | 85 (89.5%) | 38 (59.4%) | |
| Występuje | 36 (2.6%) | 10 (10.5%) | 26 (40.6%) | |
| Wiek ciążowy (tygodnie) | 38.1 ± 2.02 | 38.1 ± 1.9 | 37.9 ± 2.06 | 0.5156 |
Tabela 2: Porównanie charakterystyki demograficznej i klinicznej między grupami. Zmienne ciągłe przedstawiono jako średnią ± odchylenie standardowe (SD) lub medianę [rozstęp międzykwartylny (IQR)], w zależności od potrzeb. Zmienne kategoryczne przedstawiono jako liczbę (%). Przyrost wskaźnika masy ciała (BMI) zdefiniowano jako różnicę między BMI przed ciążą a BMI przed porodem. Grupa z prawidłowym drugim okresem porodu obejmowała uczestniczki z czasem trwania drugiego okresu ≤120 min, natomiast grupa z przedłużonym drugim okresem porodu obejmowała uczestniczki z czasem trwania drugiego okresu >120 min. Wartości P reprezentują porównania między dwiema grupami.
Parametry kliniczne i ultrasonograficzne w aktywnym pierwszym okresie
Porównano parametry ultrasonograficzne przezbrzuszne i przezkroczne w pierwszym okresie między grupą z prawidłowym a grupą z przedłużonym drugim okresem (Tabela 3). Odsetek pozycji płodu OA była istotnie niższa w grupie z przedłużonym okresem niż w grupie z prawidłowym okresem. Ponadto AOP, HSD i HPD były istotnie większe u rodzących, u których wystąpił przedłużony drugi okres (wszystkie p < 0,05).
| Zmienna | Wszystkie
(N = 159) | Prawidłowy drugi okres porodu
(N = 95) | Przedłużony drugi okres porodu
(N = 64) | wartość p |
| Pozycja embrionalna | | | | <0.0001 |
| zwyrodnienie stawów | 74 (46.5%) | 60 (63.2%) | 14 (21.9%) | |
| Nie-OA | 85 (53.5%) | 35 (36.8%) | 50 (78.1%) | |
| AOP (°) | 124 [64.3; 158] | 111 [64.3; 142] | 138 [123; 158] | <0.0001 |
| HSD (cm) | 1.97 ± 0.39 | 1.81 ± 0.33 | 2.21 ± 0.34 | <0.0001 |
| HPD (cm) | 4.61 ± 0.71 | 4.23 ± 0.56 | 5.17 ± 0.51 | <0.0001 |
Tabela 3: Porównanie parametrów ultrasonografii transabdominalnej i transkroczywej podczas pierwszego okresu porodu. Zmienne ciągłe przedstawiono jako średnią ± odchylenie standardowe (SD) lub medianę [rozstęp międzykwartylowy (IQR)], w zależności od potrzeb. Zmienne kategoryczne przedstawiono jako liczbę (%). OA – potylica z przodu; non-OA – potylica nie z przodu; AOP – kąt progresji; HSD – odległość głowa–łóneczna; HPD – odległość głowa–krocze. AOP mierzono w stopniach (°), natomiast HSD i HPD w centymetrach (cm). Grupa z prawidłowym drugim okresem porodu obejmowała uczestniczki, u których czas trwania drugiego okresu wynosił ≤120 min, natomiast grupa z przedłużonym drugim okresem porodu obejmowała uczestniczki, u których czas trwania drugiego okresu wynosił >120 min. Wartości P reprezentują porównania między obiema grupami.
Analiza jednowymiarowa potencjalnych predyktorów
Przeprowadzono analizy jednowymiarowe w celu zidentyfikowania zmiennych związanych z przedłużonym drugim okresem porodu (Tabela 4). Dziesięć zmiennych różniło się istotnie pomiędzy grupą normalną a grupą z przedłużonym okresem porodu (p < 0,10): wiek matki, przyrost BMI w czasie ciąży, cukrzyca ciążowa, masa urodzeniowa noworodka, stan psychiczny, bezwładność macicy, pozycja płodu, AOP, HSD oraz HPD. Nie zaobserwowano istotnych różnic międzygrupowych w przypadku nadciśnienia ciążowego ani wieku ciążowego (w obu przypadkach p > 0,10). Zmienne z wartością p < 0,10 zostały następnie wprowadzone do modelu regresji LASSO w celu selekcji zmiennych.
| Zmienna | Wszyscy (N = 159) | Prawidłowy drugi okres porodu (N = 95) | Przedłużony drugi okres porodu (N = 64) | Statystyka | Wartość P |
| Charakterystyka demograficzna i kliniczna | | | | | |
| Wiek (lata) | 28.8 [19.4; 4.9] | 26.3 [19.4; 3.2] | 32.4 [2.3; 4.9] | Z = −5.891 | <0.01 |
| Wzrost BMI w czasie ciąży (kg/m²) | 1.4 ± 3.62 | 10.1 ± 3.3 | 13.3 ± 3.17 | t = −5.912 | <0.01 |
| Cukrzyca ciążowa, n (%) | 18 (1.3) | 3 (3.2) | 15 (23.4) | χ² = 15.207 | <0.001 |
| Nadciśnienie ciążowe, n (%) | 35 (2.0) | 20 (21.1) | 15 (23.4) | χ² = 0.135 | 0.872 |
| Masa urodzeniowa noworodka (kg) | 2.61 [1.68; 4.60] | 2.46 [1.68; 2.85] | 3.39 [1.71; 4.60] | Z = −6.234 | <0.01 |
| Zaburzenia stanu psychicznego, n (%) | 139 (87.4) | 90 (94.7) | 49 (76.6) | χ² = 10.912 | 0.01 |
| Inercja macicy, n (%) | 36 (22.6) | 10 (10.5) | 26 (40.6) | χ² = 20.156 | <0.01 |
| Wiek ciążowy (tygodnie) | 38.1 ± 2.02 | 38.1 ± 1.9 | 37.9 ± 2.06 | t = 0.654 | 0.516 |
| Parametry ultrasonograficzne podczas pierwszego okresu porodu | | | | | |
| Położenie płodu (inne niż potyliczne przednie), n (%) | 85 (53.5) | 35 (36.8) | 50 (78.1) | χ² = 26.543 | <0.01 |
| AOP (°) | 124 [64.3; 158] | 1 [64.3; 142] | 138 [123; 158] | Z = −6.781 | <0.01 |
| HSD (cm) | 1.97 ± 0.39 | 1.81 ± 0.3 | 2.21 ± 0.34 | t = −7.45 | <0.01 |
| HPD (cm) | 4.61 ± 0.71 | 4.23 ± 0.56 | 5.17 ± 0.51 | t = −9.32 | <0.01 |
Tabela 4: Analiza jednowymiarowa czynników związanych z przedłużonym drugim okresem porodu. Zmienne ciągłe przedstawiono jako średnia ± odchylenie standardowe (SD) lub mediana [rozstęp międzykwartylny (IQR)], w zależności od rozkładu. Zmienne kategoryczne przedstawiono jako liczba (%). BMI, wskaźnik masy ciała; OA, potylica z przodu; non-OA, potylica nie z przodu; AOP, kąt progresji; HSD, odległość głowa–łóneczna; HPD, odległość głowa–krocze. Zmienne ciągłe o rozkładzie normalnym porównywano za pomocą testu t dla prób niezależnych, zmienne ciągłe o rozkładzie nienormalnym porównywano za pomocą testu U Manna–Whitneya (raportowanego jako statystyka Z), a zmienne kategoryczne porównywano za pomocą testu chi-kwadrat (χ2). Zmienne z P <0.10 uznano za kwalifikujące się do włączenia do analizy regresji metodą LASSO (least absolute shrinkage and selection operator).
Wybór zmiennych za pomocą regresji LASSO (Least Absolute Shrinkage and Selection Operator)
W celu skonstruowania oszczędnego modelu predykcyjnego przy jednoczesnym zmniejszeniu ryzyka przeuczenia oraz wpływu współliniowości, 10 potencjalnych predyktorów zidentyfikowanych w analizie jednowymiarowej (p < 0,10; Tabela 4) wprowadzono do modelu regresji LASSO w celu wyboru zmiennych. Optymalny parametr regularyzacji (λ) określono za pomocą 10-krotnej walidacji krzyżowej. Wykres ścieżki współczynników (Rysunek 2A) przedstawia stopniowe kurczenie współczynników regresji w stronę zera wraz ze wzrostem λ. Odpowiadająca mu krzywa walidacji krzyżowej (Rysunek 2B) pokazuje zależność między błędem walidacji krzyżowej a różnymi wartościami λ. Lewa pionowa linia przerywana wskazuje wartość λ, która zminimalizowała błąd walidacji krzyżowej, natomiast prawa pionowa linia przerywana wskazuje największą wartość λ mieszczącą się w granicach jednego błędu standardowego od minimum (kryterium jednego błędu standardowego). To drugie wybrano w celu uzyskania bardziej oszczędnego modelu. Stosując kryterium jednego błędu standardowego, zachowano cztery predyktory o niezerowych współczynnikach: masę urodzeniową noworodka (współczynnik = 0,50), przyrost BMI w czasie ciąży (współczynnik = 0,12), wiek matki (współczynnik = 0,05) oraz AOP (współczynnik = 0,03). Współczynniki te reprezentują ustandaryzowane współczynniki z penalizacją z końcowego modelu LASSO dla wybranej wartości λ. Pozostałe zmienne kandydackie — HSD, HPD, cukrzyca ciążowa, stan psychiczny, bezwładność macicy oraz pozycja płodu — miały współczynniki zredukowane do zera i zostały zatem wykluczone z ostatecznego zestawu predyktorów. Zachowane zmienne wprowadzono następnie do wieloczynnikowego modelu regresji logistycznej w celu oszacowania skorygowanych wartości OR wraz z odpowiadającymi im 95% CI.

Rycina 2Analiza regresji metodą LASSO (least absolute shrinkage and selection operator) w celu wyboru predyktorów. (A) Ścieżki współczynników wygenerowane podczas regresji LASSO (least absolute shrinkage and selection operator). Oś x reprezentuje logarytm parametru regularyzacji (log λ), a oś y reprezentuje wystandaryzowane współczynniki regresji. Każda krzywa ilustruje stopniowe kurczenie się współczynnika poszczególnych predyktorów w stronę zera wraz ze wzrostem kary regularyzacyjnej. Liczby widoczne u góry wskazują liczbę zmiennych zachowanych dla każdej wartości λ. Ponumerowane ścieżki współczynników odpowiadają następującym potencjalnym predyktorom: (1) wiek matki, (2) przyrost wskaźnika masy ciała (BMI) w czasie ciąży, (3) cukrzyca ciążowa, (4) masa urodzeniowa noworodka, (5) stan psychiczny, (6) bezwładność macicy, (7) pozycja płodu, (8) kąt progresji (AOP), (9) odległość głowa–łóneczko (HSD) oraz (10) odległość głowa–krocze (HPD). (B) Wykres dziesięciokrotnej walidacji krzyżowej użyty do wyznaczenia optymalnego parametru regularyzacji. Oś x przedstawia log(λ), a oś y przedstawia dewiancję dwumianową. Punkty wskazują średni błąd walidacji krzyżowej, a słupki błędów reprezentują ±1 błąd standardowy. Lewa pionowa linia przerywana wskazuje wartość λ, która minimalizuje błąd walidacji krzyżowej, natomiast prawa pionowa linia przerywana wskazuje największą wartość λ w granicach jednego błędu standardowego od minimum (kryterium jednego błędu standardowego). Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.
Wieloczynnikowa regresja logistyczna i niezależne predyktory
Cztery predyktory zachowane w modelu regresji LASSO — wiek matki, przyrost BMI w czasie ciąży, masę urodzeniową noworodka oraz AOP — wprowadzono do modelu wieloczynnikowej regresji logistycznej (Tabela 5). Wszystkie cztery zmienne pozostały niezależnie powiązane z przedłużonym drugim okresem porodu po wzajemnej korekcie (wszystkie p < 0,05). W dopasowanym modelu wieloczynnikowej regresji logistycznej masa urodzeniowa noworodka miała najwyższy skorygowany współczynnik OR (aOR = 8,45, 95% CI: 3,21–2,24). Przyrost BMI w czasie ciąży był również niezależnie powiązany z przedłużonym drugim okresem porodu (aOR = 1,38, 95% CI: 1,15–1,6). Wiek matki (aOR = 1,18, 95% CI: 1,05–1,32) oraz AOP (aOR = 1,07, 95% CI: 1,02–1,12) były w podobny sposób niezależnie powiązane z przedłużonym drugim okresem porodu. Model wieloczynnikowej regresji logistycznej był istotny statystycznie zgodnie z testem ilorazu wiarygodności (p < 0,01). Ponieważ masę urodzeniową noworodka zmierzono po porodzie, obecny model ma charakter eksploracyjny i wymagałby ponownego dopasowania z wykorzystaniem szacowanej masy płodu przed zastosowaniem klinicznym.
| Zmienna | β | SE | Wald χ²² | wartość p | dOR (dostosowany iloraz szans) | 95% CI (przedział ufności) |
| Punkt przecięcia | −15.92 | 3.01 | 27.98 | <0.001 | — | — |
| Wiek matki (lata) | 0.16 | 0.06 | 7.11 | 0.008 | 1.18 | 1.05–1.32 |
| Przyrost BMI w czasie ciąży (kg/m²²) | 0.32 | 0.09 | 12.65 | <0.001 | 1.38 | 1.15–1.66 |
| Masa urodzeniowa noworodka (kg) | 4.26 | 0.9 | 22.41 | <0.001 | 8.45 | 3.21–22.24 |
| AOP (°) | 0.07 | 0.02 | 10.76 | 0.001 | 1.07 | 1.02–1.12 |
Tabela 5: Wieloczynnikowa analiza regresji logistycznej identyfikująca niezależne predyktory przedłużonego drugiego okresu porodu. Tabela przedstawia wyniki wieloczynnikowego modelu regresji logistycznej zbudowanego z wykorzystaniem zmiennych wybranych za pomocą regresji LASSO (least absolute shrinkage and selection operator). β, współczynnik regresji; SE, błąd standardowy; Wald χ2, statystyka chi-kwadrat Walda; aOR, skorygowany iloraz szans; CI, przedział ufności; BMI, wskaźnik masy ciała; AOP, kąt progresji. Punkt przecięcia (intercept) reprezentuje stałą modelu. Ogólną istotność modelu oceniono za pomocą testu ilorazu wiarygodności (χ2 = 85.34, P <0.01).
Konstrukcja modelu predykcyjnego i nomogram
Na podstawie niezależnych czynników predykcyjnych oraz ich współczynników regresji zidentyfikowanych w wieloczynnikowej analizie regresji logistycznej (Tabela 5), skonstruowano specyficzny dla pacjentki model predykcyjny w celu oszacowania prawdopodobieństwa przedłużonego drugiego okresu porodu u pierworódek. Model wyrażono w następujący sposób:
Logit(P) = —15,92 + 0,16 ( wiek matki ) +. 0,32 ( przyrost BMI w czasie ciąży ) + 4,26 ( masa urodzeniowa noworodka ) + 0,07( AOP )
Tutaj P oznacza przewidywane prawdopodobieństwo przedłużonego drugiego okresu porodu.
Aby ułatwić wizualizację modelu, równanie regresji przedstawiono w formie nomogramu (Rycina 3). Każdy predyktor odpowiada indywidualnej skali punktowej. Na przykład, w przypadku rodzącej kobiety w wieku 28 lat, ze wzrostem BMI w czasie ciąży o 12 kg/m2, masą urodzeniową noworodka (rzeczywistą masą zmierzoną po porodzie) wynoszącą 3,6 kg oraz kątem AOP wynoszącym 15°, punktację można określić poprzez odnalezienie każdej z tych wartości na odpowiadającej im osi, rzutowanie pionowe na oś punktową (Points), zsumowanie poszczególnych wartości punktowych w celu uzyskania całkowitej liczby punktów (Total Points), a następnie rzutowanie wyniku całkowitego w dół na oś ryzyka przedłużonego drugiego okresu porodu (Risk of Prolonged Second Stage) w celu oszacowania przewidywanego prawdopodobieństwa. Nomogram ten został przedstawiony wyłącznie w celach ilustracyjnych, ponieważ uwzględnia on rzeczywistą masę urodzeniową noworodka. Zastosowanie kliniczne wymagałoby ponownego dopasowania modelu z wykorzystaniem szacowanej masy płodu.

Rycina 3Nomogram do przewidywania prawdopodobieństwa przedłużonego drugiego okresu porodu. Nomogram został opracowany na podstawie wieloczynnikowego modelu regresji logistycznej z wykorzystaniem czterech predyktorów: wieku matki, przyrostu BMI w trakcie ciąży, masy urodzeniowej noworodka (faktycznej masy zmierzonej po porodzie) oraz AOP. Aby oszacować przewidywane prawdopodobieństwo przedłużonego drugiego okresu porodu, należy odnaleźć wartość każdego predyktora na odpowiadającej mu osi i rzutować ją pionowo na oś punktów (Points), aby ustalić przypisaną liczbę punktów. Następnie należy zsumować poszczególne wyniki, aby uzyskać łączną liczbę punktów (Total Points), a następnie rzutować tę sumę na oś ryzyka (Risk), aby oszacować przewidywane prawdopodobieństwo przedłużonego drugiego okresu porodu. Nomogram został przedstawiony wyłącznie w celach ilustracyjnych, ponieważ uwzględnia faktyczną masę urodzeniową noworodka, która nie jest znana przed porodem. Zastosowanie kliniczne wymagałoby ponownego dopasowania modelu z wykorzystaniem EFW. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Walidacja modelu: dyskryminacja, kalibracja i użyteczność kliniczna
Wydajność predykcyjna modelu została oceniona wewnętrznie pod kątem dyskryminacji, kalibracji i użyteczności klinicznej. Dyskryminację modelu oceniono za pomocą krzywej ROC (Rysunek 4). Model osiągnął wartość AUC wynoszącą 0,927 (95% CI: 0,879–0,961), co wskazuje na doskonałą dyskryminację. Przy zastosowaniu maksymalnego indeksu Youdena optymalny próg prawdopodobieństwa wyniósł 0,48. Przy tym progu model wykazał czułość na poziomie 85,9% oraz swoistość na poziomie 84,2% (Tabela 6).

Rysunek 4Krzywa charakterystyki operacyjnej odbiornika (ROC) oceniająca wydajność modelu predykcyjnego. Krzywa charakterystyki operacyjnej odbiornika (ROC) obrazuje zdolność modelu predykcyjnego do rozróżnienia między prawidłowym a przedłużonym drugim okresem porodu. Niebieska krzywa reprezentuje model predykcyjny z polem pod krzywą (AUC) wynoszącym 0,927 (95% przedział ufności [CI]: 0,879–0,961). Optymalna wartość odcięcia prawdopodobieństwa wyniosła 0,48, co odpowiada czułości na poziomie 85,9% i swoistości na poziomie 84,2%. Szara przerywana linia diagonalna reprezentuje klasyfikację losową (AUC = 0,50). Oś x przedstawia 1 − swoistość (wskaźnik wyników fałszywie dodatnich [FPR]), a oś y przedstawia czułość (wskaźnik wyników prawdziwie dodatnich [TPR]). Kliknij tutaj, aby wyświetlić powiększoną wersję tej figury.
| Parametr | Wartość | 95% CI |
| AUC (pole pod krzywą) | 0.927 | 0.879–0.961 |
| Optymalny punkt odcięcia | 0.480 | — |
| Czułość | 85.90% | 75.0%–93.4% |
| Swoistość | 84.20% | 75.3%–90.9% |
| PPV | 78.60% | 67.5%–87.3% |
| NPV | 89.80% | 81.9%–95.0% |
| Dokładność | 84.90% | 78.3%–90.1% |
Tabela 6: Efektywność modelu predykcyjnego przy optymalnym progu prawdopodobieństwa. Efektywność predykcyjna modelu wieloczynnikowej regresji logistycznej została oceniona za pomocą pola pod krzywą charakterystyki operacyjnej odbiornika (ROC) (AUC). Optymalny próg prawdopodobieństwa wynoszący 0,48 został wyznaczony poprzez maksymalizację indeksu Youdena. Czułość reprezentuje proporcję prawidłowo zidentyfikowanych przypadków prawdziwie dodatnich, swoistość reprezentuje proporcję prawidłowo zidentyfikowanych przypadków prawdziwie ujemnych, dodatnia wartość predykcyjna (PPV) reprezentuje prawdopodobieństwo, że uczestnicy sklasyfikowani jako dodatni rzeczywiście mieli przedłużony drugi okres porodu, a ujemna wartość predykcyjna (NPV) reprezentuje prawdopodobieństwo, że uczestnicy sklasyfikowani jako ujemni rzeczywiście nie mieli przedłużonego drugiego okresu porodu. CI, przedział ufności.
Kalibrację modelu oceniono za pomocą bootstrapowej krzywej kalibracji wygenerowanej z 10 reprobowań oraz testu dopasowania Hosmera–Lemeshow (Rysunek 5A). Krzywa kalibracji ściśle pokrywała się z idealną linią referencyjną, co wskazuje na dobrą zgodność między prawdopodobieństwem przewidywanym a obserwowanym. Test dopasowania Hosmera–Lemeshow nie był istotny statystycznie (χ2 = 6.15, p = 0.630), co było spójne z prawidłową kalibracją w tym wewnętrznie zwalidowanym modelu. Użyteczność kliniczną oceniono za pomocą DCA (Rysunek 5B). W zakresie prawdopodobieństwa progowego od około 10% do 70%, model predykcyjny zapewniał większą korzyść netto niż strategia „lecz wszystkich” lub „nie lecz nikogo”. Przy prawdopodobieństwie progowym poniżej około 10%, większą korzyść netto zapewniała strategia „nie lecz nikogo”, natomiast przy prawdopodobieństwie progowym powyżej około 70%, strategia „lecz wszystkich” wykazywała podobne lub lepsze wyniki. Wyniki te wskazują, że model może przynieść korzyść kliniczną w średnim zakresie prawdopodobieństw progowych. Ogólnie rzecz biorąc, model predykcyjny wykazał doskonałą dyskryminację, prawidłową kalibrację i korzystną użyteczność kliniczną podczas walidacji wewnętrznej. Jednak ze względu na to, że walidacja była ograniczona do zbioru danych badania, a model uwzględniał rzeczywistą masę urodzeniową noworodka, przed wdrożeniem klinicznym wymagana jest walidacja zewnętrzna oraz ponowne dopasowanie modelu z wykorzystaniem szacowanej masy płodu.

Rycina 5Analizy kalibracji i krzywej decyzyjnej modelu predykcyjnego. (A) Krzywa kalibracyjna modelu predykcyjnego wygenerowana przy użyciu 1 0 prób bootstrapowych. Krzywa jawna (apparent curve) przedstawia wydajność modelu w zbiorze danych badania, krzywa skorygowana o błąd (bias-corrected curve) przedstawia kalibrację dostosowaną metodą bootstrapową, a ukośna linia przerywana oznacza idealną zgodność między prawdopodobieństwem przewidywanym a obserwowanym. Test dopasowania Hosmera–Lemeshowa dał wynik χ2 = 6,15 i p = 0,630. (B) Analiza krzywej decyzyjnej (DCA) oceniająca użyteczność kliniczną modelu predykcyjnego w zakresie różnych prawdopodobieństw progowych. Zielona krzywa reprezentuje model predykcyjny, pomarańczowa linia reprezentuje strategię leczenia wszystkich pacjentów, a niebieska linia reprezentuje strategię nieleczenia żadnego z pacjentów. Zacieniowany obszar wskazuje zakres prawdopodobieństwa progowego od 10% do 70% oceniany w niniejszym badaniu, a pionowe linie przerywane wyznaczają dolny (10%) i górny (70%) próg. Oś x reprezentuje prawdopodobieństwo progowe, a oś y reprezentuje korzyść netto. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Dostępność danych:
Zanonimizowany zbiór danych na poziomie uczestników, stanowiący podstawę wniosków z niniejszego badania, został przedstawiony w Tabeli uzupełniającej 1. Zbiór ten zawiera surowe dane wykorzystane do wszystkich analiz statystycznych, opracowania modeli, walidacji wewnętrznej oraz przygotowania przedstawionych tabel i rycin.
Rysunek uzupełniający 1. Ogólny schemat badania w zakresie opracowania i wewnętrznej walidacji modelu predykcyjnego. Schemat przedstawia proces rekrutacji uczestników i ocenę kwalifikowalności, gromadzenie klinicznych charakterystyk matek oraz pomiary ultrasonograficzne w trakcie aktywnej pierwszej fazy porodu, klasyfikację do grup z normalnym i przedłużonym drugim okresem porodu zgodnie z 120-minutowym kryterium Światowej Organizacji Zdrowia (WHO), wybór predyktorów za pomocą regresji LASSO (least absolute shrinkage and selection operator), konstrukcję wieloczynnikowego modelu regresji logistycznej oraz nomogramu, a także wewnętrzną ewaluację modelu z wykorzystaniem analizy charakterystyki operacyjnej odbiornika (ROC), analizy kalibracji i analizy krzywej decyzyjnej (DCA). Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca 1. Zbiór danych na poziomie uczestnika wykorzystany do opracowania modelu i walidacji wewnętrznej. Tabela ta zawiera zanonimizowane dane na poziomie uczestników wykorzystane we wszystkich analizach statystycznych, w tym przy opracowywaniu i wewnętrznej walidacji modelu predykcyjnego. Zmienne obejmują wynik drugiego okresu porodu, czas trwania drugiego okresu, wyniki neonatologiczne, dane demograficzne i kliniczne matki oraz pomiary ultrasonograficzne transabdominalne i transperinealne wykonane podczas pierwszego okresu porodu. Body mass index (BMI) = wskaźnik masy ciała; occiput anterior (OA) = potyliczne przodowe położenie płodu; angle of progression (AOP) = kąt progresji; head–symphysis distance (HSD) = odległość głowa–łóneczko; oraz head–perineum distance (HPD) = odległość głowa–krocze. Prosimy kliknąć tutaj, aby pobrać ten plik.