R 如何返回间隔至少10天的每年的月最小值

我有一个每日时间序列，我想得到每年每个月的最小值，但我想确保结果至少相隔10天。更具体地说，让我们在下面的示例数据框中进行解释

>Data
        Years   Months     Days Date        A   B
    1   2003    December    1   2003-12-01  10  10
    2   2003    December    2   2003-12-02  10  10
    3   2003    December    3   2003-12-03  10  10
    4   2003    December    4   2003-12-04  10  10
    5   2003    December    5   2003-12-05  10  10
    6   2003    December    6   2003-12-06  10  10
    7   2003    December    7   2003-12-07  10  10
    8   2003    December    8   2003-12-08  3   10
    9   2003    December    9   2003-12-09  10  10
    10  2003    December    10  2003-12-10  10  10
    11  2003    December    11  2003-12-11  10  10
    12  2003    December    12  2003-12-12  10  4
    13  2003    December    13  2003-12-13  10  10
    14  2003    December    14  2003-12-14  10  10
    15  2003    December    15  2003-12-15  10  10
    16  2003    December    16  2003-12-16  10  10
    17  2003    December    17  2003-12-17  10  10
    18  2003    December    18  2003-12-18  10  10
    19  2003    December    19  2003-12-19  10  10
    20  2003    December    20  2003-12-20  10  10
    21  2003    December    21  2003-12-21  10  10
    22  2003    December    22  2003-12-22  10  10
    23  2003    December    23  2003-12-23  10  10
    24  2003    December    24  2003-12-24  10  10
    25  2003    December    25  2003-12-25  10  10
    26  2003    December    26  2003-12-26  10  10
    27  2003    December    27  2003-12-27  10  10
    28  2003    December    28  2003-12-28  10  10
    29  2003    December    29  2003-12-29  10  10
    30  2003    December    30  2003-12-30  10  10
    31  2003    December    31  2003-12-31  10  10
    32  2004    January     1   2004-01-01  10  10
    33  2004    January     2   2004-01-02  10  10
    34  2004    January     3   2004-01-03  10  10
    35  2004    January     4   2004-01-04  10  10
    36  2004    January     5   2004-01-05  10  10
    37  2004    January     6   2004-01-06  10  10
    38  2004    January     7   2004-01-07  10  10
    39  2004    January     8   2004-01-08  10  10
    40  2004    January     9   2004-01-09  10  10
    41  2004    January     10  2004-01-10  10  10
    42  2004    January     11  2004-01-11  10  10
    43  2004    January     12  2004-01-12  10  10
    44  2004    January     13  2004-01-13  10  10
    45  2004    January     14  2004-01-14  10  10
    46  2004    January     15  2004-01-15  10  10
    47  2004    January     16  2004-01-16  10  10
    48  2004    January     17  2004-01-17  10  10
    49  2004    January     18  2004-01-18  10  10
    50  2004    January     19  2004-01-19  10  10
    51  2004    January     20  2004-01-20  10  10
    52  2004    January     21  2004-01-21  10  10
    53  2004    January     22  2004-01-22  10  10
    54  2004    January     23  2004-01-23  10  10
    55  2004    January     24  2004-01-24  10  10
    56  2004    January     25  2004-01-25  5   4
    57  2004    January     26  2004-01-26  10  10
    58  2004    January     27  2004-01-27  10  10
    59  2004    January     28  2004-01-28  10  10
    60  2004    January     29  2004-01-29  10  10
    61  2004    January     30  2004-01-30  10  10
    62  2004    January     31  2004-01-31  10  10
    63  2004    February    1   2004-02-01  10  10
    64  2004    February    2   2004-02-02  5   4
    65  2004    February    3   2004-02-03  10  10
    66  2004    February    4   2004-02-04  10  10
    67  2004    February    5   2004-02-05  10  10
    68  2004    February    6   2004-02-06  10  10
    69  2004    February    7   2004-02-07  10  10
    70  2004    February    8   2004-02-08  10  10
    71  2004    February    9   2004-02-09  7   6
    72  2004    February    10  2004-02-10  10  10
    73  2004    February    11  2004-02-11  10  10
    74  2004    February    12  2004-02-12  10  10
    75  2004    February    13  2004-02-13  10  10
    76  2004    February    14  2004-02-14  10  10
    77  2004    February    15  2004-02-15  10  10
    78  2004    February    16  2004-02-16  10  10
    79  2004    February    17  2004-02-17  10  10
    80  2004    February    18  2004-02-18  10  10
    81  2004    February    19  2004-02-19  10  10
    82  2004    February    20  2004-02-20  10  10
    83  2004    February    21  2004-02-21  10  10
    84  2004    February    22  2004-02-22  10  10
    85  2004    February    23  2004-02-23  10  10
    86  2004    February    24  2004-02-24  10  10
    87  2004    February    25  2004-02-25  10  10
    88  2004    February    26  2004-02-26  10  10
    89  2004    February    27  2004-02-27  10  10
    90  2004    February    28  2004-02-28  10  10
    91  2004    February    29  2004-02-29  10  10

