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]
Re: st: Can -bootstrap- handle vectors
From
Alex Olssen <[email protected]>
To
[email protected]
Subject
Re: st: Can -bootstrap- handle vectors
Date
Sat, 13 Aug 2011 22:35:45 +1200
Helps lots! Thanks!
On 13 August 2011 21:12, Matthew White <[email protected]> wrote:
> 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 <[email protected]> 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
> [email protected]
>
> *
> * 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/
>
*
* 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/