R/calendarprop.R

In frousseu/rNest: Estimating Nesting Chronology and Breeding Phenology in Birds

calendarprop <-
function(x,jul=1:365,column){
	dates<-x[,column] # change quantile !!!
	r<-range(dates,na.rm=T)
	p<-sapply(jul,function(i){
		sum(dates<=i,na.rm=T)/sum(!is.na(dates))
	})
	p
}