# Re: AW: st: -if- and r() in -histogram-

 From Antoine Terracol To statalist@hsphsun2.harvard.edu Subject Re: AW: st: -if- and r() in -histogram- Date Thu, 04 Feb 2010 13:41:31 +0100

<>

True, but I find it strange that references to r() work with -kdensity- (or commands like -summarize-), but hot with -histogram- (nor with -graph box-) where `r()' seems to be required.

Antoine

Martin Weiss wrote:
<>
You do get the correct graph with:

*************
clear
cap graph drop _all
set obs 1000
g x=exp(rnormal())
su x, de
hist x if x<`r(p95)', name(ha)
*************

HTH
Martin

-----Ursprüngliche Nachricht-----
Von: owner-statalist@hsphsun2.harvard.edu
[mailto:owner-statalist@hsphsun2.harvard.edu] Im Auftrag von Antoine
Terracol
Gesendet: Donnerstag, 4. Februar 2010 13:24
An: statalist@hsphsun2.harvard.edu
Betreff: Re: st: -if- and r() in -histogram-

<>

of course, all the "hist" below should be "kdensity", and the correct piece of code is

/**********************/
clear
cap graph drop _all
set obs 1000
g x=exp(rnormal())
su x, de
local p95=r(p95)
di r(p95)
kdensity x if x<r(p95), name(ka)
kdensity x if x<`p95', name(kb)
kdensity x, name(kc)
/************************/

Antoine Terracol wrote:
Now, I have not tried all types of graphs, but -kdensity- seems to behave as expected:

/**********************/
clear
cap graph drop _all
set obs 1000
g x=exp(rnormal())
su x, de
local p95=r(p95)
di r(p95)
hist x if x<r(p95), name(ka)
hist x if x<`p95', name(kb)
hist x, name(kc)
/************************/

Best,
Antoine

--
Ce message a ete verifie par MailScanner
pour des virus ou des polluriels et rien de
suspect n'a ete trouve.

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/