怎样用function来检验是否是闰年

夜筱巫

我要写个function来检验末年是否是闰年。开头这样：

1. leapyear<-function(x){}

复制代码

闰年这里的定义是:
, |3 W3 z$ T8 Q2 Q如果那一年可以整除4（比如2004），2月就有29天
如果那一年不能整除100
如果那一年可以整除400
0 G" Z& d$ ?2 j* ]$ E. t) i
( Q7 Z5 w* z# y+ G; q1 U如果x<1，这个function就要return “error”。我不大懂这个条件的意思，可否请大师解释下？
" z2 ?- W2 }. `  N9 _1 V/ o
我写的function是这样的：

1. leapyear<- function(dates, years_to_add)
   * S) N$ D% C4 ~& _3 X: F* V# v! s
2. {: ^) U4 S+ ^% p! h
   
3.   if(!require(lubridate)) stop("install lubridate")
   * a' k# e, [% d' c
4.   if(years_to_add > floor(years_to_add)) stop("error")
   4 [. R\" p- D& h; s7 d( }8 @3 F+ o, U
5.   
   5 E; F% w- Y( D& D5 Y
6.   new <- dates + years(years_to_add)) }  E2 T2 P% H3 l0 c7 F
   
7.   new[: K: h3 X! a1 s& z$ C% H
   
8.     leap_year(year(dates)) &
   ! L5 k; J( ^$ }+ R
9.       !leap_year(year(dates) + years_to_add) &
   7 P$ e+ f0 ]5 T7 y; A$ Q
10.       month(dates) == 2 & 5 C% {0 Y5 l5 g/ i4 J
    
11.       day(dates) == 29# ?+ B5 b0 p4 P6 j/ q& h8 p, K1 ~
    
12.     ] <- NA* C/ a/ ^6 I. B
    
13.   new_dates5 _# A$ @- z8 ]4 E/ O# I0 a+ U2 w$ F
    
14. }
    0 \: X/ Z$ h5 a% k: N; M
15. 8 A( d2 y. M- J* \5 @
    
16. dates <- ymd(c("2004-2-28", "2004-2-29"))

复制代码

不知道对不对，但是检验的时候显示ymd找不到。。。。

, s8 p( m( O# n请各位大师帮帮忙！小妹感激不尽！