基于縱向研究數(shù)據(jù)的非線性混合效應(yīng)模型在中醫(yī)證候演變規(guī)律研究中的應(yīng)用

伊明洋，關(guān)建新，閆麗禎，倪園園，王至婉

1.河南中醫(yī)藥大學(xué)，河南 鄭州 450046；2.河南中醫(yī)藥大學(xué)呼吸疾病中醫(yī)藥防治省部共建協(xié)同創(chuàng)新中心，河南 鄭州 450046；3.河南省中醫(yī)藥防治呼吸病重點實驗室，河南 鄭州 450046；4.河南中醫(yī)藥大學(xué)第一附屬醫(yī)院，河南 鄭州 450000

證候是中醫(yī)理論體系的核心內(nèi)容，證候演變規(guī)律研究是中醫(yī)證候研究的重要組成部分[1]，病證結(jié)合深入研究證候演變規(guī)律對構(gòu)建新的辨證論治體系并提高辨證論治水平具有重要意義[2]?？v向研究作為臨床試驗研究設(shè)計方法之一，具有較強的時間序列性及因果邏輯關(guān)系，被廣泛應(yīng)用于證候演變規(guī)律研究領(lǐng)域。然而個體化試驗設(shè)計下產(chǎn)生的大量非線性重復(fù)測量數(shù)據(jù)往往無法得到合理有效的處理[3]，基于此，本文提出使用非線性混合效應(yīng)模型處理證候演變規(guī)律縱向研究數(shù)據(jù)，以期為解決該領(lǐng)域方法學(xué)難題提供思路。

1 證候演變規(guī)律縱向研究數(shù)據(jù)特征

1.1 縱向研究特征及優(yōu)勢

縱向研究是一種前瞻性研究，因在前、果在后，易得出基于邏輯的因果判斷，在證候演變規(guī)律研究中選擇疾病發(fā)展過程中的關(guān)鍵時點作為觀察節(jié)點，采用多時點重復(fù)測量或隊列研究等方式對證候變化進行縱向研究，可描述事物的連續(xù)性變化過程，研究結(jié)果更具可靠性和科學(xué)性，是證候演變規(guī)律研究的最佳方法[4]。其優(yōu)勢在于：①可反映樣本個體間的差異，觀察影響不同個體間差異的因素，如地域、年齡、職業(yè)等。②可系統(tǒng)觀察人從未發(fā)病到發(fā)病及好轉(zhuǎn)乃至痊愈的全過程，更加貼近疾病發(fā)生發(fā)展的真實過程，研究結(jié)果可靠性更強。③較橫斷面研究，縱向研究更有效能，即若達到相同的統(tǒng)計效能，縱向研究所需的樣本量往往少于橫斷面研究，究其原因是縱向研究中同一受試者的不同重復(fù)測量值并不完全相關(guān)[5]。④在縱向研究中，每個受試者都能與自身對照，一般來說，內(nèi)部變異要小于個體間的變異，所得到的研究結(jié)果會更加敏感[6]。張明雪等[7]通過縱向研究分析冠心病合并高血壓患者4個階段證候變化，0～1期向肝腎陰虛、氣火失調(diào)證轉(zhuǎn)變，1～2期逐漸轉(zhuǎn)變?yōu)樘禑釘_心、絡(luò)虛陽升證，2～3期發(fā)展為土虛木亢證及肝陽化風(fēng)證。李偉珂等[8]將慢性阻塞性肺疾病分為急性加重期、危險窗期和穩(wěn)定期，縱向研究證候要素演變規(guī)律發(fā)現(xiàn)，從急性加重期到穩(wěn)定期總體呈現(xiàn)由實到虛的趨勢，痰熱逐漸演變?yōu)樘禎裰敝镣耆?，氣虛貫穿疾病始終，在發(fā)展過程中逐漸變?yōu)闅怅巸商?，兼有血瘀?/p>

