Re: [R] Estimating survival times with glmnet and coxph

2012-05-07 Thread Terry Therneau
What you are doing is correct. However‹ because it is not informed about the penalty parameters used in creating the fit, the coxph command is not able to correctly reconstruct the variance matrix; this is the source of all of the ³singular² messages. The confidence bands for the survival curves

[R] estimating survival times with glmnet and coxph

2012-05-07 Thread Damjan Krstajic
Dear all, I am using glmnet (Coxnet) for building a Cox Model and to make actual prediction, i.e. to estimate the survival function S(t,Xn) for a new subject Xn. If I am not mistaken, glmnet (coxnet) returns beta, beta*X and exp(beta*X), which on its own cannot generate S(t,Xn). We miss baseline