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]

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


From   Nick Cox <[email protected]>
To   "[email protected]" <[email protected]>
Subject   programming syntax problem [was: Re: st: From: "Christopher L. Aberson" ...]
Date   Sat, 15 Jun 2013 00:48:56 +0100

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
[email protected]


On 14 June 2013 23:18,  <[email protected]> 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–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index