1.2 縱向研究數(shù)據(jù)特點

證候演變規(guī)律研究中所采集的縱向數(shù)據(jù)有以下特點：①時間序列性。主要研究隨時間變化證候的演變情況，因此臨床所采集數(shù)據(jù)會隨時間的推移發(fā)生變化。②時間間隔的不平衡性。在臨床研究中，各個指標采集時間可能不盡相同，每個指標采集的間隔時間也會不同。③相關(guān)性。在同一觀測時點觀測指標內(nèi)部存在相關(guān)性，不同觀測時點獲得的觀測指標間也存在相關(guān)關(guān)系。④復(fù)雜多態(tài)性。臨床所采集數(shù)據(jù)既有分類數(shù)據(jù)也有連續(xù)數(shù)據(jù)，反應(yīng)變量類型既有連續(xù)型也有離散型。⑤非線性?；谧C候本身的復(fù)雜性，臨床采集數(shù)據(jù)重要特點即自變量與反應(yīng)變量呈非線性變化趨勢。⑥數(shù)據(jù)缺失性。縱向研究是基于長時間隨訪調(diào)查的研究方法，在長期隨訪過程中無法避免數(shù)據(jù)的缺失。⑦非正態(tài)性?？v向研究數(shù)據(jù)變化的影響因素較多、變量復(fù)雜多樣，數(shù)據(jù)難以形成嚴格意義的正態(tài)分布[9]。

1.3 縱向研究數(shù)據(jù)常用分析方法

目前在證候演變規(guī)律研究領(lǐng)域，針對縱向研究的數(shù)據(jù)處理方法主要有頻率分析、轉(zhuǎn)移概率矩陣和非線性混合效應(yīng)模型等。

1.3.1 頻率分析

頻率分析是觀察某種證型在某個時點或階段所占比例，通過數(shù)量值的變化闡釋證候在不同時間的變化規(guī)律。郭思佳等[10]使用頻率分析方法觀察特發(fā)性肺纖維化不同臨床分期的證素演變規(guī)律，發(fā)現(xiàn)氣虛貫穿疾病全過程，其中臨床1期和2期以實證為主，主要證素為痰濁和血瘀，隨著疾病進展，痰濁變化不明顯，而血瘀、陰虛等逐漸加重，總體呈現(xiàn)由實轉(zhuǎn)虛的趨勢。頻率分析法操作簡單，對數(shù)據(jù)類型要求不高，在臨床研究中應(yīng)用廣泛。但僅使用頻率分析法無法深入挖掘臨床數(shù)據(jù)蘊含的豐富信息，較為單一和片面，因此多將頻率分析法與轉(zhuǎn)移概率矩陣和非線性混合效應(yīng)模型等數(shù)據(jù)挖掘技術(shù)結(jié)合使用。

1.3.2 轉(zhuǎn)移概率矩陣

轉(zhuǎn)移概率矩陣在20世紀初由俄國數(shù)學(xué)家馬爾可夫提出，用于研究在同一系統(tǒng)中一個因素向另一個因素轉(zhuǎn)移的概率，在轉(zhuǎn)移中第n次結(jié)果只受第n-1次結(jié)果的影響，具有無后效性[11]。在證候演變規(guī)律研究中多用于研究證型由當(dāng)前時點向下一時點轉(zhuǎn)變的概率，可清晰描述各證候隨時間變化關(guān)系。傅凱麗等[12]使用轉(zhuǎn)移概率矩陣研究發(fā)現(xiàn)，后循環(huán)梗死患者初期以風(fēng)痰瘀阻證候為主，在疾病發(fā)展中期逐步轉(zhuǎn)變?yōu)闅馓撗?、痰濁?nèi)阻等證，后期以氣血虧虛、腎精不足為主，即初期以實證為主，中期表現(xiàn)為虛實夾雜，后期多虛證。使用該方法需同時滿足齊次馬爾可夫性假設(shè)和觀測獨立性假設(shè)，同時該方法無法研究癥狀與證候間的映射關(guān)系，在無公認的證候診斷標準前提下，大多數(shù)研究中的證候由研究者根據(jù)既往診斷標準判定，具有一定主觀性，導(dǎo)致證候演變規(guī)律研究結(jié)果也具有一定偏倚。

