[R] Identifying peak periods of observations in circular yearly data

[Daisy Englert Duursma]

Greetings,

I am trying to identify at which point during the year 80% of bird breeding
observations are. typically I would answer a question like this by finding
the median or quartiles but how do I deal with situations where the 80% of
the is from day 285 through day 366 (leap year) and extends to day 30?

The data is circular and and day 365 is as close to day 366 as day 1.

I am reading the manual for CircStats and circular but I could really use
some help on this.


Here is some dummy data:


obsDay<-c(rep(1:30,10),rep(45:65,2),65:180,rep(181:265,2),rep(266:330,4),rep(331:366,6))

plot(density(obsDay))



-- 
Daisy Englert Duursma
Department of Biological Sciences
Room W19F 135
Macquarie University, North Ryde, NSW 2109
Australia

Tel +61 2 9850 1302

	[[alternative HTML version deleted]]