## estimate gamma density
set.seed(1977)
x <- rgamma(200, 2, 1)
res <- logConDens(x, smoothed = TRUE, print = FALSE)
## compute function values at an arbitrary point
xs <- (res$x[100] + res$x[101]) / 2
evaluateLogConDens(xs, res)
## only compute function values for non-smooth estimates
evaluateLogConDens(xs, res, which = 1:3)
Run the code above in your browser using DataLab