Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down at the end of May, and its replacement, statalist.org is already up and running.


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

Re: programming syntax problem [was: Re: st: From: "Christopher L. Aberson" ...]


From   Nick Cox <njcoxstata@gmail.com>
To   "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu>
Subject   Re: programming syntax problem [was: Re: st: From: "Christopher L. Aberson" ...]
Date   Mon, 17 Jun 2013 22:37:56 +0100

You are using -simulate- to call your -powersim3-, but not feeding it anything.
Nick
njcoxstata@gmail.com


On 17 June 2013 21:56, Christopher L. Aberson
<Christopher.Aberson@humboldt.edu> wrote:
> Thank you for the feedback. The approaches suggested made the code
> more elegant but did not solve the core issue.
>
> To recap, I wand to run the code noted below (encased in *****) using
> this line powersim3, ry1(.30) ry2(.30) r12 (.30) n (282)
> that defines the required characteristics.
>
> It appears the issue is that the simulate command isn't recognizing
> the powersim3 command above. I expect there is a simple solution here
> -- but I'm not proficient enough in Stata to know what the simple
> answer is. Any additional suggestions would be very helpful.
>
> ******
>  program define powersim3, rclass
>  syntax, ry1(real)ry2(real)r12(real)n(
> real) [my(real 0)m1(real
> 0)m2(real 0)m3(real 0)sy(real 1) ///
>     s1(real 1)s2(real 1)s3(real 1)alpha(real .05)help]
> matrix m = (`my', `m1', `m2')
> matrix sd = (`sy', `s1', `s2')
> matrix r = (1, `ry1', `ry2'\ `ry1' , 1 , `r12'\ `ry2' , `r12', 1)
> drawnorm y x1 x2, n(`n') corr(r) means(m) sds(sd)
> regress y x1 x2
> end
>
> simulate _b _se, reps(10000): powersim3
> *****
>
>
>> Please use meaningful titles for your postings.
>>
>> Try editing out the spaces between option names and option details in
>> the cases of the options
>>
>> ry1 (real) ry2 (real) r12 (real) sy (real 1)   s1 (real 1) s2 (real 1)
>> s3 (real 1)
>>
>> Incidentally, your program does nothing with the options
>>
>> alpha(real .05) help
>>
>> although you may be planning to add details later.
>>
>> Nick
>> njcoxstata@gmail.com
>>
>>
>> On 14 June 2013 23:18,  <owner-statalist@hsphsun2.harvard.edu> wrote:
>>> I've written a tiny program for simulation of multiple regression
>>> results.  I can run this on a line by line basis (using actual values
>>> instead of the syntax command), but when I put it in an ado, it
>>> doesn't work.
>>>
>>> The basic idea is use drawnorm to generate a dataset, run regression
>>> on it, then do it over and over using the simulate option.
>>>
>>>
>>> program define powersim3, rclass
>>> syntax, ry1 (real) ry2 (real) r12 (real) [my(real 0) m1(real 0)
>>> m2(real 0) m3(real 0)] [sy (real 1) ///
>>>     s1 (real 1) s2 (real 1) s3 (real 1)] n(real) [alpha(real .05) help ]
>>> matrix m = (`my', `m1', `m2')
>>> matrix sd = (`sy', `s1', `s2')
>>> matrix r = (1, `ry1', `ry2'\ `ry1' , 1 , `r12'\ `ry2' , `r12', 1)
>>> drawnorm y x1 x2, n(`n') corr(r) means(m) sds(sd)
>>> regress y x1 x2
>>> end
>>>
>>> simulate _b _se, reps(10000): powersim3
>>>
>>> .powersim3, ry1(.30) ry2(.30) r12 (.30) n (282)
>>>
>>> This yields an error saying "ry1" required. If I take out the simulate
>>> command I get a single regression but not the repeated
>>> sampling/simulation I want.
*
*   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–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index