I have a competing risk data where a patient may die from either AIDS or
Cancer. I want to compare the cox model for each of the event of interest
with a competing risk model. In the competing risk model the cumulative
incidence function is used directly. I used the jackknife (pseudovalue) of
the cumulative incidence function for each cause (AIDS or Cancer) in a
generalized estimating equation. I used the following code:

>library(pseudo)
>cutoffs<-c(0.5,5.5,10.5,15.5)
>pseudo<-pseudoci(time,event,tmax=cutoffs)
>library(geepack)
>fit<-geese(pseudo~as.factor(Age)+as.factor(Sex),data=b,id=id, jack = TRUE,
scale.fix=TRUE,
mean.link = "logit", corstr="independence").



I want to know whether I am doing the right thing?

	[[alternative HTML version deleted]]

