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
Gen. new variables based on values for multiple existing variables, TSCS dataDear members, I have put together a data-set with information about military alliances btw states i…
Panel-Data Logistic regression Probability of LR Test is 1Dear Statalist, I have a question concerning the a logistic regression with panel data. I have a pa…
Counts for multiple variablesHello, I am analyzing a hospital dataset which provides demographic, diagnostic, and procedural data…
Drop observations based on value of previous observationHi all, I have a panel dataset on productivity. The variables are as follows: activity = the activi…
Epimodels package for Stata has been updated at SSCThanks to Kit Baum the epimodels package for Stata has been updated at SSC. Epimodels (written by S…
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