The mean residual life model for the right-censored data in the presence of covariate measurement errors.

In this article, we consider the mean residual life regression model in the presence of covariate measurement errors. In the whole cohort, the surrogate variable of the error-prone covariate is available for each subject, while the instrumental variable (IV), which is related to the underlying true covariates, is measured only for some subjects, the calibration sample. Without specifying distributions of measurement errors but assuming that the IV is missing at random, we develop two estimation methods, the IV calibration and cohort estimators, for the regression parameters by solving estimation equations (EEs) based on the calibration sample and cohort sample, respectively. To improve estimation efficiency, a synthetic estimator is derived by applying the generalized method of moment for all EEs. The large sample properties of the proposed estimators are established and their finite sample performance are evaluated via simulation studies. Simulation results show that the cohort and synthetic estimators outperform the IV calibration estimator and the relative efficiency of the cohort and synthetic estimators mainly depends on the missing rate of IV. In the case of low missing rate, the synthetic estimator is more efficient than the cohort estimator, while the result can be reversed when the missing rate is high. We illustrate the proposed method by application to data from the patients with stage 5 chronic kidney disease in Taiwan.

Mean residual life cure models for right-censored data with and without length-biased sampling.

We propose a semiparametric mean residual life mixture cure model for right-censored survival data with a cured fraction. The model employs the proportional mean residual life model to describe the effects of covariates on the mean residual time of uncured subjects and the logistic regression model to describe the effects of covariates on the cure rate. We develop estimating equations to estimate the proposed cure model for the right-censored data with and without length-biased sampling, the latter is often found in prevalent cohort studies. In particular, we propose two estimating equations to estimate the effects of covariates in the cure rate and a method to combine them to improve the estimation efficiency. The consistency and asymptotic normality of the proposed estimates are established. The finite sample performance of the estimates is confirmed with simulations. The proposed estimation methods are applied to a clinical trial study on melanoma and a prevalent cohort study on early-onset type 2 diabetes mellitus.

Maximum likelihood estimation for length-biased and interval-censored data with a nonsusceptible fraction.

Left-truncated data are often encountered in epidemiological cohort studies, where individuals are recruited according to a certain cross-sectional sampling criterion. Length-biased data, a special case of left-truncated data, assume that the incidence of the initial event follows a homogeneous Poisson process. In this article, we consider an analysis of length-biased and interval-censored data with a nonsusceptible fraction. We first point out the importance of a well-defined target population, which depends on the prior knowledge for the support of the failure times of susceptible individuals. Given the target population, we proceed with a length-biased sampling and draw valid inferences from a length-biased sample. When there is no covariate, we show that it suffices to consider a discrete version of the survival function for the susceptible individuals with jump points at the left endpoints of the censoring intervals when maximizing the full likelihood function, and propose an EM algorithm to obtain the nonparametric maximum likelihood estimates of nonsusceptible rate and the survival function of the susceptible individuals. We also develop a novel graphical method for assessing the stationarity assumption. When covariates are present, we consider the Cox proportional hazards model for the survival time of the susceptible individuals and the logistic regression model for the probability of being susceptible. We construct the full likelihood function and obtain the nonparametric maximum likelihood estimates of the regression parameters by employing the EM algorithm. The large sample properties of the estimates are established. The performance of the method is assessed by simulations. The proposed model and method are applied to data from an early-onset diabetes mellitus study.

Trajectory of Estimated Glomerular Filtration Rate and Malnourishment Predict Mortality and Kidney Failure in Older Adults With Chronic Kidney Disease.

