Dear colleagues,
I hope someone can help me with my problem. I have fitted a cox model with the following syntax: # cox01def <-cph(Surv(TEVENT,EVENT) ~ ifelse(AGE>50, (AGE-50)^2,0) + BMI + # HDL+DIABETES +HISTCAR2 + log(CREAT)+ as.factor(ALBUMIN)+STENOSIS+IMT,data # = XC, x=T, y=T, surv=T) *1 Furthermore I have estimated my beta's also with a Lasso method - Coxpath ( from Glmpath pckage) , namely: # LASSOdata <- list(x=cox01$x, time=XC$TEVENT, status=XC$EVENT) # summary(LASSOdata) # LASSOpath <- coxpath(data=SMARTdata) My problem is the following: I want to make a nomogram with the beta's from the lasso. Is that possible? I know of the nomogram-function from the "design"-library. I thought that when I would replace the coefficients of the cph-fit*1 by the lasso-coefs, the nomogram would change; However it doesn't # Lasso.coefs <- LASSOpath$b.predictor # Cox01def$coef <- Lasso.coefs If it isn't possible to solve this problem with the nomogram-function from the design library, Then I would like to know if there is an alternative function available? I hope to get a reply soon, Martijn Gondrie ( Utrecht, The Netherlands) ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
