次の 2 つの場合data.frames
:
histdevs <- data.frame(a=c(3,8,2), b=c(1,4,6), c=c(6,1,5), d=c(2,4,8), e=c(5,1,5))
newdatadevs <- data.frame(a=c(3), b=c(4), c=c(32), d=c(1), e=c(6))
時間の経過に伴う標準偏差の変動を計算したいので、たとえば 3 ブロックで計算したいので、次のループで rollapply を使用します。
for (i in 1:5) {
Over <- mean(histdevs[,i])+3*sd(histdevs[,i])
Under <- mean(histdevs[,i])-3*sd(histdevs[,i])
checkover<-rollapply(newdatadevs[,i],3,function(x)(ifelse(all(x>Over),1,0)))
checkunder<-rollapply(newdatadevs[,i],3,function(x)(ifelse(all(x<Under),1,0)))
}
次のエラーが表示されます。
Error in seq.default(start.at, NROW(data), by = by) :
wrong sign in 'by' argument
もちろん、同じことを実行できる代替機能も大歓迎です。