Estimate directional statistics for one-sigma events
Estimate directional statistics for one-sigma events
Uses the package "CircStats" to find the mean direction, rho, and kappa of the von Mises distribution summarizing the ordinal day of high- and low-flow events.
circ.s(x)
Arguments
x: Output from "compare.periods" function (of class "compflows").
Returns
circstats: data.frame with rows corresponding to the high- and low-flows for both periods. The columns list n, mu, rho, and kappa as calculated using the CircStats package.
Examples
# load datadata("sycamore")# compare for periods from 1960 to 1979 and 1980 to 1999comp = compare.periods(c("1960-01-01","1979-12-31"),c("1980-01-01","1999-12-31"), sycamore, plot=FALSE)circ.s(comp)