Spelling suggestions: "subject:"gibbs samplings"" "subject:"hibbs samplings""
1 |
競爭風險下長期存活資料之貝氏分析 / Bayesian analysis for long-term survival data蔡佳蓉 Unknown Date (has links)
當造成失敗的原因不只一種時,若各對象同一時間最多只經歷一種失敗原因,則這些失敗原因稱為競爭風險。然而,有些個體不會失敗或者經過治療之後已痊癒,我們稱這部分的群體為治癒群。本文考慮同時處理競爭風險及治癒率的混合模式,即競爭風險的治癒率模式,亦將解釋變數結合到治癒率、競爭風險的條件失敗機率,或未治癒下競爭風險的條件存活函數中,並以建立在完整資料上之擴充的概似函數為貝氏分析的架構。對於右設限對象則以插補方式決定是否會治癒或會因何種風險而失敗,並推導各參數的完全條件後驗分配及其性質。由於邊際後驗分配的數學形式無法明確呈現,再加上需對右設限者判斷其狀態,所以採用屬於馬可夫鏈蒙地卡羅法的Gibbs抽樣法及適應性拒絕抽樣法(adaptive rejection sampling) ,執行參數之模擬抽樣及設算右設限者之治癒或失敗狀態。實證部分,我們分析Klein and Moeschberger (1997)書中骨髓移植後的血癌病患的資料,並用不同模式之下的參數模擬值計算各對象之條件預測指標(CPO),換算成各模式的對數擬邊際概似函數值(LPML),比較不同模式的優劣。 / In case that there are more than one possible failure types, if each subject experiences at most one failure type at one time, then these failure types are called competing risks. Moreover, some subjects have been cured or are immune so they never fail, then they are called the cured ones. This dissertation discusses several mixture models containing competing risks and cure rate. Furthermore, covariates are associated with cure rate, conditional failure rate of each risk, or conditional survival function of each risk, and we propose the Bayesian procedure based on the augmented likelihood function of complete data. For right censored subjects, we make use of imputation to determine whether they were cured or failed by which risk and derive full conditional posterior distributions. Since all marginal posterior distributions don’t have closed forms and right censored subjects need to be identified their statuses, we take Gibbs sampling and adaptive rejection sampling of Markov chain Monte Carlo method to simulate parameter values. We illustrate how to conduct Bayesian analysis by using the bone marrow transplant data from the book written by Klein and Moeschberger (1997). To do model selection, we compute the conditional predictive ordinate(CPO) for every subject under each model, then the goodness is determined by the comparing the value of log of pseudo marginal likelihood (LMPL) of each model.
|
Page generated in 0.0715 seconds