Как вычислять периоды с момента 200-периодного максимума акции
-
28-10-2019 - |
Вопрос
Я хотел бы рассчитать количество периодов, которые прошли после 200 периода периода одномерного временного ряда. Например, вот заключительная цена шпиона:
require(quantmod)
getSymbols("SPY",from='01-01-1900')
Data <- Cl(SPY)
Теперь я могу найти 200-периодические максимумы этой серии, используя Lag
функция в QuantMod:
periodHigh <- function(x,n) {
Lags <- Lag(x,1:n)
High <- x == apply(Lags,1,max)
x[High]
}
periodHigh(Data, 200)
Но теперь я застрял. Как слияние этого обратно на оригинальную серию (Data
) и рассчитать, для каждой точки в серии, сколько периодов прошло с момента предыдущего максимума N-периода?
Решение 3
Я отредактировал код из предыдущих ответов, так что это функции, которые принимают одинаковые входы (одномерный временной ряд) и возвращали тот же выход (вектор дней с момента последнего N-Day High):
daysSinceHigh1 <- function(x,n) {
as.vector(n-rollapply(x, n, which.max))
}
daysSinceHigh2 <- function(x, n){
apply(embed(x, n), 1, which.max)-1
}
Вторая функция кажется самой быстрой, но они дают немного разные результаты:
> getSymbols("^GSPC",from='01-01-1900')
[1] "GSPC"
> system.time(x <- daysSinceHigh1(Cl(GSPC), 200))
user system elapsed
0.42 0.00 0.42
> system.time(y <- daysSinceHigh2(Cl(GSPC), 200))
user system elapsed
0.24 0.00 0.24
> all.equal(x,y)
[1] "Mean relative difference: 0.005025126"
При ближайшем рассмотрении, кажется, что в 1 -й функции есть некоторые странные края:
data <- c(1,2,3,4,5,6,7,7,6,5,6,7,8,5,4,3,2,1)
answer <- c(0,0,0,0,1,2,3,0,0,1,2,3,4,4)
x <- daysSinceHigh1(data, 5)
y <- daysSinceHigh2(data, 5)
> x
[1] 0 0 0 1 2 3 4 4 0 1 2 3 4 4
> y
[1] 0 0 0 0 1 2 3 0 0 1 2 3 4 4
> answer
[1] 0 0 0 0 1 2 3 0 0 1 2 3 4 4
> all.equal(x,answer)
[1] "Mean relative difference: 0.5714286"
> all.equal(y,answer)
[1] TRUE
Следовательно, кажется, что вторая функция (основанная на коде Андри) лучше.
Другие советы
Эта небольшая функция возвращает список с:
high
индекс номер высокой датыrecentHigh
Индекс номер самого последнего высокого дняdaysSince
количество дней с момента последнего максимумаdata
XTS объект только с высокими днями. Полезно для построения.
Код:
daysSinceHigh <- function(data, days){
highs <- days-1+which(apply(embed(data, days), 1, which.max)==1)
recentHigh <- max(highs)
daysSince <- nrow(data) - recentHigh
list(
highs=highs,
recentHigh = recentHigh,
daysSince = daysSince,
data=data[highs, ])
}
Результаты, достижения:
daysSinceHigh(Data, 200)$daysSince
[1] 90
plot(Data)
points(daysSinceHigh(Data, 200)$data, col="red")
Ответ на ваш пересмотренный вопрос:
require(zoo)
x <- sample(300:500, 1000, replace=TRUE)
str(rollapply(x, 200, function(x) which.max(x)))
# int [1:801] 14 13 12 11 10 9 8 7 6 5 ...
plot(x)
plot(200:1000, rollapply(x, 200, function(x) 200-which.max(x)))
Итак, для серии XTS:
plot( rollapply(coredata(Data), 200, function(x) 200-which.max(x)))
Вы, очевидно, не можете ничего объединить в первые 200 дат, если вы не примените более слабое определение максимума катания. (В другом сеансе, включающем данные «сдвига», я показал, как использовать Enced, чтобы подумать о том, что «следующий» периоды: Преобразование данных в R Но я не знаю, хотите ли вы построить матрицы, которые в 200 раз больше ваших входных данных.)