jeudi 23 août 2018

Trouble with Rayleigh test in R - not getting significant results

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