R/dvm.R

Defines functions dvm

Documented in dvm

dvm <- function(x, m, k, rads = FALSE, logden = FALSE) {
  if ( !rads )  x <- x * pi/180
  den <- k * cos(x - m) - log(2 * pi) - log( besselI(k, 0, expon.scaled = TRUE) ) - k
  if ( !logden )  den <- exp(den)
  den
}

Try the Directional package in your browser

Any scripts or data that you put into this service are public.

Directional documentation built on Oct. 12, 2023, 1:07 a.m.