R/llfcn.vec.R

In episplineDensity: Density Estimation with Soft Information by Exponential Epi-splines

llfcn.vec <- function (x, epiparameters, softinfo) 
{
# llfcn -- vector style
#
# Compute the reference density. The name of the function to compute
# that (e.g. "dnorm") had better be found in softinfo$KLDensity, with
# parameters softinfo$KLDensityParams
#
if (!any (names (softinfo) == "KLDensity"))
    stop ("No KL density supplied. How did we get here?")
if (!exists (softinfo$KLDensity))
    stop (paste ("KL density", softinfo$KLDensity, "not found"))

param.list <- softinfo$KLDensityParams
param.list$x <- x

c.out <- coeff (x, epiparameters)
# Shorten!
c.out.short <- c.out[(epiparameters$Ndiscr + 2):length(c.out)]


h <- do.call (softinfo$KLDensity, param.list)
#
# h holds the density; now compute and return the integrand 
#

return (h * c.out.short)
}