Displaying 1 result from an estimated 1 matches for "propdichtegam".
2005 Jul 19
1
initial points for arms in package HI
...that my target density is not logconcave.
Thanks for all comments or ideas.
Christoph Buser
## R Code:
library(HI)
## parameter for the distribution
para <- 0.1
## logdensity
logDichteGam <- function(x, u = para, v = para) {
-(u*x + v*1/x) - log(x)
}
## density except for the constant
propDichteGam <- function(x, u = para, v = para) {
exp(-(u*x + v*1/x) - log(x))
}
## calculating the constant
(c <- integrate(propDichteGam, 0, 1000, rel.tol = 10^(-12))$value)
## density
DichteGam <- function(x, u = para, v = para) {
exp(-(u*x + v*1/x) - log(x))/c
}
## calculating 1000 values by...