1.3.3 非線性混合效應(yīng)模型

非線性混合效應(yīng)模型于1977年由Sheiner等[13]提出，是線性混合效應(yīng)模型到非線性模型的推廣。非線性混合效應(yīng)模型可同時適用于連續(xù)型和離散型變量，其變量無須滿足正態(tài)分布，主要用于處理定量非線性資料，分析非線性增長趨勢差異，可很好地解決研究中出現(xiàn)缺失值和不平衡數(shù)據(jù)等問題[9]。目前該研究方法在證候演變規(guī)律研究領(lǐng)域尚未得到推廣，主要應(yīng)用于藥代動力學(xué)、林業(yè)、保險金融業(yè)等。

靳曉娟[14]使用非線性混合效應(yīng)模型研究不同樹齡的人工長白落葉松枝條生長屬性的分布規(guī)律，以樹齡和單木為兩水平構(gòu)建非線性混合效應(yīng)模型，可對枝條基徑、長度等進行合理預(yù)測，更符合樹木的生理生長特性。李春明[15]使用非線性混合效應(yīng)模型處理具有時間序列相關(guān)性的縱向數(shù)據(jù)具有較高的估計精度，且能夠較好地表現(xiàn)連續(xù)觀測數(shù)據(jù)間的誤差分布。葉楠等[16]對患者使用不同劑量的腎康注射液進行療效評價，并使用非線性混合效應(yīng)模型的擬合功能，進行上市后市場監(jiān)測。王明高等[17]使用非線性混合效應(yīng)模型處理保險行業(yè)中非壽險損失預(yù)測和費率厘定問題，解決該行業(yè)實際操作中保險損失數(shù)據(jù)的問題，對該行業(yè)的發(fā)展具有重要意義。不論是樹木生長、藥物的療效評價及保險的損失預(yù)測及費率問題，在研究過程中所采集到的數(shù)據(jù)均為重復(fù)測量數(shù)據(jù)且具有時間的連續(xù)性，這與證候演變縱向研究中獲得的數(shù)據(jù)特性相同，目前該方法雖未在證候演變規(guī)律研究領(lǐng)域推廣使用，但在其他領(lǐng)域內(nèi)應(yīng)用較為成熟，具有穩(wěn)定性，可進一步應(yīng)用于中醫(yī)證候研究領(lǐng)域。

非線性混合效應(yīng)模型在醫(yī)學(xué)領(lǐng)域中的應(yīng)用主要集中于藥代動力學(xué)及其他非線性生長曲線研究。在藥代動力學(xué)研究中，藥物在個體體內(nèi)吸收、分布及排泄的過程中采集到的血藥濃度數(shù)據(jù)與證候演變規(guī)律縱向研究中獲得的重復(fù)測量數(shù)據(jù)特征相似，都表現(xiàn)為非線性特征。如凌靜等[18]使用非線性混合效應(yīng)模型研究他克莫司在成人腎病綜合征中的群體藥代學(xué)特征，建模中采用可體現(xiàn)個體間變異和殘差變異的一級速率條件算法，同時將患者體質(zhì)、年齡及藥物劑量等作為協(xié)變量觀察是否對個體血藥濃度產(chǎn)生影響。田婷婷等[19]應(yīng)用非線性混合效應(yīng)模型構(gòu)建替考拉寧在老年患者中的群體藥動學(xué)模型，研究結(jié)果顯示該模型擬合優(yōu)度好，具有良好的穩(wěn)定性及預(yù)測性。陶玉瑛[20]基于文獻研究結(jié)果及醫(yī)院癲癇兒童服用托吡酯后體內(nèi)血藥濃度變化數(shù)據(jù)構(gòu)建非線性混合效應(yīng)模型并對其進行驗證發(fā)現(xiàn)，該模型具有穩(wěn)定的預(yù)測性能，可根據(jù)預(yù)測結(jié)果制定臨床個體化給藥方案。闕燁等[21]基于艾滋病的縱向研究數(shù)據(jù)構(gòu)建部分非線性混合效應(yīng)模型，將log10RNA作為響應(yīng)變量，CD4細胞數(shù)作為協(xié)變量，觀察不同時間采集的CD4細胞數(shù)對log10RNA的影響。以上研究結(jié)果均表明非線性混合效應(yīng)模型在醫(yī)療領(lǐng)域具有廣泛的應(yīng)用前景，結(jié)合中醫(yī)證候演變規(guī)律研究領(lǐng)域，可解決目前非線性重復(fù)測量數(shù)據(jù)處理難題，可更好地模擬證候信息隨時間變化過程，同時也可通過設(shè)置協(xié)變量了解影響證候演變相關(guān)因素。

