Hello,
I am trying to create 10 datasets of simulated survival data and store the estimates in one table.
I have created the following code:
foreach sim of num 1/10 {
clear
set obs 1000
gen trt=rbinomial(1,0.3)
survsim time state event , hazard1(dist(weibull) lambda(0.1) gamma(1.6) covariates (trt 0.7)) hazard2(dist(weibull) lambda(0.1) gamma(1.6) covariates(trt 0.7)) maxtime(5)
stset time1, failure (state1==2)
eststo stcox1_`sim': stcox i.trt, nohr
eststo stccreg1_`sim': stcrreg i.trt, compete(state1==3) noshr
stset time1, failure (state1==3)
eststo stcox2_`sim': stcox i.trt, nohr
eststo stcrreg2_`sim': stcrreg i.trt, compete(state1==2) noshr
}
esttab _all using estimates, p scalars(N_compete N_censor N_fail)
There are 4 regression per dataset. In the end I need to store 40 estimates and their mean survival times as well. However, I cannot store mean survival times because p50 is an r-class variable and esttab only produces e-class variables.
Also when I run the above code, it only saves the estimate for the last data set.
Any help would be greatly appreciated!
Thank you.
Related Posts with Survival Analysis: help with survsim and creating a loop/program to save estimates into a table
Margins command for FE Panel Data with interaction between a quadratic variable and a factor variableHi, I am running a FE Panel regression. My code looks as follows. I would like to know how I can run…
Displaying rho after biprobit, margins in the output tableDear Statalist members, I would be very thankful for an advice. I am running a biprobit model, and u…
Oaxaca-Blinder Decomposition for ordinal dependent variablesDear Statalist, I am trying to find a way to run a detailed Oaxaca-Blinder decomposition for an ordi…
Cannot locate set preferencesI accidentally undid my windows preferences and have been able to restore them, but I can't save the…
Saving regression resultsI want to run regressions using one variable at a time and store results for each of those. The foll…
Subscribe to:
Post Comments (Atom)
0 Response to Survival Analysis: help with survsim and creating a loop/program to save estimates into a table
Post a Comment