[R] Question concerning pdfs
dimitris.rizopoulos at med.kuleuven.be
Wed Jan 16 11:18:23 CET 2008
I think you need:
dx2 <- exp(dx)
dp2 <- dp / dx2
sum(c(0, diff(dx2)) * dp2)
I hope it helps.
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
----- Original Message -----
From: "Rainer M Krug" <R.M.Krug at gmail.com>
To: "r-help" <r-help at stat.math.ethz.ch>
Sent: Wednesday, January 16, 2008 10:34 AM
Subject: [R] Question concerning pdfs
> I have a dataset of absolute growth rates g ranging from close to 0
> around whatever of which I want to calculate a pdf.
> If I use density(g) the pdf will extend to below 0 so I logtransform
> and do d <- density(log(g)).
> Now I would like to transform this pdf back, i.e. d$x <-
> exp(d$x) but what do I have to do with d$y?
> (below is a small example)
> Thanks for your help,
> g <- exp(rnorm(100)) ## Just to generate the example data
> d <- density(log(g))
> dx <- d$x
> dp <- d$y
> sum(c(0, (diff(dx))) * dp) ## this is equal to one
> plot(dx, dp)
> dx2 <- exp(dx)
> dp2 <- ??????? ## what should I do here?
> plot(dx2, dp2) ## what should I do here?
> sum(c(0, (diff(exp(dx2)))) * dp2) ## this should be one
> sum( c(0, (diff(dx2))) * ( dp1 / c(1, (diff(exp(dx2)))) * c(0,
> (diff(dx1)))) ) ##This is obviously also one, but can
> use this to define my dp2? i.e.
> dp2 <- dp1 / c(1, (diff(exp(dx2)))) * c(0, (diff(dx1)))) ###????
> R-help at r-project.org mailing list
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help