[R] Fw: Hist & kernel density estimates
marammagdysalem at yahoo.com
Wed Aug 19 18:02:08 CEST 2009
For the hist estimate
+ right=FALSE,xlim=c(9000,16000),ylim=ylim,main="Histogram of q(scott)")
For the kernel estimate>options(scipen=4)
> d <- density(q, bw = "nrd0",kernel="gaussian")
In fact the variable q is a vector of 1000 simulated values; that is I generated 1000 samples from the pareto distribution, from each sample I calculated the value of q ( a certain fn in the sample observations), and thus I was left with 1000 values of q and I don't know the distribution of q.
Hence, I used the attached codes for histogram and kernel density estimation toestimate the density of q.
But what I'm really intersed in is to estimate the probability that q is greater than a certain value , for ex.,P(q>11000), using the density estimates I obtained.
Could u help me with a fn or some document to do this?
Thank u so much
Attached are the codes of a histogram & a kernel density estimate and the output they produced. I'll copy the codes here in case there's something wrong with the attachement
-------------- next part --------------
R-help at r-project.org mailing list
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