Statalist

test post

2006-05-11T12:36:00.000-04:00

Testing

Re: st: RE: RE: RE: list subjects with a similar value

2006-03-22T20:44:00.000-05:00

>Personally if you want to be able to revert to older versions of files >I'd recommend simply creating a copy before doing major revisions and >simply append a date in numeric format at some point to the filename >(before '.' would be most appropriate).

Thank you Neil. Are there resources for systematic and efficient "best practices" in do-file colding, that beginning Stata users can emulate and integrate into their own work? Michael

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: simultaneous equation with qualitative variables

2006-03-22T20:43:00.000-05:00

how can i use stata to estimate simultaneous equations with qualitative variables as follows: y1=a1+b1y2+c1x1+e1 y2=a2+b2y1+c2x2+e2 y1*=1 if y1>0 y1*=0 if y1<=0 y2*=1 if y2>0 y2*=0 if y2<=0

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

RE: st: RE: RE: RE: list subjects with a similar value

2006-03-22T20:23:00.000-05:00

MS-Word I understand to be a word processor. Recall that many members of Statalist do not use and certainly are not expert in Windows or Microsoft products generally.

In Vim, the way that seems to be most natural is a Unix way. You can have two files open in two windows and set it so that differences are highlighted. So, one could be a previous version and the other a working version.

I gather that Word behaves differently.

In general, good text editors will have something loosely similar. None that I know of regards it as a virtue to emulate Word.

Nick n.j.cox@durham.ac.uk

Michael McCulloch > Sorry; I was referring to changes in coding that one writes > in a do-file. > >Changes to what?

> > > Jennifer response brings to mind a question that recently > > > occurred to me: > > > Is there a Stata-compatible text editor that, like MS-WORD, > > > can highlight changes? > > > > > >Whoops, of course you're right. When I play with things to > > > work out code > > > >I usually don't keep the changes, so I mangle working datasets > > > >willy-nilly and didn't think to change the conditioned keep to a > > > >conditioned list. > > > >The use of bysorting and _N is much neater and more flexible.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

Re: st: RE: RE: RE: list subjects with a similar value

2006-03-22T20:16:00.000-05:00

On 3/23/06, Michael McCulloch <mm@pinest.org> wrote: > Jennifer response brings to mind a question that recently occurred to me: > Is there a Stata-compatible text editor that, like MS-WORD, can highlight > changes?

Although I don't use it myself there are version control features in Emacs. These work with programming control systems such as CVS or SVN, and I'm not sure what your mileage would be writing do/ado files under such schemas, but I suspect it is possible (don't quote me on that though :).

This isn't exactly the same as M$-words document tracking (which personally I find hideously hard to follow, particularly when there are multiple authors making revisions, I've seen some docs that end up looking like the old TV-test screens :-), but it does allow you to track the changes that you are making. One of the main problems (as I see it) is that to write do/ado-files you need a _text_ editor, and word is not a text-editor, but a word-processor, so all the colour changes that you see are essentially mark-ups of the original text, and such mark-ups would render the Stata code uninterpretable.

See http://www.gnu.org/software/emacs/manual/html_node/Version-Control.html for more on Emacs' VC system.

Personally if you want to be able to revert to older versions of files I'd recommend simply creating a copy before doing major revisions and simply append a date in numeric format at some point to the filename (before '.' would be most appropriate).

HTH's

Neil -- "The best safety device in climbing is always situated between your ears" - Ross Weiter, Perth Rock Climbing Guide (2002)

Email - nshephard@gmail.com / neilshep@cyllene.uwa.edu.au Website - http://slack.ser.man.ac.uk/ Blog - http://slack---line.blogspot.com/ Flickr - http://www.flickr.com/photos/slackline/

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

RE: st: RE: RE: RE: list subjects with a similar value

2006-03-22T20:05:00.000-05:00

Sorry; I was referring to changes in coding that one writes in a do-file.

