I'm trying to analyze some circular data using R. I'm running a rayleigh test on a subset of a data set, using this code:
CC_062818_5AS <- subset(Assay_Data, select = T1:T17,
Nest_ID == "CC-062818-5AS")
That produces a vector that looks like this:
T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 T11 T12 T13 T14 T15 T16 T17
1 135 105 99 95 86 82 73 121 98 88 86 86 83 113 109 104 99
And I run the Rayleigh test like this:
rayleigh.test(CC_062818_5AS, mu = 86.8)
And get this result:
Rayleigh Test of Uniformity
Alternative with Specified Mean Direction: 5.118591
Test Statistic: 0.1542
P-value: 0.1863
Warning messages:
1: In as.circular(x) :
an object is coerced to the class 'circular' using default
value for the following components:
type: 'angles'
units: 'radians'
template: 'none'
modulo: 'asis'
zero: 0
rotation: 'counter'
conversion.circularxradians0counter2pi
I'm a complete stats beginner, and R beginner. The data I have should be significantly oriented - they're about as oriented as it gets in my field. The data are in degrees. I tried coercing the data to circular myself so that it would be in degrees and not the default radians, but it didn't make a difference. What am I doing wrong?
Aucun commentaire:
Enregistrer un commentaire