Objective: The trajectory patterns of estimated glomerular filtration rates (eGFR) in chronic kidney disease (CKD) older adults with malnourishment and their association with subsequent patient outcomes have not been elucidated. We aimed to assess the eGFR trajectory patterns for predicting patient survival and kidney failure in the elderly without or with malnourishment. Materials and Methods: Based on a prospective longitudinal cohort, CKD patients aged 65 years or older were enrolled from 2001 to 2013. Among the 3,948 patients whose eGFR trajectory patterns were analyzed, 1,872 patients were stratified by the absence or presence of malnourishment, and 765 patients were identified and categorized as having malnourishment. Four eGFR trajectory patterns [gradual decline (T0), early non-decline and then persistent decline (T1), persistent increase (T2), and low baseline and then progressive increase (T3)] were classified by utilizing a linear mixed-effect model with a quadratic term in time. The malnourishment was defined as body mass index < 22 kg/m2, serum albumin < 3.0 mg/dL, or Geriatric Nutritional Risk Index (GNRI) < 98. This study assessed the effectiveness of eGFR trajectory patterns in a median follow-up of 2.27 years for predicting all-cause mortality and kidney failure. Results: The mean age was 76.9 ± 6.7 years, and a total of 82 (10.7%) patients with malnourishment and 57 (5.1%) patients without malnourishment died at the end of the study. Compared with the reference trajectory T0, the overall mortality of T1 was markedly reduced [adjusted hazard ratio (aHR) = 0.52, 95% confidence interval (CI) 0.32-0.83]. In patients with trajectory, T3 was associated with a high risk for kidney failure (aHR = 5.68, 95% CI 3.12-10.4) compared with the reference, especially higher risk in the presence of malnourishment. Patients with high GNRI values were significantly associated with a lower risk of death and kidney failure, but patients with malnourishment and concomitant alcohol consumption had a higher risk of kidney failure. Conclusions: Low baseline eGFR and progressively increasing eGFR trajectory were high risks for kidney failure in CKD patients. These findings may be attributed to multimorbidity, malnourishment, and decompensation of renal function.

The TVGH-NYCU Thal-Classifier: Development of a Machine-Learning Classifier for Differentiating Thalassemia and Non-Thalassemia Patients.

Thalassemia and iron deficiency are the most common etiologies for microcytic anemia and there are indices discriminating both from common laboratory simple automatic counters. In this study a new classifier for discriminating thalassemia and non-thalassemia microcytic anemia was generated via combination of exciting indices with machine-learning techniques. A total of 350 Taiwanese adult patients whose anemia diagnosis, complete blood cell counts, and hemoglobin gene profiles were retrospectively reviewed. Thirteen prior established indices were applied to current cohort and the sensitivity, specificity, positive and negative predictive values were calculated. A support vector machine (SVM) with Monte-Carlo cross-validation procedure was adopted to generate the classifier. The performance of our classifier was compared with original indices by calculating the average classification error rate and area under the curve (AUC) for the sampled datasets. The performance of this SVM model showed average AUC of 0.76 and average error rate of 0.26, which surpassed all other indices. In conclusion, we developed a convenient tool for primary-care physicians when deferential diagnosis contains thalassemia for the Taiwanese adult population. This approach needs to be validated in other studies or bigger database.

The Effect of Far-Infrared Therapy on the Peritoneal Membrane Transport Characteristics of Uremic Patients Undergoing Peritoneal Dialysis: An Open-Prospective Proof-of-Concept Study.

Long-term peritoneal dialysis (PD) can lead to detrimental changes in peritoneal membrane function, which may be related to the accumulation of glucose degradation products. A previous study demonstrated that 6 months of far-infrared (FIR) therapy may decrease glucose degradation products in PD dialysate. Due to limited literature on this matter, this study aims to investigate the effect of FIR therapy on the peritoneal membrane transport characteristics of PD patients. Patients were grouped according to baseline peritoneal transport status: lower transporters (low and low-average) and higher transporters (high-average and high). Both groups underwent 40 min of FIR therapy twice daily for 1 year. In lower transporters, FIR therapy increased weekly dialysate creatinine clearance (6.91 L/wk/1.73 m2; p = 0.04) and D/P creatinine (0.05; p = 0.01). In higher transporters, FIR therapy decreased D/P creatinine (-0.05; p = 0.01) and increased D/D0 glucose (0.05; p = 0.006). Fifty percent of high transporter patients shifted to high-average status after FIR therapy. FIR therapy may decrease D/P creatinine for patients in the higher transporter group and cause high transporters to shift to high-average status, which suggests the potential of FIR therapy in improving peritoneal membrane function in PD patients.

Influence of Traditional Chinese Medicine on Medical Adherence and Outcome in Estrogen Receptor (+) Breast Cancer Patients in Taiwan: A Real-World Population-Based Cohort Study.

