## The function is currently defined as
function (v, v2, x, distr.k, sigma.k, distr.p0, mu.p0, sigma.p0)
{
alpha <- p0(v, distr = distr.p0, mu = mu.p0, sigma = sigma.p0)/
p0(v2, distr = distr.p0, mu = mu.p0, sigma = sigma.p0)
Prod <- 1
for (i in seq(length(x))) {
fac <- dk(x[i], distr = distr.k, mu = v, sigma = sigma.k)/
dk(x[i], distr = distr.k, mu = v2, sigma = sigma.k)
Prod <- Prod * fac
}
f <- alpha * Prod
return(f)
}
Run the code above in your browser using DataLab