R/inv.gauss.R

In flexCWM: Flexible Cluster-Weighted Modeling

####################################################
## Re-parameterized Inverse Gaussian distribution ##
####################################################

dig <- function(x,mu,var,log=FALSE){
  
  statmod::dinvgauss(x, mean=mu, dispersion=var/mu^3, log=log)
  
}

rig <- function(n,mu,var){
  
  statmod::rinvgauss(n, mean=mu, dispersion=var/mu^3)
  
}