<>
Nick: How is -kdensity- different in this respect? Maybe it does not create
r-results internally?
HTH
Martin
-----Ursprüngliche Nachricht-----
Von: [email protected]
[mailto:[email protected]] Im Auftrag von Nick Cox
Gesendet: Donnerstag, 4. Februar 2010 14:05
An: [email protected]
Betreff: RE: AW: st: -if- and r() in -histogram-
`r(95)' is evaluated before -histogram- even opens: -histogram- sees only
the numeric result.
That's not true of r(95). In the first instance the text "r(95)" is used and
only later it is evaluated.
That is, you want -histogram- to work with a pre-existing r(95), but inside
-histogram- other r-class results are produced and that undermines your aim.
Nick
[email protected]
Antoine Terracol
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.
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: [email protected]
> [mailto:[email protected]] Im Auftrag von Antoine
> Terracol
> Gesendet: Donnerstag, 4. Februar 2010 13:24
> An: [email protected]
> 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)
>> /************************/
*
* 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/
*
* 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/