我想做aggregate（）所做的事情

---

有什么想法吗？

如果你有兴趣的话，我有一个混乱的解决方案：）

首先，让我们确保月份已正确排序，并为月/年组合创建一个因子

data$Months<-factor(data$Months, levels=month.name)
data$MY<-interaction(data$Months, data$Years, drop=T)

一旦我们有了这些值，我们就可以找到每个月不重叠的最小值

rr<-Map(function(d,off) {
    d<-as.matrix(d)
    stopifnot(ncol(d)==length(off))
    for(i in seq_along(off)) {
         if(off[i]>0)
             d[1:off[i],i]<-Inf
    }
    apply(d,2,min)
}, ss, offsets)
do.call(rbind,rr)

---

我不确定您需要如何格式化结果，但这至少可以提取您想要的值。

这里有一个不同的策略，创建一个

nonoverlappmin

函数。在这里，我们假设数据已经在每个组中正确排序。我将确保数据正确排序，并将创建一个组合因子，以在一个变量中跟踪月/年

data$Months <- factor(data$Months, levels=month.name)
data$MY <- interaction(data$Months, data$Years, drop=T)

---

该方法使用循环逐步查找最小值，然后将下一个

dist

值替换为

Inf

，这样它们就不会被选为最小值。循环在一组一组的值列表中逐步运行。

我认为有时最好回到基础，而不是试图找到最高效的向量实现。请记住，开发人员时间比CPU时间更重要：P

一个简单的for循环就可以了

read.table("Data.txt", header=T, sep="\t", stringsAsFactors=F) -> Data
result = matrix(ncol=4, nrow=0)
min_indA = -100; min_indB = -100; minA = 100; minB = 100
curMonth = "December"
curYear = 2003
for(i in 1:nrow(Data)) {
    if(curMonth == Data[i,"Months"] & curYear == Data[i,"Years"]) {
        if(Data[i,"A"] < minA & i - min_indA >= 10) {
            minA = Data[i,"A"]
            cur_indA = i
        }
        if(Data[i,"B"] < minB & i - min_indB >= 10) {
            minB = Data[i,"B"]
            cur_indB = i
        }
    } else {
        result = rbind(result, c(curYear, curMonth, minA, minB))
        minA = Data[i,"A"]; minB = Data[i,"B"]; min_indA = cur_indA; min_indB = cur_indB;
        curMonth = Data[i,"Months"]
        curYear = Data[i,"Years"]
    }
}
result = rbind(result, c(curYear, curMonth, minA, minB))

---

我认为解决这个问题最有效的方法是递归函数

#Load data
require("data.table")
Data <- fread("min10.csv")
Data <- data.table(Data)
Data[,Date:=as.Date(Date)]

---

这个解决方案只有十几行。我们首先将输入数据框拆分为数据框列表

ym

，每个数据框代表一年/一个月。然后我们对我们希望计算最小值的列进行分析。对于每一列，我们对

ym

组件进行迭代，以便对于每个组件，即对于每个data.frame，我们将其子集为

s

，一个在前一个

minDate

之后至少10天的行数据帧，计算最小行

ix

，更新

minDate

并返回

结果

：

ym <- split(DF, format(DF$Date, "%Y-%m"))
sapply(c("A", "B"), function(col) {
   minDate <- min(DF$Date) - 10
   result <- vector(length = length(ym)) 
   for(i in seq_along(ym)) {
       s <- subset(ym[[i]], Date >= minDate + 10)
       ix <- which.min(s[[col]])
       minDate <- s$Date[ix]
       result[i] <- min(s[[col]][ix])
   }
   setNames(result, names(ym))
})

