AdapEnetClass (version 1.2)

WEnetCC.aft: Weighted elastic net with censoring constraints for censored data based on AFT models

Description

This function is used to fit AFT models by using weighted elastic net with censorig constraints.

Usage

WEnetCC.aft(X, Y, delta, weight, C, s, lambda2)

Arguments

X

covariate matrix under study, particularly for AFT modelling. The order of matrix covariate is typically n by p.

Y

typically the logarithmic of the survival time under AFT models. Otherwise survival time.

delta

status. it includes value 1 for uncensored and value 0 for censored subject.

weight

vector of observation weights.

C

this is a positive value that accounts for the penalties of violations of constraints. C is typically allowed to take values in a grid such as (0, 0.5, 1, 1.5, ..., 10).

s

this is the optimal equivalent specification for lasso in terms of fraction of the L1 norm. This is obtained from the AEnet.aft function.

lambda2

regularization parameter for the L2 norm of the coefficients. This is typically assumed to take values in a relatively small grid, say (0, 0.5, 1.0, 1.5, 2.0, ..., 5).

Value

Details

This function is used to fit AFT models using weighted elastic net with censorig constraints (Khan and Shaw, 2015). This is an extension of the weighted elastic approach that allows the censoring constraints to be implemented into the optimization framework. The method works for both cases of datasets high-dimensional where the number of variables (p) is greater than the number of subjects (n) and low-dimensional.

This method use the same optimal pair of (lambda1, lambda2) as found in AEnet.aft. Then C is typically allowed to take values in a grid such as (0, 0.5, 1.0, 1.5, 2.0,...,5), and the optimal value for C obtained by 5-fold cross-validation. Here C typically depends upon how stringently one wants the model to satisfy the censoring constraints compared to how good is the prediction for uncensored data.

References

Hong, D. and Zhang, F. (2010). Weighted Elastic Net Model for Mass Spectrometry Imaging Processing. Mathematical Modelling of Natural Phenomena, 5, 115-133.

Khan and Shaw (2015) imputeYn: Imputing the last largest censored observation/observations under weighted least squares. R package version 1.3, https://cran.r-project.org/package=imputeYn.

Khan and Shaw (2015). Variable Selection for Survival Data with a Class of Adaptive Elastic Net Techniques. Statistics and Computing (published online; DOI: 10.1007/s11222-015-9555-8). Also available in http://arxiv.org/abs/1312.2079.

See Also

cv.AWEnetCC

Examples

#For full data typically used for AFT models (using imputeYn (2015) package)
dat<-data(n=100, p=10, r=0, b1=c(rep(5,5),rep(0,5)), sig=1, Cper=0)

#This needs to run for generating weights of the observations
l<-mrbj(cbind(dat$y, dat$delta) ~ dat$x, mcsize=10, trace=FALSE, gehanonly=TRUE)

#WEnetCC.aft: weighted elastic net with censoring constraints
wt<-round(l$gehansd)
ft<-WEnetCC.aft(dat$x, dat$y, dat$delta, weight=wt, C=1, s=959596, lambda2=0.5)
ft