Tengo un marco de datos de stock diariamente desde el cual trato de obtener el Ratio de Sharpe
> head(df)
Date .SXQR .SXTR .SXNR .SXMR .SXAR .SX3R .SX6R .SXFR .SXOR .SXDR
1 2000-01-03 364.94 223.93 489.04 586.38 306.56 246.81 385.36 403.82 283.78 455.39
2 2000-01-04 345.04 218.90 474.05 566.15 301.13 239.24 374.64 390.41 275.93 434.92
3 2000-01-05 338.22 215.88 464.20 542.29 298.22 239.55 373.26 383.48 272.54 430.05
4 2000-01-06 343.13 218.18 470.82 529.33 300.69 249.75 377.26 383.48 272.47 434.15
5 2000-01-07 349.46 220.10 478.87 531.65 306.50 255.17 381.19 390.23 273.76 447.02
6 2000-01-10 356.20 223.01 484.07 581.82 310.84 252.75 387.74 393.75 278.76 453.80
Con SharpeRatio.annualized
También probé con una de las bibliotecas ya construidas SharpeRatio.annualized
en df2 cuyos nombres son formatos de fecha estándar:
>df2 <- df
>rownames(df2) <- df2$Date
>SharpeRatio.annualized(df2, Rf = 0, scale = NA, geometric=TRUE)
Error in checkData(R) :
The data cannot be converted into a time series. If you are trying to pass in names from a data object with one column, you should use the form 'data[rows, columns, drop = FALSE]'. Rownames should have standard date formats, such as '1985-03-15'.
Sin embargo, parece que no se tienen en cuenta los nombres como fecha.
SharpeRatio a la mano
Intenté hacerlo usando el método de esta publicación de blog que calcula la función Return.excess de PerformanceAnalytics, luego calculo la proporción de Sharpe manualmente. Ya tengo las devoluciones mensuales de:
percent_change2 <- function(x)last(x)/first(x) - 1
monthly_return <- df %>%
group_by(gr = floor_date(Date, unit = "month")) %>%
summarize_at(vars(-Date, -gr), percent_change2) %>%
ungroup() %>%
select(-gr) %>%
as.matrix()
Sin embargo, cuando llegamos a la Ratio de Sharpe usando:
portfolio_excess_returns <- Return.excess(portfolio_monthly_returns,
Rf = .0003)
sharpe_ratio_manual <- round(
mean(portfolio_excess_returns) / StdDev(portfolio_excess_returns), 4
)
Tengo el error siguiente
> portfolio_excess_returns <- Return.excess(monthly_return,
+ Rf = .0003)
Error in checkData(R) :
The data cannot be converted into a time series. If you are trying to pass in names from a data object with one column, you should use the form 'data[rows, columns, drop = FALSE]'. Rownames should have standard date formats, such as '1985-03-15'.
De hecho, ya no tengo las fechas:
.SXQR .SXTR .SXNR .SXMR .SXAR .SX3R
[1,] -0.09048063 -0.070602420 -0.05688696 0.138527917 -0.079429802 -0.07370042
[2,] -0.02547240 -0.017752597 0.08995378 0.227683167 -0.007488398 0.02070630
[3,] 0.04344800 0.072023472 -0.04419912 -0.118953102 0.042191076 0.05908376
[4,] 0.07704088 -0.001580195 0.03846154 -0.009596822 -0.030854538 0.01882887
[5,] -0.02398056 -0.058094333 -0.05729563 -0.193718178 -0.034808588 0.02313015
[6,] 0.05252934 0.007013127 -0.04359662 -0.016661155 -0.051900140 0.02085395
¿Cómo puedo agregarlos a la función month_returns para calcular la SharpeRatio?
Reproducir el error
El archivo de datos es aquí.
Para crear df solo necesitas:
library(dplyr)
library(lubridate)
df <- read.xlsx("Data.xlsx", sheet = "Sector-STOXX600", startRow = 2,colNames = TRUE, detectDates = TRUE, skipEmptyRows = FALSE)
df[2:19] <- data.matrix(df[2:19])