R/Hf.R

Defines functions InverseMillsRatio

InverseMillsRatio <- function(q, mean, sd) {
  x <- (q - mean) / sd
  pdf <- dnorm(x, log = TRUE)
  cdf <- pnorm(x, log = TRUE)
  d <- exp(pdf - cdf)
  d[is.na(d)] <- 0
  return(d)
}
zy26/LTMGSCA documentation built on Jan. 24, 2020, 11:39 p.m.