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]

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/

- Prev by Date:
**Re: st: converting a SAS program to Stata code for use with the HCUP NIS** - Next by Date:
**Re: st: Saving bootstrap samples** - Previous by thread:
**programming syntax problem [was: Re: st: From: "Christopher L. Aberson" ...]** - Next by thread:
**Re: programming syntax problem [was: Re: st: From: "Christopher L. Aberson" ...]** - Index(es):