Bookmark and Share

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


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

st: Paste values of scalar in observation (bootstrap)


From   Camila Mendes <[email protected]>
To   [email protected]
Subject   st: Paste values of scalar in observation (bootstrap)
Date   Wed, 16 Oct 2013 10:02:05 -0300

Hi everyone:

I need help to paste bootstraped coefficients and their lower and
upper 95% CI in new variables to produce graphs.

I tried the code below. It looks like scalars are retrieved with
sucess, but I am not able to paste them into specific observations
when using replace.

Can anyone help me?

Thanks,

Camila

---- begin of do-file ----
clear all
forvalues j=0/1{
foreach l in P2 P3 P4 {
use `l'_per_`j', clear
set seed 1
reg estimate xaxis
predict yhat
gen yhat_bs=yhat
gen ll_95per=.
gen ul_95per=.
local N = _N
foreach i in 1/11 { //number of time periods to be predicted
bootstrap pred=(_b[_cons] + _b[xaxis]*`i'), reps(1000) seed(1): reg
estimate xaxis
matrix b=e(b)
scalar sb=b[1,1]
matrix ci_normal=e(ci_normal)
local par=b[1,1]
local ll=ci_normal[1,1]
local ul=ci_normal[2,1]
replace yhat_bs=`par' if xaxis==`i'
replace ll_95per=`ll' if xaxis==`i'
replace ul_95per=`ul' if xaxis==`i'
}
save `l'_per_`j'_lin_trend, replace
}
}
---- end of do-file ----
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/


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