（我们只使用

DF

的

“Date”

、“

A

”和“

B”

列，因此我们可以先将

DF

减少到该列。）

注意：我们假设此数据框为输入：

DF <-
structure(list(Years = c(2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L), Months = structure(c(1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("December", "February", 
"January"), class = "factor"), Days = c(1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 
20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 1L, 
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 
16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 
29L, 30L, 31L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 
25L, 26L, 27L, 28L, 29L), Date = structure(c(12387, 12388, 12389, 
12390, 12391, 12392, 12393, 12394, 12395, 12396, 12397, 12398, 
12399, 12400, 12401, 12402, 12403, 12404, 12405, 12406, 12407, 
12408, 12409, 12410, 12411, 12412, 12413, 12414, 12415, 12416, 
12417, 12418, 12419, 12420, 12421, 12422, 12423, 12424, 12425, 
12426, 12427, 12428, 12429, 12430, 12431, 12432, 12433, 12434, 
12435, 12436, 12437, 12438, 12439, 12440, 12441, 12442, 12443, 
12444, 12445, 12446, 12447, 12448, 12449, 12450, 12451, 12452, 
12453, 12454, 12455, 12456, 12457, 12458, 12459, 12460, 12461, 
12462, 12463, 12464, 12465, 12466, 12467, 12468, 12469, 12470, 
12471, 12472, 12473, 12474, 12475, 12476, 12477), class = "Date"), 
    A = c(10L, 10L, 10L, 10L, 10L, 10L, 10L, 3L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 5L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 5L, 10L, 10L, 10L, 10L, 10L, 10L, 7L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L), B = c(10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 4L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 4L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 4L, 10L, 10L, 10L, 10L, 10L, 10L, 6L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L)), .Names = c("Years", "Months", 
"Days", "Date", "A", "B"), row.names = c("1", "2", "3", "4", 
"5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", 
"16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", 
"27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", 
"38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", 
"49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", 
"60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", 
"71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", 
"82", "83", "84", "85", "86", "87", "88", "89", "90", "91"), class = "data.frame")

---

DF谢谢，但它给了我--Error:length（vals）==length（groups）不是TRUE@PanagiotisO. 这是一个防御性检查，我这样做是为了确保您传递了正确的值。你输入的向量名称正确吗？您是否如我所述创建了
MY
变量？您尝试传递给函数的向量长度是多少？是否有办法同时获得所有列的矩阵/列表？我这么问是因为我不仅仅有“A”和“B”。我有150列数据。@PanagiotisO。如果您有很多很多列，我使用split的方法可能是一个更有效的选择。否则，您也可以在data.frame列的列表上使用此函数。
data$Months <- factor(data$Months, levels=month.name)
data$MY <- interaction(data$Months, data$Years, drop=T)

nonoverlapmin <- function(vals, groups, dist) {
    stopifnot(length(vals)==length(groups))
    groups<-ordered(groups)
    r <- numeric(nlevels(groups))
    names(r) <- levels(groups)
    for (v in levels(groups)) {
        i <- which.min(vals[groups<=v])
        r[v] <- vals[i]
        vals[ 1:min(max(i+dist, max(which(groups==v))),length(vals))]<-Inf
    }
    r
}

nonoverlapmin(data$A, data$MY, 10)
# December.2003  January.2004 February.2004 
#             3             5             7 

nonoverlapmin(data$B, data$MY, 10)
# December.2003  January.2004 February.2004 
#             4             4             6 

read.table("Data.txt", header=T, sep="\t", stringsAsFactors=F) -> Data
result = matrix(ncol=4, nrow=0)
min_indA = -100; min_indB = -100; minA = 100; minB = 100
curMonth = "December"
curYear = 2003
for(i in 1:nrow(Data)) {
    if(curMonth == Data[i,"Months"] & curYear == Data[i,"Years"]) {
        if(Data[i,"A"] < minA & i - min_indA >= 10) {
            minA = Data[i,"A"]
            cur_indA = i
        }
        if(Data[i,"B"] < minB & i - min_indB >= 10) {
            minB = Data[i,"B"]
            cur_indB = i
        }
    } else {
        result = rbind(result, c(curYear, curMonth, minA, minB))
        minA = Data[i,"A"]; minB = Data[i,"B"]; min_indA = cur_indA; min_indB = cur_indB;
        curMonth = Data[i,"Months"]
        curYear = Data[i,"Years"]
    }
}
result = rbind(result, c(curYear, curMonth, minA, minB))

     [,1]   [,2]       [,3] [,4]
[1,] "2003" "December" "3"  "4" 
[2,] "2004" "January"  "5"  "4" 
[3,] "2004" "February" "7"  "6"

#Load data
require("data.table")
Data <- fread("min10.csv")
Data <- data.table(Data)
Data[,Date:=as.Date(Date)]

#Build recursive function
findmin10 <- function(Data,Var){

  Data$Var1 <- get(Var,Data)

  #Find min date for value A
  Data[,minVar:=min(Var1),by=c("Years","Months")]
  Data[,minVarDate:=(Var1==minVar)*1]
  Summ <- Data[minVarDate==1][,ord:=.I]
  Summ[,Date.Diff:=c(NA,head(as.numeric(Date[ord+1]-Date[ord]),-1))]
  To.Delete.Date <- Summ[Date.Diff<10]$Date

  #Utilize recursion until 10 day spacing requirement is met
  if (length(To.Delete.Date)!=0){
    Data <- Data[!Date%in%To.Delete.Date]
    findmin10(Data,Var=Var)
  } else {
    return(Summ[,list(Years,Months,Var1,VarName=Var)])
  }
}

#Run through multiple variables you want to find the min 10 for
outtable <- rbindlist(lapply(c("A","B"),FUN=function(x) findmin10(Data=Data,Var=x)))

#Cast it out to make it look like desired result
library("reshape2")
dcast.data.table(outtable,Years+Months~VarName,value.var="Var1")

#    Years   Months A B
# 1:  2003 December 3 4
# 2:  2004 February 7 6
# 3:  2004  January 5 4

ym <- split(DF, format(DF$Date, "%Y-%m"))
sapply(c("A", "B"), function(col) {
   minDate <- min(DF$Date) - 10
   result <- vector(length = length(ym)) 
   for(i in seq_along(ym)) {
       s <- subset(ym[[i]], Date >= minDate + 10)
       ix <- which.min(s[[col]])
       minDate <- s$Date[ix]
       result[i] <- min(s[[col]][ix])
   }
   setNames(result, names(ym))
})

        A B
2003-12 3 4
2004-01 5 4
2004-02 7 6

DF <-
structure(list(Years = c(2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 
2004L, 2004L, 2004L, 2004L), Months = structure(c(1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("December", "February", 
"January"), class = "factor"), Days = c(1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 
20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 1L, 
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 
16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 
29L, 30L, 31L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 
25L, 26L, 27L, 28L, 29L), Date = structure(c(12387, 12388, 12389, 
12390, 12391, 12392, 12393, 12394, 12395, 12396, 12397, 12398, 
12399, 12400, 12401, 12402, 12403, 12404, 12405, 12406, 12407, 
12408, 12409, 12410, 12411, 12412, 12413, 12414, 12415, 12416, 
12417, 12418, 12419, 12420, 12421, 12422, 12423, 12424, 12425, 
12426, 12427, 12428, 12429, 12430, 12431, 12432, 12433, 12434, 
12435, 12436, 12437, 12438, 12439, 12440, 12441, 12442, 12443, 
12444, 12445, 12446, 12447, 12448, 12449, 12450, 12451, 12452, 
12453, 12454, 12455, 12456, 12457, 12458, 12459, 12460, 12461, 
12462, 12463, 12464, 12465, 12466, 12467, 12468, 12469, 12470, 
12471, 12472, 12473, 12474, 12475, 12476, 12477), class = "Date"), 
    A = c(10L, 10L, 10L, 10L, 10L, 10L, 10L, 3L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 5L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 5L, 10L, 10L, 10L, 10L, 10L, 10L, 7L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L), B = c(10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 4L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 4L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 4L, 10L, 10L, 10L, 10L, 10L, 10L, 6L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
    10L, 10L, 10L, 10L, 10L)), .Names = c("Years", "Months", 
"Days", "Date", "A", "B"), row.names = c("1", "2", "3", "4", 
"5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", 
"16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", 
"27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", 
"38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", 
"49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", 
"60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", 
"71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", 
"82", "83", "84", "85", "86", "87", "88", "89", "90", "91"), class = "data.frame")