On Nov 11, 2010, at 2:50 PM, David Winsemius wrote:
On Nov 11, 2010, at 2:09 PM, Emily wrote:
I'm having the same problem
(???: from a three year-old posting for which you didn't copy any
context.)
and was wondering whether you ever found a
solution? It gives me the error "Error in Survr(id, time, event) :
Data
doesn't match. Every subject must have a censored time" even though
all my
subjects are right-censored, and to be sure, I've even used the
addCenTime
function. Any input appreciated!
Your data has a lot of 0 events at the end of calendar years. That
does not
seem to be the expected format for the Survr records. It appears to
define an
invalid record as one where the only censoring event is at the time
of the
^valid^
last observation. Here's the first line in Survr that is throwing
the error:
if (length(unique(id)) != length(event[event == 0])) {
stop("Data doesn't match. Every subject must have a censored
time")
I suspect you need to collapse your single-year entries with 0
events into multiple year entries with an event.
--
David.
Here's my sample data:
id=c(rep(1,4),rep(2,4),rep(3,4),rep(4,5))
start
=
c
("1996
-01
-01
","1997
-01
-01
","1998
-01
-01
","1998
-03
-15
","1996
-01
-01
","1996
-04
-15
","1997
-01
-01
","1998
-01
-01
","1996
-01
-01
","1997
-01
-01
","1998
-01
-01
","1998
-09
-30
","1996-01-01","1997-01-01","1997-12-15","1998-01-01","1998-06-14")
stop
=
c
("1997
-01
-01
","1998
-01
-01
","1998
-03
-15
","1999
-01
-01
","1996
-04
-15
","1997
-01
-01
","1998
-01
-01
","1999
-01
-01
","1997
-01
-01
","1998
-01
-01
","1998
-09
-30
","1999
-01
-01
","1997-01-01","1997-12-15","1998-01-01","1998-06-14","1999-01-01")
time
=c(366,365,73,292,105,261,365,365,366,365,272,93,366,348,17,164,201)
event=c(0,0,1,0,1,0,0,0,0,0,1,0,0,1,0,1,0)
enum=c(rep(seq(1,4,1),4),5)
var
=
c
(21312,21869,22441,22441,3105,3105,3086,3075,130610,133147,135692,135692,11686,11976,11976,12251,12251
)
data=data.frame(id,start,stop,time,event,enum,var)
dataOK=addCenTime(data)
m<-gcmrec(Survr(id,time,event)~var, data=dataOK)
--
View this message in context:
http://r.789695.n4.nabble.com/R-error-using-Survr-function-with-gcmrec-tp858931p3038374.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
David Winsemius, MD
West Hartford, CT
______________________________________________
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.
David Winsemius, MD
West Hartford, CT
______________________________________________
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.