Re: [R] Nagelkerkes R2N

2009-07-15 Thread North, Bernard V
I am interested Andrea is whether you ever established why your R2 was 1. I have had a similar situation previously. My main issue though, which I'd be v grateful for advice on, is why I am obtaining such negative values -0.3 for Somers Dxy using validate.cph from the Design package given m

Re: [R] Nagelkerkes R2N

2009-05-13 Thread Frank E Harrell Jr
A new version of Design will be posted to CRAN in the next 2 days. After than, update your system, including an update to the survival package. Then re-try. Your formula is wrong as it can't be negative. LR should be the likelihood ratio chi-square stat : -2 times the difference in the two

[R] Nagelkerkes R2N

2009-05-13 Thread Andrea Weidacher
Hello All, as I´m new to R and survival analysis, I´ve got a question about the Design::validate function: My Code: cox <- cph(Surv(t,status) ~ var1 + var2 + var3, data=data, x=TRUE, y=TRUE, surv=TRUE) cox.val <- validate(cox, B=10, dxy=TRUE, pr=TRUE); My output (cox.val): inde