基于R軟件單組率Meta分析研究9種中醫(yī)體質(zhì)的人群分布

姚實林

(安徽中醫(yī)藥大學(xué)科技處，安徽 合肥 230012)

王琦提出9種基本的中醫(yī)體質(zhì)分類，并設(shè)計了含有60個條目的中醫(yī)體質(zhì)量表，制定了《中醫(yī)體質(zhì)分類與判定標準》。目前，已有很多學(xué)者采用9種中醫(yī)體質(zhì)分類量表對一般人群和疾病人群進行中醫(yī)體質(zhì)辨識，但不同的研究得出的9種體質(zhì)類型構(gòu)成比存在一定差異。本研究擬采用R軟件單組率Meta分析對一般人群的體質(zhì)類型構(gòu)成比進行合并，為中國人群的中醫(yī)體質(zhì)類型分布規(guī)律提供大樣本流行病學(xué)依據(jù)。

1 資料與方法

1.1 納入標準 ①采用《中醫(yī)體質(zhì)分類與判定標準》對中醫(yī)體質(zhì)進行分類；②受試人群為一般人群或體檢人群；③年齡>15歲。

1.2 排除標準 ①采用同一次調(diào)查的全部數(shù)據(jù)或部分數(shù)據(jù)發(fā)表的論文；②體質(zhì)類型不足9種；③有一部分受試者未被劃分為某種體質(zhì)類型，如以“傾向性體質(zhì)”代之；④數(shù)據(jù)統(tǒng)計存在錯誤，如9種體質(zhì)例數(shù)之和與樣本量存在矛盾；⑤調(diào)查對象為疾病人群；⑥樣本量小于500；⑦學(xué)位論文和會議論文。

1.3 檢索策略 在中國知網(wǎng)、萬方數(shù)據(jù)知識服務(wù)平臺、維普網(wǎng)、PubMed數(shù)據(jù)庫、Cochrane數(shù)據(jù)庫中檢索中醫(yī)體質(zhì)類型流行學(xué)調(diào)查的橫斷面研究文獻，檢索起止時間為2004年1月至2020年12月。中文文獻檢索策略：題名=體質(zhì) AND 摘要=(流調(diào) OR 流行病學(xué) OR 流行病學(xué)調(diào)查 OR 調(diào)查 OR 問卷 OR 量表)AND 全文=((一般人群 OR 自然人群 OR 體檢人群)AND 中醫(yī))；英文文獻檢索策略：constitution[Title] AND (epidemiological[Abstract] OR Survey[Abstract] OR investigation[Abstract] OR research[Abstract] OR epidemiological survey[Abstract] OR scale[Abstract] OR questionnaire[Abstract])AND population[All Fields] AND (Chinese medicine[All Fields] OR TCM[All Fields])。

1.4 資料的提取 認真閱讀所納入的文獻，采用Excel軟件對文獻的基本特征進行提取，基本字段包括著者及出版年、每種體質(zhì)的例數(shù)和總樣本量、可能的異質(zhì)性來源，對每種體質(zhì)類型建立一個表格，從每篇文獻提取的基本信息均被核對3次。

1.5 納入文獻的方法學(xué)質(zhì)量評價 采用美國衛(wèi)生保健研究和質(zhì)量機構(gòu)(Agency for Healthcare Research and Quality，AHRQ)推薦的橫斷面研究質(zhì)量評價標準對納入的文獻質(zhì)量進行評價，該標準包括11個條目，從資料來源、納入排除標準、研究時間段、研究是否連續(xù)、評價者的主觀因素、質(zhì)量控制、排除分析受試者的原因、混雜因素控制措施、缺失數(shù)據(jù)處理方法、量表應(yīng)答率及數(shù)據(jù)完整性、隨訪結(jié)果共11個方面對研究質(zhì)量進行評價，總分為11分。

1.6 Meta分析方法 參照羅美玲等方法，采用R軟件(3.6.3版)進行單組率Meta分析。單組率Meta分析的異質(zhì)性大于兩組率Meta分析的異質(zhì)性，故采用隨機效應(yīng)模型進行Meta分析。對原始率(某類型體質(zhì)例數(shù)與總例數(shù)比值)采用對數(shù)轉(zhuǎn)換、Logit轉(zhuǎn)換、反正弦轉(zhuǎn)換和雙反正弦轉(zhuǎn)換，采用Shapiro-Wilk檢驗對原始率和4種轉(zhuǎn)換率進行正態(tài)性檢驗，選擇

P

值最大的率進行Meta分析。將“男女性別比”、“地區(qū)”(研究所在地區(qū))、“是否兼夾”(是否報告兼夾體質(zhì))、“人群”(是體檢人群還是一般人群)作為影響中醫(yī)體質(zhì)類型分布的潛在因素，先進行Meta回歸分析，若某因素對回歸模型有統(tǒng)計學(xué)意義(

P

<0

.

05)，則對該因素進行亞組分析。

