$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Przedstawiona metoda określa ilościowo sprawność motoryczną niemowląt, klasyfikując typy pozycji i ruchów w każdej sekundzie sesji nagraniowej. W związku z tym pakiet wyników z potoku automatycznej analizy zawiera pełną matrycę klasyfikacyjną (plik uzupełniający 3) oraz graficzne podsumowanie (plik uzupełniający 4) z całej sesji nagraniowej. W zależności od konkretnego pytania badawczego, wyniki te mogą być sprawdzane na różnych poziomach.
Inspekcja wyników została wykorzystana do opracowania i walidacji tej metody.
Poniżej przedstawiamy cztery poziomy kontroli wyników wykorzystywane do opracowania i walidacji tej metodologii. Plik uzupełniający 5 przedstawia reprezentatywne przykłady eksperymentów walidacji kluczy, które zostały wcześniej opublikowane w pełnym zakresie10,11,12.
Po pierwsze, automatyczne algorytmy wyszkolone do wykrywania ruchu i postawy zostały zweryfikowane pod kątem obserwacji zachowań ruchowych niemowląt na poziomie sekunda po sekundzie. Korzystaliśmy z usług kilku równolegle wyszkolonych ekspertów, którzy na ślepo przeglądali zsynchronizowane nagrania wideo z nagraniami ubieralnymi. Wszystkie różne kategorie postawy i ruchu zostały porównane oddzielnie z indywidualnymi adnotacjami ludzkimi i znaleźliśmy bardzo wysoką zgodność między algorytmem a człowiekiem w odniesieniu do postaw (średnia kappa 0,93); Stwierdzono zasadniczą zgodność w odniesieniu do kategorii przemieszczania (kappa specyficzna dla podkategorii, głównie w przedziale 60-80%). Zobacz plik uzupełniający 5A, aby zapoznać się z przykładową macierzą pomyłek10. Oceniono również poziomy zgodności między oceniającymi, aby potwierdzić, że algorytmy działają na poziomie zbliżonym do ludzkiego równoważnego10,11.
Po drugie, oceniliśmy, jak dobrze kwantyfikacja kategorii ruchu i postawy oparta na klasyfikatorze będzie pasować do odpowiadającej jej kwantyfikacji z adnotacji drugiego poziomu dla ludzi. Przykładowe wykresy punktowe pokazane w pliku uzupełniającym 5B10 pokazują, że kilka kluczowych kategorii ma prawie idealne dopasowanie (współczynnik korelacji >0,96) między kwantyfikacją algorytmiczną a wizualną przez człowieka. To bezpośrednio potwierdza ideę, że specyficzne dla wieku rozkłady ilościowe ruchliwości (Rysunek 3A,B) są wiarygodne12.
Po trzecie, idea holistycznej oceny dojrzałości motorycznej została potwierdzona przez trening przewidywania wieku rozwojowego na podstawie kombinacji wielkości ruchu (patrz wyżej), które bardzo ściśle korelowały z rzeczywistym wiekiem u typowo rozwijających się niemowląt (r=0,89; Plik uzupełniający 5C). Następnie przewidywanie wieku zostało przeskalowane do 0-100 jako miara bez jednostek BIMS (Baba Infant Motor Score10), a jego użyteczność do budowania wykresów wzrostu motorycznego (plik uzupełniający 5D) została zweryfikowana przy użyciu typowo rozwijającej się kohorty niemowląt, która wykazuje zależne od wieku i przewidywalne trajektorie wzrostu. Potwierdziliśmy również jego względną dokładność, pokazując, że dokładność wykresów wzrostu motorycznego dobrze wypada w porównaniu z dobrze znanymi wykresami wzrostu fizycznego12.
Po czwarte, możliwość wykrywania nieprawidłowości za pomocą podanych wskaźników została zweryfikowana w eksperymencie potwierdzającym zasadność, w którym wykazano, że poszczególne miary motoryczne wyraźnie odróżniają niemowlęta odpowiednio od słabych i dobrych sprawności motorycznych (plik uzupełniający 5E)10.
Potencjalne pytania do dalszego badania na różnych poziomach analiz
Rysunek 3 pokazuje przykłady dalszych zastosowań informacji dostarczanych przez kombinezon do noszenia i jego zautomatyzowany potok analizy. Po pierwsze, wczesny rozwój postawy i umiejętności ruchowych można wykreślić jako funkcję wieku i porównać z rozkładami zależnymi od wieku (Rysunek 3A, "wykresy wzrostu"12), lub rozwój można śledzić w czasie dla każdej osoby (Rysunek 3B). Po drugie, gdy pytanie badawcze wymaga bardziej holistycznej oceny motoryki dużej, można użyć kombinacji postawy i proporcji ruchowych danej osoby (jak pokazano w pliku uzupełniającym 5D, obliczonym na podstawie szeregów czasowych w pliku uzupełniającym 5C) w celu wygenerowania indeksu, takiego jak BIMS (Baba Infant Motor Score). Takie miary wspierają bezpośrednie zastosowanie w wykresach wzrostu motorycznego (plik uzupełniający 5C,D) i obliczaniu pochodnych statystycznych, takich jak wynik z (Rysunek 3A). Wykorzystanie pełnej serii detekcji algorytmu (Rysunek 3C i plik uzupełniający 3) pozwala na badanie szczegółowej struktury czasowej ruchliwości niemowląt z pytaniami takimi jak "Ile zmian postawy wykonuje niemowlę w jednostce czasu?" lub "Jaki jest rozkład epok stania podczas spontanicznej zabawy?".

