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/