The name is absent



24


Stata Technical Bulletin


STB-4


display (_b [XX] -7._orb) / (_b [XX] ∕sqrt (.result (6) ) )

* Previous four lines display the original-sample regression slope

* (.orb), original-sample standard error (.orSE), bootstrap

* regression slope, and bootstrap studentized slope,
log off

macro define .bsample=%.bsample+l
>

log close
drop .all
infile orb orSE bootb stub using bstemp.log
save mcit%.mcit, replace

* Each Monte Carlo iteration saves a dataset containing

* results from 2,000 resamplings,
macro define .mcit=%.mcit+l

>

macro define .mcit=l

capture erase mctemp.log

log using mctemp.log

log off

while %.mcit<101 -(

quietly use mcit%.mcit, clear
log on

display %.mcit

quietly summ orb

display .result(3)

quietly summ orSE

display .result(3)

quietly summ bootb, detail

display .result(3)

display sqrt(.result(4))

display .result(7)

display .result(13)

quietly summ stub, detail

display .result(7)

display .result(13)

display
log off

* Preceding lines display the Monte Carlo sample number, original-sample

* slope and standard error, bootstrap mean and standard deviation of

* regression slope; and bootstrap percentiles and studentized percentiles,
macro define .mcit=%.mcit+l

>

log close

drop .all

infile iterate orb orSE bootb bootSE p05 p95 t05 t95 ∕*

*/ using mctemp.log;

label variable orb "original sample b coefficient"

label variable orSE "original sample SE of b"

label variable bootb "mean bootstrap b coefficient"

label variable bootSE "SD of bootstrap b coefficient"

label define correct 0 "wrong" 1 "correct"
*

generate stanlo=orb-l.665*orSE

generate stanhi=orb+l.665*orSE

generate Stanwide=Stanhi-Stanlo

label variable Stanwide "width standard t interval"

generate Stancor=O

replace Stancor=I if stanlo<3 & 3<stanhi

label variable Stancor "standard 90% c.i. correct?"
label values Stancor correct

*

generate perwide=p95-p05

label variable perwide "width 5th-95th percentile"

generate percor=0

replace percor=l if p05<3 & 3<p95

label variable percor "percentile 90% c.i. correct?"
label values percor correct

*

generate hyblo=orb-(p95-orb)

generate hybhi=orb-(p05-orb)

generate hybwide=hybhi-hyblo

label variable hybwide "width hybrid percentile interval"

generate hybcor=0

replace hybcor=l if hyblo<3 & 3<hybhi



More intriguing information

1. The name is absent
2. Output Effects of Agri-environmental Programs of the EU
3. Ruptures in the probability scale. Calculation of ruptures’ values
4. The name is absent
5. Analyzing the Agricultural Trade Impacts of the Canada-Chile Free Trade Agreement
6. Placenta ingestion by rats enhances y- and n-opioid antinociception, but suppresses A-opioid antinociception
7. Olive Tree Farming in Jaen: Situation With the New Cap and Comparison With the Province Income Per Capita.
8. The name is absent
9. Foreign Direct Investment and Unequal Regional Economic Growth in China
10. The name is absent
11. Modelling the health related benefits of environmental policies - a CGE analysis for the eu countries with gem-e3
12. The urban sprawl dynamics: does a neural network understand the spatial logic better than a cellular automata?
13. Moffett and rhetoric
14. 5th and 8th grade pupils’ and teachers’ perceptions of the relationships between teaching methods, classroom ethos, and positive affective attitudes towards learning mathematics in Japan
15. The name is absent
16. Exchange Rate Uncertainty and Trade Growth - A Comparison of Linear and Nonlinear (Forecasting) Models
17. Automatic Dream Sentiment Analysis
18. Equity Markets and Economic Development: What Do We Know
19. The name is absent
20. NEW DEVELOPMENTS IN FARM PRICE AND INCOME POLICY PROGRAMS: PART I. SITUATION AND PROBLEM