File: rao.test.Rd

package info (click to toggle)
r-cran-circular 0.4-93-2
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 1,492 kB
  • sloc: ansic: 463; fortran: 69; sh: 13; makefile: 2
file content (56 lines) | stat: -rw-r--r-- 1,688 bytes parent folder | download | duplicates (3)
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
49
50
51
52
53
54
55
56
\name{rao.test}
\title{Rao's Tests for Homogeneity}
\alias{rao.test}
\alias{print.rao.test}

\description{
Performs Rao's test for homogeneity on k populations of angular data. 
}

\usage{
rao.test(\dots, alpha=0)
\method{print}{rao.test}(x, digits = 4, \dots)
}

\arguments{
  \item{\dots}{a sequence of \code{\link{circular}} for the
	\code{rao.test} and further arguments passed to or from other
	methods for the \code{print.rao.test} function.}
  \item{alpha}{numeric value specifying the significance level of the
  test. Default is 0, in which case p-values for the test statistic is
  printed.}
  \item{x}{an object from the \code{rao.test}.}
  \item{digits}{integer indicating the precision to be used.}
}

\value{
A list with the statistic and p.value for the mean and the dispersion
and the value of alpha. 
}

\note{
The test is performed, and the results are written to the screen. Test results are given for both the test of equality of polar vectors, and of dispersions.  If alpha is specified, the test statistic is printed, along with the level critical value.  If alpha is not specified, a p-value for the test is printed.  
}

\details{
Critical values and p-values are determined according to the chi-squared approximation of the test statistic.
}

\author{Claudio Agostinelli and Ulric Lund}

\references{
Jammalamadaka, S. Rao and SenGupta, A. (2001). Topics in Circular Statistics, Section 7.6.1, World Scientific Press, Singapore.

Rao, J.S. (1967). Large sample tests for the homogeneity of angular data, Sankhya, Ser, B., 28.
}

\examples{
x <- rvonmises(100, circular(0), kappa=10)
y <- rvonmises(100, circular(0), kappa=10)

rao.test(x, y)

}


\keyword{htest}