BACKGROUND: Medical adherence is often higher in clinical trials than in real world practice. The aim of this study was to investigate the effects of traditional Chinese medicine (TCM) on medical adherence to hormonal therapy (HT) and survival outcome in ER (+) breast cancer patients in Taiwan. SUBJECTS AND METHODS: Using a nationwide longitudinal population-based database, we enrolled patients with newly diagnosed ER-positive breast cancer who had received HT, and followed for up to 5 years (N = 872). Medication adherence in terms of medication possession ratios (MPR) and patient outcome were evaluated with or without TCM exposure. We applied logistic regression and Cox proportional hazards (PH) analysis to identify factors, including TCM exposure, associated with adherence to HT and mortality. RESULTS: MPR to HT in general decreased over the 5-year period post breast cancer diagnosis. Both TCM and MPR to HT ≥ 80% were significantly associated with reduced risk of breast cancer-associated mortality. Subgroup analysis revealed that TCM annual visits ≥ 3 times with CHP prescription 1~90 days per year affected mortality reduction most significantly (HR: 0.26; 95% CI = 0.08-0.83; p < 0.05) compared to other TCM use. In contrast, using TCM (either short-term or long-term) was not associated with MPR in HT. CONCLUSIONS: Our results supported the potential advantage of TCM on breast cancer-associated mortality, whereas TCM use does not compromise medical adherence to HT. This study offers important insights in integrative therapy for HT in patients with estrogen receptor (+) breast cancer.

Impact of time-varying center volume on technique failure and mortality in peritoneal dialysis.

BACKGROUND: Most studies on volume-outcome association used the number of patients at a particular period as the independent variable. However, peritoneal dialysis (PD) is a chronic treatment, and center volume usually changes over a patient's treatment period. Accordingly, this study used the time-varying center volume to explore the volume-outcome association in PD. METHODS: We conducted a nationwide population-based retrospective cohort study, which included patients who began chronic PD between 2001 and 2010. The risk factors of 5-year technique failure and mortality were analyzed using cause-specific and subdistribution hazard models, respectively. The annual number of patients initiating PD in each patient's treatment center was modeled as a time-varying variable with four categories. RESULTS: We included 9071 patients who started PD in 100 centers where the number of incident patients ranged from 1 to 107 patients per year (median, 25; interquartile range, 13-42). The estimated 5-year patient and technique survival rates were 64.7% and 66.6%, respectively. Being treated in centers in the largest volume category (the number of incident PD patients ≥43 per year) was associated with significantly lower cause-specific and cumulative hazards for technique failure. No association was found between facility volume and hazards of mortality. CONCLUSIONS: Receiving PD in high-volume facilities was associated with a lower risk in technique failure. No association was found between facility volume and mortality risk.

Semiparametric mixture cure model analysis with competing risks data: Application to vascular access thrombosis data.

The article is motivated by a nephrology study in Taiwan, which enrolled hemodialysis patients who suffered from vascular access thrombosis. After treatment, some patients were cured of thrombosis, while some may experience recurrence of either type (acute or nonacute) of vascular access thrombosis. Our major interest is to estimate the cumulative incidence probability of time to the first recurrence of acute thrombosis after therapy. Since the occurrence of one type of vascular access thrombosis precludes occurrence of the other type, patients are subject to competing risks. To account for the presence of competing risks and cured patients, we develop a mixture model approach to the regression analysis of competing-risks data with a cure fraction. We make inference about the effects of factors on both the cure rate and cumulative incidence function (CIF) for a failure of interest, which are separately specified in the logistic regression model and semiparametric regression model with time-varying and time-invariant effects. Based on two-stage method, we develop novel estimation equations using the inverse probability censoring weight techniques. The asymptotic properties of the estimators are rigorously studied and the plug-in variance estimators can be obtained for constructing interval estimators. We also propose a lack-of-fit test for assessing the adequacy of the proposed model and several tests for time-varying effects. The simulation studies and vascular access thrombosis data analysis are conducted to illustrate the proposed method.

The association of sociodemographic factors and needs of haemodialysis patients according to Maslow's hierarchy of needs.