2 結(jié)果

2.1 文獻檢索結(jié)果 在中國知網(wǎng)檢索到文獻718篇，閱讀文題后選取62篇，閱讀全文后剔除37篇，剩余25篇；在萬方數(shù)據(jù)知識服務(wù)平臺檢索到文獻73篇，閱讀文題后剩余23篇，其中14篇與中國知網(wǎng)重復(fù)，閱讀全文后排除7篇，剩余2篇；在維普網(wǎng)中檢索到文獻13篇，閱讀文題后剩余10篇，均與知網(wǎng)或萬方數(shù)據(jù)知識服務(wù)平臺重復(fù)；在Pubmed數(shù)據(jù)庫中檢索到文獻16篇，在Cochrane數(shù)據(jù)庫中檢索到文獻8篇，閱讀文題后均予以排除。在剩下的25篇文獻中，其中6篇文獻因采用同一次流行病學(xué)調(diào)查的部分數(shù)據(jù)而被排除，最終獲得19篇文獻。

2.2 納入文獻的基本特征 19篇文獻的總樣本量為175 331。男女性別比為0.654～2.567(25%、50%、75%位數(shù)分別為0.697、0.878、1.011)。按白明華等方法對研究所在區(qū)域進行劃分，華東地區(qū)4篇，華南地區(qū)7篇，華北地區(qū)2篇(其中1篇文獻包括華南和華北地區(qū)的調(diào)查數(shù)據(jù))，西南地區(qū)1篇，西北地區(qū)3篇，地區(qū)均衡分布(包括華南、華北、華東、西北、西南地區(qū)的調(diào)查數(shù)據(jù))3篇。有2篇文獻報告了兼夾體質(zhì)，17篇文獻未報告兼夾體質(zhì)，10篇文獻以體檢人群為調(diào)查對象，9篇文獻以一般人群為調(diào)查對象。AHRQ評分為5～8分。見表1。

表1 納入文獻的基本特征

2.3 平和質(zhì)構(gòu)成比的Meta分析結(jié)果 分別對“男女性別比”“地區(qū)”“人群”“是否兼夾”進行Meta回歸分析，結(jié)果顯示“地區(qū)”為平和質(zhì)構(gòu)成比的異質(zhì)性來源，見表2。敏感性分析無法降低

I

值。按地區(qū)進行亞組分析，結(jié)果顯示，華北地區(qū)平和質(zhì)構(gòu)成比(0.59，95% CI為0.48～0.69)最高，西南地區(qū)平和質(zhì)構(gòu)成比(0.51，95% CI為0.48～0.54)次之，華南地區(qū)平和質(zhì)構(gòu)成比(0.18，95% CI為0.12～0.25)最低，合并的平和質(zhì)構(gòu)成比為0.29(95% CI為0.23～0.36)。Meta分析結(jié)果見圖1。

表2 人群所在地區(qū)對平和質(zhì)構(gòu)成比影響的Meta回歸分析結(jié)果

圖1 平和質(zhì)構(gòu)成比的Meta分析結(jié)果

2.4 發(fā)表性偏倚分析 對漏斗圖的對稱性進行線性回歸分析，結(jié)果顯示漏斗圖對稱(

t

=-0

.

350，

ν

=18，

P

=0

.

731)，提示納入的文獻無發(fā)表性偏倚。見圖2。

圖2 平和質(zhì)Meta分析漏斗圖對稱性的Egger回歸分析

2.5 8種偏頗體質(zhì)的Meta分析結(jié)果 結(jié)果顯示，氣虛質(zhì)、陽虛質(zhì)、陰虛質(zhì)、痰濕質(zhì)、濕熱質(zhì)、血瘀質(zhì)、氣郁質(zhì)、特稟質(zhì)的構(gòu)成比分別為0.15(95% CI：0.13～0.18)、0.15(95% CI：0.13～0.18)、0.07(95% CI：0.05～0.10)、0.09(95% CI：0.07～0.12)、0.08(95% CI：0.06～0.10)、0.04(95% CI：0.03～0.05)、0.06(95% CI：0.05～0.08)、0.03(95% CI：0.02～0.04)，“男女性別比”是陰虛質(zhì)構(gòu)成比的影響因素，“是否兼夾”是陰虛質(zhì)、痰濕質(zhì)、濕熱質(zhì)、血瘀質(zhì)構(gòu)成比的影響因素，“地區(qū)”是血瘀質(zhì)、氣郁質(zhì)、特稟質(zhì)構(gòu)成比的影響因素。見表3。

