 Research article
 Open Access
 Published:
Cancer survival analysis using semisupervised learning method based on Cox and AFT models with L_{1/2} regularization
BMC Medical Genomics volume 9, Article number: 11 (2016)
Abstract
Background
One of the most important objectives of the clinical cancer research is to diagnose cancer more accurately based on the patients’ gene expression profiles. Both Cox proportional hazards model (Cox) and accelerated failure time model (AFT) have been widely adopted to the high risk and low risk classification or survival time prediction for the patients’ clinical treatment. Nevertheless, two main dilemmas limit the accuracy of these prediction methods. One is that the small sample size and censored data remain a bottleneck for training robust and accurate Cox classification model. In addition to that, similar phenotype tumours and prognoses are actually completely different diseases at the genotype and molecular level. Thus, the utility of the AFT model for the survival time prediction is limited when such biological differences of the diseases have not been previously identified.
Methods
To try to overcome these two main dilemmas, we proposed a novel semisupervised learning method based on the Cox and AFT models to accurately predict the treatment risk and the survival time of the patients. Moreover, we adopted the efficient L_{1/2} regularization approach in the semisupervised learning method to select the relevant genes, which are significantly associated with the disease.
Results
The results of the simulation experiments show that the semisupervised learning model can significant improve the predictive performance of Cox and AFT models in survival analysis. The proposed procedures have been successfully applied to four real microarray gene expression and artificial evaluation datasets.
Conclusions
The advantages of our proposed semisupervised learning method include: 1) significantly increase the available training samples from censored data; 2) high capability for identifying the survival risk classes of patient in Cox model; 3) high predictive accuracy for patients’ survival time in AFT model; 4) strong capability of the relevant biomarker selection. Consequently, our proposed semisupervised learning model is one more appropriate tool for survival analysis in clinical cancer research.
Background
An important objective of clinical cancer research is to develop tools to accurately predict the survival time and risk profile of patients based on the DNA microarray data and various clinical parameters. There are several existing techniques in the literature for performing this type of survival analysis. Among of them, both Cox proportional hazards model (Cox) [1] and the accelerated failure time model (AFT) [2] have been widely used. Cox model is the most popular approach by far in survival analysis to assess the significance of various genes in the survival risk of patients through the hazard function. On the other hand, the requirement for analyzing failure time data arises in investigating the relationship between a censored survival outcome and highdimensional microarray gene expression profiles. Therefore, AFT model has been studied extensively in recent years. However, various current cancer survival analysis mechanisms have not demonstrated themselves to be very accurate as expected. The accuracy problems, in essence, are related to some fundamental dilemmas in cancer survival analysis. We believe any attempt to improve the accuracy of survival analysis method has to compromise between these two dilemmas:

The small sample size and censored survival data versus high dimensional covariates dilemma in Cox model
Highdimensional survival analysis in particular has attracted much interest due to the popularity of microarray studies involving survival data. This is statistically challenging because the number of genes, p, is typically hundreds of times larger than the number of microarray samples, n (p> > n). For survival analysis, sample size is reduced significantly by the availability of followup data for the analyzed samples. In fact, in publicly available gene expression databases, only a small fraction of humantumor microarray datasets provides clinical followup data. A “lowrisk” or “highrisk” classification based on Cox model usually relies on traditional supervised learning techniques, in which only completed data (i.e., data from samples with clinical followup) can be used for learning, while censored data (i.e., data from samples without clinical followup) are disregarded. Thus, the small sample size and censored survival data remain a bottleneck in obtaining robust and accurate classifiers with Cox model. Recently a technique called semisupervised learning [3] in machine learning suggests that censored data, when used in conjunction with limited amount of completed data, can produce considerable improvement in learning accuracy. Indeed, semisupervised learning has been proved to be effective in solving different biological problems, such as protein classification [4, 5], drugprotein interaction prediction [6] and prediction of interactions between disease and human proteins [7]. Moreover, there are some semisupervised learning approaches worked on the gene expression data. For example, “corrected” Cox scores were used for semisupervised prediction using principal component regression by Bair and Tibshirani [8] and the semisupervised classification using nearestneighbor shrunken centroid clustering by Tibshirani et al. [9].