AIM: To explore the association between the sociodemographic factors and the needs of patients undergoing haemodialysis in Taiwan. BACKGROUND: Concomitant discomfort, including physical and mental aspects, affects the patients' quality of life and their willingness to undergo haemodialysis. Maslow's hierarchy of needs is a well-known tool to assess different levels of human needs. METHOD: We conducted a small-scale cross-sectional observational study using a structured needs assessment questionnaire on 159 patients from the Taipei Veterans General Hospital haemodialysis unit. RESULTS: The overall mean scores of physical, mental, spiritual, other needs and needs in relation to medical staff care were 4.0 ± 0.8, 3.2 ± 0.8, 2.7 ± 1.0, 3.1 ± 0.9 and 4.1 ± 0.7, respectively. The results showed that the patients' highest need was in relation to medical staff care, followed by physical needs. Further analysis showed that patients who are still employed during the treatment process have higher mental, spiritual and other needs. Patient who is financially supported by their family has higher physical needs. Patients taken care of by paid caregivers have lower spiritual needs and other needs. This is also the same with patients who are religious as opposed to those who are nonreligious. Patients who have attained tertiary education have higher other needs compared with patients who have only achieved up to primary or secondary education. CONCLUSION: The study is the first in Taiwan to identify and quantify the needs of patients undergoing haemodialysis. When the needs of the patients are identified in relation to their sociodemographic factors, the medical staff can give the appropriate treatment in order to meet the needs and improve the patients' well-being. RELEVANCE TO CLINICAL PRACTICE: Healthcare providers should not only focus on the patients' physiological needs, but should determine and address their other needs in various aspects in order to improve the quality and efficacy of the dialysis care process.

Semiparametric transformation models for interval-censored data in the presence of a cure fraction.

Mixed case interval-censored data arise when the event of interest is known only to occur within an interval induced by a sequence of random examination times. Such data are commonly encountered in disease research with longitudinal follow-up. Furthermore, the medical treatment has progressed over the last decade with an increasing proportion of patients being cured for many types of diseases. Thus, interest has grown in cure models for survival data which hypothesize a certain proportion of subjects in the population are not expected to experience the events of interest. In this article, we consider a two-component mixture cure model for regression analysis of mixed case interval-censored data. The first component is a logistic regression model that describes the cure rate, and the second component is a semiparametric transformation model that describes the distribution of event time for the uncured subjects. We propose semiparametric maximum likelihood estimation for the considered model. We develop an EM type algorithm for obtaining the semiparametric maximum likelihood estimators (SPMLE) of regression parameters and establish their consistency, efficiency, and asymptotic normality. Extensive simulation studies indicate that the SPMLE performs satisfactorily in a wide variety of settings. The proposed method is illustrated by the analysis of the hypobaric decompression sickness data from National Aeronautics and Space Administration.

Methotrexate might reduce ischemic stroke in patients with rheumatoid arthritis: a population-based retrospective cohort study.

AIM: To investigate the effects of hydroxychloroquine, sulfasalazine and methotrexate on ischemic stroke in patients with rheumatoid arthritis (RA). METHODS: This population-based retrospective cohort study included 7904 RA patients and 15 808 non-RA patients between 2000 and 2010. All of the participants were sampled from the National Health Insurance Research Database (NHIRD) of Taiwan. Using univariate analyses, these two groups of patients were compared to evaluate the differences in disease-modifying anti-rheumatic drugs usage and demographic variables. Cox proportional hazard models and Schoenfeld residuals test were performed to estimate the hazard ratios for ischemic stroke and proportional hazard assumptions of these drugs, respectively. RESULTS: The mean age of participants was about 53 years old, and about 70% of RA patients were women. The hazard ratio for ischemic stroke was 1.21 (95% CI: 1.10-1.34; P < 0.01) in the case group compared with the control group, and this significant difference persisted throughout the 10-year period. With respect to RA patients, while hydroxychloroquine showed an insignificant protective effect on ischemic stroke, sulfasalazine and methotrexate were found out to have inconsistent effects during these 10 years. The proportional hazard assumption test of methotrexate at > 0.5 defined daily dose (8.75 mg/week) was violated at a significant level after adjustment (P = 0.0002). CONCLUSIONS: At a dosage of > 0.5 defined daily dose, short-term methotrexate might decrease ischemic stroke risk in RA patients, while hydroxychloroquine and sulfasalazine were neutral.

Conditional maximum likelihood estimation in semiparametric transformation model with LTRC data.

