Looking at cumulative incidence of CMV and EBV infections at 12 months. I can compile a plot of cuminc based on first event to occur, but I have two patients who had BOTH CMV and EBV events during the 0-12 month period. How do I account for those?
I have been using x <- cuminc(timetoevent, status) then plot(x). Here my status is coded as 0 for no event, 1 for CMV, or 2 for EBV.