福建農(nóng)林大學(xué)-數(shù)學(xué)建模培訓(xùn)-mcm2006B艾滋病療法的評(píng)價(jià).ppt_第1頁(yè)
福建農(nóng)林大學(xué)-數(shù)學(xué)建模培訓(xùn)-mcm2006B艾滋病療法的評(píng)價(jià).ppt_第2頁(yè)
福建農(nóng)林大學(xué)-數(shù)學(xué)建模培訓(xùn)-mcm2006B艾滋病療法的評(píng)價(jià).ppt_第3頁(yè)
福建農(nóng)林大學(xué)-數(shù)學(xué)建模培訓(xùn)-mcm2006B艾滋病療法的評(píng)價(jià).ppt_第4頁(yè)
福建農(nóng)林大學(xué)-數(shù)學(xué)建模培訓(xùn)-mcm2006B艾滋病療法的評(píng)價(jià).ppt_第5頁(yè)
已閱讀5頁(yè),還剩18頁(yè)未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

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é)建模題目

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論