Re: [R] val.surv

2011-08-21 Thread Frank Harrell
bit: > > w <- structure(list(harefit = f, p = est.surv, actual = actual, > pseq = pseq, actualseq = actualseq, u = u, fun = fun, > n = nrow(S), d = sum(S[, 2]), units = units), class = > "val.survh") > return(w) > > But a much

Re: [R] val.surv

2011-08-21 Thread David Winsemius
, d = sum(S[, 2]), units = units), class = "val.survh") return(w) But a much better option would be to report the error to Frank Harrell. I'm copying him since I think your .txt files probably reached the list as well as my mailbox. -- David. ___

Re: [R] val.surv

2011-08-21 Thread Salvo Mac
ot;G:\\test.txt", header=T, sep="\t")) f.1<-cph(Surv(time,event)~age, x=T, y=T,surv=T, data=train) val.surv(f.1, newdata=test, u=10) #plot(calibrate(f.1, u=30, B=20)) From: David Winsemius To: Salvo Mac Cc: "r-help@R-project.org" Se

Re: [R] val.surv

2011-08-20 Thread David Winsemius
ay without data and code. -- David. From: David Winsemius Cc: "r-help@R-project.org" Sent: Sunday, August 21, 2011 3:19 AM Subject: Re: [R] val.surv On Aug 20, 2011, at 8:08 PM, Salvo Mac wrote: Thanks David However, I tried your trick on val.sur

Re: [R] val.surv

2011-08-20 Thread Salvo Mac
? From: David Winsemius Cc: "r-help@R-project.org" Sent: Sunday, August 21, 2011 3:19 AM Subject: Re: [R] val.surv On Aug 20, 2011, at 8:08 PM, Salvo Mac wrote: > Thanks David > > However, I tried your trick on val.surv with newdata=test['age'] but still

Re: [R] val.surv

2011-08-20 Thread David Winsemius
object. I said it was a guess. Now stop wasting our time and offer what is needed. -- david. Salvo From: David Winsemius Cc: "r-help@R-project.org" Sent: Sunday, August 21, 2011 12:55 AM Subject: Re: [R] val.surv On Aug 20, 2011, at 3:32 PM

Re: [R] val.surv

2011-08-20 Thread Salvo Mac
tion: Warning message: In est.surv + S[, 1] :   longer object length is not a multiple of shorter object length Salvo From: David Winsemius Cc: "r-help@R-project.org" Sent: Sunday, August 21, 2011 12:55 AM Subject: Re: [R] val.surv On Aug 20, 2011, at

Re: [R] val.surv

2011-08-20 Thread David Winsemius
On Aug 20, 2011, at 3:32 PM, Salvo Mac wrote: Dear R-users, I have two questions regarding validation and calibration of Survival regression models. 1. I am trying to calibrate and validate a cox model using val.surv. here is my code: f.1<-cph(Surv(time,event)~age, x=T, y=T, data=train

[R] val.surv

2011-08-20 Thread Salvo Mac
 Dear R-users,   I  have two questions regarding validation and calibration of Survival regression models.   1.  I am trying to calibrate and validate a cox model using val.surv. here is my code:  f.1<-cph(Surv(time,event)~age, x=T, y=T, data=train)  test1<-test[,"age"]  val.surv(f.1, newdata=data