At 04:47 PM 3/22/2006, you wrote: >Changes to what? > >Nick >n.j.cox@durham.ac.uk > >Michael McCulloch > > > Jennifer response brings to mind a question that recently > > occurred to me: > > Is there a Stata-compatible text editor that, like MS-WORD, > > can highlight changes? > > > >Whoops, of course you're right. When I play with things to > > work out code > > >I usually don't keep the changes, so I mangle working datasets > > >willy-nilly and didn't think to change the conditioned keep to a > > >conditioned list. > > >The use of bysorting and _N is much neater and more flexible. > >* >* For searches and help try: >* http://www.stata.com/support/faqs/res/findit.html >* http://www.stata.com/support/statalist/faq >* http://www.ats.ucla.edu/stat/stata/

Best wishes, Michael

____________________________________

Michael McCulloch Pine Street Clinic 124 Pine Street, San Anselmo, CA 94960-2674 tel 415.407.1357 fax 415.485.1065 email: mm@pinest.org web: www.pinest.org www.pinestreetfoundation.org

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

RE: st: RE: RE: RE: list subjects with a similar value

2006-03-22T19:49:00.000-05:00

Changes to what?

Nick n.j.cox@durham.ac.uk

Michael McCulloch > Jennifer response brings to mind a question that recently > occurred to me: > Is there a Stata-compatible text editor that, like MS-WORD, > can highlight changes? > >Whoops, of course you're right. When I play with things to > work out code > >I usually don't keep the changes, so I mangle working datasets > >willy-nilly and didn't think to change the conditioned keep to a > >conditioned list. > >The use of bysorting and _N is much neater and more flexible.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: RE: RE: RE: Macro display format

2006-03-22T18:11:00.000-05:00

It's a good question.

I guess: -scatter- calls -graph- calls ... something that clears your r-class results.

-graph- needs all sorts of little calculations to work out what to show. It tends to do this on the fly, but either way, you can lose your r-class results.

Whatever it is, it is quite deep down.

Nick n.j.cox@durham.ac.uk

Alex Ogan

> Sorry if this is a silly question. > > I did the following sequence of commands. I closed the > scatter as soon > as it opened. No other commands. > > Why does r(mean) go away after you use it in the scatter with the > formatting extended function? > > . sysuse auto, clear > (1978 Automobile Data) > > . quietly summ length > > . di `r(mean)' > 187.93243 > > . di `: di %9.1f `r(mean)'' > 187.9 > > . di `: di %9.1f `r(mean)'' > 187.9 > > . twoway scatter length mpg, text(200 35 "The mean of length is:`: di > %2.1f `r(mean)''") > > . di `: di %9.1f `r(mean)'' > > > . di `r(mean)' > > > . >

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: "sureg" with long-format data

2006-03-22T18:04:00.000-05:00

Hello all,

I have searched the web and available Stata faqs and helpfiles for an answer to an sureg question and am unable to find one on my own. Much obliged if anyone can help me through this particular puzzle.

I have a dataset oriented lengthwise where observations are grouped by country and year. I would like to run sureg using this vertical orientation. For example, in the sureg syntax, I would like to write "sureg (yvar xvar)" as equation 1 representing the first 35 observations, followed by "(yvar xvar)" as equation 2 representing the next 35 observations, and so on. My difficulty is that without re-orienting the data width-wise (which creates a cumbersome heap of newly-named X and Y variables and stops me from being able to tweak the regressions around the edges) I have no way to distinguish for Stata that instead of each equation having different Y and X variable names, each equation has the same variable names but should contain a different set of 35 observations. My abbreviated dataset looks like:

Country Year Yvar Xvar Australia 2004 .5 14 Australia 2003 .6 17 .... .. Australia 1970 .2 10

Austria 2004 .9 35 Austria 2003 .8 14 .... .. Austria 1970 .15 14

Belgium 2004 .13 7 .... (and so on for approximately 100 countries)

Ideally, the sureg command would allow something like "by country: sureg (Yvar Xvar)" but instead, I have to reorient the data widthwise and write "sureg (YvarAustralia XvarAustralia) (YvarAustria XvarAustria)..." for about one-hundred different equations (and change all 100 when I make adjustments). Requiring this horizontal orientation to make sureg work seems remarkably inefficient in the command line, but I can't seem to find another way around it. Many thanks for any and all ideas.

