Hi, rvm is used to generate pseudo-random numbers from a von Mises
distribution
Usage: rvm(n, mean, k)
Maybe this is the code you want:
require(CircStats)
require(circular)
data <- circular(c(0.41, 5.02, 2.41, 1.57, 4.03, 6.14, 5.46, 0.79, 1.55,
1.18, 5.54, 4.21, 4.47, 0.28, 1.70, 5.31, 0.68, 0.88, 5.48, 5.16, 2.47, 2.72,
6.00, 5.68, 5.50, 6.19, 2.75, 2.97, 3.02, 4.80, 1.89, 2.06, 1.80))
r.test(data)
HTH
Xiaohua
On 6/28/06, Pascual Lopez <[EMAIL PROTECTED]> wrote:
>
> Hello everybody, I´m new on this list and I´m looking for help with
> circular statistics.
>
> I would like to perform a Rayleigh test for a vector of radians.
>
> According to the example I type in the workspace of R:
>
> data <- rvm(0.41, 5.02, 2.41, 1.57, 4.03, 6.14, 5.46, 0.79, 1.55, 1.18,
> 5.54, 4.21, 4.47, 0.28, 1.70, 5.31, 0.68, 0.88, 5.48, 5.16, 2.47, 2.72,
> 6.00, 5.68, 5.50, 6.19, 2.75, 2.97, 3.02, 4.80, 1.89, 2.06, 1.80)
> r.test(data)
>
> And the programm returns the following:
>
> data <- rvm(0.41, 5.02, 2.41, 1.57, 4.03, 6.14, 5.46, 0.79, 1.55, 1.18,
> 5.54, 4.21, 4.47, 0.28, 1.70, 5.31, 0.68, 0.88, 5.48, 5.16, 2.47, 2.72,
> 6.00, 5.68, 5.50, 6.19, 2.75, 2.97, 3.02, 4.80, 1.89, 2.06, 1.80)
> Error in rvm(0.41, 5.02, 2.41, 1.57, 4.03, 6.14, 5.46, 0.79, 1.55, 1.18
> , :
> unused argument(s) ( ...)
> > r.test(data)
> Error in sin(x) : Non-numeric argument to mathematical function
>
> Data are in radians and I have loaded the CircStat package previously. How
> could I perform this analysis?
> Why it returns "unused argument(s) (...)? How can I entry my data vector?
>
> Thank you in advance.
>
> Pascual
>
>
>
> [[alternative HTML version deleted]]
>
>
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
>
[[alternative HTML version deleted]]
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