Left-truncated data often arise in epidemiology and individual follow-up studies due to a biased sampling plan since subjects with shorter survival times tend to be excluded from the sample. Moreover, the survival time of recruited subjects are often subject to right censoring. In this article, a general class of semiparametric transformation models that include proportional hazards model and proportional odds model as special cases is studied for the analysis of left-truncated and right-censored data. We propose a conditional likelihood approach and develop the conditional maximum likelihood estimators (cMLE) for the regression parameters and cumulative hazard function of these models. The derived score equations for regression parameter and infinite-dimensional function suggest an iterative algorithm for cMLE. The cMLE is shown to be consistent and asymptotically normal. The limiting variances for the estimators can be consistently estimated using the inverse of negative Hessian matrix. Intensive simulation studies are conducted to investigate the performance of the cMLE. An application to the Channing House data is given to illustrate the methodology.

Impact of Sjogren's syndrome on Parkinson's disease: A nationwide case-control study.

OBJECTIVE: To investigate whether Sjogren's syndrome would have an influence on the development of Parkinson's disease. METHODS: A population-based case-control study was conducted. Participants consisted of 7716 subjects with newly diagnosed Parkinson's disease and a population of 75129 matched control subjects between 2000 and 2010. We measured the risk of Parkinson's disease in association with Sjogren's syndrome by using adjusted odds ratios. RESULTS: A total of 143 Parkinson's disease subjects (1.9%) and 893 control subjects (1.2%) suffered from Sjogren's syndrome (p < 0.001). The crude odds ratio for Parkinson's disease among subjects with Sjogren's syndrome was 1.56 (95% CI 1.30-1.86; p < 0.01). After adjustment for potential confounders which have been proposed that would increase the risk of development of Parkinson's disease, Sjogren's syndrome was found to be significantly associated with the risk of Parkinson's disease with an odds ratio of 1.37 (95% CI 1.15-1.65; p < 0.01). CONCLUSION: This study preliminarily proposed that Sjogren's syndrome was significant associated with an increased risk of Parkinson's disease.

Semiparametric regression analysis of failure time data with dependent interval censoring.

Interval-censored failure-time data arise when subjects are examined or observed periodically such that the failure time of interest is not examined exactly but only known to be bracketed between two adjacent observation times. The commonly used approaches assume that the examination times and the failure time are independent or conditionally independent given covariates. In many practical applications, patients who are already in poor health or have a weak immune system before treatment usually tend to visit physicians more often after treatment than those with better health or immune system. In this situation, the visiting rate is positively correlated with the risk of failure due to the health status, which results in dependent interval-censored data. While some measurable factors affecting health status such as age, gender, and physical symptom can be included in the covariates, some health-related latent variables cannot be observed or measured. To deal with dependent interval censoring involving unobserved latent variable, we characterize the visiting/examination process as recurrent event process and propose a joint frailty model to account for the association of the failure time and visiting process. A shared gamma frailty is incorporated into the Cox model and proportional intensity model for the failure time and visiting process, respectively, in a multiplicative way. We propose a semiparametric maximum likelihood approach for estimating model parameters and show the asymptotic properties, including consistency and weak convergence. Extensive simulation studies are conducted and a data set of bladder cancer is analyzed for illustrative purposes. Copyright © 2017 John Wiley & Sons, Ltd.

The effect of anti-rheumatic medications for coronary artery diseases risk in patients with rheumatoid arthritis might be changed over time: A nationwide population-based cohort study.

OBJECTIVES: To determine whether anti-rheumatic drug usage is associated with risk of coronary artery diseases (CAD) in incident Rheumatoid Arthritis (RA) patients. METHODS: Data were obtained from the Taiwan National Health Insurance Research Database. The study cohort comprised 6260 patients who were newly diagnosed with RA between 2001-2010. The study endpoint was occurrence of CAD according to the ICD-9-CM codes. We used the WHO Defined Daily Dose (DDD) as a tool to assess the drugs exposure. The Cox proportional hazards regression model was used to estimate the hazard ratio (HR) of disease after controlling for demographic and other co-morbidities. When the proportionality assumption is violated, a spline curve of the Scaled Schoenfeld residuals is fitted to demonstrate the estimated effect on CAD over time for drug usage. RESULTS: Among RA patients, use of celecoxib, and etoricoxib was associated with significantly decreased incidence of CAD. The adjusted HR(95% CI) of CAD for low-dose celecoxib (DDD≦1) and high-dose user were 0.47(0.34, 0.65) and 0.37(0.24, 0.58) during the 4 year follow-up time; however, it became 0.98(0.70, 1.37) and1.29(0.85, 1.95). Adjusted HR(95% CI) of CAD for etoricoxib users remained 0.47(0.26, 0.84). CONCLUSIONS: This study revealed association of decreased CAD risk in RA patients taking 2 different kinds of COX-2i in comparison with nonusers. The effect might be changed over time, after about 4 years.

