0

我有每日数据（4011 天）和指标（1 个工作日，2 个周末）。我想用相应的指标找到每周最大值。例如（假设）数据是：

mydat <- matrix(c(0.027,0.034,0.019,0.021,0.026,0.024,0.058,0.026,0.064,
     0.066,0.026,0.101,0.069,0.054,rep(2,2),rep(1,5),rep(2,2),rep(1,5)), ncol=2)

我尝试使用以下代码。我设法获得最大序列（在这种情况下，每周最大值），但我不希望指标中的最大序列。这是代码

week.max <- function(vec){    
   if(length(vec[is.na(vec)]) == 7){       
   return(NA)     
   } 
   else{
     return(max(vec, na.rm = T))
   } 
}

max.week.dat <- apply(mydat, 2, function(x) tapply(x, rep(1:(length(x)/7), 
  each=7, len=length(x)), week.max))

和结果

matrix(c(0.058,0.101,2,2),ncol=2)

我想要这样的输出：

matrix(c(0.058,0.101,1,1),ncol=2)

提前谢谢了。

4

1 回答 1

1

这是数据（第三周多出一天）

mydat <- data.frame(value = c(0.027,0.034,0.019,0.021,0.026,0.024,0.058,0.026,0.064,
     0.066,0.026,0.101,0.069,0.054,0.95), ind = c(rep(2,2),rep(1,5),rep(2,2),rep(1,5),2))

你的功能

week.max <- function(vec){    
  if(length(vec[is.na(vec)]) == 7){       
return(NA)     
 } 
 else{
return(max(vec, na.rm = T))
 } 
}

添加周信息

mydat$week <- c(rep(1:2,each=7),3)

使用与此处相同的解决方案

library(plyr)
ddply(mydat, .(week), subset, subset = value==week.max(value), select = -week)
于 2012-05-23T05:06:46.380 回答