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

夜筱巫

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

1. leapyear<-function(x){}

复制代码

闰年这里的定义是:
3 l! b7 q* T, D8 U% m9 _; `/ C% w如果那一年可以整除4（比如2004），2月就有29天
如果那一年不能整除100
如果那一年可以整除400
3 v" {* b8 v' Y$ s
如果x<1，这个function就要return “error”。我不大懂这个条件的意思，可否请大师解释下？
, g) }+ y8 K% l: ~% S
我写的function是这样的：

1. leapyear<- function(dates, years_to_add)
   ; p6 V$ B0 C7 I% v
2. {
   + C8 D+ |( p) Z5 p+ w
3.   if(!require(lubridate)) stop("install lubridate")
   ! q$ E, c1 s. ^- y) E
4.   if(years_to_add > floor(years_to_add)) stop("error")
   ; e% F8 y/ |5 ^9 \! ^5 t/ |' t
5.   / t9 a+ Y+ v8 F3 I# R+ ^
   
6.   new <- dates + years(years_to_add)
   ! z* l) ^3 e\" R; C+ U  z' X5 N
7.   new[2 [. U: B9 B& y) }$ c! {
   
8.     leap_year(year(dates)) & / P2 \4 B+ y( V\" R. k& C! |
   
9.       !leap_year(year(dates) + years_to_add) & 8 p# z% u& P* B# Z7 t. o0 d
   
10.       month(dates) == 2 &
    # v. O; c8 Y\" {7 A3 u4 X
11.       day(dates) == 29
    , R\" q! C9 ]4 W6 t
12.     ] <- NA
    4 b( m' F: X/ a3 L3 [' z
13.   new_dates
    \" [- N, m5 P/ k! U+ }. d
14. }+ Z( \2 E- F  ^7 [
    
15. 
    # D3 T# A1 j' K- L3 P
16. dates <- ymd(c("2004-2-28", "2004-2-29"))

复制代码

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

请各位大师帮帮忙！小妹感激不尽！