R/polynom.R

In LuchTiarna/PsyMetFuns: Package for extrapolating psychometric information

#'@name polynom
#'@rdname polynom
#'@aliases
#'poly
#'poly.cdf
#'poly.pdf
#'
#'@title Polynom
#'@description
#'Polynomial function is a core type function.
#'\cr
#'It's CDF formula is: y = (x / a) ^ b
#'\cr
#'It's PDF formula is: y = (x / a) ^ (b- 1) * b / a
#'@param x Vector of x parametres
#'@param params polynom function has two parameters \n a \n b
#'@param b adjusts the curve slope
#'
#'@return Vector of result vaues
NULL

#'@rdname polynom
polynom <- function(x, params) {
  UseMethod("polynom")
}
#'@rdname polynom
polynom.orig <- function(x, params) {
  UseMethod("polynom.orig")
}
#'@rdname polynom
polynom.inverse <- function(x, params) {
  UseMethod("polynom.inverse")
}
#'@rdname polynom
polynom.orig.cdf <- function(x, params) {
  if(length(params) != 2) { warning("Polynomial function has two parameters: a and b"); return(rep(NaN, length(x))) }
  return((x / params[1]) ^ params[2])
}
#'@rdname polynom
polynom.orig.pdf_x <- function(x, params) {
  if(length(params) != 2) { warning("Polynomial function has two parameters: a and b"); return(rep(NaN, length(x))) }
  return((x / params[1]) ^ (params[2]- 1) * params[2] / params[1])
}
#'@rdname polynom
polynom.orig.pdf_p1 <- function(x, params) {
  if(length(params) != 2) { warning("Polynomial function has two parameters: a and b"); return(rep(NaN, length(x))) }
  return(params[2] * (x/params[1]) ^ params[2] / params[1])
}
#'@rdname polynom
polynom.orig.pdf_p2 <- function(x, params) {
  if(length(params) != 2) { warning("Polynomial function has two parameters: a and b"); return(rep(NaN, length(x))) }
  return(((x / params[1]) ^ params[2]) * log(x / params[1]))
}
#'@rdname polynom
polynom.inverse_x.cdf <- function(x, params) {
  if(length(params) != 2) { warning("Polynomial function has two parameters: a and b"); return(rep(NaN, length(x))) }
  return(params[1] * (x ^ (1/params[2])))
}
#'@rdname polynom
polynom.inverse_x.pdf_x <- function(x, params) {
  if(length(params) != 2) { warning("Polynomial function has two parameters: a and b"); return(rep(NaN, length(x))) }
  return(params[1]/params[2] * x ^ (1/params[2]-1))
}