aft.qp: Fiting Regularized Weighted Least Squares Method by Quadratic...

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/imputeYn.R

Description

The AFT model is solved by quadratic programming.

Usage

1
aft.qp(X, Y, delta)

Arguments

X

the covariate matrix of size n by p.

Y

survival time.

delta

status.

Details

The AFT model is solved using weighted least squares with ridge penalty and censoring constraints (Khan and Shaw, 2013b). The optimization is based on quadratic programming that facilitates incorporating additional censoring constraints into the model. The ridge penalty used here is 0.01 sqrt(2log(p))

Value

Only gives the parameter estimates as below

beta0

the intercept of the AFT model

beta

the coefficients of the covariates

Author(s)

Hasinur Rahaman Khan and Ewart Shaw

References

Khan and Shaw. (2013a). On Dealing with Censored Largest Observations under Weighted Least Squares. CRiSM working paper, Department of Statistics, University of Warwick, UK, No. 13-07. Also available in http://arxiv.org/abs/1312.2533.

Khan and Shaw (2013b). Variable Selection with The Modified Buckley-James Method and The Dantzig Selector for High-dimensional Survival Data. Proceedings 59th ISI World Statistics Congress, 25-30 August 2013, Hong Kong, p. 4239-4244.

See Also

imputeYn

Examples

1
2
3
4
5
6
7
8
9
# For uncorrelated dataset
data1<-data(n=100, p=2, r=0, b1=c(2,4), sig=1, Cper=0)
fit<-aft.qp(data1$x, data1$y, data1$delta)
fit

# For correlated dataset
data2<-data(n=100, p=2, r=0.5, b1=c(2,4), sig=1, Cper=0)
fit2<-aft.qp(data2$x, data2$y, data2$delta)
fit2

imputeYn documentation built on May 29, 2017, 2:18 p.m.