On Nov 20, 2010, at 9:34 PM, Shant Ch wrote:
David, I did look at ?density many times. I think I didn't explain
what I have to find.
Suppose I have a data.
x<- c(rnorm(40,5,3),rcauchy(30,0,4),rexp(30,6))
Suppose I don't have information about the mixture, I have been
given only the data.
density(x) will give the 6 number summary of the data, given as x
and also the 6 number summary of the density of density given as y.
I am not sure what the number six (6) represents in the two palces it
occurs:
> str(density(x))
List of 7
$ x : num [1:512] -59.2 -59 -58.8 -58.5 -58.3 ...
$ y : num [1:512] 3.30e-05 5.27e-05 8.19e-05 1.24e-04
1.82e-04 ...
$ bw : num 1.37
$ n : int 100
$ call : language density.default(x = x)
$ data.name: chr "x"
$ has.na : logi FALSE
- attr(*, "class")= chr "density"
Perhaps you want to look at:
?approxfun
It would let you interpolate at points that are not part of the set
density(x)$x
__
David.
I want to find the density of the given data at x=1. I basically
want the value of y(=density) for x=1 i.e. kernel density at x=1.
Shant
From: David Winsemius <dwinsem...@comcast.net>
To: Shant Ch <sha1...@yahoo.com>
Cc: r-help@r-project.org
Sent: Sat, November 20, 2010 8:54:32 PM
Subject: Re: [R] density at particular values
On Nov 20, 2010, at 8:07 PM, Shant Ch wrote:
> Hello everyone!
>
> I want to use density function of R to compute the density at
x(=0, say). But it
> is giving me the 5-number summary and mean of the data and
densities at that
> point.
> I just want the densities at different values specified by me. Can
anyone let me
> know how can I find that?
Here's what you should have done (even before posting):
?density
<Read the help page to see the structure of what density() returns.>
"Value
x the n coordinates of the points where the density is estimated.
y the estimated density values. These will be non-negative, but can
be zero."
Realize that the "specified by me" part is either going to be
modified to "pick an existing estimate near my specification" or
that you will need to approximate the value. So what is the actual
problem (and the actual data setup) ?
--David.
>
> For example
>
>
> Thanks in advance for your help.
>
>
> Shant
>
>
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@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.
David Winsemius, MD
West Hartford, CT
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@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.