• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    體檢指標(biāo)健康預(yù)警的灰色-時(shí)序組合模型

    2023-07-27 22:33:44朱人杰葉春明
    關(guān)鍵詞:灰色

    朱人杰 葉春明

    摘要:對于個(gè)體健康體檢數(shù)據(jù)而言,傳統(tǒng)的以大樣本為基礎(chǔ)的數(shù)學(xué)模型無法滿足體檢數(shù)據(jù)的建模 需求?;趥€(gè)體體檢數(shù)據(jù)特征分析,首先構(gòu)建適用于個(gè)體體檢指標(biāo)健康預(yù)警的近似非齊次指數(shù)序 列的改進(jìn)離散灰色模型。其次,為降低單個(gè)模型預(yù)測精度的有限性,利用方差倒數(shù)法為離散灰色 模型和差分自回歸移動(dòng)平均模型賦權(quán)重,在模型誤差平方和達(dá)到最小時(shí)取得最佳的權(quán)重值。從而 將兩個(gè)模型的預(yù)測結(jié)果進(jìn)行組合, 實(shí)現(xiàn)對健康指標(biāo)的建模與趨勢分析,及時(shí)掌握個(gè)體健康指標(biāo)的 變化并發(fā)現(xiàn)潛在的疾病隱患。預(yù)測模型在實(shí)驗(yàn)數(shù)據(jù)集上的相對模擬誤差與最優(yōu)基準(zhǔn)模型相比有所 下降,表明灰色– 時(shí)序組合模型具有更高的模擬精度,解決了傳統(tǒng)的依據(jù)單次體檢指標(biāo)進(jìn)行靜態(tài)分 析的弊端以及單個(gè)模型預(yù)測結(jié)果的局限性,更加關(guān)注個(gè)體差異,能有效提升健康預(yù)警的效果。

    關(guān)鍵詞:灰色– 時(shí)序組合模型 ;體檢指標(biāo) ; 離散灰色模型 ;差分自回歸移動(dòng)平均模型 ;健康預(yù)警

    中圖分類號(hào):? TP 931???????????? 文獻(xiàn)標(biāo)志碼:?? A

    Grey time series combination model for health warning of physical examination indexes

    ZHU Renjie1,2, YE Chunming1

    (1. Business School, University of Shanghai for Science and Technology, Shanghai 200093, China;2. Shanghai Easthospital Affiliated to Shanghai Tongji University, Shanghai 200120, China)

    Abstract: For individual health examination data, the traditional mathematical model based on large samples can not meet the modeling requirements of physical examination data. Based on the analysis of the characteristics of individual physical examination data, an improved discrete grey model of approximately non-homogeneous index series suitable for individual physical examination indicator health warning was first constructed. Secondly, in order to reduce the limitation of the prediction accuracy of a single model, the inverse variance method was used to assign weights to the discrete grey model and the differential autoregressive moving average model, and the best weight value was obtained when the sum of squares of the model errors reached the minimum. Thus, the prediction results of the two models were combined to achieve the modeling and trend analysis of health indicators, timely grasp the changes of individual health indicators and discover potential disease hazards. The relative simulation error of the prediction model on the experimental data set decreases in comparison with the optimal benchmark model, which indicates that the grey time series combination model has higher simulation? accuracy.? The?? shortcomings?? of? traditional?? static?? analysis?? based?? on ??single?? physical examination indicators? and the? limitations? of single model prediction results? are? solved.? Individual differences are emphasized, and the effect of health warning can be effectively improved.

    Keywords: grey time series combination model; physical examination index; discrete grey model; differential autoregressive integrated moving average model; health warning

    隨著時(shí)代的發(fā)展、社會(huì)需求和疾病譜的改變,以預(yù)防為主的大健康理念逐漸深入人心,民眾健康預(yù)防的意識(shí)也逐漸增強(qiáng),定期體檢演變?yōu)橐环N健康生活習(xí)慣和社會(huì)趨勢。健康體檢產(chǎn)生的個(gè)體體檢數(shù)據(jù),可以幫助醫(yī)生和體檢者分析個(gè)體自身生理狀況和潛在的疾病隱患。但是,醫(yī)生對體檢者身體狀況的診斷,或者體檢者對自身健康情況的判斷大多是基于單次體檢結(jié)果的高低對比,無法發(fā)現(xiàn)個(gè)體體檢指標(biāo)在不同時(shí)期的動(dòng)態(tài)變化趨勢。因此,分析個(gè)體體檢指標(biāo)的發(fā)展規(guī)律和變化趨勢,發(fā)現(xiàn)體檢者的潛在疾病隱患,從而提前采取預(yù)防和治療措施,降低個(gè)體未來的患病風(fēng)險(xiǎn),對于保障個(gè)體健康具有重大的現(xiàn)實(shí)意義。

    灰色模型對于“少數(shù)據(jù)”、“貧信息”的樣本具有較高的預(yù)測精度,能夠通過研究對象有限的數(shù)據(jù),挖掘出數(shù)據(jù)發(fā)展規(guī)律和新信息,從而實(shí)現(xiàn)對序列未來值的預(yù)測[1-2]。在疾病預(yù)測方面,灰色預(yù)測模型主要用于發(fā)病率、死亡率的預(yù)測[3-4]。而其他典型預(yù)測方法雖然在疾病預(yù)測中發(fā)揮了重要的作用,但是各類模型的適用范圍有所差異。時(shí)間序列模型通過將疾病數(shù)據(jù)隨時(shí)間推移形成的序列視為一個(gè)隨機(jī)序列,并用一定的數(shù)學(xué)模型來近似擬合這個(gè)序列,常用的時(shí)間序列模型為 ARIMA( autoregressive integrated moving average model)模型[5-6]?;诟怕收摰鸟R爾可夫鏈模型通常是基于系統(tǒng)現(xiàn)在的狀態(tài)來預(yù)測系統(tǒng)未來可能存在的狀態(tài),例如劉瓊等[2]利用隱馬爾科夫模型對乙肝發(fā)病數(shù)量時(shí)間序列進(jìn)行預(yù)測[7]。隨著神經(jīng)網(wǎng)絡(luò)的發(fā)展, BP 神經(jīng)網(wǎng)絡(luò)模型也被大量應(yīng)用于疾病預(yù)測中,并且在疾病預(yù)測中具有較好的識(shí)別效果[8-9]。多元回歸模型常用于傳染病發(fā)病率的趨勢預(yù)測,建模過程中應(yīng)用直線或曲線擬合原始傳染病數(shù)據(jù),用數(shù)字和等式來表達(dá)傳染病的流行規(guī)律[10-13]。近年來,國內(nèi)外學(xué)者將灰色模型與其他模型進(jìn)行組合,融合多個(gè)模型的優(yōu)勢,開展疾病預(yù)測研究。王永斌等[14]將灰色模型和廣義回歸神經(jīng)網(wǎng)絡(luò)模型相結(jié)合,預(yù)測我國塵肺病發(fā)病人數(shù)。嚴(yán)薇榮等[15]在進(jìn)行傷寒副傷寒發(fā)病率預(yù)測時(shí),將 GM(1,1)模型和 Markov模型進(jìn)行組合得到新的預(yù)測模型,提高了傳染病發(fā)病率的預(yù)測精度。時(shí)冬青等[16]綜合 GM(1,1)模型和馬爾可夫鏈進(jìn)行預(yù)測,實(shí)驗(yàn)結(jié)果表明組合模型在職業(yè)病預(yù)測中的高預(yù)測精度。

    目前,對于個(gè)體體檢指標(biāo)的研究主要集中于兩個(gè)方面:一方面是分析個(gè)體體檢指標(biāo)對于疾病診斷的影響或?qū)膊〉念A(yù)測價(jià)值[17-18];另一方面是在疾病風(fēng)險(xiǎn)預(yù)測中,將多個(gè)或群體健康體檢指標(biāo)作為預(yù)測特征來預(yù)測疾病發(fā)病率或患病情況[19-20]。然而,針對個(gè)體體檢指標(biāo)未來發(fā)展趨勢預(yù)測的研究還較少。通過上述分析可以發(fā)現(xiàn),以上研究大多采用群體健康指標(biāo)數(shù)據(jù)集開展疾病預(yù)測,而對于個(gè)體健康體檢指標(biāo)的預(yù)測較少,并且個(gè)體健康體檢數(shù)據(jù)的特征也增加了個(gè)體健康指標(biāo)預(yù)測的難度。為此,需要構(gòu)建有效的個(gè)體健康指標(biāo)預(yù)測模型,以期準(zhǔn)確預(yù)測體檢指標(biāo)未來變化趨勢或范圍,實(shí)現(xiàn)個(gè)體健康狀況的有效預(yù)警管理?;谏鲜龇治?,考慮體檢指標(biāo)數(shù)據(jù)為小樣本數(shù)據(jù),并且更偏向于是一個(gè)非齊次指數(shù)序列,為提高模型的泛化性和準(zhǔn)確性,本文構(gòu)建了一個(gè)離散灰色模型。同時(shí),為提高預(yù)測精度,將 ARIMA 模型和灰色模型進(jìn)行組合預(yù)測,從而充分利用各個(gè)模型的優(yōu)勢。

    1 個(gè)體健康體檢指標(biāo)特征分析

    隨著人們對于健康和自我保健追求的愈加強(qiáng)烈,健康消費(fèi)市場迅猛發(fā)展,個(gè)人定期健康體檢已成為常態(tài)。個(gè)人在醫(yī)療機(jī)構(gòu)進(jìn)行體檢,得到各類身體指標(biāo)檢查數(shù)據(jù)。這些具有時(shí)間間隔的數(shù)據(jù)匯總后形成了時(shí)間序列,對這些時(shí)間序列數(shù)據(jù)進(jìn)行數(shù)據(jù)分析和預(yù)測,可以有效地輔助醫(yī)生和患者了解當(dāng)前身體狀況和指標(biāo)的未來變化趨勢,幫助人們提前采取應(yīng)對策略,做好疾病預(yù)防。

    由于體檢指標(biāo)時(shí)間序列數(shù)據(jù)有其獨(dú)有的特征,在構(gòu)建時(shí)間序列預(yù)測模型時(shí)有必要基于其特征進(jìn)行設(shè)計(jì)。以單個(gè)體檢指標(biāo)m為例,指標(biāo)m在時(shí)間跨度1? n之間的檢查結(jié)果構(gòu)成一個(gè)時(shí)間序列 Xm =(xm(t1), xm(t2), ··· , xm(tn))。單個(gè)體檢指標(biāo)時(shí)間序列具有如下特征:

    a.數(shù)據(jù)量小。

    隨著時(shí)間變化,個(gè)人健康狀況受年齡變化、外界環(huán)境等因素影響,使得體檢指標(biāo)具有階段性和時(shí)效性。通常來說,極早期的體檢指標(biāo)對于分析個(gè)人當(dāng)前身體健康狀況的可用價(jià)值較低,許多體檢數(shù)據(jù)集中僅保留體檢者最近6~8年的體檢指標(biāo)數(shù)據(jù)。因此,體檢指標(biāo)時(shí)間序列 Xm =(xm(t1), xm(t2), ··· , xm(tn))的樣本數(shù)量非常有限,一般取樣本個(gè)數(shù)介于6~10之間。

    b.數(shù)據(jù)的不確定性。

    個(gè)人體檢指標(biāo)數(shù)值常受到生理狀況、心理變化、外界環(huán)境等多方面因素的影響,甚至由于測量儀器、檢測技術(shù)水平的參差不齊也會(huì)導(dǎo)致指標(biāo)數(shù)據(jù)的不準(zhǔn)確。所以個(gè)體在進(jìn)行體檢時(shí),總會(huì)對異常指標(biāo)進(jìn)行多次“復(fù)查”,將多次體檢結(jié)果的可能值或取值范圍作為最終檢查結(jié)果。這導(dǎo)致了體檢指標(biāo)序列的區(qū)間出現(xiàn)不確定或離散不確定的情況。

    c.時(shí)間間隔不一致。

    時(shí)間序列Xm =(xm(t1), xm(t2), ··· , xm(tn))的時(shí)間間隔計(jì)算公式一般為Δt = tk+1? tk ,k =1, 2, ··· , n?1,當(dāng)Δt *常數(shù)時(shí),將時(shí)間序列Xm稱為非等時(shí)距序列?,F(xiàn)實(shí)生活中,由于各種因素導(dǎo)致個(gè)體未能按期進(jìn)行健康體檢,從而導(dǎo)致體檢時(shí)間序列數(shù)據(jù)集中缺失某一時(shí)間段的數(shù)據(jù),出現(xiàn)時(shí)間“斷層”問題。

    d.數(shù)據(jù)類型異構(gòu)。

    體檢指標(biāo)數(shù)據(jù)類型異構(gòu)是指時(shí)間序列Xm中不同體檢指標(biāo)具有不同的數(shù)據(jù)類型。舉例來說,時(shí)間序列Xm中可能存在某一元素?cái)?shù)據(jù)類型是一個(gè)區(qū)間值,某一元素?cái)?shù)據(jù)類型為離散灰數(shù),還有元素?cái)?shù)據(jù)類型為實(shí)數(shù),這就使得Xm具有數(shù)據(jù)類型異構(gòu)的特征。

    e.數(shù)據(jù)具有上下波動(dòng)性。

    體檢指標(biāo)受到自身以及外部等多個(gè)因素制約,從而使得單個(gè)個(gè)體體檢時(shí)間序列并非呈現(xiàn)明顯的單調(diào)遞變或恒定不變的規(guī)律,通常是在一定數(shù)值范圍內(nèi)表現(xiàn)出反復(fù)的上下波動(dòng)的特征。

    2 灰色–時(shí)序組合預(yù)測模型 NDGM- ARIMA

    2.1 改進(jìn) GM(1,1)模型——NDGM(1,1)

    由于體檢指標(biāo)數(shù)據(jù)是一個(gè)數(shù)據(jù)量少的小樣本數(shù)據(jù)集,通常數(shù)據(jù)量級(jí)在幾至幾十。而灰色模型 GM(1,1)對于“少數(shù)據(jù)”、“貧信息”的樣本具有較高的預(yù)測精度。因此,本文考慮使用灰色模型 GM(1,1)。GM(1,1)模型是灰色系統(tǒng)理論中經(jīng)典的預(yù)測模型,模型的基本思路是利用原始數(shù)據(jù)得到一組原始數(shù)據(jù)序列,對原始數(shù)據(jù)序列進(jìn)行累加生成新的數(shù)據(jù)序列,以此來削弱原始數(shù)據(jù)的隨機(jī)性,突出和增強(qiáng)原始數(shù)據(jù)的規(guī)律性,實(shí)現(xiàn)對原始數(shù)據(jù)未來變化規(guī)律的模糊預(yù)測。

    GM(1,1)具體實(shí)現(xiàn)步驟如下:

    步驟1設(shè)原始數(shù)據(jù)構(gòu)成的序列為X(0),對原始序列進(jìn)行一次累加生成(1-AGO)得到新的數(shù)據(jù)序列X(1)。

    步驟2構(gòu)建新生成序列X(1)的緊鄰均值生成序列,記為Z(1)。由此得到 GM(1,1)模型的灰色微分方程x(0)(k)+az(1)(k)= b 。

    步驟3基 于最小二乘原理,可得到參數(shù)a, b滿足的條件為h(?)=(a, b)T =(BT B)?1 BTY,矩陣B 是構(gòu)造累加矩陣,向量 Y 為常數(shù)項(xiàng)向量。

    步驟4由 序列 X(0) , X(1) , Z(1)可得到 GM(1,1)模型的白化微分方程,將 GM(1,1)模型白化方程的解稱為時(shí)間響應(yīng)函數(shù)。

    步驟5求解得到白化微分方程的時(shí)間響應(yīng)序列后,通過累減生成還原得到原始序列為?(x)(0)(k+1)=?(x)(1)(k+1)??(x)(1)(k),即灰色 GM(1,1)的預(yù)測方程表達(dá)式,對其進(jìn)行求導(dǎo)還原就可得到序列還原值。

    傳統(tǒng)的 GM(1,1)模型是用一階微分方程對單個(gè)變量實(shí)現(xiàn)預(yù)測的模型,其建模過程主要是利用齊次指數(shù)序列來擬合原始數(shù)據(jù)。因此, GM(1,1)模型對于具有近齊次指數(shù)的原始序列具有較好的擬合與預(yù)測性能。但是,現(xiàn)實(shí)生活中存在許多不確定因素,絕大部分的時(shí)間序列都不符合指數(shù)增長規(guī)律。對于體檢指標(biāo)序列,這類序列由于數(shù)值結(jié)果不確定性大、時(shí)間間隔不統(tǒng)一導(dǎo)致的數(shù)值缺失,以及數(shù)據(jù)上下波動(dòng)等原因,使得體檢指標(biāo)序列更符合近似非齊次指數(shù)序列變化特征。同時(shí),傳統(tǒng)的 GM(1,1)模型中參數(shù)估計(jì)方程是離散的,模型預(yù)測方程是連續(xù)的,為了解決離散參數(shù)估計(jì)和連續(xù)預(yù)測表示之間跳躍所產(chǎn)生的模擬誤差,本文借鑒了謝乃明等[21]提出的離散灰色模型 DGM(1,1)基本思想,使改進(jìn)灰色模型的參數(shù)估計(jì)和模型預(yù)測都是離散形式。

    結(jié)合上述體檢序列特征分析和預(yù)測模型性能分析,為了構(gòu)建適用于體檢指標(biāo)序列的預(yù)測模型,本文構(gòu)建一個(gè)近似非齊次指數(shù)序列的離散 GM(1,1)模型(non-homogenous discrete grey model),簡寫為 NDGM(1,1)模型。

    同樣地 ,設(shè)原始非負(fù)序列為 X(0): X(0)= (x(0)(1), x(0)(2), ··· , x(0)(n))。其中,x(0)(i)>0, i =1, 2, ··· , n。經(jīng)過一次累加生成得到新序列 X(1):X(1)=(x(1)(1), x(1)(2), ··· , x(1)(n)),從而得到離散灰色模型 NDGM(1,1)的表達(dá)式為 x(1)(t +1)+ ax(1)(t)= bt+ c,則模型的白化微分方程表達(dá)式為

    式中,參數(shù)列 h(?)=(a, b, c)T 為 NDGM(1,1)模型待求解參數(shù)。

    求解 NDGM(1,1)模型白化方程的時(shí)間響應(yīng)序列,首先公式對應(yīng)的齊次方程為

    解出齊次方程的通解為x(1)(t)= C1e?at。利用常數(shù)變易法,令C1= f(t),則x(1)(t)= f(t)e?at。對x(1)(t)=f(t)e?at兩端同時(shí)求導(dǎo)后代入式(2)可得

    將式(5)代入x(1)(t)= C1e?at中,可知

    當(dāng) t =1時(shí),可得 x(1)(1)= t ? + +Ce?a,解出 C 的表達(dá)式為

    將式(7)代入式(6)得到 NDGM(1,1)模型的時(shí)間響應(yīng)序列表達(dá)式為

    則式(8)經(jīng)過累減還原得到還原式為

    當(dāng)t =2, 3, 4, ··· , n時(shí),?(x)(0)(t)為模型所得擬合值;當(dāng)t = n+1, n+2, ···時(shí),?(x)(0)(t)為模型所得預(yù)測值。

    令α= e?a,β=(1? e?a ),γ=(1? e?a )? +,則式(8)可表示為

    式(10)的參數(shù)列C(?)=(α,β, γ)T ,由最小二乘法得到參數(shù)的估計(jì)值,當(dāng)式(11)所示的誤差平方和達(dá)到最小時(shí)可求解出參數(shù)α,β, γ。

    參數(shù)列 C(?)=(α,β, γ)T應(yīng)滿足條件(α,β, γ)T =(BT B)?1 BTY,其中

    a, b, c的估計(jì)值分別為

    將參數(shù)估計(jì)值?(a), b(?), ?(c)代入式(9)所得的還原式,即可求出原始數(shù)據(jù)序列的模擬值和預(yù)測值。

    NDGM(1,1)模型建立后,為了評價(jià)模型運(yùn)行的可行性,需要對模型進(jìn)行精度檢驗(yàn),本文利用后殘差檢驗(yàn)法進(jìn)行檢驗(yàn)。記原始序列X(0): X(0)=(x(0)(1), x(0)(2), ··· , x(0)(n))和殘差序列ε(0)=(ε(1), ε(2),··· , ε(n))=(x(0)(1)??(x)(0)(1), x(0)(2)??(x)(0)(2), ··· , x(0)(n)??(x)(0)(n))的方差分別為 S1(2), S 2(2),計(jì)算公式分別為

    式中:?(x)(0)表示原始序列的均值,計(jì)算公式為?(x)(0)=x(0)(k);ε(?)(0)為序列殘差均值,且ε(?)(0)= ε(0)(k)。

    后驗(yàn)殘差檢驗(yàn)法是利用后驗(yàn)差比值c和小概率誤差p進(jìn)行檢驗(yàn),二者計(jì)算方法為

    若 NDGM(1,1)模型滿足表1所示的模型精度標(biāo)準(zhǔn),則說明構(gòu)建的 NDGM(1,1)模型合格。

    2.2 ARIMA 模型的構(gòu)建

    將時(shí)間序列定義為一組按時(shí)間先后順序排列的數(shù)據(jù)集合,時(shí)間序列預(yù)測就是指利用模型分析和處理時(shí)間序列,根據(jù)時(shí)間序列呈現(xiàn)出的規(guī)律,構(gòu)建有效的模型對數(shù)據(jù)未來發(fā)展趨勢進(jìn)行預(yù)測。常用于預(yù)測平穩(wěn)時(shí)間序列的時(shí)間序列模型包括自回歸模型 AR(n)、自回歸移動(dòng)平均模型 ARMA(p,q)、差分自回歸移動(dòng)平均模型 ARIMA(p,d ,q)。

    ARIMA(p,d ,q)模型的建模過程為,首先將非平穩(wěn)時(shí)間序列經(jīng)處理后轉(zhuǎn)化為平穩(wěn)時(shí)間序列,然后將因變量只對其滯后值(階數(shù))以及隨機(jī)誤差項(xiàng)的現(xiàn)值和滯后值進(jìn)行回歸分析。 ARIMA(p,d ,q)模型對于短期時(shí)間序列預(yù)測具有較高的預(yù)測精度。其中: AR 表示自回歸; MA 表示移動(dòng)平均; p, q分別表示 AR ,MA 的階數(shù); d表示差分的階數(shù),一般取值為1階或2階序列達(dá)到穩(wěn)定。ARIMA(p,d ,q)模型為

    式中:?d =(1? B)d為差分運(yùn)算;{εt }表示零均值白噪聲序列;Φ(B)=1??1B?···??p Bp ,Θ(B)=1?θ1B?···?θp Bp分別表示模型 ARIMA(p,d ,q)的自回歸系數(shù)多項(xiàng)式和移動(dòng)平均系數(shù)多項(xiàng)式;B 表示延遲算子,并且滿足Bnxt = xt?n 。

    ARIMA(p,d ,q)模型的建模包括時(shí)間序列預(yù)處理、模型識(shí)別和定階、模型檢驗(yàn)、模型驗(yàn)證及優(yōu)化和模型預(yù)測5個(gè)步驟。本文利用 Eviews 軟件進(jìn)行 ARIMA 模型確定和指標(biāo)預(yù)測,具體過程如下:

    a.時(shí)間序列預(yù)處理。適用于 ARIMA(p,d ,q)模型的時(shí)間序列必須為平穩(wěn)非白噪聲時(shí)間序列,對于非平穩(wěn)時(shí)間序列,需進(jìn)行數(shù)據(jù)預(yù)處理使原始序列滿足平穩(wěn)化和零均值的條件。將實(shí)驗(yàn)序列數(shù)據(jù)錄入 Eviews 軟件后,通過繪制原始序列的時(shí)序圖來判斷序列的平穩(wěn)性。若序列是非平穩(wěn)狀態(tài),采用取對數(shù)或差分處理等操作進(jìn)行處理,處理完后進(jìn)行 ADF 單位根檢驗(yàn)序列平穩(wěn)性。

    b.模型識(shí)別和定階。對于模型的識(shí)別和定階本質(zhì)上就是確定參數(shù)p ,q的值,基于數(shù)據(jù)預(yù)處理后的平穩(wěn)時(shí)間序列,計(jì)算出實(shí)驗(yàn)數(shù)據(jù)集的自相關(guān)系數(shù) ACF 和偏自相關(guān)系數(shù) PACF。對預(yù)處理后的序列通過 Eviews 軟件的 Correlogram 得到序列自相關(guān)圖和數(shù)值,采用 AIC 準(zhǔn)則為預(yù)測模型的階數(shù)p和 q取合適的值。

    c.模型檢驗(yàn)。對識(shí)別和定階后的 ARIMA 模型進(jìn)行參數(shù)估計(jì),模型的檢驗(yàn)包括參數(shù)估計(jì)的顯著性檢驗(yàn)和殘差序列的隨機(jī)性檢驗(yàn),即驗(yàn)證殘差之間的獨(dú)立性。確定 ARIMA 模型各項(xiàng)階數(shù)后,在 Eviews 中創(chuàng)建估計(jì)方程式得到 Prob.值, Prob.值若小于5%則模型是顯著的,可靠性較高。

    d.模型的驗(yàn)證和優(yōu)化。根據(jù)模型檢驗(yàn)結(jié)果對模型的階數(shù)進(jìn)行調(diào)整和優(yōu)化,使構(gòu)建出的模型滿足顯著性檢驗(yàn)要求。即若步驟 c 中得到的模型估計(jì)結(jié)果未通過檢驗(yàn),則返回修改模型階數(shù)p和q,重新進(jìn)行檢驗(yàn)。

    e.模型擬合和預(yù)測。利用構(gòu)建好的 ARIMA 模型對實(shí)驗(yàn)時(shí)間序列進(jìn)行擬合,并預(yù)測數(shù)據(jù)未來的趨勢。對于檢驗(yàn)通過的 ARIMA 模型利用 Eviews 中的 Forecast 模塊,在 sample 欄中選擇需預(yù)測的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行逐步向前預(yù)測。

    2.3 組合預(yù)測模型 NDGM-ARIMA

    各類預(yù)測模型的研究重點(diǎn)和關(guān)注方向都有所不同,因此,對同一個(gè)實(shí)驗(yàn)數(shù)據(jù)集進(jìn)行預(yù)測,不同的模型會(huì)產(chǎn)生不同的結(jié)果。為了提高預(yù)測模型的預(yù)測精度以及模型的適用性,本文將 NDGM(1,1)模型和 ARIMA(p,d ,q)模型進(jìn)行組合,簡稱 NDGM- ARIMA 模型。組合預(yù)測模型綜合考慮兩個(gè)模型的預(yù)測結(jié)果,通過為單個(gè)模型的預(yù)測結(jié)果賦予最佳的權(quán)重系數(shù),最大限度地利用多個(gè)模型的樣本信息。構(gòu)建組合模型,也在一定程度上減少了單個(gè)預(yù)測模型受外界因素的干擾,考慮問題更加全面系統(tǒng),從而提高模型預(yù)測的精度。

    本文構(gòu)建的 NDGM-ARIMA 組合預(yù)測模型用于實(shí)現(xiàn)個(gè)人體檢指標(biāo)序列的預(yù)測,模型具體的表達(dá)式為: X(?)(t)= wG(?)(t)+(1? w)A(?)(t)。其中: G(?)(t)表示NDGM(1,1)模型 t 時(shí)刻的預(yù)測值;A(?)(t)表示 ARIMA模型 t 時(shí)刻的預(yù)測值; w為組合模型權(quán)重值,取值范圍為w ∈[0, 1],表示單個(gè)模型預(yù)測結(jié)果的重要程度。

    在組合預(yù)測模型中,如何恰當(dāng)?shù)厍蠼獬鰴?quán)重系數(shù)是關(guān)鍵。確定權(quán)重系數(shù)常用方法包括算術(shù)平均法、最優(yōu)加權(quán)法、方差倒數(shù)法等。算術(shù)平均法是在對模型重要性缺乏了解時(shí)常用的權(quán)重選定方法,但是該方法缺乏對單個(gè)模型重要性的掌握,對每個(gè)模型賦予相同的權(quán)重,不分優(yōu)先順序使得預(yù)測效果不佳。最優(yōu)加權(quán)法需要求解線性或非線性規(guī)劃,計(jì)算復(fù)雜并且計(jì)算結(jié)果有可能為負(fù),在實(shí)際應(yīng)用中具有較大的局限性。方差倒數(shù)法則是通過預(yù)測模型的誤差平方和的計(jì)算來反映預(yù)測精度,相較于算數(shù)平均法和最優(yōu)加權(quán)法,直接應(yīng)用預(yù)測誤差平方和更能反映各個(gè)模型在組合預(yù)測中的重要程度,賦予的權(quán)重?cái)?shù)值更為合理有效。而且方差倒數(shù)法易操作,獲得的預(yù)測效果好。因此,為求解預(yù)測模型最佳的組合權(quán)重大小,本文采用方差倒數(shù)這一方法。方差倒數(shù)的目的是使組合預(yù)測模型的誤差平方和盡可能小。因此,需要對組合模型中誤差平方和大的模型賦較小的權(quán)重值,對誤差平方和小的模型賦較大的權(quán)重值。

    采用方差倒數(shù)進(jìn)行組合權(quán)重賦值,首先計(jì)算出單個(gè)預(yù)測模型的預(yù)測誤差平方和。用ei表示第 i個(gè)模型的誤差平方和,其計(jì)算方式如式(18)所示。

    式中: xi為原始數(shù)據(jù);?(x)ti 為其對應(yīng)的預(yù)測值;(xi ??(x)ti)為預(yù)測誤差。

    計(jì)算出單個(gè)模型的誤差平方和在全部模型中的占比,這一占比即該模型的權(quán)重值大小。利用模型的預(yù)測誤差得到權(quán)重系數(shù)的計(jì)算公式為

    式中,Σ(m)wi ?1=1, j =1, 2, ··· , m。

    由式(19)可以發(fā)現(xiàn),當(dāng)單個(gè)模型的誤差平方和越大時(shí),獲得的權(quán)重越小,則模型預(yù)測精度越低,預(yù)測結(jié)果的價(jià)值度越低。

    3 個(gè)體體檢指標(biāo)預(yù)測實(shí)驗(yàn)及結(jié)果分析

    3.1 實(shí)驗(yàn)數(shù)據(jù)集描述

    心血管疾病已成為當(dāng)前社會(huì)的一種高發(fā)疾病,該類疾病的高危致病因素眾多,包括高血壓、糖尿病、肥胖、血脂異常、吸煙和過度飲酒等。由相關(guān)統(tǒng)計(jì)數(shù)據(jù)可發(fā)現(xiàn),近年來,心血管疾病患者死亡率極高,所以人們必須對此類疾病引起重視,加強(qiáng)自身健康管理。患者通過定期健康體檢,可以幫助醫(yī)生和患者及時(shí)了解當(dāng)前身體狀況,發(fā)現(xiàn)關(guān)鍵病因信號(hào),提前進(jìn)行預(yù)防和治療,降低患病的風(fēng)險(xiǎn)。因此,構(gòu)建適當(dāng)?shù)念A(yù)測模型,實(shí)現(xiàn)對人體主要健康指標(biāo)序列的有效預(yù)測,具有重要的現(xiàn)實(shí)意義。

    本文采用天池公開數(shù)據(jù)集中的心臟病體檢數(shù)據(jù)集進(jìn)行分析,數(shù)據(jù)集中包含多名體檢者連續(xù)多年的體檢數(shù)值,例如血脂水平中甘油三酯、總膽固醇、高密度脂蛋白膽固醇、低密度脂蛋白膽固醇4項(xiàng)指標(biāo)和空腹血糖指標(biāo)等數(shù)值。實(shí)驗(yàn)選擇空腹血糖指標(biāo)作為實(shí)驗(yàn)數(shù)據(jù)序列,血糖指標(biāo)是檢測心血管疾病和糖尿病的關(guān)鍵指標(biāo),同時(shí)也是人體健康管理中重要的體檢指標(biāo),關(guān)注血糖值的變化可以有效監(jiān)測到心血管類疾病??崭寡侵笜?biāo)的正常取值為3.9~6.1 mmol/L 。在4個(gè)不同年齡段(20~30歲,30~40歲,40~50歲,50~60歲)中隨機(jī)選擇一名體檢者,對4名體檢者的空腹血糖指標(biāo)進(jìn)行擬合和預(yù)測。4名體檢者2005—2014年指標(biāo)的空腹血糖體檢時(shí)間序列為表2,將4名實(shí)驗(yàn)對象樣本分別用X1,X2,X3,X4表示。

    3.2 模型預(yù)測結(jié)果及分析

    為了更加直觀地分析組合預(yù)測模型的性能,利用 ARIMA(p,d,q),GM(1,1),NDGM(1,1),NDGM- ARIMA 組合預(yù)測模型4個(gè)模型對血糖體檢時(shí)間序列進(jìn)行擬合和預(yù)測,通過分析各模型的預(yù)測值和相對模擬誤差Δ(t)來分析組合預(yù)測模型的預(yù)測性能。相對模擬誤差計(jì)算公式如下:

    3.2.1 ARIMA(2,2,1)模型預(yù)測

    實(shí)驗(yàn)分別對4名體檢者空腹血糖時(shí)間序列建立相應(yīng)的 ARIMA 預(yù)測模型。例如45歲體檢者 的時(shí)間序列由原始序列 X(0)=(x(0)(1), x(0)(2), ··· , x(0)(10))=(5.69, 5.03, ··· , 5.78)可知,原始序列是一 個(gè)非平穩(wěn)時(shí)間序列,首先進(jìn)行差分處理轉(zhuǎn)化為平穩(wěn)序列。將數(shù)據(jù)輸入 Eviews 軟件中,對原始序列進(jìn)行 ADF 檢驗(yàn), ADF 檢驗(yàn)結(jié)果如圖1所示。可發(fā)現(xiàn)當(dāng)二階差分時(shí),所有 t 值的絕對值均小于 ADF 檢驗(yàn)統(tǒng)計(jì)量的絕對值,且p 值為0.0090,小于0.05,說明原序列已轉(zhuǎn)化為平穩(wěn)時(shí)間序列,則 ARIMA 模型的差分階數(shù)為d =2。

    接著,對模型進(jìn)行識(shí)別,確定模型的 ACF 和 PACF。利用 Eviews 軟件 Correlogram 相關(guān)圖查看序列二階差分的 ACF 和 PACF 值,得到如圖2所示的自相關(guān)圖。

    由圖2可知,時(shí)間序列的自相關(guān)系數(shù) ACF 在1階截尾,偏自相關(guān)系數(shù) PACF 在2階截尾。因此,構(gòu)建 ARIMA(2,2,1)模型對空腹血糖體檢序列進(jìn)行預(yù)測。之后,在 Eviews 軟件中進(jìn)行建模,采用列表法對 ARIMA 方程進(jìn)行定義: data c ar(1) ar(2) ma(1),根據(jù)定義后的模型得到 ARIMA(2,2,1)模型具體的表達(dá)式為

    對于27,35,57歲體檢者血糖序列,同樣利用 Eviews 軟件建立最優(yōu)的 ARIMA 模型,得到 27歲體檢者血糖序列的時(shí)間序列預(yù)測模型為 ARIMA(2,1,1),35歲對應(yīng)模型為 ARIMA(3,1,2),57歲對應(yīng)模型為 ARIMA(3,2,2)。

    3.2.2 NDGM(1,1)模型預(yù)測

    同樣地,對于4個(gè)時(shí)間序列建立對應(yīng)的 GM(1, 1)模型和離散灰色模型 NDGM(1,1)。以表2所示的45歲體檢者的血糖數(shù)據(jù)為具體例子進(jìn)行建模,可知該體檢者空腹血糖指標(biāo)原始序列為 X(0)=(x(0)(1), x(0)(2), ··· , x(0)(10))=(5.09, 5.03, ··· , 5.78),利 用 python 代碼建立序列X(0)的 GM(1,1)模型,得到模型參數(shù)a =?0.01026,b =4.946395,則關(guān)于空腹血糖指標(biāo)預(yù)測的 GM(1,1)模型的時(shí)間響應(yīng)表達(dá)式為

    進(jìn)一步對優(yōu)化后的 NDGM(1,1)模型的參數(shù)α(?) ,β(?) , γ(?)及a ,b ,c進(jìn)行參數(shù)估計(jì),計(jì)算出具體的數(shù)值結(jié)果,α(?)=0.1286,β(?)=4.5564,γ(?)=5.3875,a =2.0507, b =10.7234,c =5.6018,得到 NDGM(1,1)模型為 ?(x)(0)(t)=0.2791(1?e2.0507)e?2.0507(t?1)+5.2291(23)

    同理可得:27歲體檢者的 GM(1,1)模型和 NDGM(1,1)模型表達(dá)式分別如式(24)和(25)所示;35歲體檢者的 GM(1,1)模型和 NDGM(1,1)模型表達(dá)式分別如式(26)和(27)所示;57歲體檢者的 GM(1,1)模型和 NDGM(1,1)模型表達(dá)式分別如式(28)和(29)所示。

    3.2.3灰色–時(shí)間序列組合模型 NDGM-ARIMA 預(yù)測

    將4個(gè)預(yù)測模型 ARIMA(2,2,1), GM(1,1), NDGM(1,1)和 NDGM-ARIMA 組合模型分別對 4名體檢者2005—2014年空腹血糖體檢序列進(jìn)行預(yù)測,各個(gè)模型對45歲體檢者血糖的預(yù)測結(jié)果如表3所示,4名體檢者的整體預(yù)測結(jié)果如圖3所 示。利用式(19)的權(quán)重系數(shù)計(jì)算方法確定組合模型的權(quán)重系數(shù),得到在對45歲體檢者進(jìn)行預(yù)測時(shí), NDGM(1,1)模型和 ARIMA(2,2,1)模型的權(quán)重系數(shù)分別為0.6286,0.3714。

    由圖3所示的4名體檢者的預(yù)測結(jié)果曲線和實(shí)際數(shù)據(jù)曲線對比分析可知,論文對于35歲體檢者的空腹血糖指標(biāo)預(yù)測結(jié)果并非是對比模型中最佳的。這有可能是因?yàn)樵跀?shù)據(jù)集中,該體檢者初始體檢年份血糖指標(biāo)與最終體檢年份的指標(biāo)數(shù)值相差較大。由于存在各種外界因素導(dǎo)致的兩個(gè)體檢數(shù)據(jù)的不準(zhǔn)確和差距較大,使得模型的誤差較大,從而導(dǎo)致預(yù)測精度下降。但是,通過進(jìn)一步分析35歲體檢者空腹血糖指標(biāo)預(yù)測值可以發(fā)現(xiàn),構(gòu)建的組合模型與最優(yōu)預(yù)測模型二者間的預(yù)測值相差極小。同時(shí),組合預(yù)測模型在其余3個(gè)年齡段的體檢者的血糖值擬合上都更接近真實(shí)數(shù)值,說明了組合模型對于絕大多數(shù)的體檢數(shù)據(jù)預(yù)測是有效的,也證明了組合模型預(yù)測結(jié)果的真實(shí)性、高可信度。

    進(jìn)一步對45歲體檢者血糖指標(biāo)預(yù)測具體數(shù)值進(jìn)行分析。與傳統(tǒng)的 GM(1,1)模型對比,改進(jìn)的灰色模型 NDGM(1,1)在實(shí)驗(yàn)序列上的擬合值雖然存在部分預(yù)測值差于 GM(1,1)模型,但是從兩個(gè)模型的平均相對誤差來看, NDGM(1,1)模型的平均相對誤差為0.0501,GM(1,1)模型的平均相對誤差為0.0528,NDGM(1,1)模型的平均相對誤差小于 GM(1,1)。這一實(shí)驗(yàn)結(jié)果顯示,構(gòu)建的 NDGM(1,1)模型在體檢指標(biāo)序列預(yù)測上整體的預(yù)測效果優(yōu)于 GM(1,1)模型,說明構(gòu)建的改進(jìn)灰色預(yù)測模型在預(yù)測精度上得到了一定程度的提升。

    通過 NDGM-ARIMA 模型與3個(gè)單個(gè)預(yù)測模型的對比,組合模型的擬合值和相對誤差都優(yōu)于單個(gè)灰色預(yù)測和時(shí)間序列模型,這說明組合模型確實(shí)適用于健康體檢指標(biāo)的預(yù)測,模型的擬合值更加接近實(shí)際體檢數(shù)據(jù)。另外,這也證明了組合模型能夠更好地結(jié)合單個(gè) NDGM(1,1)模型和時(shí)間序列模型的優(yōu)點(diǎn),在一定程度上克服了單一預(yù)測模型的局限性,提高了模型的預(yù)測精度。

    利用 NDGM-ARIMA 組合模型預(yù)測45歲體檢者2015—2018年的血糖數(shù)值,預(yù)測結(jié)果如表4所示。

    已知空腹血糖的正常范圍為3.6~6.1 mmol/L,結(jié)合2015—2018年的預(yù)測值來分析該體檢者身體狀況變化趨勢。由預(yù)測數(shù)值可發(fā)現(xiàn),該體檢者的空腹血糖指標(biāo)數(shù)值呈現(xiàn)緩慢上升的趨勢,預(yù)計(jì)到2017年血糖指標(biāo)數(shù)值將達(dá)到6.12 mmol/L,已經(jīng)突破人體空腹血糖正常值最大臨界值,體檢者極有可能患糖尿病等疾病,危害身體健康。因此,由預(yù)測結(jié)果可以得出,體檢者未來幾年患糖尿病和心血管疾病的潛在風(fēng)險(xiǎn)很大,必須注意自身糖分的攝入,加強(qiáng)身體日常管理,提前做好預(yù)防措施或采取及時(shí)的治療手段。

    4 結(jié) 論

    傳統(tǒng)體檢指標(biāo)分析僅局限于單次指標(biāo)數(shù)值高低的靜態(tài)分析,忽略了因個(gè)體差異導(dǎo)致的體檢數(shù)據(jù)的動(dòng)態(tài)變化趨勢。因此,構(gòu)建合理有效的數(shù)據(jù)模型來挖掘體檢指標(biāo)的發(fā)展規(guī)律,準(zhǔn)確預(yù)測體檢數(shù)值的變化趨勢和未來取值范圍,并基于預(yù)測結(jié)果對個(gè)體健康狀況實(shí)施預(yù)警管理,通過監(jiān)測人體主要健康指標(biāo)的變化,及時(shí)發(fā)現(xiàn)潛在的患病因子或風(fēng)險(xiǎn)因素,進(jìn)一步采取有效的預(yù)防和治療措施,對于實(shí)現(xiàn)個(gè)體健康管理具有重要的現(xiàn)實(shí)意義。

    為了構(gòu)建適用于個(gè)體主要健康體檢指標(biāo)的預(yù)測模型,加強(qiáng)模型在體檢指標(biāo)上的預(yù)測性能,本文提出一個(gè)改進(jìn)灰色模型和時(shí)間序列模型相結(jié)合的組合預(yù)測模型。首先分析體檢指標(biāo)序列的特征,考慮到體檢指標(biāo)序列是一個(gè)近似非齊次指數(shù)序列,以及 GM(1,1)模型中的離散和連續(xù)之間的誤差,構(gòu)建了一個(gè)近似非齊次指數(shù)序列的離散灰色模型 NDGM(1,1)。其次,為了將單個(gè)預(yù)測模型的優(yōu)勢結(jié)合在一起,論文將時(shí)間序列預(yù)測模型 ARIMA(p,d ,q)和 NDGM(1,1)模型進(jìn)行組合得到 NDGM-ARIMA 模型。在盡可能保證組合模型誤差平方和達(dá)到最小的情況下,為兩個(gè)模型的預(yù)測結(jié)果賦予最佳權(quán)重系數(shù),并將加權(quán)后的結(jié)果作為最終的模型擬合結(jié)果和預(yù)測結(jié)果。NDGM-ARIMA 組合模型在血糖體檢指標(biāo)數(shù)據(jù)集上的預(yù)測結(jié)果表明,組合模型在體檢指標(biāo)序列上的預(yù)測精度有所提高,保證了預(yù)測結(jié)果的有效性和準(zhǔn)確性,從而可以利用預(yù)測結(jié)果有效地分析出個(gè)人主要健康體檢指標(biāo)的變化趨勢,實(shí)現(xiàn)人們健康管理的目標(biāo)。

    但是,本文模型存在一定的局限性。首先,本文研究數(shù)據(jù)集為等時(shí)距的近似非齊次指數(shù)序列,然而,實(shí)際應(yīng)用中存在大量的非等間距的近似非齊次指數(shù)序列,容易導(dǎo)致因數(shù)據(jù)序列類型不符合預(yù)測模型而出現(xiàn)較大的建模誤差。因此,如何進(jìn)一步拓展灰色預(yù)測模型的適用范圍將成為未來的研究方向。其次,本文組合模型中僅使用了方差倒數(shù)法求解各預(yù)測模型權(quán)重,但是單一賦權(quán)的方式可能存在較大的權(quán)重求解誤差,因此在對多種賦權(quán)方法研究的基礎(chǔ)上,是否可通過將兩種及以上賦權(quán)方法結(jié)合起來進(jìn)行求權(quán),從而提高預(yù)測模型建模精度,同樣是本文進(jìn)一步的研究方向。

    參考文獻(xiàn):

    [1]周偉杰, 張宏如, 黨耀國, 等.新息優(yōu)先累加灰色離散模型的構(gòu)建及應(yīng)用[J].中國管理科學(xué) , 2017, 25(8):140–148.

    [2]曾波, 劉思峰, 曲學(xué)鑫.一種強(qiáng)兼容性的灰色通用預(yù)測模型及其性質(zhì)研究[J].中國管理科學(xué) , 2017, 25(5):150–156.

    [3]陸冬磊, 吳春峰, 段勝剛, 等.應(yīng)用灰色模型 GM(1, 1)預(yù)測上海市副溶血性弧菌引起的食源性疾病發(fā)病率[J].環(huán)境與職業(yè)醫(yī)學(xué), 2015, 32(8):728–730.

    [4]曾波, 劉思峰, 白云, 等.基于灰色系統(tǒng)建模技術(shù)的人體疾病早期預(yù)測預(yù)警研究[J].中國管理科學(xué), 2020, 28(1):144–152.

    [5]華來慶, 申廣榮, 熊林平, 等. ARIMA 模型在黃瓜霜霉病疾病指數(shù)時(shí)間序列建模中的應(yīng)用研究[J].第二軍醫(yī)大學(xué)學(xué)報(bào), 2006, 27(7):729–732.

    [6] VALIPOUR M, BANIHABIB M E, BEHBAHANI S M R. Comparison? of the? ARMA,? ARIMA,? and? the autoregressive artificial neural network models inforecasting the monthly inflow of Dez dam reservoir[J]. Journal of Hydrology, 2013, 476:433–441.

    [7]劉瓊, 楊建華.隱馬爾科夫模型在乙肝發(fā)病預(yù)測中的應(yīng)用[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí), 2017, 47(19):203–210.

    [8]聶雄, 陳華, 伍思霖.基于灰度共生矩陣和 BP 神經(jīng)網(wǎng)絡(luò)的乳腺腫瘤識(shí)別[J].電子技術(shù)應(yīng)用 , 2019, 45(7):97–101,116.

    [9]王振飛, 陳金磊, 鄭志蘊(yùn), 等.面向心血管疾病的自適應(yīng)模塊化神經(jīng)網(wǎng)絡(luò)預(yù)測模型[J].小型微型計(jì)算機(jī)系統(tǒng) ,2019, 40(1):232–235.

    [10] ASILT?RK ?, ?UNKA? M. Modeling and prediction of surface roughness in turning operations using artificial neural network and multiple regression method[J]. Expert Systems With Applications, 2011, 38(5):5826–5832.

    [11] MCCLELLAND G H, IRWIN J R, DISATNIK D, et al. Multicollinearity is a red herring in the search for moderator variables: a guide to interpreting moderated multiple regression models and a critique of Iacobucci, Schneider, Popovich, and Bakamitsos (2016)[J]. Behavior Research Methods, 2017, 49(1):394–402.

    [12] SHIOTA S, OKAMOTO Y, OKADA G, et al. The neural correlates of the metacognitive function of other perspective: a multiple regression analysis study[J]. Neuroreport, 2017, 28(11):671–676.

    [13] YOO H Y, LEE J H, KIM D S, et al. Enhancement of glucose yield from canola agricultural residue by alkali pretreatment based on multi-regression models[J]. Journal of Industrial and EngineeringChemistry, 2017, 51:303–311.

    [14]王永斌, 李向文, 柴峰, 等.采用灰色–廣義回歸神經(jīng)網(wǎng)絡(luò)組合模型預(yù)測我國塵肺病發(fā)病人數(shù)的方法探討[J].環(huán)境與職業(yè)醫(yī)學(xué), 2016, 33(10):984–987.

    [15]嚴(yán)薇榮, 徐勇, 楊小兵, 等.基于 ARIMA-GRNN 組合模型的傳染病發(fā)病率預(yù)測[J].中國衛(wèi)生統(tǒng)計(jì), 2008, 25(1):82–83.

    [16]時(shí)冬青, 宋文華, 張桂釧, 等.基于灰色 GM(1, 1)–馬爾科夫模型的職業(yè)病預(yù)測研究[J].中國安全生產(chǎn)科學(xué)技術(shù), 2017, 13(4):176–180.

    [17]范文俊, 劉靜怡, 史菲, 等.新型炎性指標(biāo)對冠狀動(dòng)脈疾病的診斷預(yù)測價(jià)值[J].醫(yī)學(xué)研究雜志 , 2021, 50(1):80–85.

    [18]章玫紅.血清 Hcy 水平在心腦血管類疾病風(fēng)險(xiǎn)預(yù)測中的價(jià)值分析[J].中國社區(qū)醫(yī)師, 2021, 37(16):103–104.

    [19]劉丹, 趙森, 顏志良, 等.基于堆疊自動(dòng)編碼器的 miRNA-疾病關(guān)聯(lián)預(yù)測方法[J].計(jì)算機(jī)科學(xué), 2021, 48(10):114–120.

    [20]苗立志, 白瑞思蒙, 劉成良, 等.面向非平衡數(shù)據(jù)的癌癥患者生存預(yù)測分析[J].計(jì)算機(jī)工程, 2021, 47(12):316–320.

    [21]謝乃明 , 劉思峰.一類離散灰色模型及其預(yù)測效果研究[J].系統(tǒng)工程學(xué)報(bào), 2006, 21(5):520–523.

    (編輯:丁紅藝)

    猜你喜歡
    灰色
    瓊瑤 從灰色成長中生發(fā)粉色浪漫
    海峽姐妹(2020年5期)2020-06-22 08:26:08
    淺灰色的小豬
    灰色層染
    灰色時(shí)代
    Coco薇(2017年2期)2017-04-25 17:59:38
    她、它的灰色時(shí)髦觀
    Coco薇(2017年2期)2017-04-25 17:57:49
    在線直播如何走出灰色地帶
    新聞傳播(2016年22期)2016-07-12 10:09:54
    感覺
    中國衛(wèi)生(2015年3期)2015-11-19 02:53:28
    抽油井間抽周期的灰色預(yù)測
    斷塊油氣田(2012年5期)2012-03-25 09:53:49
    灰色
    文苑·感悟(2009年8期)2009-02-11 11:07:44
    黄色视频,在线免费观看| 麻豆久久精品国产亚洲av| 一个人看视频在线观看www免费| 人妻丰满熟妇av一区二区三区| 中出人妻视频一区二区| 成熟少妇高潮喷水视频| 国产精品永久免费网站| 国产又黄又爽又无遮挡在线| 欧美极品一区二区三区四区| 成人高潮视频无遮挡免费网站| 成人国产一区最新在线观看| a级毛片a级免费在线| 在线观看午夜福利视频| 久久国产精品人妻蜜桃| 亚洲国产高清在线一区二区三| 国产不卡一卡二| 国产精品爽爽va在线观看网站| 一级a爱片免费观看的视频| 久久久国产成人精品二区| x7x7x7水蜜桃| 欧美3d第一页| 国产精品嫩草影院av在线观看 | 国产白丝娇喘喷水9色精品| 国产精品一区二区性色av| 日韩欧美国产一区二区入口| 美女免费视频网站| 亚洲精品色激情综合| 长腿黑丝高跟| 精品熟女少妇八av免费久了| 精品国内亚洲2022精品成人| 亚洲精品在线观看二区| 制服丝袜大香蕉在线| 狂野欧美白嫩少妇大欣赏| 亚洲av成人精品一区久久| 波野结衣二区三区在线| 免费av毛片视频| 免费搜索国产男女视频| 精品99又大又爽又粗少妇毛片 | 精品久久久久久,| 最新中文字幕久久久久| 国产精品女同一区二区软件 | 亚洲第一电影网av| 欧美精品国产亚洲| 亚洲五月婷婷丁香| 99精品在免费线老司机午夜| 精品一区二区三区视频在线观看免费| 国产av麻豆久久久久久久| 九九久久精品国产亚洲av麻豆| 欧美精品啪啪一区二区三区| 中文亚洲av片在线观看爽| 在线a可以看的网站| 又紧又爽又黄一区二区| 亚洲人成网站在线播放欧美日韩| 国产高潮美女av| 91av网一区二区| 国产免费av片在线观看野外av| 99久久无色码亚洲精品果冻| 亚洲第一电影网av| 永久网站在线| 国产三级黄色录像| 热99re8久久精品国产| 免费在线观看日本一区| 大型黄色视频在线免费观看| 国产精品精品国产色婷婷| 日韩有码中文字幕| 熟妇人妻久久中文字幕3abv| 久久精品国产自在天天线| 欧美成人性av电影在线观看| 成年女人毛片免费观看观看9| 9191精品国产免费久久| 亚洲成人精品中文字幕电影| 亚洲七黄色美女视频| 99久国产av精品| 成人永久免费在线观看视频| 少妇丰满av| 色播亚洲综合网| 国产久久久一区二区三区| 免费av不卡在线播放| 国产伦精品一区二区三区视频9| 亚洲av熟女| 欧美极品一区二区三区四区| 午夜两性在线视频| 久久精品影院6| 99精品久久久久人妻精品| 麻豆一二三区av精品| 午夜福利高清视频| 深爱激情五月婷婷| 成年女人毛片免费观看观看9| 最近在线观看免费完整版| 成年女人看的毛片在线观看| 亚洲在线自拍视频| 国产伦精品一区二区三区四那| 亚洲av成人精品一区久久| 天天一区二区日本电影三级| 久久久久久九九精品二区国产| 热99在线观看视频| 最近视频中文字幕2019在线8| 桃红色精品国产亚洲av| 国产单亲对白刺激| 99热精品在线国产| 久久国产精品影院| 97人妻精品一区二区三区麻豆| 免费高清视频大片| 我要看日韩黄色一级片| 国模一区二区三区四区视频| 夜夜夜夜夜久久久久| 精品久久久久久久久久免费视频| 人妻制服诱惑在线中文字幕| 12—13女人毛片做爰片一| 婷婷色综合大香蕉| 久久精品国产自在天天线| 国产精品一区二区三区四区久久| 偷拍熟女少妇极品色| 午夜福利视频1000在线观看| 最新在线观看一区二区三区| 啦啦啦观看免费观看视频高清| 757午夜福利合集在线观看| 亚洲成人中文字幕在线播放| 桃红色精品国产亚洲av| 可以在线观看毛片的网站| 嫩草影视91久久| 99国产综合亚洲精品| 精品午夜福利视频在线观看一区| 亚洲电影在线观看av| 国产av不卡久久| 久久久久性生活片| 特级一级黄色大片| 亚洲激情在线av| 欧美乱妇无乱码| 成人精品一区二区免费| 搡女人真爽免费视频火全软件 | 午夜福利高清视频| 99国产精品一区二区蜜桃av| 国产色爽女视频免费观看| 亚洲成人久久性| a级毛片a级免费在线| 国产成人欧美在线观看| 国语自产精品视频在线第100页| 香蕉av资源在线| 欧美乱妇无乱码| 在线国产一区二区在线| 久久国产乱子免费精品| 免费观看的影片在线观看| 成人欧美大片| 一区福利在线观看| 国产真实乱freesex| 中文资源天堂在线| 国产精品一及| 久久久色成人| 激情在线观看视频在线高清| 成年女人永久免费观看视频| 国产成人啪精品午夜网站| 国产黄片美女视频| 日韩成人在线观看一区二区三区| 他把我摸到了高潮在线观看| 国内毛片毛片毛片毛片毛片| 青草久久国产| 精品久久久久久成人av| 久久久久免费精品人妻一区二区| 成人av一区二区三区在线看| 有码 亚洲区| 免费观看的影片在线观看| 赤兔流量卡办理| 精品久久久久久久久亚洲 | 国产大屁股一区二区在线视频| 亚洲国产高清在线一区二区三| 国产成年人精品一区二区| 男人的好看免费观看在线视频| 中文字幕精品亚洲无线码一区| 99精品在免费线老司机午夜| 韩国av一区二区三区四区| 成年女人看的毛片在线观看| 亚洲真实伦在线观看| 麻豆国产97在线/欧美| 男人的好看免费观看在线视频| 又爽又黄a免费视频| 看十八女毛片水多多多| 久久婷婷人人爽人人干人人爱| 91在线精品国自产拍蜜月| 欧美3d第一页| 国产精品99久久久久久久久| 波多野结衣高清无吗| 久久久久久久午夜电影| 成年人黄色毛片网站| 国产在视频线在精品| 白带黄色成豆腐渣| 一级a爱片免费观看的视频| 偷拍熟女少妇极品色| 免费在线观看亚洲国产| 午夜视频国产福利| 精品人妻偷拍中文字幕| 亚洲国产欧洲综合997久久,| 欧美极品一区二区三区四区| 午夜福利在线在线| 国产精品嫩草影院av在线观看 | 亚洲经典国产精华液单 | 久久久久免费精品人妻一区二区| 国产高清视频在线观看网站| 亚洲欧美精品综合久久99| 夜夜看夜夜爽夜夜摸| 国产一区二区三区在线臀色熟女| 精品人妻一区二区三区麻豆 | 国产午夜精品论理片| 九九久久精品国产亚洲av麻豆| 日韩 亚洲 欧美在线| 亚洲欧美日韩东京热| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 12—13女人毛片做爰片一| av视频在线观看入口| 亚洲欧美日韩高清在线视频| 亚洲美女视频黄频| 又爽又黄无遮挡网站| 久久精品91蜜桃| 美女xxoo啪啪120秒动态图 | 亚洲av熟女| 窝窝影院91人妻| 最新中文字幕久久久久| 国产亚洲欧美98| 九九在线视频观看精品| 99国产综合亚洲精品| 十八禁国产超污无遮挡网站| 一区福利在线观看| 欧美黄色淫秽网站| 亚洲精品亚洲一区二区| 欧美乱妇无乱码| 国内精品一区二区在线观看| or卡值多少钱| 成人一区二区视频在线观看| 国产欧美日韩一区二区三| 久久99热6这里只有精品| 久久精品91蜜桃| 91在线精品国自产拍蜜月| 精品国内亚洲2022精品成人| 黄片小视频在线播放| 桃色一区二区三区在线观看| 一进一出抽搐gif免费好疼| 欧美在线一区亚洲| 色5月婷婷丁香| 欧美黑人欧美精品刺激| 精品久久久久久久久久久久久| 亚洲av成人不卡在线观看播放网| 中国美女看黄片| 老女人水多毛片| 国产精品久久久久久久电影| 美女免费视频网站| 亚洲精品乱码久久久v下载方式| 久久久久久久久久成人| 深夜a级毛片| 丝袜美腿在线中文| 中文在线观看免费www的网站| 日本黄色视频三级网站网址| 久久久久久久精品吃奶| 特级一级黄色大片| 亚洲性夜色夜夜综合| 日本一本二区三区精品| 免费av毛片视频| 欧美中文日本在线观看视频| 午夜福利在线在线| 亚洲成av人片免费观看| 又紧又爽又黄一区二区| 丰满乱子伦码专区| 一卡2卡三卡四卡精品乱码亚洲| 日日夜夜操网爽| 很黄的视频免费| 人人妻人人看人人澡| 亚洲精品亚洲一区二区| 18禁在线播放成人免费| 亚洲av电影不卡..在线观看| 亚洲经典国产精华液单 | 亚洲人与动物交配视频| 国产亚洲精品综合一区在线观看| 亚洲av成人av| 国产色爽女视频免费观看| 欧美bdsm另类| 国产日本99.免费观看| 麻豆久久精品国产亚洲av| 午夜福利在线观看吧| 亚洲无线观看免费| 一个人免费在线观看的高清视频| 日韩大尺度精品在线看网址| 五月玫瑰六月丁香| 国产高清三级在线| 草草在线视频免费看| 欧美丝袜亚洲另类 | 国产一级毛片七仙女欲春2| 国产伦人伦偷精品视频| netflix在线观看网站| 宅男免费午夜| 老司机午夜十八禁免费视频| 亚洲av二区三区四区| 中文字幕熟女人妻在线| 亚洲国产精品久久男人天堂| 日本免费一区二区三区高清不卡| 99在线视频只有这里精品首页| 日本成人三级电影网站| 中亚洲国语对白在线视频| 国产午夜精品论理片| 国产三级黄色录像| 日日夜夜操网爽| 美女黄网站色视频| 国产大屁股一区二区在线视频| eeuss影院久久| 色视频www国产| 大型黄色视频在线免费观看| 精品人妻熟女av久视频| 每晚都被弄得嗷嗷叫到高潮| 欧美成狂野欧美在线观看| 在线观看66精品国产| 欧美成人性av电影在线观看| 欧美日韩瑟瑟在线播放| 精品免费久久久久久久清纯| 亚洲一区二区三区不卡视频| 中文字幕熟女人妻在线| 91在线观看av| 欧美性感艳星| 成人永久免费在线观看视频| 制服丝袜大香蕉在线| 国产亚洲精品av在线| 有码 亚洲区| 淫妇啪啪啪对白视频| 久久精品国产亚洲av香蕉五月| 亚洲第一欧美日韩一区二区三区| 精品久久久久久,| 久久久久性生活片| 又粗又爽又猛毛片免费看| 精品一区二区三区视频在线| 国产高潮美女av| 午夜精品在线福利| av在线蜜桃| 免费观看的影片在线观看| 麻豆成人av在线观看| 久久久久久久精品吃奶| 看十八女毛片水多多多| 色精品久久人妻99蜜桃| 欧美激情久久久久久爽电影| 国产精品久久久久久人妻精品电影| 亚洲va日本ⅴa欧美va伊人久久| 深爱激情五月婷婷| 特级一级黄色大片| 日日干狠狠操夜夜爽| 亚洲午夜理论影院| 噜噜噜噜噜久久久久久91| 国产精品电影一区二区三区| 国产av不卡久久| 天堂网av新在线| 国产精品一及| 久久久久久久亚洲中文字幕 | 成人精品一区二区免费| 亚洲av电影在线进入| 成人精品一区二区免费| 国产精品自产拍在线观看55亚洲| 丁香六月欧美| 亚洲av日韩精品久久久久久密| 可以在线观看毛片的网站| 久久人人爽人人爽人人片va | 精品一区二区免费观看| 日韩欧美在线二视频| 亚洲片人在线观看| 国产久久久一区二区三区| 久久亚洲精品不卡| 色尼玛亚洲综合影院| 男人舔女人下体高潮全视频| 亚洲人成网站在线播| 婷婷色综合大香蕉| 欧美高清性xxxxhd video| 一区福利在线观看| 日韩欧美国产在线观看| 亚洲18禁久久av| 最近最新中文字幕大全电影3| 亚洲av.av天堂| 中文字幕人妻熟人妻熟丝袜美| 国产三级中文精品| 国产成人影院久久av| 女人十人毛片免费观看3o分钟| 精品久久久久久久久久免费视频| 两个人的视频大全免费| 国产黄色小视频在线观看| 日韩欧美在线二视频| 91字幕亚洲| 亚洲国产精品成人综合色| 免费人成在线观看视频色| 老熟妇仑乱视频hdxx| 天美传媒精品一区二区| 国产aⅴ精品一区二区三区波| 九九在线视频观看精品| 国产精品爽爽va在线观看网站| 日本在线视频免费播放| 老熟妇仑乱视频hdxx| 亚洲aⅴ乱码一区二区在线播放| 99精品在免费线老司机午夜| 少妇被粗大猛烈的视频| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利欧美成人| 精品久久久久久久久久久久久| 亚洲av.av天堂| 国产色爽女视频免费观看| 深夜精品福利| 看片在线看免费视频| 丰满乱子伦码专区| 国产成+人综合+亚洲专区| 51午夜福利影视在线观看| 久久久久国产精品人妻aⅴ院| 国产高清有码在线观看视频| 亚洲av五月六月丁香网| 高清毛片免费观看视频网站| 欧美日韩黄片免| 麻豆成人午夜福利视频| 波多野结衣高清无吗| 久久99热6这里只有精品| 无遮挡黄片免费观看| 色综合欧美亚洲国产小说| 亚洲片人在线观看| 十八禁国产超污无遮挡网站| 高清日韩中文字幕在线| 男人狂女人下面高潮的视频| 看十八女毛片水多多多| 亚洲经典国产精华液单 | 国产成人啪精品午夜网站| 我的老师免费观看完整版| 深爱激情五月婷婷| 天堂网av新在线| 欧美在线黄色| 欧美色视频一区免费| 亚洲国产日韩欧美精品在线观看| 99热只有精品国产| 校园春色视频在线观看| 宅男免费午夜| 国产伦精品一区二区三区四那| a级毛片免费高清观看在线播放| 十八禁国产超污无遮挡网站| 草草在线视频免费看| 成熟少妇高潮喷水视频| 午夜福利18| 一级黄片播放器| 变态另类成人亚洲欧美熟女| a级毛片免费高清观看在线播放| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩无卡精品| 成人性生交大片免费视频hd| 麻豆成人av在线观看| а√天堂www在线а√下载| 露出奶头的视频| 日本精品一区二区三区蜜桃| 国产单亲对白刺激| 国产男靠女视频免费网站| 99精品久久久久人妻精品| 亚洲中文字幕日韩| 久久久久久久精品吃奶| netflix在线观看网站| 亚洲熟妇熟女久久| 国产又黄又爽又无遮挡在线| 久久精品国产自在天天线| 久久99热6这里只有精品| 亚洲 国产 在线| 成人av在线播放网站| 日本黄色视频三级网站网址| 亚洲欧美日韩高清专用| 自拍偷自拍亚洲精品老妇| 欧美另类亚洲清纯唯美| 久久久久久久久久黄片| 亚洲欧美日韩高清在线视频| 国产欧美日韩一区二区三| 欧美高清性xxxxhd video| 午夜福利欧美成人| 国产精品久久久久久久久免 | 国产精品av视频在线免费观看| 国产单亲对白刺激| 婷婷精品国产亚洲av在线| 精品熟女少妇八av免费久了| 国内精品美女久久久久久| 国产精品久久视频播放| 国产白丝娇喘喷水9色精品| 亚洲七黄色美女视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品爽爽va在线观看网站| 国产精品电影一区二区三区| 亚洲自偷自拍三级| 乱码一卡2卡4卡精品| 国产蜜桃级精品一区二区三区| 亚洲av电影不卡..在线观看| 精华霜和精华液先用哪个| 性欧美人与动物交配| 久久这里只有精品中国| 99热精品在线国产| 一夜夜www| 熟女电影av网| 一级黄片播放器| 午夜a级毛片| 人妻丰满熟妇av一区二区三区| 欧美性猛交╳xxx乱大交人| 国产黄片美女视频| 欧美日韩福利视频一区二区| www日本黄色视频网| 九色国产91popny在线| 午夜福利免费观看在线| 国产黄片美女视频| 欧美日韩福利视频一区二区| 精品久久久久久久久亚洲 | 国产一级毛片七仙女欲春2| aaaaa片日本免费| 国内揄拍国产精品人妻在线| 99精品在免费线老司机午夜| 欧美成人性av电影在线观看| 精品一区二区三区人妻视频| 首页视频小说图片口味搜索| 97超视频在线观看视频| 91狼人影院| 性色avwww在线观看| 97碰自拍视频| 直男gayav资源| 亚洲内射少妇av| 在线十欧美十亚洲十日本专区| 青草久久国产| 日韩欧美精品v在线| 国产 一区 欧美 日韩| 亚洲无线在线观看| 黄色一级大片看看| 男人和女人高潮做爰伦理| 亚洲成人精品中文字幕电影| 十八禁国产超污无遮挡网站| 我要搜黄色片| 亚洲精品在线美女| 国产成人欧美在线观看| 日本精品一区二区三区蜜桃| 国产免费男女视频| 禁无遮挡网站| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区| 欧美bdsm另类| 中文在线观看免费www的网站| 亚洲18禁久久av| 天堂网av新在线| 色哟哟·www| 在线播放无遮挡| 亚洲av第一区精品v没综合| 婷婷亚洲欧美| 一区二区三区免费毛片| 日韩中文字幕欧美一区二区| 免费看美女性在线毛片视频| 91狼人影院| 两个人视频免费观看高清| 免费无遮挡裸体视频| 男女视频在线观看网站免费| 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看| 永久网站在线| 欧美日韩福利视频一区二区| 人人妻,人人澡人人爽秒播| 特大巨黑吊av在线直播| 99在线人妻在线中文字幕| 国产精品一区二区三区四区免费观看 | 91麻豆av在线| 婷婷六月久久综合丁香| 亚洲天堂国产精品一区在线| 美女黄网站色视频| 高潮久久久久久久久久久不卡| 少妇的逼好多水| 人人妻人人看人人澡| 欧美一区二区亚洲| 成人高潮视频无遮挡免费网站| 男女视频在线观看网站免费| 中文亚洲av片在线观看爽| 亚洲国产色片| av中文乱码字幕在线| 亚洲精品456在线播放app | 一本久久中文字幕| 久久国产乱子伦精品免费另类| 成人美女网站在线观看视频| 一区二区三区四区激情视频 | 性色avwww在线观看| 国产欧美日韩精品一区二区| 国内少妇人妻偷人精品xxx网站| 欧美日韩福利视频一区二区| 亚洲专区中文字幕在线| 很黄的视频免费| 亚洲自拍偷在线| 国产成人啪精品午夜网站| 88av欧美| 9191精品国产免费久久| 亚洲国产精品成人综合色| 久久精品91蜜桃| 精品熟女少妇八av免费久了| 1000部很黄的大片| 最近在线观看免费完整版| 国内精品美女久久久久久| 亚洲成av人片免费观看| 夜夜看夜夜爽夜夜摸| 十八禁人妻一区二区| 在线观看66精品国产| 男人的好看免费观看在线视频| 亚洲精品亚洲一区二区| 成年女人毛片免费观看观看9| 久久久成人免费电影| 淫妇啪啪啪对白视频| 国产精品精品国产色婷婷| 一个人看的www免费观看视频| 18禁裸乳无遮挡免费网站照片| 成年女人看的毛片在线观看| 亚洲人与动物交配视频| 国产欧美日韩一区二区三| 老女人水多毛片| 欧美精品国产亚洲| 国产探花在线观看一区二区| 国产一区二区在线av高清观看| 可以在线观看毛片的网站| 国产一区二区激情短视频| 国产av麻豆久久久久久久| 久久久久久大精品| 搡老熟女国产l中国老女人| 精品免费久久久久久久清纯| 一级黄色大片毛片| 国产成人aa在线观看| av国产免费在线观看| 国产又黄又爽又无遮挡在线|