Для цикла с функцией скользящего / скользящего среднего? - PullRequest
0 голосов
/ 17 апреля 2019

По существу (в R) я хочу применить функцию скользящего среднего за период времени (например, переменные даты и времени), чтобы увидеть, как конкретная метрика изменяется со временем.Однако метрика сама по себе является функцией.Баллы могут быть 1 (профессиональная), 0 (нейтральная) или -1 (отрицательная).Функция для метрики:

function(pro, neg, total) {
x <- (pro / total) * 100
y <- (neg / total) * 100
x - y
}

Таким образом, процентное значение 1 минус процентное значение -1 - это значение метрики.

Учитывая временные метки для каждого записанного балла, я хочу оценитьметрика как скользящее среднее по всем строкам.Я предполагал, что для цикла будет лучшим способом применить это, но я застрял в том, как это сделать.

У кого-нибудь есть какие-нибудь мысли / советы?

1 Ответ

2 голосов
/ 17 апреля 2019

Как уже упоминалось в комментариях, rollapply() из zoo - хороший вариант.Я позволил себе создать пример данных, извиняюсь, если они не похожи на ваши.

library(zoo)

f <- function(x, l) {
    p <- sum(x == 1) / l
    n <- sum(x == -1) / l
    (p - n)*100
}

# Or more efficiently
f <- function(x, l=length(x)) {
    (sum(x)/l)*100
}

set.seed(1)
N <- 25
dtf <- data.frame(time=as.Date(15000+(1:N)), score=sample(-1:1, N, rep=TRUE))

score <- read.zoo(dtf)
l <- 8
zts <- cbind(score, rolling=rollapply(score, l, f, l, fill=NA))

zts
#            score rolling
# 2011-01-27    -1      NA
# 2011-01-28     0      NA
# 2011-01-29     0      NA
# 2011-01-30     1    12.5
# 2011-01-31    -1    25.0
# 2011-02-01     1    12.5
# 2011-02-02     1     0.0
# 2011-02-03     0   -25.0
# 2011-02-04     0     0.0
# 2011-02-05    -1   -12.5
# 2011-02-06    -1   -12.5
# 2011-02-07    -1   -12.5
# 2011-02-08     1     0.0
# 2011-02-09     0    25.0
# 2011-02-10     1    37.5
# 2011-02-11     0    62.5
# 2011-02-12     1    62.5
# 2011-02-13     1    50.0
# 2011-02-14     0    37.5
# 2011-02-15     1    25.0
# 2011-02-16     1     0.0
# 2011-02-17    -1      NA
# 2011-02-18     0      NA
# 2011-02-19    -1      NA
# 2011-02-20    -1      NA
...