Cheers, Eric

Eric K. Bielke Regulatory Economics Advisor www.telecom.co.nz Level 2, Telecom House, 68-86 Jervois Quay, P O Box 570, Wellington, New Zealand

This communication, including any attachments, is confidential. If you are not the intended recipient, you should not read it - please contact me immediately, destroy it, and do not copy or use any part of this communication or disclose anything about it. Thank you. Please note that this communication does not designate an information system for the purposes of the Electronic Transactions Act 2002.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: RE: RE: Macro display format

2006-03-22T17:59:00.000-05:00

Sorry if this is a silly question.

I did the following sequence of commands. I closed the scatter as soon as it opened. No other commands.

Why does r(mean) go away after you use it in the scatter with the formatting extended function?

. sysuse auto, clear (1978 Automobile Data)

. quietly summ length

. di `r(mean)' 187.93243

. di `: di %9.1f `r(mean)'' 187.9

. twoway scatter length mpg, text(200 35 "The mean of length is:`: di %2.1f `r(mean)''")

. di `: di %9.1f `r(mean)''

. di `r(mean)'

-----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Nick Cox Sent: Wednesday, March 22, 2006 5:51 PM To: statalist@hsphsun2.harvard.edu Subject: st: RE: Macro display format

No. You should try

` : di %9.1f `r(mean)''

See help extended_fcn.

Tip: I would go `: di %2.1f `r(mean)'' even if you are sure that format is too restrictive. You're likely to be wrong, as Stata will stretch the space to avoid damage. However, with %9.1f you are likely to get the ugly spaces that are a consequence of what you asked for.

Nick n.j.cox@durham.ac.uk

Thomas Speidel > I am trying to include the content of a macro within a graph, but I'm > having problems with the display format. > > For example: > > sysuse auto, clear > qui: summ length > twoway scatter length mpg, text(200 35 "The mean of length is: > `r(mean)'") > > How do I change the format of the macro to display something > like %9.2f? > > I tried: > ... , text(200 35 "The mean of length is: `%9.1f `r(mean)''") > > Am I missing some triple compound quote? :-)

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

This message is intended solely for the designated recipient(s). It may contain confidential or proprietary information and may be subject to confidentiality protections. If you are not a designated recipient, you may not review, copy, or distribute this message. If you receive this in error, please notify the sender by reply e-mail and delete this message.

Arrowstreet Capital, L.P. is an Equal Opportunity Employer.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: RE: Macro display format

2006-03-22T17:52:00.000-05:00

No. You should try

` : di %9.1f `r(mean)''

See help extended_fcn.

Nick n.j.cox@durham.ac.uk

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: RE: RE: RE: list subjects with a similar value

2006-03-22T17:48:00.000-05:00

Whoops, of course you're right. When I play with things to work out code I usually don't keep the changes, so I mangle working datasets willy-nilly and didn't think to change the conditioned keep to a conditioned list. The use of bysorting and _N is much neater and more flexible.

N.J.Cox wrote: >Whoa! The question was just about _listing_. >You just changed Raoul's dataset by throwing >much of it away.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: Macro display format

2006-03-22T17:43:00.000-05:00

I am trying to include the content of a macro within a graph, but I'm having problems with the display format.

For example:

sysuse auto, clear qui: summ length twoway scatter length mpg, text(200 35 "The mean of length is: `r(mean)'")

How do I change the format of the macro to display something like %9.2f?

I tried: ... , text(200 35 "The mean of length is: `%9.1f `r(mean)''")

Am I missing some triple compound quote? :-)

Thanks, Thomas

-- Thomas Speidel Statistical Associate Clinical Trials Unit Tom Baker Cancer Centre 1331 - 29th Street N.W. Calgary, AB, T2N 4N4

Tel. (403) 521-3370 Email: thomassp@cancerboard.ab.ca

This e-mail and any attachments may contain confidential and privileged information. If you are not the intended recipient, please notify the sender immediately by return e-mail, delete this e-mail and destroy any copies. Any dissemination or use of this information by a person other than the intended recipient is unauthorized and may be illegal.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

