福建農(nóng)林大學(xué)-數(shù)學(xué)建模培訓(xùn)-mcm2006B艾滋病療法的評(píng)價(jià).ppt

1、艾滋病療法的評(píng)價(jià) 及療效的預(yù)測(cè),全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽 2006年B題,艾滋病是當(dāng)前人類社會(huì)最嚴(yán)重的瘟疫之一，從1981年發(fā)現(xiàn)以來(lái)的20多年間，它已經(jīng)吞噬了近3000萬(wàn)人的生命。 艾滋病的醫(yī)學(xué)全名為“獲得性免疫缺損綜合癥”，英文簡(jiǎn)稱AIDS，它是由艾滋病毒（醫(yī)學(xué)全名為“人體免疫缺損病毒”, 英文簡(jiǎn)稱HIV）引起的。這種病毒破壞人的免疫系統(tǒng)，使人體喪失抵抗各種疾病的能力，從而嚴(yán)重危害人的生命。人類免疫系統(tǒng)的CD4細(xì)胞在抵御HIV的入侵中起著重要作用，當(dāng)CD4被HIV感染而裂解時(shí)，其數(shù)量會(huì)急劇減少，HIV將迅速增加，導(dǎo)致AIDS發(fā)作。 艾滋病治療的目的，是盡量減少人體內(nèi)HIV的數(shù)量，同時(shí)產(chǎn)生更多的

2、CD4，至少要有效地降低CD4減少的速度，以提高人體免疫能力。,迄今為止人類還沒(méi)有找到能根治AIDS的療法，目前的一些AIDS療法不僅對(duì)人體有副作用，而且成本也很高。許多國(guó)家和醫(yī)療組織都在積極試驗(yàn)、尋找更好的AIDS療法。 現(xiàn)在得到了美國(guó)艾滋病醫(yī)療試驗(yàn)機(jī)構(gòu)ACTG公布的兩組數(shù)據(jù)。 ACTG320（見附件1）是同時(shí)服用zid., lam., ind. 3種藥物的300多名病人每隔幾周測(cè)試的CD4和HIV的濃度（每毫升血液里的數(shù)量）。193A（見附件2）是將1300多名病人隨機(jī)地分為4組，每組按下述4種療法中的一種服藥，大約每隔8周測(cè)試的CD4濃度（這組數(shù)據(jù)缺HIV濃度，它的測(cè)試成本很高）。4種療

3、法的日用藥分別為：600mg zid.或400mg did.，這兩種藥按月輪換使用；600 mg zid.加2.25 mg zal.；600 mg zid.加400 mg did.；600 mg zid.加400 mg did.，再加400 mg nev.。,請(qǐng)你完成以下問(wèn)題： （1）利用附件1的數(shù)據(jù)，預(yù)測(cè)繼續(xù)治療的效果，或者確定最佳治療終止時(shí)間（繼續(xù)治療指在測(cè)試終止后繼續(xù)服藥，如果認(rèn)為繼續(xù)服藥效果不好，則可選擇提前終止治療）。 （2）利用附件2的數(shù)據(jù)，評(píng)價(jià)4種療法的優(yōu)劣（僅以CD4為標(biāo)準(zhǔn)），并對(duì)較優(yōu)的療法預(yù)測(cè)繼續(xù)治療的效果，或者確定最佳治療終止時(shí)間。 (3) 艾滋病藥品的主要供給商對(duì)不發(fā)達(dá)國(guó)

4、家提供的藥品價(jià)格如下：600mg zid. 1.60美元，400mg did. 0.85美元，2.25 mg zal. 1.85美元，400 mg nev. 1.20美元。如果病人需要考慮4種療法的費(fèi)用，對(duì)（2）中的評(píng)價(jià)和預(yù)測(cè)（或者提前終止）有什么改變。,附件1 ACTG320數(shù)據(jù) 同時(shí)服用3種藥物的300多名病人每隔幾周測(cè)試的CD4和HIV的濃度。 第1列是病人編號(hào)，第2列是測(cè)試CD4的時(shí)刻（周），第3列是測(cè)得的CD4（乘以0.2個(gè)/ml），第4列是測(cè)試HIV的時(shí)刻（周），第5列是測(cè)得的HIV（單位不詳）。 PtIDCD4Date CD4Count RNADateVLoad 23424 01

5、780 5.5 23424 42284 3.9 23424 81268 4.7 23424 25171 25 4 23424 40 99 40 5 23425 0 140 5.3 23425 4 624 2.4 23425 91109 3.7 23425 23122 23 2.6 23425 40320,附件2 193A數(shù)據(jù) 1300多名病人按照4種療法服藥大約每隔8周測(cè)試的CD4濃度。 第1列是病人編號(hào)，第2列是4種療法的代碼： 第3列是病人年齡，第4列是測(cè)試CD4的時(shí)刻（周），第5列是測(cè)得的CD4，取值log(CD4+1). ID 療法 年齡 時(shí)間 Log(CD4 count+1) 123

