[R] R error using Survr function with gcmrec

[David Winsemius]

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 at 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 at 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