2 非線性混合效應(yīng)模型方法學(xué)概述

2.1 模型構(gòu)建

非線性混合效應(yīng)模型在使用過程中為適應(yīng)不同需求，可轉(zhuǎn)化為相應(yīng)的表達式，其基本表達式為：

式中，i表示研究對象個數(shù)，j表示重復(fù)觀察次數(shù)，yij是第i個觀測對象第j次重復(fù)測量的因變量值，xij表示自變量值，eij表示第i個研究對象第j次觀察的隨機誤差值，f(·)為非線性函數(shù)關(guān)系，?ij表示該函數(shù)關(guān)系中的參數(shù)向量。β和bi分別為固定效應(yīng)和隨機效應(yīng)參數(shù)，Aij、Bij為設(shè)計矩陣[22]。

Pinheiro等[23]依據(jù)產(chǎn)生隨機效應(yīng)因子的個數(shù)將模型區(qū)分為單水平和嵌套多水平2種不同類型。證候演變規(guī)律研究中產(chǎn)生隨機效應(yīng)的因素較多，因此建議選擇嵌套多水平模型，具體模型的構(gòu)建應(yīng)根據(jù)所收集臨床資料要求建立。

2.2 參數(shù)估計

該模型的參數(shù)估計方法主要有廣義最小二乘估計法、一階線性化算法、條件一階線性化算法、高斯-埃爾米特求積法等[24]。這些算法各有利弊，其中廣義最小二乘估計法精度低，高斯-埃爾米特求積法計算精度高但運算量大且速度較慢，一階線性化算法和條件一階線性化算法計算方法簡單且精度高，在實際應(yīng)用中被廣泛使用[25]。目前，符利勇等[26]在條件一階線性化算法的基礎(chǔ)上提出一種改進的隨機效應(yīng)參數(shù)計算方法，并利用樹高數(shù)據(jù)模擬分析，發(fā)現(xiàn)這種改進的算法可獲得更好的擬合效果并可更好地反映個體間的差異，此算法可在SAS軟件中實現(xiàn)。

2.3 初始值選擇

非線性混合效應(yīng)模型是否能夠成功運行很大程度上取決于初始值的選擇，尤其當(dāng)擬合模型較為復(fù)雜時，隨意選擇初始值會造成迭代不收斂或得到的參數(shù)估計僅為局部最優(yōu)解而非全局最優(yōu)解，易影響參數(shù)估計的有效性[27]。該模型在開始迭代時，必須賦予所有參數(shù)初始值，初期可選擇較為簡單的模型，使用簡單模型的運算結(jié)果作為初始值，再通過數(shù)值積分或附加高斯積分求出最大似然的直接估計[27]。

3 實例分析