A semiparametric mixture cure survival model for left-truncated and right-censored data.

In follow-up studies, the disease event time can be subject to left truncation and right censoring. Furthermore, medical advancements have made it possible for patients to be cured of certain types of diseases. In this article, we consider a semiparametric mixture cure model for the regression analysis of left-truncated and right-censored data. The model combines a logistic regression for the probability of event occurrence with the class of transformation models for the time of occurrence. We investigate two techniques for estimating model parameters. The first approach is based on martingale estimating equations (EEs). The second approach is based on the conditional likelihood function given truncation variables. The asymptotic properties of both proposed estimators are established. Simulation studies indicate that the conditional maximum-likelihood estimator (cMLE) performs well while the estimator based on EEs is very unstable even though it is shown to be consistent. This is a special and intriguing phenomenon for the EE approach under cure model. We provide insights into this issue and find that the EE approach can be improved significantly by assigning appropriate weights to the censored observations in the EEs. This finding is useful in overcoming the instability of the EE approach in some more complicated situations, where the likelihood approach is not feasible. We illustrate the proposed estimation procedures by analyzing the age at onset of the occiput-wall distance event for patients with ankylosing spondylitis.

Heteroscedastic transformation cure regression models.

Cure models have been applied to analyze clinical trials with cures and age-at-onset studies with nonsusceptibility. Lu and Ying (On semiparametric transformation cure model. Biometrika 2004; 91:331?-343. DOI: 10.1093/biomet/91.2.331) developed a general class of semiparametric transformation cure models, which assumes that the failure times of uncured subjects, after an unknown monotone transformation, follow a regression model with homoscedastic residuals. However, it cannot deal with frequently encountered heteroscedasticity, which may result from dispersed ranges of failure time span among uncured subjects' strata. To tackle the phenomenon, this article presents semiparametric heteroscedastic transformation cure models. The cure status and the failure time of an uncured subject are fitted by a logistic regression model and a heteroscedastic transformation model, respectively. Unlike the approach of Lu and Ying, we derive score equations from the full likelihood for estimating the regression parameters in the proposed model. The similar martingale difference function to their proposal is used to estimate the infinite-dimensional transformation function. Our proposed estimating approach is intuitively applicable and can be conveniently extended to other complicated models when the maximization of the likelihood may be too tedious to be implemented. We conduct simulation studies to validate large-sample properties of the proposed estimators and to compare with the approach of Lu and Ying via the relative efficiency. The estimating method and the two relevant goodness-of-fit graphical procedures are illustrated by using breast cancer data and melanoma data. Copyright © 2016 John Wiley & Sons, Ltd.

The partly Aalen's model for recurrent event data with a dependent terminal event.

Recurrent event data are commonly observed in biomedical longitudinal studies. In many instances, there exists a terminal event, which precludes the occurrence of additional repeated events, and usually there is also a nonignorable correlation between the terminal event and recurrent events. In this article, we propose a partly Aalen's additive model with a multiplicative frailty for the rate function of recurrent event process and assume a Cox frailty model for terminal event time. A shared gamma frailty is used to describe the correlation between the two types of events. Consequently, this joint model can provide the information of temporal influence of absolute covariate effects on the rate of recurrent event process, which is usually helpful in the decision-making process for physicians. An estimating equation approach is developed to estimate marginal and association parameters in the joint model. The consistency of the proposed estimator is established. Simulation studies demonstrate that the proposed approach is appropriate for practical use. We apply the proposed method to a peritonitis cohort data set for illustration.