Re: st: Re: programming assist, too many unique values for levels

["Michael Blasnik" <michael.blasnik@verizon.net>]

..
Woops, there's a mistake in one line:

gen bphigh=bp_systolic >= `threshold' & bp_systolic !=.

should be

gen bphigh=bp_systolic >= `threshold' if bp_systolic !=.

This change will cause it to have missing values when bp is missing, which will 
give the correct proportion (assuming you don't want to count missing values in 
the denominator)

Michael

----- Original Message ----- 
From: "Andrew O'Connor DO" <aoconnor@metrohealth.org>
To: <statalist@hsphsun2.harvard.edu>
Sent: Tuesday, May 01, 2007 4:25 PM
Subject: Re: st: Re: programming assist, too many unique values for levels



> thank you
> Andrew O'Connor, DO, MPH
> Division of Nephrology
> Center for Healthcare Research and Policy
> MetroHealth Medical Center/Case Western Reserve University
> (216)778-8484
> > > > michael.blasnik@verizon.net 05/01/07 4:19 PM >>>
> You don't need the levels command or the loops at the end, just use
> -egen-
> commands:
>
> egen time=sum(obstime), by(pt)
> egen time_out=sum(out_range), by(pt)
> gen bphigh=bp_systolic >= `threshold' & bp_systolic !=.
> egen proportion=mean(bphigh), by(pt)
<snip> 
*
*   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/