Przedstawiona metoda ilościowo ocenia sprawność motoryczną niemowląt poprzez klasyfikację rodzajów postaw i ruchów dla każdej sekundy sesji nagraniowej. W związku z tym pakiet wyników z automatycznego potoku analizy obejmuje pełną macierz klasyfikacji (Plik uzupełniający 3) oraz podsumowanie graficzne (Plik uzupełniający 4) z całej sesji nagraniowej. W zależności od konkretnego pytania badawczego, wyniki te mogą być analizowane na różnych poziomach.
Do opracowania i walidacji tej metody wykorzystano weryfikację wyników.
Poniżej przedstawiamy cztery poziomy weryfikacji wyników wykorzystanych do opracowania i walidacji tej metodologii. Plik uzupełniający 5 zawiera reprezentatywne przykłady kluczowych eksperymentów walidacyjnych, które zostały wcześniej opublikowane w pełnym szczególe10,11,12.
W pierwszej kolejności algorytmy automatyczne wytrenowane do detekcji ruchu i postawy zostały zwalidowane względem sekunda po sekundzie przeprowadzonych obserwacji ludzkich dotyczących zachowań ruchowych niemowląt. Wykorzystaliśmy kilku równolegle szkolonych ekspertów, którzy w sposób zaślepiony analizowali zsynchronizowane nagrania wideo wraz z zapisami z urządzeń ubieralnych. Wszystkie poszczególne kategorie postaw i ruchów porównywano oddzielnie z indywidualnymi adnotacjami ludzkimi; stwierdzono bardzo wysoką zgodność między algorytmem a człowiekiem w przypadku postaw (średnia kappa 0,93), natomiast w przypadku kategorii ruchów stwierdzono znaczną zgodność (kappa dla poszczególnych podkategorii mieściła się głównie w zakresie 60–80%). Patrz Plik uzupełniający 5A dla przykładowej macierzy pomyłek10Ponadto oceniono poziom zgodności między sędziami, aby potwierdzić, że algorytmy działają na poziomie zbliżonym do ludzkiego.10,11.
Po drugie, oceniliśmy, w jakim stopniu kwantyfikacja kategorii ruchu i postawy oparta na klasyfikatorze pokrywa się z odpowiadającą jej kwantyfikacją z adnotacji ludzkich drugiego poziomu. Przykładowe wykresy rozrzutu przedstawione w Pliku uzupełniającym 5B10 wykazują, że dla kilku kluczowych kategorii występuje niemal idealna zgodność (współczynnik korelacji >0,96) pomiędzy algorytmiczną a ludzką kwantyfikacją wizualną. Bezpośrednio potwierdza to tezę, że specyficzne dla wieku rozkłady kwantyfikacji motoryki (Rysunek 3A,B) są wiarygodne12.
Po trzecie, koncepcja holistycznej oceny dojrzałości motorycznej została zwalidowana poprzez opracowanie modelu przewidywania wieku rozwojowego na podstawie kombinacji parametrów ruchu (patrz wyżej), który był bardzo silnie skorelowany z rzeczywistym wiekiem niemowląt o typowym rozwoju (r=0,89; Plik uzupełniający 5C). Następnie przewidywany wiek przeskalowano do wartości 0-100 jako bezwymiarową miarę BIMS (Baba Infant Motor Score10), a użyteczność tej miary do tworzenia wykresów wzrostu motorycznego (Plik uzupełniający 5D) zwalidowano na kohorcie niemowląt o typowym rozwoju, u których zaobserwowano wyraźne, zależne od wieku i przewidywalne trajektorie wzrostu. Potwierdziliśmy również względną dokładność tej metody, wykazując, że precyzja wykresów wzrostu motorycznego jest porównywalna z powszechnie znanymi wykresami wzrostu fizycznego12.
Po czwarte, potencjał wykrywania nieprawidłowości za pomocą podanych wskaźników został potwierdzony w eksperymencie dowodowym, w którym wykazano, że poszczególne parametry motoryczne wyraźnie różnicują niemowlęta o odpowiednio słabej i dobrej sprawności motorycznej (Plik uzupełniający 5E)10.
Potencjalne pytania do dalszych badań na różnych poziomach analizy
Rycina 3 przedstawia przykłady dalszego wykorzystania informacji dostarczanych przez kombinezon ubieralny oraz zautomatyzowany potok analizy. Po pierwsze, wczesny rozwój postawy i umiejętności ruchowych można przedstawić jako funkcję wieku i porównać z rozkładami zależnymi od wieku (Rycina 3A, „wykresy wzrostu”12) lub śledzić rozwój w czasie dla każdego osobno (Rycina 3B). Po drugie, gdy pytanie badawcze wymaga bardziej holistycznej oceny motoryki grubej, można zastosować kombinację proporcji postawy i ruchów danej osoby (jak pokazano w plikach uzupełniających5D, obliczonych z szeregów czasowych w pliku uzupełniającym 5C), aby wygenerować indeks taki jak BIMS (Baba Infant Motor Score). Takie miary umożliwiają bezpośrednie zastosowanie w wykresach wzrostu motorycznego (pliki uzupełniające 5C,D) oraz obliczanie pochodnych statystycznych, takich jak z-score (Rycina 3A). Wykorzystanie pełnych szeregów czasowych detekcji algorytmu (Rycina 3C oraz plik uzupełniający 3) pozwala na prowadzenie badań nad szczegółową strukturą czasową motoryki niemowląt, odpowiadając na pytania takie jak „Ile przejść między postawami wykonuje niemowlę w jednostce czasu?” lub „Jaki jest rozkład epok stania podczas spontanicznej zabawy?”.