如何使用非線性混合效應(yīng)模型解決證候演變規(guī)律研究中的實際問題是眾多研究者較關(guān)心的問題，以下通過具體的研究范例分析該方法在實際研究中的應(yīng)用。

王麗穎[28]使用非線性混合效應(yīng)模型等多種數(shù)據(jù)分析方法研究不同藥物干預(yù)下冠心病心絞痛血瘀證的動態(tài)演變規(guī)律，探索隨時間變化藥物對證候演變過程的影響，明確證候改變的拐點，進一步評價中醫(yī)證候臨床療效。該研究首先使用2種藥物開展隨機、雙盲多中心臨床試驗，以證候積分量表作為證候評價工具；同時使用潛在類別模型對血瘀證進行藥效學(xué)解構(gòu)，尋找證候的“內(nèi)實”（最能表現(xiàn)某病機的關(guān)鍵內(nèi)容，即群體在某一病變過程中的共性規(guī)律，是選擇干預(yù)措施的主要依據(jù)）及潛在類別；最后使用非線性混合效應(yīng)模型結(jié)合向量相似度分析法探索證候整體及證候“內(nèi)實”的演變規(guī)律，以此作為證候療效評價的切入點。該研究以每位入組患者為不同的個體，以最能反映該證候的主要癥狀胸痛、胸悶、心悸、唇色紫黯、舌質(zhì)紫黯和脈細澀作為組成觀察向量（EO），以這些癥狀的最大積分值作為對照向量（EMAX），計算得到每位患者在不同時間的向量相似度值。將相似度值設(shè)為因變量，時間設(shè)為自變量（TIME），構(gòu)建療效（EFF）隨時間變化的非線性混合效應(yīng)模型：

EFF=EO-EMAX*(TIME**GAMM)/(ET50**GAMM+TIME**GAMM)

式中GAMM為時間指數(shù)，ET50為藥效下降1/2的時間。據(jù)用藥前的證候相似度值及相關(guān)文獻支持選擇相似度值0.8作為臨界點，代入非線性混合效應(yīng)模型中計算得出2組證候轉(zhuǎn)變時間分別為7.28周和10.72周，說明2種藥物對緩解冠心病血瘀證均有效，根據(jù)方程擬合結(jié)果，A藥療效略優(yōu)于B藥。同時將心絞痛積分、年齡、試驗中心和組別等作為協(xié)變量加入模型發(fā)現(xiàn)，心絞痛積分、試驗中心對2組證候變化均有影響。該研究使用非線性混合效應(yīng)模型研究血瘀證演變規(guī)律，既可發(fā)現(xiàn)證候轉(zhuǎn)變的具體時間，又能進一步探究影響證候演變的相關(guān)因素，表明該模型在證候演變規(guī)律研究中具有可行性，后續(xù)研究可進一步推廣使用。

4 討論

近年來，隨著證候演變規(guī)律研究的不斷深入，對數(shù)據(jù)處理的要求不斷提高。基于此，將非線性混合效應(yīng)模型應(yīng)用于該研究領(lǐng)域，其優(yōu)勢主要有以下2個方面：①證候本身是一個巨大的復(fù)雜非線性系統(tǒng)，其各個子系統(tǒng)不能呈現(xiàn)規(guī)律性成比例的線性運動[29]。傳統(tǒng)的計算方法多要求數(shù)據(jù)呈現(xiàn)線性規(guī)律，然而多時點采集的中醫(yī)四診信息難以滿足該要求。如針對慢性阻塞性肺疾病患者，不同時點的咳嗽積分難以呈現(xiàn)嚴格的線性規(guī)律，同時針對證候本身而言，在疾病的不同階段也難以形成單一的線性變化。因此，針對研究數(shù)據(jù)的非線性特征，非線性混合效應(yīng)模型可很好地擬合不同時間點的信息資料，滿足其非線性要求，更加科學(xué)準確地發(fā)掘證候動態(tài)演變規(guī)律。②證候演變規(guī)律研究基于復(fù)雜的中醫(yī)四診信息，多時點的重復(fù)測量資料既有單個反應(yīng)變量重復(fù)測量值間的相關(guān)性，也有多個反應(yīng)變量間的相關(guān)性。以慢性阻塞性肺疾病為例，急性加重期第1日咳嗽與咳痰積分間具有相關(guān)性，第7日的咳嗽積分與第1日也具有相關(guān)性，且由于具有不同的相關(guān)性，會產(chǎn)生個體間及個體內(nèi)的隨機誤差。針對以上情況，非線性混合效應(yīng)模型既可識別個體間差異也可兼顧個體內(nèi)差異，同時能計算不同層次的誤差估計值。該模型允許固定效應(yīng)和隨機效應(yīng)同時進入非線性部分，通過計算得到固定和隨機效應(yīng)的估計值及單個反應(yīng)變量多次重復(fù)測量的相關(guān)系數(shù)和多個反應(yīng)變量間的相關(guān)系數(shù)，使數(shù)據(jù)間的相關(guān)性更加細化[23]。

