Badanie to zostało przejrzane i zatwierdzone przez Komisję Etyczną Shanghai Pudong New Area Beicai Community Health Service Center (numer zatwierdzenia 20250425AR; zatwierdzone 25 kwietnia 2025 r.). Wszyscy uczestnicy wyrazili pisemną świadomą zgodę przed włączeniem do badania.
Projekt badania i populacja badawcza
Celem tego badania przekrojowego była identyfikacja obiektywnych wskaźników związanych z konstytucjami wilgotnymi (flegmowo-wilgotną i wilgotno-gorączkową) przy użyciu BIA. Uczestników zrekrutowano w ramach społecznego programu przesiewowego zdrowia przeprowadzonego w Beicai Community Health Service Center (Pudong New Area, Szanghaj) w okresie od 1 lipca 2025 r. do 5 stycznia 2026 r. Do badania włączono kolejno wszystkie osoby, które w tym czasie przeszły zarówno ocenę konstytucji według TCM, jak i pomiar BIA. Wykluczono uczestników z brakującymi danymi w jakiejkolwiek kluczowej zmiennej. Ostatecznie włączono łącznie 232 uczestników, w tym 129 w grupie konstytucji harmonijnej i 103 w grupie konstytucji wilgotnej.
Kryteria diagnostyczne dla konstytucji wilgotnej
Ocenę konstytucji przeprowadzono przy użyciu opatentowanego przenośnego inteligentnego lustra do analizy konstytucji w tradycyjnej medycynie chińskiej (TCM) (Hi-Face22; chiński patent na wzór użytkowy nr ZL 2021 2 1286496.5), które wykorzystuje zestandaryzowany system klasyfikacji dziewięciu konstytucji TCM i generuje ilościową diagnozę konstytucjonalną na podstawie cech twarzy oraz odpowiedzi w zintegrowanym kwestionariuszu. Uczestników zidentyfikowanych jako posiadający przeważający typ konstytucji śluzowo-wilgotnej lub wilgotno-gorącej włączono do grupy z konstytucją wilgotną; osoby spełniające kryteria konstytucji harmonijnej i niewykazujące innych zaburzeń konstytucjonalnych włączono do grupy kontrolnej2. Z naszej wiedzy nie opublikowano recenzowanej walidacji tego urządzenia w odniesieniu do niezależnego standardu referencyjnego (np. oceny eksperckiej TCM lub zwalidowanego kwestionariusza), co zostało uznane za ograniczenie badania.
Pomiar BIA i składu ciała
Pomiary BIA przeprowadzono przy użyciu wieloczęstotliwościowego, ośmiopunktowego segmentowego analizatora impedancji bioelektrycznej z elektrodami dotykowymi (InBody 770, patrz Tabela materiałów), w ścisłej zgodności ze standardowymi procedurami operacyjnymi producenta. Urządzenie przykłada prądy przemienne o wielu częstotliwościach (1, 5, 50, 250, 500 i 1 000 kHz) i jest czułe na fazę, bezpośrednio mierząc oporność (R), reaktancję (Xc) oraz kąt fazowy (PhA) dla całego ciała i pięciu segmentów (prawe ramię, lewe ramię, tułów, prawa noga, lewa noga) za pomocą tetrapolarnych konfiguracji elektrod. Oporność, reaktancja i kąt fazowy tułowia były mierzone między elektrodami prawej dłoni i prawej stopy i rejestrowane przy częstotliwości 50 kHz. Należy zauważyć, że urządzenia BIA bezpośrednio mierzą parametry bioelektryczne, na które wpływ ma stan płynów ustrojowych i elektrolitów; raportowane objętości płynów i wartości składu ciała są szacunkami wyprowadzonymi z równań predykcyjnych urządzenia. Przed pomiarem badani byli zobowiązani do bycia na czczo przez co najmniej 2 h, unikania intensywnego wysiłku fizycznego przez 12 h oraz opróżnienia pęcherza moczowego. Podczas pomiaru badani stali boso na elektrodach stóp i trzymali elektrody dłoni obiema rękami. Zarejestrowano następujące parametry składu ciała: całkowitą wodę ustrojową (TBW), wodę wewnątrzkomórkową (ICW), wodę pozakomórkową (ECW), stosunek ECW/TBW (wartość liczbowa i procentowa), stosunek TBW do beztłuszczowej masy ciała (TBW/FFM), procent tkanki tłuszczowej (PBF), masę tkanki tłuszczowej (BFM), powierzchnię tłuszczu trzewnego (VFA), wskaźnik masy tłuszczowej (FMI), stopień otyłości, masę mięśni szkieletowych (SMM), beztłuszczową masę ciała (SLM), masę beztłuszczową (FFM), masę komórek ciała (BCM), białko, obwód ramienia (AMC) oraz podstawową przemianę materii (BMR). Zarejestrowano następujące segmentowe parametry bioelektryczne: kąt fazowy całego ciała (PA_Whole), kąt fazowy prawej kończyny górnej (PA_RA), kąt fazowy lewej kończyny górnej (PA_LA), kąt fazowy tułowia (PA_TR), kąt fazowy prawej kończyny dolnej (PA_RL), kąt fazowy lewej kończyny dolnej (PA_LL), impedancję tułowia (Impedance_TR) oraz reaktancję tułowia (Reactance_TR). Jednocześnie zmierzono następujące parametry antropometryczne: wzrost, wagę, wskaźnik masy ciała (BMI), stosunek obwodu talii do obwodu bioder (WHR), obwód brzucha oraz obwody lewej i prawej łydki.
Pomiar siły ścisku dłoni
Siłę uścisku dłoni mierzono za pomocą cyfrowego dynamometru ręcznego (patrz Tabela materiałów). Po wykonaniu dwóch maksymalnych skurczów izometrycznych każdą dłonią, naprzemiennie lewą i prawą, z zachowaniem jednej minuty przerwy między skurczami, obliczono średnią z dwóch prób dla dłoni dominującej, która została wykorzystana do analizy (Mean_Grip)
Analiza statystyczna
Ogólna struktura analityczna (przesiewanie cech w oparciu o metodę LASSO, a następnie wieloczynnikowa regresja logistyczna, wraz z analizami stratyfikowanymi ze względu na płeć oraz analizami wrażliwości) została zdefiniowana a priori. Jednakże konkretne cechy włączone do modelu końcowego zostały określone na podstawie selekcji LASSO sterowanej danymi, a analizy stratyfikowane według płci oraz analizy wrażliwości należy uznać za eksploracyjne. W przypadku żadnej z zmiennych uwzględnionych w analizie nie odnotowano brakujących wartości. Zmienne ciągłe o rozkładzie normalnym przedstawiono jako średnia ± odchylenie standardowe; zmienne o rozkładzie innym niż normalny przedstawiono jako mediana (rozstęp międzykwartylny); zmienne kategoryczne przedstawiono jako częstość (procent). Porównania między grupami przeprowadzono za pomocą testu t-Studenta dla prób niezależnych, testu U Manna-Whitneya lub testu chi-kwadrat, przyjmując, że P < 0,05 oznacza istotność statystyczną.
Następnie przesiewowo wyselekcjonowano cechy BIA za pomocą regresji LASSO; optymalny parametr kary λ określono poprzez 10-krotną walidację krzyżową, a wybrano cechy o niezerowych współczynnikach przy wartości λ, która minimalizowała błąd walidacji krzyżowej (λ.min). Zdolność dyskryminacyjną modelu oceniono za pomocą pola pod krzywą charakterystyki operacyjnej odbiornika (AUC). Procedura 10-krotnej walidacji krzyżowej została zastosowana wyłącznie do wyboru parametru kary LASSO λ i nie służyła jako walidacja końcowego modelu; wewnętrzną walidację końcowego modelu przeprowadzono przy użyciu korekty optymizmu metodą bootstrap (1 000 reprobowań) oraz 5-krotnej walidacji krzyżowej, a współliniowość oceniono za pomocą współczynnika inflacji wariancji (VIF), przyjmując, że wartości poniżej pięciu są akceptowalne. Kalibrację modelu oceniono za pomocą testu dopasowania Hosmera–Lemeshow oraz skorygowanego metodą bootstrap nachylenia kalibracji. Korzyść kliniczną netto oceniono za pomocą analizy krzywych decyzyjnych (DCA). Analizę DCA przeprowadzono dla pełnego zakresu prawdopodobieństwa progowego od 0 do 1,0, a korzyść netto wyrażono w stosunku do strategii referencyjnych „lecz wszystkich” oraz „nie lecz nikogo”. Scenariusz decyzji klinicznej dotyczył tego, czy zalecić dalszą ocenę lub interwencję w przypadku konstytucji wilgotnej (dampness constitution) wśród uczestników przesiewowych z populacji ogólnej.
Następnie cechy zidentyfikowane w ramach przesiewu LASSO, wraz z wiekiem i płcią, włączono do wieloczynnikowego modelu regresji logistycznej, aby ocenić niezależny związek każdego wskaźnika z konstytucją wilgotną, co pozwoliło na wyznaczenie ilorazów szans (OR) oraz 95% przedziałów ufności (CI).
Analizy w podgrupach przeprowadzono z podziałem na płeć, aby ocenić stabilność powiązań. Analizy wrażliwości obejmowały: (1) zastosowanie winsoryzacji na poziomie 1% i 99% dla zmiennych ciągłych oraz (2) dodatkową korektę o stosunek obwodu talii do bioder lub obwód brzucha. Wszystkie analizy statystyczne wykonano przy użyciu programu R w wersji 4.5.2.