Rysunek 1: Przegląd wielosensorowego systemu ubieralnego oraz typowy przebieg badania od rejestracji do analizy. Rysunek 1 dostosowany na podstawie pracy Airaksinen et al.12 opublikowany na licencji CC_BY. Zdjęcie niemowlęcia opublikowano za zgodą rodziców. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Kategoryzacja postawy i ruchu oraz przykładowa wizualizacja surowych danych i wyników analizy (A) Schemat kategoryzacji postawy i ruchu stosowany przez klasyfikatory algorytmiczne w pełni zautomatyzowanego rurociągu analitycznego dla urządzenia ubieralnego z wieloma czujnikami. Rysunek ten został przedrukowany z pracy Airaksinen et al.10 (B) Przykład 10 min surowych danych z czujnika akcelerometru z każdego z czterech ramion, pochodzących z zapisów MAIJU. Poziome paski poniżej przedstawiają wyniki automatycznych klasyfikatorów dla detekcji postawy (górny pasek) i ruchu (dolny pasek) dla tego samego 10 min interwału. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3Przykładowe wyniki uzyskane z surowych danych wyjściowych klasyfikatora. (A(A) Przykład porównania danych dotyczących wieku niemowlęcia oraz skali Baba Infant Motor Score (BIMS). Zielona krzywa w kształcie litery S przedstawia trajektorię rozwojową BIMS w większej populacji. Kropka przedstawia przykład pojedynczej osoby z pomiarem w 14. miesiącu życia, z wynikiem BIMS ok. 74, co odpowiada poziomowi nieco poniżej średniej typowej dla danego wieku (zielona linia). (B) Przykład kohorty z indywidualnymi trajektoriami rozwoju motorycznego z zastosowaniem skali Baba Infant Motor Score (BIMS). Każda linia reprezentuje niemowlę badane w kilku punktach czasowych w zależności od wieku (kropki na linii). Linie są zabarwione zgodnie z przeciętnym odchyleniem względem średniej typowej dla wieku (niebieska krzywa w kształcie litery S w tle; patrz również panel A).C) Reprezentatywna macierz wyjściowa z automatycznego klasyfikatora pochodząca z potoku przetwarzania Babacloud. Pierwsza kolumna przedstawia czas upływający od początku nagrania (w sekundach) dla każdego okna analizy w klasyfikacji (czas trwania okna 2,3 s, z 50% zakładką). Druga i trzecia kolumna pokazują detekcje klasyfikatora odpowiednio dla postawy i ruchu. Trzecia i czwarta kolumna to klasyfikatory pomocnicze przedstawiające epoki, w których niemowlę było przenoszone przez inną osobę oraz gdy niemowlę zajmowało się samodzielną zabawą. (D) Ilustracja raportu podsumowującego. Panele A i B zostały opracowane na podstawie pracy Airaksinen i wsp.10. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Plik uzupełniający 1: Krótki przewodnik po rejestracji za pomocą urządzenia ubieralnego MAIJU. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Krótki przewodnik po transferze danych z Maijulogger na serwer analizy w Babacloud. Prosimy kliknąć tutaj, aby pobrać ten plik.
Plik uzupełniający 3: Przykład szczegółowej macierzy klasyfikacji zawierającej wszystkie klasyfikacje sekunda po sekundzie dla postawy, ruchu, noszenia niemowlęcia oraz swobodnej zabawy. Prosimy kliknąć tutaj, aby pobrać ten plik.
Plik uzupełniający 4: Pełny przykład pliku PDF z graficznym raportem podsumowującym. Plik ten pochodzi z pliku algorytmu i zawiera (A) podstawowe informacje wstępne o sesji nagraniowej (ID i wiek badanego, datę nagrania, czas trwania oraz całkowitą liczbę epok wykorzystanych do końcowych analiz). (B) Graficzną prezentację pełnego nagrania, wskazującą rozkład postaw podczas całego zapisu oraz epoki wykluczone z oceny ilościowej. (C) Typowe rozkłady przedstawione na wykresach skrzypcowych (violin plots) dla wszystkich sześciu postaw (po lewej stronie) i 12 typów ruchów (po prawej stronie). Podobnie, rysunki po prawej stronie przedstawiają typy ruchów odpowiadające każdej postawie, ukazując również sekwencję stopniowo rozwijanych sprawności ruchowych (kropki wskazują wyniki z pojedynczego nagrania, a kolorowe wykresy skrzypcowe wskazują rozkłady postaw/ruchów zależne od wieku w odpowiednim zbiorze danych). Należy zaznaczyć, że surowe wartości wskazane przez kropki pokazują rzeczywiste ilości danej sprawności ruchowej i mogą być bezpośrednio wykorzystywane w innych kontekstach. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 5: Eksperymenty walidacyjne na różnych poziomach analizy. (A) Macierze pomyłek pokazujące zgodność między adnotacjami ludzkimi (klasa docelowa) a detekcjami algorytmu (klasa przewidywana) dla kategorii postawy oraz ruchu. (B) Porównanie ilościowych parametrów motorycznych w pełnych sesjach nagraniowych pomiędzy adnotacjami ludzkimi a detekcjami pochodzącymi z automatycznego klasyfikatora. (C) Korelacja między przewidywaniem wieku rozwojowego na podstawie danych z urządzeń ubieralnych (lewa oś Y) a jego przeliczeniem w celu wygenerowania wyniku BIMS (prawa oś Y). Rzeczywisty wiek niemowlęcia w czasie nagrania przedstawiono na osi X. (D) Korelacja przewidywanego wieku z wiekiem rzeczywistym przy zastosowaniu dopasowanej funkcji. Wartości obrazują dopasowanie do modelu przy użyciu średniej grupowej w przedstawionych oknach czasowych (kolor niebieski), wszystkich wartości surowych (kolor czarny) lub przy uwzględnieniu pomiarów powtórzonych dla każdej osoby (kolor zielony). (E) Porównanie indywidualnych miar motorycznych pomiędzy grupami niemowląt o dobrej i słabej sprawności sugeruje, że kilka automatycznie wykrytych wskaźników motorycznych może różnicować te grupy niemowląt. Panele A, B, C zostały opracowane na podstawie Airaksinen et al.10. Panel D został opracowany na podstawie Airaksinen et al.12. Panel E został opracowany na podstawie Airaksinen et al.11. Kliknij tutaj, aby pobrać ten plik.