# R/medianHL.circular.R In circular: Circular Statistics

#### Documented in medianHL.circular

```#############################################################
#
#   median.circular function
#   Author: Claudio Agostinelli and Alessandro Gagliardi
#   E-mail: claudio@unive.it
#   Date: May, 12, 2015
#   Version: 0.4-1
#
#   Copyright (C) 2015 Claudio Agostinelli and Alessandro Gagliardi
#
#############################################################

medianHL.circular <- function(x, na.rm=FALSE, method=c("HL1","HL2","HL3"), prop=NULL) {
method <- match.arg(method)
if (!is.null(prop))
if (prop <= 0 | prop >=1)
stop("'prop' is outside (0,1)")
if (na.rm)
x <- x[!is.na(x)]
if (length(x)==0) {
warning("No observations (at least after removing missing values)")
return(circular(NA))
}
if (is.circular(x)) {
dc <- circularp(x)
} else {
dc <- list(type="angles", units="radians", template="none", modulo="asis", zero=0, rotation="counter")
}
x <- conversion.circular(x, units="radians", zero=0, rotation="counter")
attr(x, "class") <- attr(x, "circularp") <-  NULL
circmedian <- conversion.circular(circular(drop(circmedian)), dc\$units, dc\$type, dc\$template, dc\$modulo, dc\$zero, dc\$rotation)
attr(circmedian, "medians") <- conversion.circular(circular(drop(MinusPiPlusPiRad(attr(circmedian, "medians")))), dc\$units, dc\$type, dc\$template, dc\$modulo, dc\$zero, dc\$rotation)
attr(attr(circmedian, "medians"), "class") <- attr(attr(circmedian, "medians"), "circularp") <-  NULL
return(circmedian)
}

MedianHLCircularRad <- function(x, method, prop) {
x <- x%%(2*pi)
n <- length(x)
mediancirc <- NA
methods <- c("HL2","HL1","HL3")
if (is.null(prop))
{