除上述優(yōu)點外，非線性混合效應(yīng)模型在證候演變規(guī)律研究中的應(yīng)用還有以下方面值得關(guān)注：①臨床研究不可避免出現(xiàn)數(shù)據(jù)缺失，且無法保證各觀測時間點間隔完全均衡。非線性混合效應(yīng)模型不要求每個個體測量時間相同，同時允許缺失值的出現(xiàn)，也不要求測量間隔時間完全相等，可同時估算個體曲線和全部樣本曲線的平均值[20]。②證候演變過程極為復(fù)雜，易受生存環(huán)境、體質(zhì)、性別、心情等諸多因素的影響。非線性混合效應(yīng)模型在計算過程中可加入解釋變量，如地域、年齡、疾病階段、輕重程度等，分析這些因素對證候演變過程的影響。③使用非線性混合效應(yīng)模型可定量把握證候演變的轉(zhuǎn)折點，為臨床辨證論治提供較為客觀依據(jù)。

非線性混合效應(yīng)模型在證候演變規(guī)律研究領(lǐng)域尚處于探索階段，基于數(shù)理統(tǒng)計方面的復(fù)雜性，該模型在具體應(yīng)用中仍存在亟待解決的問題：①方程模型的構(gòu)建。在進行模型擬合前可先根據(jù)原始數(shù)據(jù)做出散點圖，再依據(jù)相關(guān)專業(yè)知識、研究經(jīng)驗及實測點的分布情況構(gòu)建符合數(shù)據(jù)要求的方程模型[30]。②初始值的選擇?？蓞⒖嫉慕鉀Q方法是，同時選擇幾個初始值代入運算，觀察得到結(jié)果是否一致，若不一致，則通過選取運算所得圖形中的點，計算各參數(shù)的近似值，再反向作為初始值帶入模型計算擬合，觀察該初始值是否為最優(yōu)初始值[23]。③收斂性評價和模型擬合優(yōu)度。模型的構(gòu)建及初始值的選擇是否滿足收斂條件和模型的擬合優(yōu)度對非線性分析過程至關(guān)重要。

將非線性混合效應(yīng)模型應(yīng)用于證候演變規(guī)律縱向研究數(shù)據(jù)處理有利于實現(xiàn)對證候動態(tài)變化時間的把握，也可明確影響證候演變的相關(guān)因素，對指導(dǎo)臨床辨證論治具有重要意義。但該模型的構(gòu)建較為復(fù)雜，在實際應(yīng)用過程中存在較多方法學(xué)難題。在實際應(yīng)用中除解決上述困難外，還應(yīng)注重證候診斷標準及術(shù)語的規(guī)范，可通過建立電子病歷解決數(shù)據(jù)完整性等問題，為后續(xù)研究提供堅實基礎(chǔ)[31]。