R/loglikelihood_global_admkr.R

In bbefkr: Bayesian bandwidth estimation and semi-metric selection for the functional kernel regression with unknown error density

loglikelihood_global_admkr <- function(h, resid)
{
    b = h[2]
    epsilon = scale(resid)
    std = sd(resid)
    cont = (2.0*pi)^(-0.5)
    logf = vector(,length(resid))
    for(i in 1:length(resid))
    {
        temp = epsilon[i] - epsilon[-i]
        res = sum(cont * exp(-0.5*((temp/b)^2))/b)
        logf[i] = log(res/length(temp)/std)
    }
    return(sum(logf))
}