表3 8種偏頗體質(zhì)構(gòu)成比的Meta分析結(jié)果(隨機效應(yīng)模型)

3 討論

中醫(yī)體質(zhì)是采用中醫(yī)陰陽氣血津液、病因病機等理論對人體形態(tài)結(jié)構(gòu)、生理功能和心理特征等方面相對穩(wěn)定特質(zhì)的中醫(yī)學(xué)描述。中醫(yī)體質(zhì)是可以分類的，在2萬例中醫(yī)體質(zhì)類型流行病學(xué)調(diào)查的基礎(chǔ)上，形成了《中醫(yī)體質(zhì)分類與判定》標準，基于該標準的9種中醫(yī)體質(zhì)量表被廣泛應(yīng)用于一般人群和疾病人群的中醫(yī)體質(zhì)類型判定。由于各研究的中醫(yī)體質(zhì)類型構(gòu)成比存在差異，有必要采用Meta分析方法對各個研究數(shù)據(jù)進行合并，為中醫(yī)體質(zhì)類型分布全國常模的建立提供參考。

體質(zhì)的形成與自然環(huán)境因素相關(guān)，《素問·異法方宜論》論及東方之人、西方之人、北方之人、南方之人、中央之人的體質(zhì)差異。不同性別者的體質(zhì)類型存在一定差異：女性為陰柔之體，陰盛陽衰；男性為陽剛之體，氣多血少，陰弱陽旺。兼夾體質(zhì)在人群中廣泛存在，常表現(xiàn)為以一種體質(zhì)為主，兼夾其他一種或數(shù)種體質(zhì)。兼夾體質(zhì)常被歸類為其主要體質(zhì)類型，即按照偏頗體質(zhì)評分最大值法劃分體質(zhì)類型。筆者在提取某體質(zhì)構(gòu)成比時，包含該體質(zhì)在兼夾體質(zhì)的構(gòu)成比和該單一體質(zhì)的構(gòu)成比，因此報告兼夾體質(zhì)的研究，其偏頗體質(zhì)的構(gòu)成比應(yīng)高于未報告兼夾體質(zhì)的研究。所納入的研究中，既有以一般人群為研究對象者，也有以體檢人群為研究對象者，故將人群作為影響體質(zhì)類型構(gòu)成比的潛在影響因素。在進行Meta分析時，并未對“男女性別比”“是否兼夾”“地區(qū)”“人群”均進行亞組分析，而是先對這4種潛在的影響因素進行Meta回歸分析，僅對進入回歸模型的影響因素進行亞組分析，以減少亞組分析的盲目性。結(jié)果只有“人群”對9種中醫(yī)體質(zhì)類型構(gòu)成比的影響均無統(tǒng)計學(xué)意義(

P

>0

.

05)，而“男女性別比”“是否兼夾”“地區(qū)”僅是部分中醫(yī)體質(zhì)類型構(gòu)成比的影響因素。結(jié)果既驗證了環(huán)境因素是體質(zhì)的后天成因和男女體質(zhì)有別的中醫(yī)體質(zhì)理論，又提示有些研究在流行病學(xué)調(diào)查的質(zhì)量控制和中醫(yī)體質(zhì)判定標準的掌握程度方面尚需提高。合并后的9種中醫(yī)體質(zhì)構(gòu)成比之和為0.96，非常接近1，說明Meta分析結(jié)果具有較好的可信度。

續(xù)表3

本研究尚存在一些不足。其一，大多數(shù)研究的方法學(xué)質(zhì)量不高。其二，包含兼夾體質(zhì)的研究較少，僅有2篇文獻報告了兼夾體質(zhì)。Meta回歸分析結(jié)果中，“是否兼夾”應(yīng)是8種偏頗體質(zhì)的影響因素，但“是否兼夾”僅是部分偏頗體質(zhì)的影響因素，其原因可能與報告兼夾體質(zhì)的研究過少(僅2篇文獻)及兼夾體質(zhì)過于復(fù)雜而未被充分報告有關(guān)。其三，對人群的年齡分布情況未按固定的格式進行介紹，有僅介紹最小值者，有僅介紹平均值者；大部分研究中受試人群的年齡分布范圍過寬，而白明華等采用老年版中醫(yī)體質(zhì)量表對65歲以上人群進行體質(zhì)辨識，筆者僅提取其研究中15～64歲人群的數(shù)據(jù)。因此，無法按照年齡段進行Meta回歸和亞組分析。

雖然存在以上不足，但本研究得出的9種中醫(yī)體質(zhì)的總體構(gòu)成比對全國常模的建立具有借鑒意義。期待更多一般人群中醫(yī)體質(zhì)流行病學(xué)調(diào)查結(jié)果對今后Meta分析結(jié)果進行修正和完善。