function [ llv ] = nlln(mu) % computes the non-parametric loglikelihood value % global resids data sigmahat; llv=0; n=size(resids,1); for i=1:n; llv=llv+log(kdensity((data(i)-mu),resids)); end llv=-llv; end