R/dsiniw.R

In lgallindo/SinIW: The SinIW Distribution

#' The density function of the SinInverseWeibull probability distribution.
#' @export
#'
#' @param x vector of quantiles.
#' @param alpha Alpha parameter.
#' @param theta Theta parameter.
#' @return A vector with n observations of the SinInverseWeibull distribution.
#' @examples
#' dsiniw(0.5,1,1)
#' dsiniw(0.5,0.5,0.7)

dsiniw <- function(x,alpha,theta){
  alpha*theta*(pi/2)*x^(-theta-1)*exp(-alpha*x^(-theta))*cos((pi/2)*exp(-alpha*x^(-theta)))
}