The similar phenotype disease versus different genotype cancer dilemma in the AFT model
In the accelerated failure time model, to increase the available sample size and get the more accurate result, each censored observation time is replaced with the imputed value using some estimators, such as the inverse probability weighting (IPW) [10] method, mean imputation method, BuckleyJames method [11] and rankbased method. In fact, these estimation methods assume that the AFT model was used for the patients with similar phenotype cancer, and the survival times should satisfy the same unspecified common probability distribution. Nevertheless, the disparity we see in disease progression and treatment response can be attributed to that the similar phenotype cancer may be completely different diseases on the molecular genotype level. So we need to identify different cancer genotypes. Can we do it based exclusively on the clinical data? For example, patients can be assigned to a “lowrisk” or a “highrisk” subgroup based on whether they were still alive or whether their tumour had metastasized after a certain amount of time. This approach has also been used to develop procedures to diagnose patients [12]. However, by dividing the patients into subgroups just based on their survival times, the resulting subgroups may not be biologically meaningful. Suppose, for example, the underlying cell types of each patient are unknown. If we were to assign patients to “lowrisk” and “highrisk” subgroups based on their survival times, many patients would be assigned to the wrong subgroup, and any future predictions based on this model would be suspect. Therefore, we need propose more accurate classification methods by identifying these underlying cancer subtypes based on microarray data and clinical data together, and build a model that can determine which subtype is present in future patients.
Our idea in this study is to strike a tactical balance between the two contradictory dilemmas. We propose a novel semisupervised learning method based on the combination of Cox and AFT models with L_{1/2} regularization for highdimensional and low sample size biological data. In our semisupervised learning framework, the Cox model can classify the “lowrisk” or a “highrisk” subgroup though samples as many as possible to improve its predictive accuracy. Meanwhile, the AFT model can estimate the censored data in the subgroup, in which the samples have the same molecular genotype.
Methods
Cox proportional hazards model (Cox)
The Cox proportional hazards model is now the most widely used for survival analysis to classify the patients into “lowrisk” or “highrisk” subgroup after prognostic. Under the Cox model, the hazard function for the covariate matrix x with sample size n and the number of genes p is specified as λ(t) = λ_{0}(t)exp(β′x), where t is the survival time and the baseline hazard function λ_{0}(t) is common to all subjects, but is unspecified or unknown. Let ordered risk set at time t(r) be denoted by Rr = {j∈1,…, n:tj ≥ t(r)}. Assume that censoring is non informative and that there are no tied event times. The Cox log partial likelihood can then be defined as
Where D denotes the set of indices for observed events.
Accelerated failure time model (AFT)
The AFT model is a linear regression model for survival analysis, in which the logarithm of response t_{i} is related linearly to covariates x_{i}:
where h(.) is the log transformation or some other monotone function. In this case, the Cox assumption of multiplicative effect on hazard function is replaced with the assumption of multiplicative effect on outcome. In other words, it is assumed that the variables x_{ i } act multiplicatively on time and therefore affect the rate at which individual i proceeds along the time axis. Because censoring is present, the standard least squares approach cannot be employed to estimate the regression parameters in Eq. (2) even when p < n.
One approach for AFT model implementation entails the replacement of censored t_{ i } with imputed values. In order to simplify the method, we use KaplanMeier weight approach to estimate the censored data in the least square criterion. Since for high dimensional and low simple size data, the KaplanMeier weight estimator is more efficient than the BuckleyJames and rank based approaches. Moreover, it also has rigorously and strong theoretical justifications under reasonable conditions [13]. For each censored t_{ i } with the conditional expectation of t_{ j } given t_{ j } > t_{ i } [14], the imputed value h(t_{ i }) can then be given by
where Ŝ is the KaplanMeier estimator (Kaplan and Meier, 1958) of the survival function and ΔŜ(t_{(r)}) is the step of Ŝ at time t_{(r)} [15].
L_{1/2} regularization
In recent years, various regularization methods for survival analysis under the Cox and AFT models have been proposed, which perform both continuous shrinkage and automatic gene selection simultaneously. For example, Coxbased methods utilizing kernel transformations [16], threshold gradient descent minimization [17], and lasso penalization [18] have been proposed. Likewise, a few authors have proposed variable selection methods based on accelerated failure time models. Most of these procedures are based on L_{1} norm, however, the results of L_{1} regularization are not good enough for spartity, especially in biology research. Theoretically, the L_{q} (0 < q < 1) type regularization with the lower value of q would lead to better solutions with more sparsity. Moreover, among L_{q} regularizations with q ∈ (0, 1), only L_{1/2} and L_{2/3} regularizations permit an analytically expressive thresholding representation [19]. In the literature [19], Xu et al. investigated that when 0 < q < 1/2, there are not obvious difference in the variable selection performance of L_{q} (0 < q < 1/2) regularization, but solving the L_{1/2} regularization is much efficient compared to the L_{0} regularization. On the other hand, the L_{1/2} regularization can yield most sparse solutions among L_{q} (1/2 < q < 1) regularizations. Moreover, they also proved some attractive properties of the L_{1/2} regularization, such as unbiasedness, sparsity and oracle properties. Our previous works have also demonstrated the efficiencies of L_{1/2} regularization for Cox and AFT models respectively [20]. The sparse L_{1/2} regularization model has expressed as:
where l is loss function and λ is tuning parameter. Since the penalty function of L_{1/2} regularization is nonconvex, which raises numerical challenges in fitting the Cox and AFT models. Recently, coordinate descent algorithms [21] for solving nonconvex regularization approach (such as SCAD, MCP) have been shown significantly efficiency and convergence [22]. The algorithms optimize a target function with respect to a single parameter at a time, iteratively cycling through all parameters until reached its convergence. Since the computational burden increases only linearly with the number of the covariates p, coordinate descent algorithms can be a powerful tool for solving highdimensional problems.
Therefore, in this paper, we introduce a novel univariate half thresholding operator of the coordinate descent algorithm for the L_{1/2} regularization, which can be expressed as:
where ỹ ^{(j)}_{ i } = ∑_{k ≠ j}x_{ ik }β_{ k } as the partial residual for fitting β_{ j }, ω_{ j } = ∑ ^{n}_{i = 1} x_{ ij }(y_{ i } − ỹ_{ i }^{(j)}), and \( {\varphi}_{\lambda}\left(\omega \right)= arccos\Big(\frac{\lambda }{8}{\left(\frac{\left\omega \right}{3}\right)}^{\frac{3}{2}} \).
Remark: In our previous work [23], we used \( \frac{3}{4}{\left(\uplambda \right)}^{\frac{2}{3}} \) for represent L_{1/2} regularization thresholding operator. Here, we introduced a new half thresholding representation \( \frac{\sqrt[3]{54}}{4}{\left(\uplambda \right)}^{\frac{2}{3}} \). This new value is more precisely and effectively than the old one. Since it is known that the quantity of the solutions of a regularization problem depends seriously on the setting of the regularization parameter λ. Based on this novel thresholding operator, when λ is chosen by some efficient parameters tuning strategy, such as crossvalidation, the convergence of the algorithm is proved [24].
Our proposed semisupervised learning method
Figure 1 illustrates the overview of our proposed semisupervised learning development and evaluation workflow. Microarray gene expression data on a specific cancer type are collected, processed, and separated into completed samples and censored samples. In order to identify tumor subclasses that were both biologically meaningful and clinically relevant, we applied the L_{1/2} regularized Cox model on the completed data to select a group of outcomerelated genes firstly. Thus, all samples including completed and censored cases can be subsequently classified into “lowrisk” and “highrisk” classes. Once such classes are identified, we can evaluate the censored data using the mean imputation approach based on the completed data belonged to the same risk classes, because they are correlated to similar disease biologically meaningful at the molecular level. When the censored data replaced by the appropriate imputation values, the L_{1/2} regularized AFT model can be used to select a list of genes that correlate with the clinical variable of interest, and reevaluate the censored data based on these selected genes. A stratified Kfold crossvalidation is used for regularization parameter tuning. We repeated this semisupervised learning procedure including Cox and AFT steps multiple time with increasing number of available training data and estimating the censored data based on the similar genotype disease.
In the semisupervised learning framework, the predictive accuracy of the Cox and AFT models would be improved because the number of the training data increased and the censored data were imputed reasonably. The L_{1/2} regularization approach can select the significant relevant gene sets based on the Cox and AFT models respectively.
In our proposed semisupervised learning method, the censored data are evaluated from the same risk class to improve prediction performance. However, there are some observable errors in the imputations of the censored data. For example, the estimated survival time by AFT model was even less than the censored time. We regarded them as error estimations, and would not use them for model training.
In this paper, two parameters were used to test the performances obtained by different methods.
Integrated BrierScore (IBS)
The Brier Score (BS) [25] is defined as a function of time t > 0 by:
where Ĝ(⋅) denotes the KaplanMeier estimation of the censoring distribution and Ŝ(⋅X_{ i }) stands to estimate survival for the patient i. Note that the BS(t) is dependent on the time t, and its values are between 0 and 1. The good predictions at the time t result in small values of BS. The integrated Brier Score (IBS) is given by:
The IBS is used to assess the goodness of the predicted survival functions of all observations at every time between 0 and max(t_{ i }).
Concordance Index (CI)
The Concordance Index (CI) can be interpreted as the fraction of all pairs of subjects which predicted survival times are correctly ordered among all subjects that can actually be ordered. By the CI definition, we can determine t_{ i } > t_{ j } when f_{ i } > f_{ j } and δ_{j} = 1 where f(⋅) is survival function. The pairs for which neither t_{ i } > t_{ j } nor t_{i<}t_{ j }can be determined are excluded from the calculation of CI. Thus, the CI is defined as:
Note that the values of CI are between 0 and 1, the perfect predictions of the building model would lead to 1 while have a CI of 0.5 at random.
Results
Simulated experiment
We adopted the simulation scheme in R. Bender’s work [26]. The generation procedure of the simulated data is as follows:

Step 1:
we generate γ_{i0}, γ_{i1},…, γ_{ip} (i = 1,…,n) independently from standard normal distribution and set: \( {X}_{ij}={\gamma}_{ij}\sqrt{1\rho }+{\gamma}_{i0}\sqrt{\rho } \) (j = 1,…, p) where ρ is the correlation coefficient.

Step 2:
The survival time y_{i} is written as: \( {\mathrm{y}}_{\mathrm{i}}=\frac{1}{\alpha } log\left(1\frac{\alpha * \log (U)}{\omega * \exp \left(\beta X\right)}\right) \) which U is an uniformly distributed variable, ω is the scale parameter, α is the shape parameter.

Step 3:
Censoring time point y_{i}′ (i = 1,…n) is obtained from an random distribution E (θ), where θ is determined by specify censoring rate.

Step 4:
Here we define y_{i} = min(y_{i}, y_{i}′) and δ_{i} = I(y_{i} < y_{i}′), the observed data represented as (y_{i}, x_{i}, δ_{i}) for the model are generated.
In our simulated experiments, we build highdimensional and low sample size datasets. In every dataset, the dimension of the predictive genes is p = 1000, in which 10 prognostic genes and their corresponding coefficients are nonzero. The coefficients of the remaining 990 genes are zero. About 40 % of the data in each subgroup are right censored. We considered the training sample sizes are n = 100, 200, 300 and the correlation coefficients of genes areρ = 0 and ρ = 0.3 respectively. The simulated data were applied to the single Cox, single AFT and semisupervised learning approach with Cox and AFT models. For gene selection, we use L_{1/2} regularization approach and the regularization parameters are tuned by 5fold cross validation. To assess the variability of the experiment, each method is evaluated on a test set including 200 samples, and replicated over 50 random training and test partitions.
Figure 2 shows the percentage of data distribution processed by our semisupervised learning model with L_{1/2} regularization in different parameter settings (a: n = 100, ρ = 0.3; b: n = 100, ρ = 0; c: n = 200, ρ = 0.3; d: n = 200, ρ = 0; e: n = 300, ρ = 0.3; f: n = 300, ρ = 0;). The first cylinder represents the simulated dataset, and the cylinders af present the form of the dataset processed by our semisupervised learning model. Compared to the original dataset, the most censored data can be reasonable estimated to the available data by semisupervised learning model. For example, when the training sample n = 300 and the correlation coefficientρ = 0, just 2.41 % censored data cannot conjugate into the available samples because their imputed survival time based on the AFT model is smaller than their observed censored time. Moreover, we can see that with the sample size increases or the correction coefficient decreases, more censored data can be correctly estimated to available training data.
The classification accuracy under the correlation coefficient ρ = 0.3 with different training sample size setting was demonstrated in Fig. 3, the sum of red and blue part represent the samples which can be correctly classified by the Cox model. The first cylinder in each group represents the result obtained by Cox model, and the second one represents the result obtained by our semisupervised learning model. No matter in which group, the semisupervised learning model obtained the high improvements of the classification performance. When the training sample size n = 100, 200, 300, more than 32.23, 20.55 and 15.63 % samples were correctly classified by semiCox model when comparing with the results of the single Cox model.
The precision of our semisupervised learning model with L_{1/2} regularization was given in Table 1. The precision is got from the number of correct selected genes divided the total number of selected genes by the methods. With the sample size increase or the correction coefficients of the features decrease, the classification performances of each model become better. We found the single Cox and single AFT model is difficult to select the whole correct genes in the dataset. This means these models selected too few corrected genes and many other irrelevant genes in their results. This made their prediction precision very low. Nevertheless, our semisupervised learning model solves this problem, the precisions of the semiCox or the semiAFT were both higher than that obtained by the single Cox or single AFT model. After processed by our semi supervised learning method, the number of selected correct genes was increased, and the number of total selected genes were decreased, the semiCox achieved about 130 % improvements in precision compared to the single Cox model. Although the precision improvement of semiAFT model is smaller than that of the semiCox model, it can select most correct genes under different parameter settings. Therefore we think our semisupervised learning method can significantly improve the accuracy of prediction for survival analyses with the highdimensional and low sample size gene expression data.
Simulation analysis of real microarray datasets
In this section, the proposed semisupervised learning approach was applied to the four real gene expression datasets respectively, such as DLBCL (2002) [27], DLBCL (2003) [28], Lung cancer [29], AML [30]. The brief information of these datasets is summarized in Table 2.
In order to accurately assess the performance of the semisupervised learning approach, the real datasets were randomly divided into two pieces: two thirds of the available patient samples, which include the completed and correct imputed censored data, were put in the training set used for estimation and the remaining completed and censored patients’ data would be used to test the prediction capability. We used single Cox and single AFT with L_{1/2}regularization approaches for comparisons and for each procedure, the regularization parameters are tuned by 5fold cross validation. All results in this article are averaged over 50 repeated times respectively.
As show in Fig. 4, our proposed semisupervised learning method can significantly increase the available sample size for classification model training. Especially, in Lung cancer dataset, the available samples increase from 27.91 to 94.19 %. For other three datasets, the available sample sizes also augment from 57.50, 69.56, 57.75 to 96.67, 96.73, 94.84 % respectively. Most censored data were accurately estimated by the AFT model using samples, which belong to the same genotype disease classes, and were sequentially classified into highrisk or lowrisk classes by the Cox model respectively. In addition of that, just small part of the censored data cannot conjugate into the available samples because their imputed survival time based on the AFT model is smaller than their observed censored time. The reason may be the individual differences of the patients.
The integrated brier score (IBS) and the concordance index (CI) measurements were used to evaluate the classification and prediction performance of Cox and AFT models in the semisupervised learning approach. In the IBS measure, the lower value means the more accurate prediction result. As shown in Fig. 5, the values of IBS obtained by our semisupervised learning model with L_{1/2} penalty were smaller than that obtained by the single Cox and AFT models. For example, in the Lung cancer dataset, the IBS values of the Cox and AFT models from 0.2164 and 0.2195 improve to 0.1259 and 0.1341 respectively in the semisupervised learning approach. For the other gene expression datasets DLBCL2002, DLBCL2003 and AML, the IBS values of the Cox model improve 34, 45 and 26 %, and the IBS values of the AFT model improve 34, 36 and 28 % respectively. This means that our proposed semisupervised learning approach can significantly improve the classification and prediction accuracy of the Cox and AFT models. In Fig. 6, the values of CI measure obtained by Cox and AFT with and without the semisupervised learning approaches were given respectively. The CI values belong to the regain [0.5, 1] and its larger value means the more accurate prediction results. As shown in Fig 6, for the Lung cancer dataset, the CI values of the Cox and AFT models from 0.5738 and 0.6013 improve to 0.6620 and 0.7225 respectively in the semisupervised learning approach. The improvement rate is higher than (0.66200.5738)/(0.57380.500) = 120 %. For the other gene expression datasets DLBCL2002, DLBCL2003 and AML, the CI values of the Cox models improve 39, 45 and 25 %, and the CI values of the AFT models improve 56, 45 and 36 % respectively. These also illustrated the semisupervised learning method can significantly improve the accuracy of prediction in survival analysis with the highdimensional and low sample size gene expression data.
Figure 7 gives the number of genes selected by the L_{1/2}regularized Cox and AFT models with and without the semisupervised learning framework. The semiCox and semiAFT selected less genes compared to the single Cox and the AFT model. For example, in the lung cancer dataset, the single Cox and single AFT models select 14 and 22 genes respectively. However, the Cox and AFT models just select 10 and 17 genes in semisupervised learning model. Moreover, Combined the found in the Figs. 5 and 6, the prediction accuracy of Cox and AFT in the semisupervised learning model was significantly improved using more relevant genes.
On the other hand, we find that for these all four gene expression datasets, the selected genes from Cox and AFT models are quite different and just small parts are overlapping. We think the reason may be that the regularized Cox model selects the relevant genes for lowrisk and highrisk classification. Nerveless, the genes selected by the AFT model are high correlation for the survival time of patients. So these two models may select different genes, which have different biological function. Through our below analyses, we know that the genes selected by semisupervised learning methods are significant relevant with the cancer.
Figure 8 shows the survival curves of the Cox model with and without the semisupervised learning method for AML dataset. The xaxis represents the survival days and the yaxis is the estimated survival probability. The green and read curves represent the changes of the survival probability for the “lowrisk” and “highrisk” classes respectively. As show in Fig. 8a, these two curves intersect at the time point of 564 day, which means that the single Cox cannot efficiently classify and predict the survival rate of the patients using the AML dataset. On the contrary, in Fig. 8b, the survival probabilities of the “lowrisk” and “highrisk” patients can be efficiently estimated by the semiCox model. For other three gene expression datasets, we also got the similar results, which are the classification performance of semiCox model significantly outperforms the single Cox model.
Discussion
In this section, we introduce a brief biological discussion of the selected genes for the Lung cancer dataset to demonstrate the superiority of our proposed semisupervised learning method. The number of selected genes by semisupervised learning method is less than the single Cox and AFT model, but includes some genes which are significantly associated with cancer and cannot be selected by the two single Cox and AFT models, such as GDF15, ARHGDIB and PDGFRL. GDF15 belongs to the transforming growth factorbeta superfamily, and is one kind of bone morphogenetic proteins. It was showed that GDF15 can be seen as prognostication of cancer morbidity and mortality in men [31]. ARHGDIB is the member of the Rho (or ARH) protein family; it is involved in many different cell events such as cell secretion, proliferation. It is likely to impact on the cancer [32]. The role of PDGFRL is to encode a protein contains an important sequence which is similar to the ligand binding domain of plateletderived growth factor receptor beta. Biological research has confirmed that this gene can affect the sporadic hepatocellular carcinomas. This suggests that this gene product may get the function of the tumour inhibition.
At the same time, the Cox and AFT models with and without semisupervised learning method also selected some common genes. For example,the PTP4A2, TFAP2C, GSTT2. PTP4A2 is the member of the protein tyrosine phosphatase family, overexpression of PTP4A2 will confer a transformed phenotype in mammalian cells, which suggested its role in tumorigenic is [33]. TFAP2C can encode a protein contains a sequencespecific DNAbinding transcription factor which can activate some developmental genes [34]. GSTT2 is one kind of a member of a superfamily of proteins. It has been proved to play an important role in human carcinogenesis and shows that these genes are linked to cancer with a certain relationship [35].
Through the comparison of the biological analyses of the selected genes, we found the semisupervised method based on Cox and AFT models with L_{1/2} regularization is a competitive method compared to single regularized Cox and AFT models.
Conclusion
To overcome the limitations of fully unsupervised and fully supervised approaches for survival analysis in cancer research, we have developed a discriminative semisupervised method based on Cox and AFT models with L_{1/2} regularization. This method combines the advantages of both Cox and AFT models, and overcome the dilemma in their applications. By comparison the results of Cox and AFT modes with and without the semisupervised method in simulation experiment and real microarray datasets experiment with different regularizing method, we demonstrated that 1) the censored data could be employed after appropriate processing; 2) the semisupervised classification improved prediction accuracy as compared to the state of the art single Cox model; 3) the gene selection performance gain improved with the increase number of available samples. Therefore, for clinical applications, where the goal is often to develop an accurate predicting test using fewer genes in order to control cost, the semisupervised method based on Cox and AFT models with L_{1/2} regularization can be chosen to applied, it will be an efficient and accuracy method based on the highdimensional and lowsample size data in cancer survival analysis.
References
Cox DR. Partial likelihood. Biometrika. 1975;62:269–762.
Wei LJ. The accelerated failure time model: a useful alternative to the Cox regression model in survival analysis. Stat Med. 1992;11:1871–9.
Chapelle O, Sindhwani V, Keerthi SS. Optimization techniques for semisupervised support vector machines. J Mach Learn Res. 2008;9:203–33.
Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A. 2002;99:6567–72.
Wasito I, Veritawati I. Subtype of Cancer Identification for Patient Survival Prediction Using Semi Supervised Method. JCIT. 2012;7:14.
Xia Z, Wu LY, Zhou X, et al. Semisupervised drugprotein interaction prediction from heterogeneous biological spaces. BMC Syst Biol. 2010;4 Suppl 2:S6.
Qi Y, Tastan O, Carbonell JG, et al. Semisupervised multitask learning for predicting interactions between HIV1 and human proteins. Bioinformatics. 2010;26(18):i645–52.
Koestler DC, Marsit CJ, Christensen BC, et al. Semisupervised recursively partitioned mixture models for identifying cancer subtypes. Bioinformatics. 2010;26(20):2578–85.
Golub T, Slonim D, Tamayo P, Huard C, Gaasenbeek M, Mesirov J, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999;286:531–7.
Wang Z, Wang CY. BuckleyJames boosting for survival analysis with highdimensional biomarker data. Stat Appl Genet Mol Biol. 2010;9(1):Article 24.
Seaman SR, White IR, Copas AJ, et al. Combining Multiple Imputation and Inverse‐Probability Weighting. Biometrics. 2012;68(1):129–37.
Bair E, Tibshirani R. Semisupervised methods to predict patient survival from gene expression data. PLoS Biol. 2004;2:E108.
Huang J, Ma S, Xie H. Regularized Estimation in the Accelerated Failure Time Model with High‐Dimensional Covariates. Biometrics. 2006;62(3):813–20.
Tsiatis A. Estimatingregressionparametersusinglinearranktestsforcensored data. Ann Stat. 1996;18:305–28.
Datta S. Estimatingthemeanlifetimeusingrightcensoreddata. Stat Methodol. 2005;2:65–9.
Luan Y, Li H. Modelbased methods for identifying periodically expressed genes based on time course microarray gene expression data. Bioinformatics. 2004;20:332–9.
Gui J, Li H. Threshold gradient descent method for censored data regression, with applications in pharmacogenomics. Pac Symp Biocomput. 2005a;10:272–83.
Gui J, Li H. Penalized Cox regression analysis in the highdimensional and lowsample size settings, with applications to microarray gene expression data. Bioinformatics. 2005b;21:3001–8.
Xu ZB, et al. L1/2 regularization. Sci China. 2010;40(3):1–11. series F.
Liu C, et al. The L1/2 regularization method for variable selection in the Cox model. Appl Soft Comput. 2014;14(c):498–503.
Cox DR. Regression models and lifetables. J R Statist Soc. 1972b;34:187–220.
Ernst J, et al. A semisupervised method for predicting transcription factorgene interactions in Escherichia coli. Plos Comput Biol. 2008;4(3):e1000044.
Fan J, Li R. Variable selection for Cox’s proportional hazards model and frailty model. Ann Statist. 2002;30:74–99.
Gui J, Li H. Penalized Cox regression analysis in the high dimensional and lowsample size settings, with applications to microarray gene expression data. Bioinformatics. 2005;21(13):3001–8.
Murphy AH. A new vector partition of the probability score. J Appl Meteorol. 1973;12(4):595–600.
Bender R, Augustin T, Blettner M. Generating survival times to simulate Cox proportional hazards models. Stat Med. 2005;24:1713–23.
Rosenwald A, et al. The use of molecular profiling to predict survival after chemotherapy for diffuse largeBcell lymphoma. N Engl J Med. 2002;346:1937–46.
Rosenwald A, et al. The proliferation gene expression signature is aquantitative integrator of oncogenic events that predicts survival in mantle cell lymphoma. Cancer Cell. 2003;3:185–97.
Beer DG, et al. Geneexpression profiles predict survival of patients with lung adenocarcinoma. Nat Med. 2002;8:816–24.
Bullinger L, et al. Use of geneexpression profiling to identify prognostic subclasses in adult acute myeloid leukemia. N Engl J Med. 2004;350:1605–16.
Wallentin L, et al. GDF15 for prognostication of cardiovascular and cancer morbidity and mortality in men. PLoS One. 2013;8:12.
Hatakeyama K, et al. Placenta—Specific novel splice variants of Rho GDP dissociation inhibitor beta are highly expressed in cancerous cells. BMC Res Notes. 2012;5:666.
Riker A, et al. The gene expression profiles of primary and metastatic melanoma yields a transition point of tumor progression and metastasis. BMC Med Genomics. 2008;1:13.
Ailan H, et al. Identification of target genes of transcription factor activator protein 2 gamma in breast cancer cells. BMC Cancer. 2009;9:279.
Jang SG, Kim IJ, Kang HC, et al. GSTT2 promoter polymorphisms and colorectal cancer risk. BMC Cancer. 2007;7:16.
Acknowledgements
This research was supported by Macau Science and Technology Develop Funds (Grant No. 099/2013/A3) of Macau Special Administrative Region of the People’s Republic of China.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
YL, HC and XYL developed the semisupervised learning methodology, designed and carried out the comparative study, wrote the code, and drafted the manuscript. ZBX, HZ and KSL brought up the biological problem that prompted the method biological development and verified and provided discussion on the method, and coauthored the manuscript. The authors read and approved the manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Liang, Y., Chai, H., Liu, XY. et al. Cancer survival analysis using semisupervised learning method based on Cox and AFT models with L_{1/2} regularization. BMC Med Genomics 9, 11 (2016). https://doi.org/10.1186/s1292001601696
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1292001601696
Keywords
 Cancer survival analysis
 Semisupervised learning
 Gene selection
 Regularization
 Cox proportional hazards model
 Accelerated failure time model