RE: st: list subjects with a similar value

2006-03-22T17:38:00.000-05:00

You are avoiding the command -duplicates- and doing it from first principles. That is a very good idea. -duplicates- is just a wrapper for stuff like this.

But the three steps here can be cut to two.

bysort date_of_birth : gen dob_duplicate = _N list id date_of_birth if dob_duplicate >= 2

Nick n.j.cox@durham.ac.uk

clinton.thompson@summitllc.us > there may be a more elegant way to do this, albeit this is > but one attempt: > > * obtain the number of duplicates w/in date of birth > bysort date_of_birth: gen dob_duplicate = _N > * tag each combination of DOB & the duplicates therein > egen dob_tag = tag(date_of_birth dob_duplicate) > * list the ID & DOB associated w/ each repeated DOB... > list id date_of_birth dob_tag if dob_tag > > I have a large database and would like to list the idnumber of all > > subjects with the same date of birth. How do I do this? I have tried > > .duplicate, but can figure out how to do it. Thanks.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: RE: RE: list subjects with a similar value

2006-03-22T17:35:00.000-05:00

Whoa! The question was just about _listing_. You just changed Raoul's dataset by throwing much of it away.

duplicates tag dateofbirth, gen(tag) sort dateofbirth list id dateofbirth if tag

Nick n.j.cox@durham.ac.uk

Marino, Jennifer > Maybe try something along these lines: > > duplicates tag dateofbirth, gen(tag) > drop if tag==0 > sort dateofbirth > by dateofbirth: list id

Raoul C Reulen > I have a large database and would like to list the idnumber of all > subjects with the same date of birth. How do I do this? I have tried > .duplicate, but can figure out how to do it. Thanks.

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

Re: st: list subjects with a similar value

2006-03-22T17:30:00.000-05:00

there may be a more elegant way to do this, albeit this is but one attempt:

* obtain the number of duplicates w/in date of birth bysort date_of_birth: gen dob_duplicate = _N * tag each combination of DOB & the duplicates therein egen dob_tag = tag(date_of_birth dob_duplicate) * list the ID & DOB associated w/ each repeated DOB... list id date_of_birth dob_tag if dob_tag

note that i didn't subject this to a rigorous test...but i think it works, nonetheless. --clint

> Hi, > > I have a large database and would like to list the idnumber of all > subjects with the same date of birth. How do I do this? I have tried > .duplicate, but can figure out how to do it. Thanks. > > Raoul > > * > * For searches and help try: > * http://www.stata.com/support/faqs/res/findit.html > * http://www.stata.com/support/statalist/faq > * http://www.ats.ucla.edu/stat/stata/ >

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: RE: list subjects with a similar value

2006-03-22T17:25:00.000-05:00

Maybe try something along these lines:

duplicates tag dateofbirth, gen(tag) drop if tag==0 sort dateofbirth by dateofbirth: list id

Jen Marino

-----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Raoul C Reulen Sent: Wednesday, March 22, 2006 2:14 PM To: statalist@hsphsun2.harvard.edu Subject: st: list subjects with a similar value

Hi, I have a large database and would like to list the idnumber of all subjects with the same date of birth. How do I do this? I have tried .duplicate, but can figure out how to do it. Thanks. Raoul

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: list subjects with a similar value

2006-03-22T17:15:00.000-05:00

Hi, I have a large database and would like to list the idnumber of all subjects with the same date of birth. How do I do this? I have tried .duplicate, but can figure out how to do it. Thanks. Raoul

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

Re: st: monte carlo study

2006-03-22T15:01:00.000-05:00

Rudy,

here is an example based on the help file for -simulate-

