Badanie zostało zatwierdzone przez Komisję Etyczną Szpitala Centralnego w Ziyang, filii szpitala West China Hospital w Ziyang przy Uniwersytecie w Sichuan (numer zatwierdzenia: 2026(115)). Przed włączeniem do badania od wszystkich uczestników lub ich ustawowych przedstawicieli uzyskano pisemną zgodę na udział. Wszystkie procedury przeprowadzono zgodnie z Deklaracją Helsińską oraz instytucjonalnymi wytycznymi dotyczącymi badań klinicznych z udziałem ludzi.
Projekt i warunki badania
Niniejsze jednocentrowe, nierandomizowane, kontrolowane badanie implementacyjne z grupami równoległymi oceniało ścieżkę opieki przeciwosteoporotycznej dla osób starszych po chirurgicznym leczeniu złamania biodra, prowadzoną przez pielęgniarki i dostosowaną do tolerancji żołądkowo-jelitowej (GI). Badanie przeprowadzono w Ziyang Central Hospital, West China Hospital Ziyang Branch, Sichuan University, w okresie od stycznia 2021 do grudnia 2024 roku. Uczestnicy byli obserwowani przez 24 miesiące po wypisie ze szpitala. Badanie zostało opracowane zgodnie z wytycznymi Strengthening the Reporting of Observational Studies in Epidemiology (STROBE)14. Przesiew uczestników, przydział do grup, harmonogram obserwacji oraz konstrukcja zbioru analizowego zostały podsumowane na Rysunku 1. Wstępnie określone definicje wyników oraz harmonogram oceny w 24. miesiącu zostały podsumowane w Tabeli uzupełniającej 1.
Uczestników przydzielono do grup na podstawie wdrożenia ścieżki prowadzonej przez pielęgniarki, a nie drogą randomizacji. Pacjenci przyjęci przed wdrożeniem ścieżki (styczeń 2021–grudzień 2022) otrzymali standardową opiekę, natomiast pacjenci przyjęci po wdrożeniu (styczeń 2023–grudzień 2024) byli prowadzeni zgodnie ze ścieżką prowadzoną przez pielęgniarki. Aby zminimalizować błąd czasowy, w obu grupach zastosowano identyczne kryteria kwalifikacji, harmonogramy obserwacji, definicje wyników, okna oceny DXA, obliczenia przestrzegania zaleceń oraz procedury rozstrzygania w kwestii ponownych złamań. Okres kalendarzowy został uwzględniony w analizach skorygowanych oraz analizach wrażliwości.
Uczestnicy
Kryteria kwalifikacji analizowano u kolejno przyjmowanych pacjentów hospitalizowanych w wieku ≥65 lat z niskoenergetycznym złamaniem szyjki kości udowej lub złamaniem międzykrotarowym biodra. Złamania niskoenergetyczne zdefiniowano jako złamania będące wynikiem upadku z wysokości własnego wzrostu lub niższej.
Kryteria kwalifikacji uczestników obejmowały: wiek ≥65 lat, przejście leczenia chirurgicznego z powodu niskoenergetycznego złamania szyjki kości udowej lub złamania międzykłykciowego, przeżycie do momentu wypisu ze szpitala, zakwalifikowanie do terapii przeciwosteoporotycznej po ocenie funkcji nerek, poziomu wapnia skorygowanego oraz statusu witaminy D, posiadanie kompletnych klinicznych i lekowych danych wyjściowych, zdolność do uczestnictwa w obserwacji bezpośrednio lub za pośrednictwem opiekuna oraz wyrażenie pisemnej świadomej zgody.
Z badania wykluczono uczestników, u których wystąpiła trauma wysokenergetyczna, złamania patologiczne wtórne do nowotworów lub metabolicznych chorób kości innych niż osteoporoza, przewidywany czas przeżycia wynosił <6 miesięcy, stwierdzono ciężką niewydolność nerek (szacowany współczynnik filtracji kłębuszkowej <30 mL/min/1.73 m2), nieskorygowaną hipokalcemię, aktywne krwawienie z górnego odcinka przewodu pokarmowego, ciężkie zwężenie przełyku, niestabilną chorobę przewodu pokarmowego wymagającą pilnego leczenia, stosowano terapię przeciwosteoporotyczną przez >3 miesiące w ciągu poprzedniego roku, brakowało możliwości niezawodnej obserwacji lub wsparcia opiekuna lub w przypadku odmowy udziału.
Ocena wyjściowa
Oceny wyjściowe przeprowadzono podczas indeksowej hospitalizacji, najlepiej w ciągu 48 h przed wypisem. Zmienne demograficzne i kliniczne obejmowały wiek, płeć, wskaźnik masy ciała, rodzaj złamania, procedurę chirurgiczną, czas od przyjęcia do operacji, czas pobytu w szpitalu, status palenia tytoniu, spożycie alkoholu, wcześniejsze złamania osteoporotyczne, choroby współistniejące, historię przyjmowanych leków, dostępność opiekuna oraz wyjściowe czynniki ryzyka GI. Obciążenie chorobami współistniejącymi oceniano za pomocą wskaźnika Charlson Comorbidity Index, a poziom kruchości oceniano w skali Clinical Frailty Scale.
Podstawowe badania laboratoryjne obejmowały oznaczanie stężenia wapnia w surowicy, albumin, wapnia skorygowanego, fosforanów, kreatyniny, szacowanego wskaźnika filtracji kłębuszkowej, fosfatazy alkalicznej, 25-hydroksywitaminy D oraz parathormonu, w zależności od dostępności. Przegląd przyjmowanych leków koncentrował się na niesteroidowych lekach przeciwzapalnych, lekach przeciwpłytkowych, przeciwzakrzepowych, inhibitorach pompy protonowej, glikokortykosteroidach, lekach uspokajających oraz wcześniej stosowanych lekach przeciwosteoporotycznych.
Gęstość mineralna kości (BMD) była oceniana za pomocą dwuenergetycznej absorpcjometrii promieni rentgenowskich (DXA) przed wypisem ze szpitala, o ile było to możliwe. W przypadku, gdy badanie DXA przed wypisem nie mogło zostać wykonane, skanowanie przeprowadzano w ciągu 4–6 tygodni po wypisie i uznawano je za ocenę wyjściową. Pomiary wykonywano w odcinku lędźwiowym kręgosłupa, szyjce kości udowej oraz w obrębie całej stawu biodrowego, przy czym do pomiarów biodra preferencyjnie wybierano stronę przeciwległą do złamania.
Standardowa opieka
Uczestnicy przydzieleni do grupy opieki standardowej otrzymali rutynowe postępowanie okołoperacyjne, standardowe zalecenia przy wypisie, prowadzoną przez lekarza terapię przeciwosteoporotyczną oraz rutynową kontrolę ambulatoryjną zgodnie z praktyką instytucji. Ustrukturyzowana stratyfikacja ryzyka GI, dobór leczenia prowadzony przez pielęgniarkę, proaktywne monitorowanie uzupełnień recept oraz zdefiniowane algorytmy modyfikacji leczenia nie były częścią opieki standardowej.
Ścieżka dostosowana do tolerancji układu pokarmowego prowadzona przez pielęgniarkę
Ścieżka opieki prowadzona przez pielęgniarkę została wdrożona przed wypisem ze szpitala i kontynuowana przez cały 24-miesięczny okres obserwacji. Specjalistyczna pielęgniarka ds. osteoporozy koordynowała opiekę pomiędzy chirurgami ortopedami, geriatrami lub endokrynologami, farmaceutami, personelem rehabilitacyjnym, usługami DXA i laboratoryjnymi oraz lokalnymi dostawcami opieki zdrowotnej. Schemat oceny ryzyka GI, wyboru leczenia, obserwacji oraz modyfikacji schematu terapeutycznego przedstawiono na Rysunku 2.
Przed wypisem ze szpitala tolerancję żołądkowo-jelitową (GI) sklasyfikowano jako niskie, umiarkowane lub wysokie ryzyko zgodnie z predefiniowanymi kryteriami klinicznymi (Tabela uzupełniająca 2). Klasyfikacja ryzyka uwzględniała historię schorzeń układu pokarmowego, zdolność połykania, możliwość pozostania w pozycji pionowej po podaniu leku, stosowane leki towarzyszące oraz aktualne objawy ze strony układu pokarmowego.
Wybór leczenia był indywidualnie dostosowany do ryzyka żołądkowo-jelitowego (GI), funkcji nerek, skorygowanego stężenia wapnia, poziomu 25-hydroksywitaminy D, ryzyka złamań, preferencji pacjenta oraz przewidywanego przestrzegania zaleceń w okresie obserwacji. Pacjentom z niskim ryzykiem GI preferencyjnie przepisywano bisfosfoniany doustne wraz ze zestandaryzowanymi instrukcjami podawania. Pacjenci z umiarkowanym ryzykiem otrzymywali terapię doustną tylko wtedy, gdy można było polegać na przestrzeganiu wymogów podawania leku; w przeciwnym razie zalecano terapię antyresorpcyjną niedoustną. Pacjenci zakwalifikowani do grupy wysokiego ryzyka GI otrzymywali preferencyjnie dożylny kwas zoledronowy lub podskórny denosumab po potwierdzeniu bezpieczeństwa leczenia i zatwierdzeniu przez lekarza. Zalecano suplementację wapniem i witaminą D, chyba że istnieją przeciwwskazania.
Przed wypisem z placówki pielęgniarka specjalistka wypełniła ustandaryzowaną listę kontrolną edukacji, obejmującą sposób podawania leków, harmonogram dawkowania, suplementację wapniem i witaminą D, rozpoznawanie objawów ostrzegawczych, procedury uzupełniania recept, terminy wizyt kontrolnych oraz dane kontaktowe. Pacjenci otrzymujący bisfosfoniany doustne przeszli dodatkową edukację w zakresie podawania leku na czczo z wodą, utrzymywania pozycji pionowej przez co najmniej 30 min po przyjęciu dawki, zachowania odstępu czasowego od posiłków i suplementów wapnia oraz niezwłocznego zgłaszania objawów ze strony układu pokarmowego.
Optymalizacja schematu leczenia oparta na wyzwalaczach
Optymalizacja leczenia przebiegała zgodnie z algorytmem decyzyjnym przedstawionym na Rysunku 2. Wyzwalaczami żołądkowo-jelitowymi były: nowy lub nasilający się refluks, nudności, wymioty, ból w nadbrzuszu, ból brzucha, dysfagia, zaparcie, podejrzenie krwawienia z przewodu pokarmowego lub jakikolwiek objaw żołądkowo-jelitowy utrudniający przestrzeganie zaleceń terapeutycznych. Wyzwalaczami dotyczącymi przestrzegania zaleceń były: pominięcie odnowienia recepty, przerwy w leczeniu doustnym przekraczające 30 dni, opóźnienie w podaniu iniekcji lub infuzji przekraczające 30 dni, wczesny wskaźnik PDC <0,80, zgłoszone przez pacjenta zaprzestanie leczenia lub zgłoszone przez opiekuna trudności w terapii.
Łagodne objawy stały się podstawą do wzmocnienia technik podawania leków, przeglądu leków przyjmowanych równocześnie przeprowadzonego przez farmaceutę, wdrożenia zaleceń dietetycznych i postępowania w przypadku zaparć oraz powtórnej oceny w ciągu 7–14 dni. Utrzymujące się objawy lub nawracający brak przestrzegania zaleceń terapeutycznych skutkowały wielodyscyplinarną ponowną oceną z udziałem lekarza prowadzącego, farmaceuty i pielęgniarki specjalistycznej. Postępowanie obejmowało tymczasowe przerwanie leczenia, terapię gastroprotekcyjną w odpowiednich przypadkach, przejście z terapii doustnej na iniekcyjną, uproszczenie schematu leczenia lub skierowanie do pielęgniarskiej opieki środowiskowej. Ciężkie objawy, podejrzenie krwawienia z przewodu pokarmowego, postępująca dysfagia, odwodnienie lub hospitalizacja wymagały pilnej oceny lekarskiej oraz odstawienia lub modyfikacji terapii.
Wyniki i obserwacje
Głównym punktem końcowym była ciągłość przyjmowania leków, oceniana za pomocą proporcji dni z pokryciem lekiem (PDC), osiągnięcia celu w zakresie przestrzegania zaleceń terapeutycznych oraz trwałości leczenia. Wtórne punkty końcowe obejmowały zmiany gęstości mineralnej kości (BMD) w określonych miejscach szkieletu oraz czas do wystąpienia pierwszego potwierdzonego kolejnego złamania kruchości. Punkty końcowe dotyczące bezpieczeństwa obejmowały zdarzenia niepożądane ze strony układu pokarmowego (GI), inne zdarzenia niepożądane związane z leczeniem, zmianę schematu dawkowania, tymczasowe przerwanie leczenia, trwałe zaprzestanie leczenia, hospitalizację związaną ze zdarzeniami niepożądanymi oraz zgon. Definicje punktów końcowych oraz harmonogram ocen w okresie 24 miesięcy przedstawiono w Tabeli uzupełniającej 1.
Oceny kontrolne przeprowadzono po 1 tygodniu oraz po 1, 3, 6, 12, 18 i 24 miesiącach od wypisu w obu grupach badawczych. Podczas każdej wizyty kontrolnej odnotowano ekspozycję na leki, status odnowienia recept, bariery w przestrzeganiu zaleceń terapeutycznych, objawy żołądkowo-jelitowe, inne zdarzenia niepożądane, upadki, podejrzenia ponownych złamań, hospitalizacje, zgony oraz modyfikacje leczenia. Badania DXA powtarzano w 12. i 24. miesiącu, gdy było to możliwe.
Przestrzeganie i trwałość farmakoterapii
Przestrzeganie zaleceń terapeutycznych oceniano za pomocą wskaźnika proporcji dni objętych terapią (PDC), który szacuje dostępność leków w określonym okresie obserwacji15. Okres obserwacji rozpoczynał się w dacie wypisu po indeksowej hospitalizacji z powodu złamania biodra i kończył po 24 miesiącach, w momencie zgonu, wycofania się z badania lub ostatniej potwierdzonej wizyty kontrolnej, w zależności od tego, co wystąpiło wcześniej. Wartości PDC dla konkretnych przedziałów czasowych obliczano w 1., 3., 6., 12., 18. i 24. miesiącu. PDC obliczano jako liczbę dni objętych terapią podzieloną przez całkowitą liczbę kwalifikowalnych dni obserwacji, z maksymalną wartością 1,00. Przyjęto, że PDC ≥0,80 wskazuje na odpowiednie przestrzeganie zaleceń, zgodnie z powszechnie akceptowanymi progami w badaniach nad osteoporozą16. Wskaźnik PDC interpretowano jako miarę dostępności leków, a nie potwierdzone przyjmowanie leków.
Dane dotyczące wydawania i uzupełniania leków pozyskano z dokumentacji szpitalnej apteki ambulatoryjnej, recept wypisanych przy wypisie, elektronicznej dokumentacji medycznej, zewnętrznej dokumentacji recept, paragonów z aptek ogólnodostępnych, opakowań leków przedstawionych lub sfotografowanych podczas wizyt kontrolnych oraz ustrukturyzowanych zapisów pielęgniarskich. Korzystanie z zewnętrznych aptek bez weryfikacji dokumentacyjnej zostało odnotowane, lecz wyłączone z głównej analizy PDC i uwzględnione jedynie w rozszerzonej analizie wrażliwości. W przypadku terapii doustnych liczbę dni pokrycia obliczono na podstawie wydanej ilości leku i przepisanego odstępu między dawkami. Każdy potwierdzony zastrzyk z denosumabem dostarczył 180 dni pokrycia, natomiast każda potwierdzona infuzja z kwasem zoledronicznym dostarczyła 365 dni pokrycia. Nakładające się zapasy leków przenoszono dalej bez podwójnego liczenia.
Utrzymanie leczenia zdefiniowano jako ciągłą terapię przeciwosteoporozie bez klinicznie istotnej przerwy w leczeniu. Brak utrzymania leczenia zdefiniowano jako przerwę przekraczającą 60 dni w przypadku terapii doustnej lub ponad 90 dni poza planowaną datą podania w przypadku terapii w formie wstrzyknięć lub infuzji. Odnotowano datę, czas trwania, przyczynę oraz podjęte działania naprawcze związane z każdą przerwą. Trwałe zaprzestanie leczenia zdefiniowano jako brak pokrycia lekami przeciwosteoporozowymi przez pozostałą część okresu obserwacji po przerwie w leczeniu.
Ocena BMD
Gęstość mineralna kości została zmierzona w odcinku lędźwiowym kręgosłupa, szyjce kości udowej oraz w obrębie całej głowy kości udowej przy użyciu absorpcjometrii rentgenowskiej o dwóch energiach (DXA). Bezwzględna zmiana BMD została wyrażona w g/cm2, a względną zmianę BMD obliczono następująco:

Zmiany BMD w konkretnych lokalizacjach sklasyfikowano przy użyciu metody najmniejszej istotnej zmiany (least significant change) stosowanej w seryjnej interpretacji DXA17. Wzrost BMD zdefiniowano jako relatywny przyrost >3,0% względem wartości wyjściowej, natomiast spadek zdefiniowano jako relatywną stratę >3,0%. Zmiany między −3,0% a +3,0% sklasyfikowano jako stabilny poziom BMD. Uczestników, u których brakowało poprawnych pomiarów DXA w punkcie wyjściowym i w badaniu kontrolnym w tej samej lokalizacji anatomicznej, wykluczono z analiz BMD dla konkretnych lokalizacji, lecz pozostali oni kwalifikowalni do analiz dotyczących przestrzegania zaleceń i ponownych złamań.
Stwierdzenie ponownego złamania
Kolejne złamanie fragilne zdefiniowano jako nowe złamanie niskoenergetyczne występujące po wypisie ze szpitala po pierwotnym przyjęciu z powodu złamania biodra. Podczas każdej wizyty kontrolnej uczestnicy byli badani pod kątem podejrzenia ponownych złamań poprzez rejestrację nowych upadków, bólu ogniskowego, wizyt na oddziale ratunkowym, hospitalizacji, badań obrazowych lub diagnoz złamań postawionych w innych placówkach. Ponowne złamania potwierdzano na podstawie raportów z badań obrazowych, dokumentacji szpitalnej, streszczeń wypisowych lub zewnętrznej dokumentacji medycznej. Zdarzenia zgłaszane wyłącznie podczas wywiadów, bez potwierdzenia radiologicznego lub w dokumentacji medycznej, klasyfikowano jako podejrzenia ponownych złamań i wykluczano z analizy głównej.
Datę zdarzenia zdefiniowano jako datę urazu, jeśli była dostępna, a w przeciwnym razie jako datę pierwszej wizyty lekarskiej. Do głównej analizy czasu do wystąpienia zdarzenia włączono tylko pierwsze potwierdzone ponowne złamanie. Odnotowano lokalizację złamania, źródło potwierdzenia oraz odstęp czasu między wystąpieniem objawów a potwierdzeniem diagnozy. Wszystkie przypadki ponownych złamań zostały niezależnie przeanalizowane przez dwóch klinicystów, a rozbieżności rozstrzygał starszy specjalista ortopedii.
Stopniowanie zdarzeń niepożądanych w obrębie układu pokarmowego
Niepożądane zdarzenia żołądkowo-jelitowe oceniano w skali klinicznego stopnia nasilenia obejmującej cztery zdefiniowane poziomy. Zdarzenia stopnia 1 obejmowały łagodne objawy, które nie wymagały przerwania leczenia ani dodatkowych leków. Zdarzenia stopnia 2 obejmowały objawy umiarkowane, wymagające dodatkowej opieki medycznej, konsultacji z farmaceutą lub lekarzem, zintensyfikowanego monitorowania lub przerwania leczenia na okres ≤14 dni. Zdarzenia stopnia 3 obejmowały objawy skutkujące zaprzestaniem stosowania obecnego schematu leczenia, przerwaniem leczenia powyżej 14 dni, badaniem ambulatoryjnym lub w trybie pilnym bądź zmianą drogi podania. Zdarzenia stopnia 4 stanowiły poważne powikłania żołądkowo-jelitowe, w tym podejrzenie lub potwierdzone krwawienie z przewodu pokarmowego, perforację przewodu pokarmowego, ciężkie odwodnienie lub hospitalizację. Związek każdego zdarzenia z terapią przeciwosteoporotyczną był klasyfikowany przez lekarza prowadzącego jako brak związku, możliwy związek, prawdopodobny związek lub definite związek.
Wierność implementacji i progi KPI
Wierność implementacji oceniono za pomocą zdefiniowanych wcześniej wskaźników, obejmujących: przeprowadzenie stratyfikacji ryzyka GI, wypełnienie listy kontrolnej edukacji przy wypisie, przegląd wyjściowych badań laboratoryjnych w zakresie bezpieczeństwa, rozpoczęcie leczenia w ciągu 1 miesiąca oraz w ciągu 90 dni, suplementację wapniem i witaminą D w razie wskazań, realizację zaplanowanych wizyt kontrolnych, wykonanie badania DXA, dokumentację ocen zdarzeń niepożądanych GI, realizację interwencji opartych na wyzwalaczach oraz kompletność oceny ponownych złamań. Wcześniej określone progi kluczowych wskaźników efektywności (KPI) oraz kryteria działań korygujących przedstawiono w Rysunku 3 i Tabeli uzupełniającej 3.
Predefiniowane cele KPI były następujące: przeprowadzenie stratyfikacji ryzyka GI w ≥90% przypadków do momentu wypisu lub w ciągu 1 tygodnia; wypełnienie listy kontrolnej edukacji przy wypisie w ≥85% przypadków; przegląd wyjściowych badań laboratoryjnych bezpieczeństwa w ≥85% przypadków przed wyborem leczenia; rozpoczęcie leczenia w ≥70% przypadków w ciągu 1 miesiąca oraz w ≥85% przypadków w ciągu 90 dni; ukończenie obserwacji w ≥80% przypadków do 6 miesięcy i w ≥75% przypadków w 24 miesiącach; wykonanie badania DXA w ≥70% przypadków w 24 miesiącach; udokumentowana ocena działań niepożądanych GI w ≥85% przypadków; oraz potwierdzona adjudykacja ponownych złamań w ≥95% przypadków. Nieosiągnięcie predefiniowanych celów KPI skutkowało dodatkowym kontaktem z pacjentem, zaangażowaniem opiekuna, ponownym przeszkoleniem personelu, współpracą z farmaceutą, wsparciem pielęgniarstwa środowiskowego lub ponowną oceną schematu leczenia.
Wielkość próby
Obliczenie wielkości próby oparto na głównym punkcie końcowym dotyczącym przestrzegania zaleceń terapeutycznych. Dane pilotażowe z ośrodka badawczego oraz opublikowane badania nad przestrzeganiem zaleceń w leczeniu osteoporozy sugerowały średnią wartość PDC w okresie 24 miesięcy wynoszącą 0,55 w grupie standardowej opieki oraz 0,68 w grupie ścieżki terapeutycznej, co odpowiadało przewidywanej różnicy międzygrupowej na poziomie 0,13. Przy założeniu połączonego odchylenia standardowego wynoszącego 0,25, dwustronnej wartości α = 0,05 oraz mocy statystycznej 80%, wymagano minimum 59 uczestników w każdej grupie. Uwzględniając 20% wskaźnik rezygnacji w trakcie 24-miesięcznego okresu obserwacji, docelowa liczba rekrutowanych osób wynosiła około 74 uczestników na grupę. Końcowa kohorta badawcza składała się ze 147 uczestników, w tym 73 w grupie ścieżki terapeutycznej i 74 w grupie standardowej opieki. Badanie nie miało mocy statystycznej pozwalającej na wykrycie różnic w odniesieniu do wyników ponownych złamań; w związku z tym analizy ponownych złamań uznano za wtórne i eksploracyjne.
Analiza statystyczna
Analizy statystyczne przeprowadzono przy użyciu programu R w wersji 4.3.2 oraz SPSS w wersji 27.0. Wszystkie testy statystyczne były dwustronne, a wartość P < 0,05 uznano za statystycznie istotną. Zmienne ciągłe przedstawiono jako średnią ± odchylenie standardowe lub medianę (rozstęp międzykwartylowy), w zależności od potrzeb, natomiast zmienne kategoryczne przedstawiono jako liczebności i wartości procentowe. Porównania międzygrupowe przeprowadzono przy użyciu testu t Studenta lub testu U Manna–Whitneya dla zmiennych ciągłych oraz testu chi-kwadrat lub dokładnego testu Fishera dla zmiennych kategorycznych. Porównywalność wyjściową oceniono za pomocą wartości P oraz standaryzowanych różnic średnich.
Analiza główna została przeprowadzona zgodnie z zasadą intencji leczenia (intention-to-treat), w której wszystkich uczestników analizowano zgodnie z ich pierwotnym przydziałem do grup. Przeprowadzono również analizę per-protocol po wykluczeniu uczestników z istotnymi odstępstwami od protokołu, którzy wycofali zgodę lub dla których nie było dostępnych danych dotyczących przestrzegania zaleceń po wypisie.
Wskaźnik PDC analizowano zarówno jako wynik ciągły, jak i jako zmienną binarną, stosując zdefiniowany próg przestrzegania zaleceń PDC ≥ 0,80. Wartości ciągłe PDC oceniano za pomocą regresji liniowej, natomiast binarne wyniki przestrzegania zaleceń analizowano przy użyciu regresji log-dwumianowej lub zmodyfikowanej regresji Poissona z odpornymi błędami standardowymi. Modele wielowymiarowe skorygowano o wiek, płeć, rodzaj złamania, procedurę chirurgiczną, wskaźnik chorób współistniejących Charlsona, wcześniejsze złamania osteoporotyczne, wyjściową kategorię ryzyka GI, szacowany współczynnik filtracji kłębuszkowej, wyjściowe stężenie 25-hydroksywitaminy D, wyjściową gęstość mineralną kości (BMD) oraz okres kalendarzowy.
Zmiany BMD w czasie analizowano za pomocą liniowych modeli mieszanych z efektami stałymi dla grupy terapeutycznej, czasu oraz interakcji grupy z czasem, wraz z losowymi punktami przecięcia dla poszczególnych uczestników. Oddzielne analizy przeprowadzono dla kręgosłupa lędźwiowego, całkowitej gęstości kości biodra oraz szyjki kości udowej. Czas do pierwszego potwierdzonego ponownego złamania oceniano przy użyciu krzywych przeżycia Kaplana–Meiera, testu log-rank oraz regresji proporcjonalnych hazardów Coxa. Założenie o proporcjonalności hazardów zweryfikowano za pomocą reszt Schoenfelda. Ze względu na ograniczoną liczbę zdarzeń ponownych złamań, szacunki współczynnika hazardu zinterpretowano jako wstępne. Zgon przed wystąpieniem ponownego złamania oceniano za pomocą regresji ryzyk konkurencyjnych Fine’a–Graya18.
Brakujące dane podsumowano według grupy badawczej, zmiennej, punktu czasowego obserwacji oraz przyczyny braku danych. W przypadku, gdy założenie o losowości braków (missing-at-random) uznano za prawdopodobne, zastosowano wielokrotną imputację za pomocą równań łańcuchowych (multiple imputation by chained equations) dla brakujących współzmiennych i powtarzanych wyników zdarzeń niekorzystnych19. Ze względu na różnice w stopniu ukończenia obserwacji pomiędzy grupami badawczymi, przeprowadzono wcześniej określone analizy wrażliwości dla danych brakujących w sposób nielosowy (missing-not-at-random). W analizach przestrzegania zaleceń zastosowano konserwatywny scenariusz najgorszego przypadku, w którym uczestników utraconych z obserwacji w grupie ścieżki terapeutycznej zaklasyfikowano jako nieprzestrzegających zaleceń, a utraconych w grupie standardowej opieki jako przestrzegających zaleceń w 24. miesiącu. W analizach BMD brakujące wartości DXA w grupie ścieżki terapeutycznej imputowano do 10. percentyla obserwowanych zmian, natomiast brakujące wartości w grupie standardowej opieki imputowano do 90. percentyla. W analizach ponownych złamań konserwatywny scenariusz przypisania zdarzeń zakładał, że uczestnicy utraceni z obserwacji w grupie ścieżki terapeutycznej doznali ponownego złamania w połowie okresu między ostatnim kontaktem a 24. miesiącem, podczas gdy uczestnicy utraceni w grupie standardowej opieki zostali ocenzurowani w momencie ostatniej potwierdzonej wizyty kontrolnej. Przeprowadzono analizę punktu krytycznego (tipping-point analysis), zmieniając założone wyniki wśród uczestników utraconych z obserwacji do momentu, w którym główne wnioski dotyczące różnic między grupami przestały być istotne statystycznie lub klinicznie.
Wstępnie określone analizy w podgrupach oceniły wpływ wieku, płci, kategorii ryzyka żołądkowo-jelitowego, funkcji nerek, wyjściowej kategorii BMD oraz wcześniejszych złamań osteoporotycznych. W modelach regresji oceniono człony interakcji, a wszystkie analizy w podgrupach uznano za eksploracyjne.