$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Jeśli chodzi o wyniki rekrutacji, uczestnicy byli rekrutowani głównie poprzez wysyłanie listów rekrutacyjnych i kolejne rozmowy telefoniczne w oparciu o nakreślone przepisy Systemu Opieki Zdrowotnej w Atlancie VA. Zespół badawczy zrekrutował łącznie 50 uczestników, udowadniając skuteczność zastosowanych metod w realizacji celu rekrutacyjnego (patrz Rysunek 2). Zastosowanie nowych klinicznych kryteriów diagnostycznych fibromialgii pozwoliło zespołowi badawczemu na prawidłowe odsianie osób, które nie spełniały kryteriów fibromialgii32. Pytanie potencjalnych uczestników o diagnozę fibromialgii nie jest tak solidnym środkiem, jak dodatkowe badania przesiewowe i mogło prowadzić do niewłaściwie zaplanowanych wizyt wyjściowych lub udziału w badaniu. Czterdziestu ośmiu uczestników zostało losowo przydzielonych do grup aktywnych i pozorowanych; Dwóch uczestników zostało wykluczonych, ponieważ nie spełniali kryteriów kwalifikacyjnych opartych na testach włączenia do badania.

Rysunek 2: Schemat rekrutacji. Raport i diagram przepływu rekrutacji do badania, randomizacji i alokacji interwencji. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
W obliczeniach wielkości próby użyto trzech wyników: DVPRS (ból kliniczny), test 30-sekundowego stania na krześle (funkcja) i łączność DMN-SMN (rs-fcMRI). Wszystkie obliczenia mocy opierały się na danych wstępnych. Jako główny punkt końcowy zainteresowania wybrano kliniczne zmiany bólu za pomocą DVPRS. Test 30-sekundowego stania na krześle (30sCST) został wybrany jako reprezentatywny wynik funkcjonalny, ponieważ wykazał najmniejszą zmianę między grupami. Łączność DMN-SMN została wybrana jako drugorzędny wynik zainteresowania jako biomarker neuroobrazowania dla bólu klinicznego i odpowiedzi na leczenie. Analizy wielkości próby przeprowadzono przy założeniu istotności mocy 1% i 80% (2-próbowa, 1-stronna).
Nasiona do tej analizy zostały wybrane na podstawie wstępnych danych, a także literatury na temat fibromialgii, bólu i CES16,24,26,63,64,65,66,67,68,69. Na podstawie wstępnych danych średnia (SD) zmiany po leczeniu w łączeniu lewej pierwszorzędowej kory czuciowo-ruchowej (L-S1M1) z lewą tylną korą obręczy (L-PCC) wynosi 0,041 (0,079) dla grupy leczonej i -0,026 (0,049) dla standardowej grupy leczenia; obserwowana wielkość efektu różnicy między grupami wynosi 1,03 (patrz Tabela 216,21,26,62,63,64,70,71,72,73). Badanie wymagałoby 20 osób w każdej grupie, aby osiągnąć 80% mocy do wykrycia różnicy między grupą CES a grupą standardowego leczenia w ich zmianie po leczeniu w łączności dla L-S1M1 do L-PCC na poziomie istotności 0,01 przy użyciu dwustronnego testu t, przy założeniu, że wielkość efektu różnicy między grupami wynosi 1,03, jak zaobserwowano w danych pilotażowych. Chociaż w badaniu pilotażowym zaobserwowano 17% wskaźnik wyczerpania u 12 osób w naszym poprzednim badaniu neuromodulacji małżowiny usznej (wszyscy ukończyli kontrolny rezonans magnetyczny, ale dwóch zostało utraconych na obserwację po 8 i 12 tygodniach), w celu utrzymania konserwatywnego szacunku dla obliczeń wielkości próby, w tym badaniu założono 20% wskaźnik wyczerpania. Przy oczekiwanym 20% zniszczeniu podczas wizyty po leczeniu, badanie wymagało rekrutacji 20/0,8 = 25 osób w grupie.
| Nasiona DMN (x,y,z) | Nasiona SMN (x,y,z) | Nasiona SN (x,y,z) |
| Przyśrodkowa kora przedczołowa62,67 | Prawy skorupiak64,71 | Prawa grzbietowo-boczna kora przedczołowa62 |
| Prawy PCC70 | Lewa klasa M170 | Lewa wyspa przednia 62 |
| Lewy PCC16,67 | Prawy M170 | Prawa wyspa przednia przednia62 |
| Precuneus71 | Prawa ręka S116,72 | Lewa tylna wyspa64 |
| Lewa ręka S1<-sup class="xref">16,72 | Prawa tylna izolacja63 |
| Thalamus21 | Grzbietowa przednia kora obręczy72,73 |
| | Prawe skrzyżowanie skroniowo-ciemieniowe62 |
Tabela 2: Nasiona do analizy. Nasiona DMN, SMN i SN wybrane do analizy na podstawie hipotez a priori. Każde nasiono jest opatrzone odniesieniami do wcześniejszej literatury potwierdzającej jego testowanie w zespołach bólowych.
Na podstawie wcześniejszych badań CES u cywilnych osób z fibromialgią, łącznie 50 osób (n = 25 pozorowane i n = 25 prawdziwe) powinno osiągnąć 80% mocy wykrywania różnicy w wynikach bólu między dwiema grupami17 (patrz Tabela 3). Obliczenia wielkości próby dla tego badania przeprowadzono przy użyciu sealedenvelope.com i oparto na danych wstępnych.
pkt.
TGL
jedn.
TGL
szt.
Szt.
TGL
pkt.
punktu
pkt.
pkt.
| Kontrola | Interwencja | N na grupę |
| Średnia zmiana | Sd | Średnia zmiana | Sd |
| Ból kliniczny (DVPRS) | 0,375 | 1,493 | -1,833 | Pytanie 2,229 | 10 |
| Funkcja (30sCST) | -0,250 | 1.500 | 3.000 | 4,980 | 14 |
| rs-fcMRI (S1M1 do PCC) | -0,026 | 0,049 | 0,041 | 0,079 | 20 |
Tabela 3: Obliczanie wielkości próby. Obliczenia związane z wielkością próby badawczej.
Zaimportowane dane funkcjonalne fMRIPrep zostały wygładzone za pomocą konwolucji przestrzennej z jądrem Gaussa o pełnej szerokości 8 mm i połowie maksimum (FWHM) (zobacz Rysunek 3) pokazuje funkcjonalny wynik z fMRIPrep znormalizowany do przestrzeni szablonu MNI152NLin2009cAsym (po lewej) i wygładzony obraz funkcjonalny z CONN Toolbox (po prawej). Powoduje to zwiększenie stosunku sygnału do szumu, co z kolei poprawia wykrywanie sygnałów zależnych od poziomu tlenu we krwi (BOLD).

