R/negloglik.laplace.R

`negloglik.laplace` <-
function (xpar, xx) 
{
    a <- xpar[2]
    w <- wfromt(xpar[1], a = a)
    loglik <- sum(log(1 + w * beta.laplace(xx, a)))
	if (abs(loglik)==Inf){
		loglik<-sign(loglik)*10^6
	}   
 return(-loglik)
}

Try the binhf package in your browser

Any scripts or data that you put into this service are public.

binhf documentation built on May 2, 2019, 4:07 p.m.