[R] rpart:- regression tree for survival data

2010-02-23 Thread paaventhan jeyaganth
Dear r-users, i am developing regression tree for censored data, I have difficulty purning the tree, I choose the smallest of the cp , minimizing the predictive error (xerror) but this is not enough, i still have big tree. if anybody know about he pruning Technique, could you please

[R] regression tree for censored data

2010-02-23 Thread paaventhan jeyaganth
Dear r-users, i am developing regression tree for censored data, I have difficulty purning the tree, I choose the smallest of the cp , minimizing the predictive error (xerror) but this is not enough, i still have big tree. if anybody know about he pruning Technique, could you please sent th

[R] Cox Calibration regression test & double graphs

2010-03-09 Thread paaventhan jeyaganth
Dear All, How can i do the bootstrap and calibration plot after surival analysis , i don't know how to do it after that please help me. Thanks here is my syntex library(survival) mfit <- survfit(Surv(days, status == 1)~drug, data = melanom) #how can i do the bootstra #Cox Calibratio

[R] remove text from the calibration plot

2010-03-12 Thread paaventhan jeyaganth
Dear all, how can i get rid of the text in the calibration plot below the x-axis. thanks Paaveen _ Check your Hotmail from your phone. [[alternative HTML version deleted]]

[R] removing text form the graph(outer region)

2010-03-24 Thread paaventhan jeyaganth
Dear R communities, when i do a calibration plot i have text inside the graph bottom left and right (outer region) saying x resampling added B=200, i want to get rid of this texts please advise me how can i do it Thanks very much Paaveen

[R] bagging survival tree

2010-04-01 Thread paaventhan jeyaganth
Dear R users, I have problem with bagging survial tree after finding the final tree. f<-rpart(Surv(time ,dead )~ x1 +x2+ x3+x4+x5+x6, data=crp) f.prun<-prune(f,cp=0.036701) # final tree i have 3 endnote including #the variable x3 and x4 how can i use the bagging code i use like this but it d

[R] Brier's score for bootstrap sample (coxph)

2010-04-09 Thread paaventhan jeyaganth
Dear all, How can i get brier's score for the bootsrap sample for survival analysis. this are the code i am using for the validation. f1 <- cph(Surv(time,dead ) ~ strata(x1)+strata(x2)+strata(x3), x=TRUE, y=TRUE, surv=TRUE, time.inc=12, data=new) validate(f1,B=200,u=12,dxy=T) Tha

[R] competing risk analysis

2010-04-29 Thread paaventhan jeyaganth
Deart All, I am doing competing risk analysis, I am tring store the estimate value in bbb but it gives error, it's anybody knows how to store in estimated value. Thanks Paaveen xx <- cuminc(wkdt$fu.year, wkdt$status) bbb <- cuminc.est(xx,F) Error in cuminc.est: Argument number 2 no

[R] c-statiscs 95% CI for cox regression model

2010-07-24 Thread paaventhan jeyaganth
Dear all, how can i do the calculate the C-statistics 95% confidences interval for the cox regression model. Thanks very much for your any help. Paaveenthan _ [[elided Hotmail spam]]

[R] calibration plot labels

2010-08-15 Thread paaventhan jeyaganth
Dear all, when i do the calibration plot, i put the x label y label , there is some labels are i did not put it , like "resample optimism added ..." i want to get rid of the these label , is any body know how can i get rid of these label. these are the following command i used cal <