R/drootinvgamma.R

drootinvgamma <-
function(x,shape,rate,log=F) {
	output = log(2)+shape*log(rate)-lgamma(shape)-(2*shape+1)*log(x)-rate/x^2
	if (log==F) output = exp(output)
	return(output)
}
stevenlecomber/Rgeoprofile-1.1.0 documentation built on May 30, 2019, 4:46 p.m.