------------------------- capture program drop mcarlols program define mcarlols, rclass syntax [, obs(integer 1) ] drop _all set obs `obs' drawnorm eps x g y=1+x+eps reg y x return scalar beta = _b[x] end

simulate beta=r(beta) , reps(1000): mcarlols, obs(100)

su beta, de kdensity beta, norm --------------------------------

le 22/03/2006 20:10, Rudy Fichtenbaum a ecrit : > Stata Users: > > I am still learning my way around Stata after many years of using SAS. > In SAS it is fairly easy to write a program to do a simple Monte Carlo > study to illustrate the properties of least squares estimators. > > Is there anyone that has a simple example of a Monte Carlo Study for OLS? > > Thanks, > > Rudy > * > * For searches and help try: > * http://www.stata.com/support/faqs/res/findit.html > * http://www.stata.com/support/statalist/faq > * http://www.ats.ucla.edu/stat/stata/ >

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

Re: st: monte carlo study

2006-03-22T14:54:00.000-05:00

In the unlikely event that you have not, do check out -findit monte carlo-. The example in STB reprints Vol 4, pg 207 may be helpful

TITLE STB-20 ssi6. Simplified Monte Carlo simulations.

Rudy Fichtenbaum wrote: > Stata Users: > > I am still learning my way around Stata after many years of using SAS. > In SAS it is fairly easy to write a program to do a simple Monte Carlo > study to illustrate the properties of least squares estimators. > > Is there anyone that has a simple example of a Monte Carlo Study for OLS? > > Thanks, > > Rudy > * > * For searches and help try: > * http://www.stata.com/support/faqs/res/findit.html > * http://www.stata.com/support/statalist/faq > * http://www.ats.ucla.edu/stat/stata/ > * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: Chou-Talalay method

2006-03-22T14:16:00.000-05:00

Dear all,

Does anyone have or know of a Stata program to compute the various Chou-Talalay statistics of dose-effect relationship?

Thank you, Ricardo

Ricardo Ovaldia, MS Statistician Oklahoma City, OK

__________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

st: monte carlo study

2006-03-22T14:11:00.000-05:00

Stata Users:

I am still learning my way around Stata after many years of using SAS. In SAS it is fairly easy to write a program to do a simple Monte Carlo study to illustrate the properties of least squares estimators.

Is there anyone that has a simple example of a Monte Carlo Study for OLS?

Thanks,

Rudy * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

RE: st: functions for computing prob

2006-03-22T13:58:00.000-05:00

Thanks a lot, Gary!

Lei

-----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Gary Longton Sent: Wednesday, March 22, 2006 12:47 PM To: statalist@hsphsun2.harvard.edu Subject: Re: st: functions for computing prob

It appears that I was mistaken about -tprob()-, that it does exist and appear to work, though undocumented and apparently obsolete. It looks to be the 2-tailed version of -ttail()-

. di tprob(30,1.8) .08192507

. di ttail(30,1.8) .04096253

- GL

Gary Longton wrote: > Lei Xuan wrote: > >> I am computing probabilities for t-test and z-test. >> I want to know if the functions -tprob- and -normprob- are out-of-date >> since no help files explain these functions. Are -ttail(n,t)- and >> -normal(z) >> >> right functions to compute probs? > > > see -help density functions- > or -help functions- > > Yes, the cumulative normal density function, normprob(), still works but > is out of date, no longer documented, and has been replaced by > normal(). Am not sure whether tprob() ever existed? - doesn't seem to > work and is not documented in any case. ttail(n,t) is documented under > help for functions. > > - Gary * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

Re: st: functions for computing prob

2006-03-22T13:48:00.000-05:00

It appears that I was mistaken about -tprob()-, that it does exist and appear to work, though undocumented and apparently obsolete. It looks to be the 2-tailed version of -ttail()-

. di tprob(30,1.8) .08192507

. di ttail(30,1.8) .04096253

- GL

Re: st: functions for computing prob

2006-03-22T13:32:00.000-05:00

Lei Xuan wrote:

> I am computing probabilities for t-test and z-test. > I want to know if the functions -tprob- and -normprob- are out-of-date > since no help files explain these functions. Are -ttail(n,t)- and -normal(z) > > right functions to compute probs?

see -help density functions- or -help functions-

Yes, the cumulative normal density function, normprob(), still works but is out of date, no longer documented, and has been replaced by normal(). Am not sure whether tprob() ever existed? - doesn't seem to work and is not documented in any case. ttail(n,t) is documented under help for functions.

- Gary

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/