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 on April 23, and its replacement, statalist.org is already up and running.


[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 <njcoxstata@gmail.com>
To   "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu>
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
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