6、6.427103.1355 1236.42717.57143.0445 1236.427115.57142.7726 1236.427123.57142.8332 1236.427132.57143.2189 1236.427140 3.0445 2447.84670 3.0681 2447.84678 3.8918 2447.846716 3.9703 2447.846723 3.6109 2447.846730.71433.3322 2447.846739 3.0910,CD4大致有先增后減的趨勢(shì)，HIV有先減后增的趨勢(shì)，啟示應(yīng)建立時(shí)間的二次函數(shù)模型,問(wèn)題（1） 利用附件1的數(shù)據(jù)，預(yù)測(cè)繼續(xù)

7、治療的效果, 或者確定最佳治療終止時(shí)間。,分析數(shù)據(jù),如隨機(jī)取20個(gè)病人，畫出他們CD4和HIV濃度隨時(shí)間變化的圖形（折線），,若先用一次模型，應(yīng)與二次模型做統(tǒng)計(jì)分析比較,CD4,HIV,1） 總體回歸模型 用全部數(shù)據(jù)擬合一個(gè)模型，如 yij=b0+b1tij+b2tij2，tij，yij分別為第i病人第j次測(cè)量的時(shí)間和CD4（或HIV），估計(jì)b0, b1, b2,建立模型,有以下形式的回歸模型（一次與二次模型比較，二次較優(yōu)）：,對(duì)HIV，b20, b10, t=-b1/2b2 達(dá)到最小,對(duì)CD4，b20, t=-b1/2b2 達(dá)到最大,平均地應(yīng)在2530（周）結(jié)束治療。,可對(duì)CD4統(tǒng)計(jì)b2i0

8、（存在正最大點(diǎn)）及b2i0（不存在最大點(diǎn)）的頻率，分別作為及時(shí)結(jié)束治療與繼續(xù)治療的概率；也可用它們的均值和均方差在確定分布下計(jì)算這些概率。,2) 個(gè)人回歸模型,用每個(gè)病人的數(shù)據(jù)擬合一個(gè)模型，如yij=b0i+b1itij+b2itij2，計(jì)算b0i ,b1i ,b2i的均值和均方差，用均值可得CD4的最大點(diǎn)和HIV的最小點(diǎn)。,將模型 yij=b0i+b1itij+b2itij2 的系數(shù)分解: bki=bk+ik（k=0,1,2），設(shè)ik服從零均值正態(tài)分布，利用極大似然估計(jì)可以得到bk的估計(jì)值及ik的方差。,3) 縱向數(shù)據(jù)回歸模型,bk的估計(jì)值用于計(jì)算CD4的最大點(diǎn)和HIV的最小點(diǎn)，它與ik的方

9、差一起可以在正態(tài)分布假定下計(jì)算及時(shí)結(jié)束治療與繼續(xù)治療的概率。,如果考慮病人初始狀態(tài)（t=0時(shí)的CD4和HIV）的不同對(duì)模型的影響，可以將模型中的yij定義為第i病人第j次測(cè)量的CD4（或HIV）與初始值之差或之比。,或者先按照病人初始狀態(tài)分類（如輕度、中度、重度），然后對(duì)于每一類建立回歸模型。,建立模型,問(wèn)題（2） 利用附件2的數(shù)據(jù)，評(píng)價(jià)4種療法的優(yōu)劣，并對(duì)較優(yōu)的療法預(yù)測(cè)繼續(xù)治療的效果，或者確定最佳治療終止時(shí)間。,對(duì)于每種療法隨機(jī)取20個(gè)病人，畫出他們CD4隨時(shí)間變化的圖形（折線），可以看出療法13的CD4基本上水平，略有下降，而療法4有先增后減的趨勢(shì)。啟示應(yīng)建立時(shí)間的一次與二次函數(shù)模型，經(jīng)統(tǒng)

10、計(jì)分析比較，確定哪種較優(yōu)。,分析數(shù)據(jù),可以引入4（或3）個(gè)0-1變量表示4種療法建立統(tǒng)一模型，或者對(duì)每種療法各建立一個(gè)模型（一般來(lái)說(shuō)前者較優(yōu)）。,1）回歸模型,以總體回歸模型為例，只需增加年齡變量，分別用一次與二次時(shí)間函數(shù)模型進(jìn)行比較，可知療法13用一次模型較優(yōu)，且一次項(xiàng)系數(shù)為負(fù)，即CD4在減少，從數(shù)值看療法3優(yōu)于療法2和1；療法4用二次模型較優(yōu)，即CD4先增后減，,建立模型,做療法有無(wú)顯著性差異的兩兩比較： 用1個(gè)0-1變量構(gòu)造兩種療法的統(tǒng)一模型，可以用t檢驗(yàn)作回歸系數(shù)是否為零的假設(shè)檢驗(yàn)。結(jié)果是療法1與2無(wú)顯著性差異，而療法1與3，2與3，3與4均有顯著性差異。,2）假設(shè)檢驗(yàn),療法是必選的因

