在本論文中，利用CD4 細胞數以及病毒承載量來預測愛滋病患的發病時間，並探討使用雞尾酒療法與沒有使用雞尾酒療法對愛滋病患的療效。此種包含長期追蹤共變數與存活時間的資料，常會因為長期追蹤資料的測量誤差或生物體本身的差異，以及共變數觀測值測量與存活有關時，導致推論產生偏差，因此，本研究利用聯合模型來解決此問題，在生物指標方面，使用線性隨機效應模型對長期追蹤資料做配適，並利用概似比檢定診斷長期追蹤模型的適合度；在事件時間方面，使用Cox比例風險模型描述共變數與存活時間之關係。結合這兩部分建構出多重長期追蹤資料的聯合概似模型且利用EM演算法對參數做估計。

In this thesis, we use AIDS patients’ CD4 counts and viral loads to predict their onset times and explore the curative effect whether patients were treated with HAART. Usually, the study data include longitudinal and survival time information, and, in general, result in inference bias due to the measurement errors on the longitudinal part, the differences among patients themselves, or the time-dependent covariates. Thus, we use the joint model to solve this problem. The approach uses a linear random effects model to characterize the longitudinal part and conducts the likelihood test to select a suitable longitudinal model, and utilizes the Cox proportional hazard model to describe the relationship between covariates and survival time information. Incorporated these two parts to build a multiple longitudinal data joint likelihood function of which EM algorithm is implemented to search for the maximum likelihood estimate.

施伊珊 (2007)。存活與長期追蹤資料之聯合模型－台灣愛滋病實例研究。國立中央大學統計研究所碩士論文。
孫怡婷 (2011)。聯合長期追蹤與存活資料分析─愛滋病病患之實例分析。國立中央大學統計研究所碩士論文。
陳瑞霙 (2008) 。聯合模型在雞尾酒療法療效評估之應用—利用CD4/CD8比值探討台灣愛滋病資料。國立中央大學統計研究所碩士論文。
張輔仁 (2008)。時間相依共變數之雙重存活時間分析—台灣愛滋病病患存活時間與 CD4 / CD8 比值關係之案例研究。國立中央大學統計研究所碩士論文。
Ciampi, A. and Etezadi-Amoli, J. (1985). A general model for testing the proportional hazards and the accelerated failure time hypothesis in the analysis of censored survival data with covariate. Communications in Statistics, 14, 651-667.
Cox, D. R. (1972). Regression Models and Life-Tables. Journal of the Royal Statistical Society, 34, 187-220.
Dubin, J. A., Müller, H. G. and Wang, J. L. (2001).Event history graphs for censored survival data. Statistics in Medicine, 20, 2951-2964.

Henderson, R., Diggle, P. and Dobson, A. (2000). Joint modeling of longitudinal measurements and event time data. Biostatistics, 4, 465-480.
Hsieh, F., Tseng, Y. K. and Wang, J. L. (2006). Joint Modeling of Survival and Longitudinal Data: Likelihood Approach Revisited. Biometrics, 62, 1037-1043.
Kaplan, E. L. and Meier, P. (1958). Nonparametric Estimation from Incomplete Observations. Journal of the American Statistical Association, 53, 457-481.
Kauermann G., Opsomer J. (2004). Generalized Cross-validation for Bandwidth Selection of Backfitting Estimates in Generalized Additive Models. Journal of Computational and Graphical Statistics, 13, 66-89.
Laird, N. M. and Ware, J. H. (1982). Random-effects models for longitudinal data. Biometrics, 38, 963-974.
Prentice, R. L. (1982). Covariate measurement errors and parameter estimation in a failure time regression model. Biometrika, 69, 331-342.
Tseng, Y. K., Hsieh , F. and Wang, J. L. (2005). Joint modeling of accelerated ailure time and longitudinal data. Biometrika, 92, 587-603.
Tsiatis, A. A. and Davidian, M. (2001). A semiparametric estimator for the proportional hazards model with longitudinal covariates measured with error. Biometrika, 88, 447-458.
Tsiatis, A. A. and Davidian, M. (2004). Joint Modeling of Longitudinal and Time-to-Event Data: An Overview. Statistica Sinica, 14, 809-834.
Verbeke, G. and Davidian, M. (2008). Joint Models for Longitudinal Data: Introduction and Overview. Longitudinal data analysis: handbooks of modern statistical methods Ed. Fitzmaurice, G., Davidian, M., Verbeke, G. and Molenberghs, G., 319-326. Chapman&Hall/CRC.
Wulfsohn, M. S. and Tsiatis, A. A. (1997). A Joint Model for Survival and Longitudinal Data Measured with Error. Biometrics, 53, 330-339.
Zeng D., Cai J. (2005). Asymptotic Results for Maximum Likelihood Estimators in Joint Analysis of Repeated Measurements and Survival Time. The Annals of Statistics, 33, 2132-2163.