[R] truehist and density plots

Cuvelier Etienne ecu at info.fundp.ac.be
Tue Apr 28 15:43:33 CEST 2009


Perhaps like this :

temp =density(b[20:50])
truehist(b, ymax=max(temp$y))
lines(temp)

Etienne

carol white a écrit :
> Consider a vector of 100 elements (attached files). then, 
> 
> truehist(b)
> lines(density(b[20:50]))
> 
> How is it possible to have density plots of all subsets like b[20:50] within histogram (without exceeding the max of historgram on y axis)?
> 
> Is it more clear?
> 
> Best,
> 
> --- On Tue, 4/28/09, Uwe Ligges <ligges at statistik.tu-dortmund.de> wrote:
> From: Uwe Ligges <ligges at statistik.tu-dortmund.de>
> Subject: Re: [R] truehist and density plots
> To: "carol white" <wht_crl at yahoo.com>
> Cc: r-help at stat.math.ethz.ch
> Date: Tuesday, April 28, 2009, 5:42 AM
> 
> 
> carol white wrote:
>> Hi,
>> I wanted to plot the histogram of a vector and then, plot the density
> function of subsets of the vector on the histogram. So I use truehist in MASS
> package and lines(density) as follows:
>> length(b) = 1000
>> truehist(b)
>> lines(density(b[1:100]))
> 
> 
> I do not undertsand what you mean. Can you please provide a *reproducible*
> example?
> 
> Uwe Ligges
> 
> 
>> however the density plot of the first 100 points exceeds the max of y axis
> (see attached). how is it possible to make a graphics so that the density plot
> of the subsets doesn't go beyond the maximum of all points in the complete
> set?
>> Cheers,
>>
>> Carol
>>
>>
>>
>>
>>
>>
>>       
>>
>> ------------------------------------------------------------------------
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
> 
> 
> 
>       
> 
> 
> ------------------------------------------------------------------------
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list