R/pkimball.R

In JohnsonBrent/kimball: Several functions to impliment a truncated normal distribution

#' Returns the probability density function
#'
#' @param x A positive numerical input vector requiring evaluation
#' @param min A parameter for the minimum of the truncated normal distribution. Must be positive. Default is zero.
#' @param L Mean parameter. Must be positive.
#' @param S A parameter for the maximum of the truncated normal distribution. Must be positive.
#' @return A numerical result of length x.
#' @importFrom stats dnorm optimize pnorm
#' @examples
#' pkimball(3, 8, 4, min=0)
#' @export

pkimball <- function(x, S, L, min=0) {
  L <- L-min  # L = mean
  S <- S-min  # S = max
  x <- x-min
  h <- hmin(S,L)
  phi.h <- dnorm(-h,0,1)
  PHI.h <- pnorm(h,0,1)
  w <- h + phi.h/PHI.h
  Ft <- pnorm(w*x/L-h,0,1) / PHI.h
  return(Ft)
}