r - 超过最大值百分比的日期

Question

我有一个如下所示的数据框：

df <- data.frame(site=paste0('site', sort(rep(1:5, 20))),
                 date=as.Date(paste0(sample(1:28, 100, replace=T), '/', 
                                     sample(1:12, 100, replace=T), '/', 
                                     2013), 
                              '%d/%m/%Y'), 
                 count=rep(seq(1, 1000, length.out=20), 10))

对于每个站点，我需要 > 500 的最早日期。count换句话说，计数大于最大计数的 50% 的第一个日期。

Answer 1

解决方案plyr。

如果您想要 count > 500 的第一个日期：

ddply(df, .(site), summarise, date=min(date[count>500]))

如果您想要 count > 50%*max(count) 的第一个日期（对于每个站点）：

ddply(df, .(site), summarise, date=min(date[count>max(count)*0.5]))

这里都给出相同的结果：

   site       date
1 site1 2013-01-15
2 site2 2013-02-04
3 site3 2013-03-13
4 site4 2013-02-04
5 site5 2013-01-07