Rysunek 3: Pojedynczy obiekt porównujący niewygładzony funkcjonalny obraz w przestrzeni MNI (po lewej) z jego wygładzonym odpowiednikiem przy 8 mm FWHM. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Dane zostały następnie odszumione za pomocą standardowego potoku odszumiania74, w tym regresji potencjalnych efektów zakłócających charakteryzujących szeregi czasowe istoty białej (5 składowych szumu CompCor), szeregi czasowe CSF (5 składowych szumu CompCor), parametry ruchu i ich pochodne pierwszego rzędu (12 czynników)75, skany wartości odstających (poniżej 295 czynników)48oraz trendy liniowe (2 czynniki) w każdym przebiegu funkcjonalnym, po których następuje filtrowanie częstotliwości pasmowoprzepustowej szeregów czasowych BOLD76 między 0,008 Hz a 0,09 Hz. CompCor49,77 składowe szumu w istocie białej i płynie mózgowo-rdzeniowym oszacowano, obliczając średni sygnał BOLD, a także największe składowe główne prostopadłe do średniej BOLD, parametry ruchu, oraz skany wartości odstających w obrębie zerodowanych masek segmentacji każdego podmiotu (patrz Rysunek 4). Na podstawie liczby warunków szumu zawartych w tej strategii odszumiania, oszacowano, że efektywne stopnie swobody sygnału BOLD po odszumieniu wahają się od 33 do 240,6 (średnio 173,4) u wszystkich badanych. Odszumianie spowodowało redukcję szumów fizjologicznych i innych zewnętrznych szumów pochodzących z danych, które mogły skutkować efektami zakłócającymi.

