File: stephens.test.R

package info (click to toggle)
r-cran-circular 0.5-1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,492 kB
  • sloc: ansic: 464; fortran: 69; sh: 13; makefile: 2
file content (11 lines) | stat: -rw-r--r-- 304 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
 
StephensTestRad <- function(x, mu=0, alpha) {
  n <- length(x)
  Muobs <- MeanCircularRad(x)
  Robs <- RhoCircularRad(x)*n
  Cobs <- Robs*cos(Muobs-mu) 
  Rcrit <- Ralpha(x=Cobs, n=n, alpha=alpha)
  test <- TRUE
  if (Robs > Rcrit) test <- FALSE
  result <- c(test, Robs, Cobs, Rcrit)
  return(result)
}