hello there,
I am trying to replicate this graph published here:

My dilemma is this the correct way on calculating the 10 year cumulative revision rate?
Revised = procedure has failed, where 1= revised
survival is the survival time variable calculated from when the procedure was done to the time the procedure failed
volgrp = three level categorical variable as <10; 10-30 or >30
I wonder if I am being too simple, and perhaps should be looking at another way of calcualting the cumulative revision rate using survival stats.
I am trying to replicate this graph published here:
My dilemma is this the correct way on calculating the 10 year cumulative revision rate?
Revised = procedure has failed, where 1= revised
survival is the survival time variable calculated from when the procedure was done to the time the procedure failed
volgrp = three level categorical variable as <10; 10-30 or >30
Code:
gen revised10=0 replace revised10 ==1 if revised ==1 and survival <=1 keep if revised10==1 graph bar(percent), over(cemented_cemntless) over(volgrp) asyvars blabel(bar, size(tiny) format(%9.1f)) ytitle("%Percentage")