tests/test-median.R

suppressMessages(library("circular"))
# ?median.circular


  med <- median.circular(circular(c(0, pi/2, pi, 3/2*pi)))

## expect:
## median is equal to NA
  stopifnot(
    is.na(med)
  )

## expect:
## all the points are minimizers of the function
  stopifnot(
    all.equal(attr(med, "medians"), c(0, pi/2, pi, 3/2*pi), tol = 2e-7)
  )

Try the circular package in your browser

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

circular documentation built on July 4, 2017, 9:03 a.m.