艾滋病療法的評(píng)價(jià)及療效的預(yù)測(cè)1.doc

艾滋病療法的評(píng)價(jià)及療效的預(yù)測(cè)摘要：本文利用灰色預(yù)測(cè)模型預(yù)測(cè)了艾滋病繼續(xù)治療的效果，并確定了最佳治療終止時(shí)間對(duì)問(wèn)題中的各種療法，按年齡進(jìn)行分類，對(duì)每一類人群評(píng)價(jià)了4種療法的優(yōu)劣，并對(duì)較優(yōu)的療法預(yù)測(cè)繼續(xù)治療的效果，確定最佳治療終止時(shí)間最后考慮每一種療法的費(fèi)用，用層次分析法對(duì)各種療法再一次進(jìn)行評(píng)價(jià)，得到病人最滿意的療法關(guān)鍵詞:灰色預(yù)測(cè)模型 層次分析法 HIV的濃度 CD4的值 殘差檢驗(yàn)1 問(wèn)題重述（略）2 模型假設(shè)1）對(duì)于題中附件1里的數(shù)據(jù)我們假設(shè)如下：將其按測(cè)試的時(shí)間排序，其中在同一時(shí)間段測(cè)試結(jié)果數(shù)量少于30個(gè)的數(shù)據(jù)，我們忽略不計(jì)（例如，對(duì)于第六周測(cè)試時(shí)只有11個(gè)人，我們將把這十一個(gè)數(shù)據(jù)忽略）；2）假設(shè)題中附件2里同一時(shí)間里的數(shù)據(jù)過(guò)少時(shí)，我們?cè)谟?jì)算時(shí)也將這些數(shù)據(jù)忽略不計(jì)；3）假設(shè)影響病人治療的只有費(fèi)用和療效兩個(gè)因素，不考慮其它因素；3 符號(hào)說(shuō)明1、為原始序列2、為生成累加序列3、為相對(duì)殘差序列4、為平均相對(duì)殘差 5、a稱為發(fā)展系數(shù)6、u稱為灰色作用量7、為時(shí)間響應(yīng)序列8、 為待估參數(shù)向量9、(0)為絕對(duì)殘差序列10、CD4（Date）為測(cè)試CD4的時(shí)刻（周）11、Log(CD4 count+1) 為測(cè)得的CD4濃度12、為最大特征值13、CI為檢驗(yàn)判斷矩陣一致性的指標(biāo)14、CR為一致性比例15、RI為修正系數(shù)16、和為本層次的要素和（i,j=1,2,n）17、為上一層次要素4 問(wèn)題分析問(wèn)題1是依據(jù)同時(shí)服用zidovudine，lamivudine和indinavir這3種藥物的300多名病人每隔幾周測(cè)試的CD4和HIV的濃度，預(yù)測(cè)繼續(xù)治療的效果，或者確定最佳治療終止時(shí)間.我們運(yùn)用灰色預(yù)測(cè)模型GM（1，1）模型來(lái)進(jìn)行假設(shè)，需要對(duì)數(shù)據(jù)進(jìn)行分析，通過(guò)灰色模型和相應(yīng)的計(jì)算來(lái)預(yù)測(cè)未來(lái)時(shí)間里CD4和HIV的濃度及相關(guān)數(shù)據(jù). 問(wèn)題2是將1300多名病人分為4組，每組按4種療法中的一種服藥，大約每隔8周測(cè)試一次CD4的濃度.最后評(píng)價(jià)4種療法的優(yōu)劣，并對(duì)較優(yōu)的療法預(yù)測(cè)繼續(xù)治療的效果，或者確定最佳治療終止時(shí)間.我們運(yùn)用灰色預(yù)測(cè)模型，來(lái)預(yù)測(cè)4種療法的治療效果，即我們要通過(guò)數(shù)據(jù)找到規(guī)律，通過(guò)相應(yīng)的計(jì)算（運(yùn)用EXCEL軟件）來(lái)預(yù)測(cè)未來(lái)時(shí)間里4種療法在不同年齡段CD4的濃度及相關(guān)數(shù)據(jù)，找到哪一種療法療效較好，與問(wèn)題1的方法相同，運(yùn)用灰色模型預(yù)測(cè)出繼續(xù)治療的效果. 問(wèn)題3中病人需要考慮4種療法的費(fèi)用，因個(gè)人經(jīng)濟(jì)狀況的不同，這4種療法費(fèi)用的多少可能會(huì)影響到病人選擇哪一種療法，對(duì)于此問(wèn)題，我們需要將這4種療法預(yù)測(cè)的治療效果和其所需費(fèi)用綜合考慮，運(yùn)用層次分析法進(jìn)行分析，決定采取何種療法.5 模型的建立與求解1、 問(wèn)題1：首先，我們通過(guò)對(duì)題中所給的附件1里的數(shù)據(jù)進(jìn)行分類排序分析.對(duì)CD4和HIV的數(shù)據(jù)按照時(shí)間（周）從小到大進(jìn)行排序，生成有較強(qiáng)規(guī)律性的數(shù)據(jù)序列，對(duì)于每個(gè)時(shí)間點(diǎn)CD4的數(shù)據(jù)，我們分別求其平均值，并對(duì)于人數(shù)較多的時(shí)間對(duì)應(yīng)的CD4的值，從中挑出來(lái)，作成圖形圖1（CD4濃度隨時(shí)間的變化圖形） 圖2（HIV濃度隨時(shí)間的變化圖形） 對(duì)于這兩個(gè)圖，我們?nèi)∷鼈兊漠惓Ｖ担磳?duì)于圖1我們?nèi)∑錁O大值，對(duì)于圖2我們?nèi)O小值.首先，用GM（1，1）模型對(duì)檢測(cè)CD4的極大值點(diǎn)進(jìn)行建模預(yù)測(cè)，分別得到第3周，第9周，第23周，第39周的原始數(shù)據(jù)，其相應(yīng)的和如下：表1K1234CD4 （Date）周392339原始數(shù)據(jù) 392339累加后數(shù)據(jù)3123547根據(jù)表1，可知=3，9，23，39， 構(gòu)造累加生成數(shù)列=3，12，35，74，構(gòu)造數(shù)據(jù)矩陣B和數(shù)據(jù)向量Y計(jì)算 由Matlab計(jì)算得 得a=-0.6217 ;u=5.9478，然后計(jì)算得 u/a= -9.5670將數(shù)據(jù)代入（2）式，得 對(duì)于CD4的預(yù)測(cè)，我們代入數(shù)據(jù)分析：當(dāng) k=4時(shí)67.5228-39（39周的序號(hào)）= 28.5228即，從開始用藥到第28.5228+40=68.5228周時(shí)停止用藥最好.另外取CD4的極小值點(diǎn)和CD4的穩(wěn)定點(diǎn)，對(duì)它們分別進(jìn)行類似于對(duì)CD4的極大值的處理方法，可以得到：CD4出現(xiàn)極小值點(diǎn)在測(cè)試用藥的40周后的第23.850周及出現(xiàn)穩(wěn)定點(diǎn)在在測(cè)試用藥的40周后的第44.55周.接下來(lái)，用GM（1，1）模型對(duì)檢測(cè)HIV的極小值點(diǎn)進(jìn)行建模預(yù)測(cè)，分別得到第3周，第9周，第25周，第40周的原始數(shù)據(jù)，其相應(yīng)的和如下：表4K1234HIV（Date）周392540原始數(shù)據(jù)392540累加后數(shù)據(jù)3123777根據(jù)上表4，按照GM(1,1) 求解過(guò)程得到a=-0.6044 ;u=6.7361，計(jì)算得 u/a= -11.1451將數(shù)據(jù)代入，得 對(duì)于HIV的預(yù)測(cè)，我們代入數(shù)據(jù)分析：當(dāng) k=4時(shí)70.5474-40（40周的序號(hào)）=30.5474即，從開始用藥到第30.5474+40=70.5474周時(shí)停止用藥最好.取HIV的極大值點(diǎn)和HIV的穩(wěn)定點(diǎn)，對(duì)它們進(jìn)行類似于對(duì)HIV的極小值的處理方法，可以得到：HIV出現(xiàn)極大值點(diǎn)在測(cè)試用藥的40周后的第23周及出現(xiàn)穩(wěn)定點(diǎn)在測(cè)試用藥的40周后的第20周.綜合上面對(duì)問(wèn)題一CD4和HIV的預(yù)測(cè)，我們得出結(jié)論：在測(cè)試用藥的40周后的大約第28周時(shí)停止用藥.2、 問(wèn)題二：對(duì)于附錄二的數(shù)據(jù)，我們先對(duì)其按療法排序，再在每種療法中按年齡排序，分別把各療法的數(shù)據(jù)分成三個(gè)階段：青年（15歲30歲），中年（30歲50歲），老年（50歲70歲）；首先繪出青年（15歲30歲）的各種療法效果的曲線圖；圖3（療法1中CD4的曲線圖） 圖4（療法2中CD4的曲線圖）圖5（療法3中CD4的曲線圖） 圖6（療法4中CD4的曲線圖）對(duì)于15歲到30歲，我們可由圖看出療法效果由好到差排序?yàn)椋?，3，2，1仿照上面相同的辦法同樣得到中年（30歲50歲）各種療法效果的曲線圖，通過(guò)比較這些曲線圖，可以看出療法效果由好到差排序?yàn)椋?，3，2，1同樣對(duì)于50歲到70歲，我們也可由圖看出療法效果由好到差排序?yàn)椋?，3，2，1；綜上，我們?cè)u(píng)價(jià)如下，不論在任何年齡段，4種療法總體療效由好到差依次為：4，3，2，1，即第4種療法最優(yōu)下面，我們確定其最佳治療終止時(shí)間.對(duì)療法4的數(shù)據(jù)不分年齡，按照時(shí)間進(jìn)行排序，對(duì)各個(gè)時(shí)間點(diǎn)對(duì)應(yīng)的CD4的數(shù)據(jù)求平均值，如果該時(shí)間點(diǎn)所對(duì)應(yīng)的數(shù)據(jù)較多，把該時(shí)間點(diǎn)和所對(duì)應(yīng)的CD4的值取出來(lái)得到下表及所對(duì)應(yīng)的圖形：表8（療法4中Log(CD4 count+1)的均值）時(shí)間Log(CD4count+1)02.83564973.02882383.24883793.317596153.301195163.299549173.360405233.024497243.064096263.194985303.176694312.893515323.0544332.919191342.783876392.736491402.957543用GM（1，1）模型對(duì)檢測(cè)CD4的濃度進(jìn)行建模預(yù)測(cè)，分別得到第9周，第17周，第24周，第32周的原始數(shù)據(jù)，其相應(yīng)的和如下：表9K1234CD4（Date）周9172432原始數(shù)據(jù)X（0）9172432累加后數(shù)據(jù)X（1）9265082根據(jù)上表9，可知=9，17，24，32， 構(gòu)造累加生成數(shù)列 =9，26，50，82，構(gòu)造數(shù)據(jù)矩陣B2和數(shù)據(jù)向量Y2計(jì)算 由Matlab軟件，計(jì)算得：得a=-0.3079 ;u=11.8639，計(jì)算得 u/a= -38.5317將數(shù)據(jù)代入（1）式，得通過(guò)分析，我們得出療法4最優(yōu).對(duì)于CD4的預(yù)測(cè)，我們代入數(shù)據(jù)分析：當(dāng) k=4時(shí)42.3456-32= 10.3456即，