11、素，如果還考慮年齡和初始狀態(tài)，用雙因素分析則比較復(fù)雜，可以先按這些因素分類，再做療法的單因素分析。,3）方差分析,論文中出現(xiàn)問(wèn)題的評(píng)析,1只做數(shù)據(jù)擬合，不做統(tǒng)計(jì)檢驗(yàn)。,用回歸分析方法做統(tǒng)計(jì)檢驗(yàn)： 得到的模型有無(wú)顯著意義，它的置信度多大，用它作預(yù)測(cè)時(shí)準(zhǔn)確程度如何。,最小二乘擬合： 已知一組數(shù)據(jù)(xk,yk), k=1,2,n 和一個(gè)形式已定、參數(shù)待定的模型（曲線）y=f(x, ), 確定參數(shù)使數(shù)據(jù)與模型在誤差平方和最小的意義下擬合得最好。,有現(xiàn)成的計(jì)算公式和軟件確定參數(shù),y=a1x+b1=0.2013x+1.0040,z=a2x+b2=0.2025x+1.0025,兩個(gè)模型的結(jié)果一樣，但是可靠性

12、和準(zhǔn)確性一樣嗎？,需要做統(tǒng)計(jì)檢驗(yàn)！,簡(jiǎn)例,最小二乘擬合,R2=0.3528在z的變化中只有35%是由x決定的,p值大于0.05在常用的置信度=5%下該模型沒(méi)有顯著意義,a2的置信區(qū)間包含零點(diǎn)a=0落在置信水平為95%的區(qū)間內(nèi),若x=5.5,預(yù)測(cè)區(qū)間 y：2.0621，2.1606 ，z：0.3899，3.8432，,預(yù)測(cè)值 y=2.1113, z=2.1165,回歸分析結(jié)果,z的預(yù)測(cè)沒(méi)有什么價(jià)值！,2對(duì)原始數(shù)據(jù)先取平均，再用平均值做擬合,數(shù)據(jù)中大多數(shù)病人的測(cè)試時(shí)間是0,4,8,（周），個(gè)別病人是5,7, （周），為了得到CD4(或HIV)的變化趨勢(shì)，先在每個(gè)時(shí)間點(diǎn)上對(duì)CD4取平均，再用平均值做

13、擬合,這樣做有什么問(wèn)題呢？,注意：有的時(shí)間點(diǎn)上有上百個(gè)病人的CD4，而有的時(shí)間點(diǎn)上只有幾個(gè)病人的CD4。,數(shù)據(jù)(x,y)：x=0, y=0:0.05:2; x=1,y=1:0.05:3; x=2,y=0,2, 在圖上共84個(gè)點(diǎn)。,簡(jiǎn)例,對(duì)數(shù)據(jù)(x,y)擬合一條直線，得y=0.7647x+1.0784,先在x=0,1,2對(duì)y取平均得到3個(gè)點(diǎn), 再擬合直線，得y=1.3333,取平均后的3個(gè)點(diǎn)在做擬合時(shí)的權(quán)重相同; 但是x=0,1的點(diǎn)都是41個(gè)數(shù)據(jù)的平均值, 而x=3的點(diǎn)只是2個(gè)數(shù)據(jù)的平均值,3擬合過(guò)度,用4次、5次甚至更高次數(shù)的多項(xiàng)式作擬合,數(shù)據(jù)的起伏是由與時(shí)間本身無(wú)關(guān)的其它隨機(jī)因素引起的，不應(yīng)該用增加時(shí)間的高次項(xiàng)來(lái)擬合,注意擬合與插值在應(yīng)用場(chǎng)合的差別,一般地說(shuō)，多項(xiàng)式擬合不要超過(guò)3次，對(duì)于本題如果用統(tǒng)計(jì)檢驗(yàn)做顯著性分析，可以發(fā)現(xiàn)多數(shù)情況是以時(shí)間的2次多項(xiàng)式為好,適用本題的顯然是擬合而非插值,4只按照?qǐng)D形做直觀的定性判斷，不做定量分析,為了比較4種療法的優(yōu)劣，對(duì)每種療法的CD4（或取平均）做散點(diǎn)圖, 或折線（散點(diǎn)連線）圖, 直觀地比較療法的優(yōu)劣,沒(méi)有普遍的指導(dǎo)意義，對(duì)于其它問(wèn)題或數(shù)據(jù)，這樣做不一定能得到結(jié)果,一種療法的CD4散點(diǎn)圖(或折線圖)在另一種療法的上面, 從統(tǒng)計(jì)意義上并不能表明其顯著性,作為數(shù)學(xué)建模題目