1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
var <- function(x, ...) UseMethod("var")
var.default <- function(x, y = NULL, na.rm = FALSE, use, ...) stats::var(x=x, y=y, na.rm=na.rm, use=use)
#var.matrix <- function(x, ...) {
# apply(x, 2, var, ...)
#}
var.data.frame <- function(x, ...) {
sapply(x, var, ...)
}
#############################################################
# #
# Original Splus: Ulric Lund #
# E-mail: ulund@calpoly.edu #
# #
#############################################################
#############################################################
# #
# var.circular function #
# Author: Claudio Agostinelli #
# Email: claudio@unive.it #
# Date: June, 24, 2011 #
# Copyright (C) 2011 Claudio Agostinelli #
# #
# Version 0.5 #
#############################################################
var.circular <- function (x, na.rm=FALSE, ...) {
if (is.matrix(x)) {
apply(x, 2, var.circular, na.rm=na.rm)
} else {
if (na.rm)
x <- x[!is.na(x)]
x <- conversion.circular(x, units="radians", zero=0, rotation="counter")
attr(x, "class") <- attr(x, "circularp") <- NULL
VarCircularRad(x=x)
}
}
VarCircularRad <- function(x) {
rbar <- RhoCircularRad(x)
circvar <- 1-rbar
return(circvar)
}
|
