Hi,
I am running the following program to obtain estimates with bias corrected bootstrapping. I get the estimate for x1 with the following program, but how to obtain the estimates for other variables (x2-x8) in the model? especially the i.x8, which is a factor variable with 8 categories (first as reference).
program myprog , rclass
mi estimate, saving(miest1, replace) eform post: glm y x1 x2 x3 x4 x5 x6 x7 i.x8 c1 c2, fam(poisson) link(log) nolog vce(robust)
matrix bb_x1= e(b_mi)
scalar b_x1=(bb_x1[1,1])
return scalar b_x1=x1[1,1]
bootstrap exp(r(b_x1)) , seed(12345) reps(50): myprog
estat bootstrap, all
Any advice is highly appreciated.
Best wishes,
Massao
I am running the following program to obtain estimates with bias corrected bootstrapping. I get the estimate for x1 with the following program, but how to obtain the estimates for other variables (x2-x8) in the model? especially the i.x8, which is a factor variable with 8 categories (first as reference).
program myprog , rclass
mi estimate, saving(miest1, replace) eform post: glm y x1 x2 x3 x4 x5 x6 x7 i.x8 c1 c2, fam(poisson) link(log) nolog vce(robust)
matrix bb_x1= e(b_mi)
scalar b_x1=(bb_x1[1,1])
return scalar b_x1=x1[1,1]
bootstrap exp(r(b_x1)) , seed(12345) reps(50): myprog
estat bootstrap, all
Any advice is highly appreciated.
Best wishes,
Massao