Rysunek 1: Przegląd wielosensorowego systemu ubieralnego i typowy przebieg badania od nagrania do analizy. Rysunek 1 zaadaptowany z Airaksinen et al.12opublikowany na licencji CC_BY. Zdjęcie niemowlęcia jest publikowane za zgodą rodziców. Kliknij tutaj, aby zobaczyć większą 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 używany przez klasyfikatory algorytmiczne w ramach w pełni zautomatyzowanego potoku analizy dla wielosensorowego urządzenia do noszenia. Ten rysunek jest przedrukowany z Airaksinen et al.10 (B) Przykład 10-minutowego surowego czujnika akcelerometru z każdego z czterech ramion, które pochodzą z nagrań MAIJU. Poziome paski poniżej przedstawiają dane wyjściowe automatycznego klasyfikatora dla wykrywania postawy (górny pasek) i ruchu (dolny pasek) dla tej samej 10-minutowej epoki. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Przykładowe wyniki pochodzące z surowych danych wyjściowych klasyfikatora. (A) Przykład porównania danych między wiekiem niemowlęcia a wynikiem motorycznym niemowląt Baba (BIMS). Zielona krzywa w kształcie litery S przedstawia trajektorię rozwoju BIMS w większej populacji. Kropka przedstawia przykładową osobę mierzoną w wieku 14 miesięcy, z BIMS ~74, co odpowiada nieco poniżej średniego poziomu typowego dla wieku (zielona linia). (B) Przykładowa kohorta z indywidualnymi trajektoriami rozwoju motorycznego z wykorzystaniem Baba Infant Motor Score over (BIMS). Każda linia reprezentuje niemowlę zarejestrowane w kilku punktach wieku (kropki w linii). Linie są pokolorowane dla średniego odchylenia w stosunku do ś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, ponieważ pochodzi z potoku Babacloud. Pierwsza kolumna przedstawia czas, jaki upłynął od rozpoczęcia rejestracji (w sekundach) dla każdego okna analizy w klasyfikacji (czas trwania okna 2,3 s, przy 50% nałożeniu). Druga i trzecia kolumna pokazują wykrywanie klasyfikatora odpowiednio dla postawy i ruchu. Trzecia i czwarta kolumna to pomocnicze klasyfikatory przedstawiające epoki, w których niemowlę było noszone przez kogoś innego i gdy niemowlę angażowało się w samodzielną zabawę. (D) Ilustracja sprawozdania podsumowującego. Panele A i B są zaadaptowane z Airaksinen et al.10. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Plik uzupełniający 1: Krótki przewodnik po nagrywaniu za pomocą urządzenia do noszenia na MAIJU. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Krótki przewodnik po przesyłaniu danych z Maijulogger do serwera analitycznego w Babacloud. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 3: Przykład szczegółowej tabeli klasyfikacji, która podaje wszystkie klasyfikacje sekunda po sekundzie dla postawy, ruchu, noszenia niemowlęcia i czasu na swobodną zabawę. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 4: Pełny plik PDF przykładowy graficznego raportu podsumowującego. Jest to pobierane z pliku algorytmu, który zawiera: (A) podstawowe informacje o sesji nagraniowej (identyfikator i wiek podmiotu, datę nagrania, czas trwania i całkowitą liczbę epok użytych do końcowych analiz). (B) Graficzne przedstawienie pełnego zapisu wskazujące rozkład postaw w pełnym nagraniu i epokach, które są wyłączone z oceny ilościowej. C) Typowe rozkłady wyświetlane z wykresami skrzypiec dla wszystkich sześciu pozycji (lewa strona) i 12 typów ruchu (prawa strona). Podobnie, rysunki po prawej stronie przedstawiają rodzaj ruchów wskazywanych przez każdą postawę, pokazując również sekwencję stopniowo rozwijającej się sprawności motorycznej (kropki wskazujące wyniki z indywidualnego nagrania i kolorowe wykresy skrzypiec wskazujące zależne od wieku rozkłady postawy/ruchu odpowiedniego zestawu danych). Warto zauważyć, że surowe wartości oznaczone kropkami pokazują rzeczywiste wielkości danej wydajności silnika i mogą być używane bezpośrednio 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) zarówno dla kategorii postawy, jak i ruchu. (B) Porównanie ilościowych parametrów motorycznych podczas pełnych sesji nagraniowych między adnotacjami ludzkimi a detekcjami pochodzącymi z automatycznego klasyfikatora. (C) Korelacja między przewidywaniem wieku rozwojowego na podstawie danych z urządzeń do noszenia (oś Y po lewej stronie) a ich ponownym skalowaniem w celu wygenerowania wyniku BIMS (oś Y po prawej stronie). Rzeczywisty wiek niemowlęcia w momencie nagrywania jest pokazany na osi X. (D) Korelacja między przewidywanym wiekiem a rzeczywistym wiekiem w przypadku korzystania z dopasowanej funkcji. Wartości przedstawiają dopasowanie do modelu przy użyciu średniej grupowej w przedstawionych oknach czasowych (kolor niebieski), wszystkich wartości surowych (kolor) lub przy uwzględnianiu powtarzających się pomiarów każdej osoby (kolor zielony). (E) Porównanie indywidualnych miar motorycznych między grupami niemowląt osiągających dobre i słabe wyniki sugeruje, że kilka automatycznie wykrywanych wskaźników motorycznych może różnicować te grupy niemowląt. Panel A,B,C został zaadaptowany z Airaksinen et al.10. Panel D jest zaadaptowany z Airaksinen et al.12. Panel E jest zaadaptowany z "Airaksinen et al.11. Kliknij tutaj, aby pobrać ten plik.