I've got data on vials of drosophila infected with a virus. Then at each time point, I record how many are still alive. I would like to do an interval censored survival analysis to determine how many are alive at the end.
----------------------- copy starting from the next line -----------------------
Code:
* Example generated by -dataex-. For more info, type help dataex clear input byte(STRAIN infected startingsample time1 time2 time3 time4 time5 time6) str1 gender byte age 1 1 20 19 19 18 18 18 15 "F" 2 2 0 15 14 14 14 14 12 12 "M" 4 3 1 30 29 29 29 29 29 29 "F" 3 4 1 15 15 15 15 13 12 10 "F" 3 5 0 10 10 9 9 8 7 7 "F" 3 6 1 20 17 16 16 16 14 14 "M" 4 7 0 9 8 5 3 0 0 0 "F" 2 8 0 12 10 10 10 10 9 7 "M" 3 end
I used this code below
[/CODE]
stintreg i.infected gender , interval(t6 startingsample) distribution(weibull)
This code doesn't allow me to enter all the timepoints in the interval. Do I only need to enter the last sample recorded regardless of other time points and the starting sample? Thank you in advance.
0 Response to Interval censored survival using aggregate data
Post a Comment