Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: bootstrap reject()

From   Shehzad Ali <>
To   Stata List <>
Subject   st: bootstrap reject()
Date   Tue, 1 Nov 2011 11:59:29 +0000 (GMT)

Dear all,

A colleague asked this question. Is there a way to specify a range (or rejection region) for the parameter of interest when using the -bootstrap- program in Stata? She is using GLM for her cost data, followed by predicting costs for the treatment and control arms, then taking the difference between mean costs in treatment and control arms. The standard error of the difference in mean costs is then estimated by bootstrapping the process. In order to execute this, she is using the following general approach:

capture program drop bootstrap1

glm model followed by predicting costs for treatment=1 and treatment=0. Then the difference between mean predicted costs for the two groups is saved as a scalar 'costdiff'.

Subsequently, Stata's -bootstrap- is used to repeat this process a 1,000 times.

bootstrap costdiff=r(costdiff), reps(1000) saving(bootrep1, replace) seed(1234): 

However, some of the predicted cost differences are very huge (>5% of the bootstrap samples) and the SE is also huge (several times the SE in the observed sample). Is there a way, the predicted costs could be restricted within a range (for instance, the cost difference should be within a range that is close to the observed sample) or would this approach be incorrect?

Thank you

*   For searches and help try:

© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index