年月日からDOYを返す関数【R】

doy<-function(yyyy,mm,dd){
dom<-c(0,31,28,31,30,31,30,31,31,30,31,30)
if(yyyy%%4==0){
dom[2]<-29
}
doy<-sum(dom[1:mm])+dd
return (doy)
}