Re: st: Can -bootstrap- handle vectors

[Matthew White <mwhite@poverty-action.org>]

I don't know of a way to use Mata to put numeric values in a local.
st_local only works with string values, I think. But you could return
each element of a Mata vector to separate Stata scalars. Say X is a
Mata column vector:
mata: st_numscalar("r(rows)", rows(X))
forvalues row = 1/`r(rows)' {
    mata: st_numscalar("r(v`row')", X[`row', 1])
}

Hope that helps!

On Sat, Aug 13, 2011 at 11:40 AM, Alex Olssen <alex.olssen@gmail.com> wrote:
> Hi Statalisters,
>
> It seems to me that -bootstrap- cannot handle vectors?  Is this
> correct?  If so, does anyone know if there is a way to use a loop to
> send each element of a Mata vector to a separate Stata local so I can
> bootstrap the vector element by element?
>
> Sample code is below.  All help greatly appreciated!
>
> ************************************************************
> capture program drop my_vector
> program define my_vector, rclass
>        syntax varlist
>        mata: m_my_vector("`varlist'")
>        return list
>        return add
> end
>
> capture mata mata drop m_my_vector()
> mata:
> void m_my_vector(string scalar varlist)
> {
>        real matrix data, data_cov
>        st_view(data, ., varlist)
>        data_cov = variance(data)
>        st_matrix("r(v1)", data_cov[., 1])
>        st_numscalar("r(v11)", data_cov[1, 1])
> }
> end
>
> sysuse auto, clear
> my_vector price length weight
> di r(v11)
> matrix list r(v1)
>
> bootstrap r(v11), reps(200) seed(13082011): my_vector price length weight
> bootstrap r(v1), reps(200) seed(13082011): my_vector price length weight
> ************************************************************
>
> Cheers,
> Alex Olssen
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>



-- 
Matthew White
Evaluation Coordinator
Urban Micro-Insurance Project
Innovations for Poverty Action

+254 (0)701 025 276
mwhite@poverty-action.org

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/