Rysunek 4: Kontrole jakości. Wykresy kontroli jakości z CONN Toolbox przedstawiające wpływ odszumiania na łączność funkcjonalną (FC), średni sygnał globalny i maksymalny ruch. Uwaga (A) górny wykres przedstawia dane dla pojedynczego tematu i sesji, wykresy (B,C) B i C są wynikiem odszumiania na poziomie grupy. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Analizy na poziomie grupy zostały przeprowadzone przy użyciu Ogólnego Modelu Liniowego (GLM)78. Dla każdego pojedynczego woksela oszacowano oddzielny GLM, z miarami łączności pierwszego poziomu na tym wokselu jako zmiennymi zależnymi (jedna niezależna próbka na obiekt i jeden pomiar na zadanie lub warunek eksperymentalny, jeśli dotyczy), a grupy lub inne identyfikatory na poziomie podmiotu jako zmienne niezależne. Hipotezy na poziomie wokseli oceniano przy użyciu wielowymiarowej statystyki parametrycznej z losowymi efektami u badanych i estymacją kowariancji próbki w wielu pomiarach. Wnioskowania przeprowadzono na poziomie poszczególnych klastrów (grup sąsiadujących ze sobą wokseli). Wnioskowanie na poziomie klastra oparto na statystykach parametrycznych z teorii losowego pola Gaussa62,79. Wyniki zostały określone przy użyciu kombinacji progu na poziomie woksela tworzącego klaster p < 0,005 oraz progu wielkości klastra z rodzinną korekcją błędów p < 0,001680. Wykonanie tych kroków powoduje porównanie wartości łączności funkcjonalnej na poziomie grupy i warunków pozornych na podstawie regionów zainteresowania (ROI). Wyniki te można wizualizować na wiele sposobów za pomocą graficznego interfejsu użytkownika eksploratora wyników w programie CONN Toolbox. Aby zobaczyć wizualizację wyświetlania głośności przykładowych wyników na poziomie grupy, gdzie czerwony obszar wskazuje regiony o większej dodatniej łączności z ROI, a niebieski to regiony z większą ujemną łącznością z ROI, zobacz Rysunek 5.

Rysunek 5: Wyświetlanie objętości nasion po zakrętie obręczy. Obraz wykazuje większą dodatnią łączność z przednim zakrętem obręczy i (czerwony) oraz większą negatywną łączność z przedklinkiem (niebieski) dla stanu True względem Sham, pFWEc < 0,05. Liczba ta przedstawia częściowe dane na poziomie grupy (n = 34). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Plik uzupełniający 1: Instrukcje dotyczące zamawiania urządzeń CES Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Dziennik urządzenia CES Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 3: fMRIPrep Boilerplate Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 4: Instrukcje CES CONN Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 5: Wykresy kodu CES R Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 6: R Code CES eddy-qc Anova Kliknij tutaj, aby pobrać ten plik.
Dodatkowe Rysunek 1: Skoroszyt CONN. (A) Współzmienne drugiego poziomu. (B) Warunki doświadczalne. (C) Odszumianie. (d) Wyniki. (E) Dane strukturalne. Kliknij tutaj, aby pobrać ten rysunek.