Skip to Content

avec circular.R

  1. library(circular)
  2. angles = c(341.0, 359.0, 334.0, 15.0, 330.0, 301.0, 299.0, 9.0, 7.0, 353.0, 28.0, 25.0, 23.0, 30.0, 350.0, 25.0, 22.0, 8.0, 356.0, 27.0)
  3. # création d'un objet circular modulo 360°
  4. anglecir =  circular(angles, type="angles", units="degrees",modulo="2pi", template='geographics')
  5. summary(anglecir)
  6.        n     Min.  1st Qu.   Median     Mean  3rd Qu.     Max.      Rho 
  7.  20.0000  30.0000  23.5000   7.5000   1.0610 -12.2500 -61.0000   0.8962 
  8.  
  9. # détails
  10. # direction moyenne angulaire
  11. mean(anglecir)
  12. Circular Data: 
  13. Type = angles 
  14. Units = degrees 
  15. Template = geographics 
  16. Modulo = 2pi
  17. Zero = 1.570796 
  18. Rotation = clock 
  19. [1] 1.060902
  20. # longueur résultante moyenne
  21. rho.circular(anglecir)
  22. 0.8938667
  23. # variance angulaire
  24. 1 - rho.circular(anglecir)
  25. 0.1037632
  26. # variance angulaire selon Batschelet (1981) = 2 * (1 - R)
  27. angular.variance(anglecir)
  28. [1] 0.2075264
  29. # déviation angulaire standard
  30. angular.deviation(anglecir)
  31. [1] 0.4555506
  • 933 lectures


Theme by Dr. Radut.