我正在使用 PerformanceAnalytics 包来分析一些月度返回。 charts.RollingRegression 应该根据某个基准绘制 n 个月的滚动回归。
数据只是 6 个返回系列,从 4 月 8 日到 2014 年 12 月,试图回归 SPY。
indexReturns <- read.table("quantIndices.csv", stringsAsFactors = FALSE, sep = ",", fill = TRUE, row.names = 1, header=TRUE)
hfIndexReturns <- read.table("quantHFIndices.csv", stringsAsFactors = FALSE, sep = ",", fill = TRUE, row.names = 1, header=TRUE)
peerReturns <- read.table("quantPeers.csv", stringsAsFactors = FALSE, sep = ",", fill = TRUE, row.names = 1, header=TRUE)
splits <- as.data.frame(strsplit(rownames(indexReturns), "/"))
rownames(indexReturns) <- unname(sapply(splits, function(x) paste0(x[3], "-", x[1], "-", x[2])))
splits <- as.data.frame(strsplit(rownames(peerReturns), "/"))
rownames(peerReturns) <- unname(sapply(splits, function(x) paste0(x[3], "-", x[1], "-", x[2])))
Ret <- xts(peerReturns, order.by = as.Date(row.names(peerReturns)))
Rb <- xts(indexReturns, order.by = as.Date(row.names(indexReturns)))
charts.RollingRegression(Ret, Rb[,2, drop = FALSE], Rf = 0.001, na.pad = TRUE)
这将产生以下图表:
我希望它省略“无意义的)前 12 个月,但没有关于如何完成的文档,我能找到的该图表的任何其他描述看起来像这样:
查看源代码,在函数的主要内容中,我看到:
for (column.a in 1:columns.a) {
for (column.b in 1:columns.b) {
merged.assets = merge(Ra.excess[, column.a, drop = FALSE],
Rb.excess[, column.b, drop = FALSE])
if (attribute == "Alpha")
column.result = rollapply(na.omit(merged.assets),
width = width, FUN = function(x) lm(x[, 1,
drop = FALSE] ~ x[, 2, drop = FALSE])$coefficients[1],
by = 1, by.column = FALSE, fill = na.pad, align = "right")
if (attribute == "Beta")
column.result = rollapply(na.omit(merged.assets),
width = width, FUN = function(x) lm(x[, 1,
drop = FALSE] ~ x[, 2, drop = FALSE])$coefficients[2],
by = 1, by.column = FALSE, fill = na.pad, align = "right")
if (attribute == "R-Squared")
column.result = rollapply(na.omit(merged.assets),
width = width, FUN = function(x) summary(lm(x[,
1, drop = FALSE] ~ x[, 2, drop = FALSE]))$r.squared,
by = 1, by.column = FALSE, align = "right")
column.result.tmp = xts(column.result)
colnames(column.result.tmp) = paste(columnnames.a[column.a],
columnnames.b[column.b], sep = " to ")
column.result = xts(column.result.tmp, order.by = time(column.result))
if (column.a == 1 & column.b == 1)
Result.calc = column.result
else Result.calc = merge(Result.calc, column.result)
}
}
我们可以看到没有 na.pad 被传递到最终的“R-Squared”函数,这导致了我希望看到的前两个图表的图表。我想解决这个问题,但我无法编辑包代码。我尝试使用“assignInNamespace”,但它不起作用。该功能似乎可以工作,但功能代码在包中没有变化。我也想删除图表中的前导空格,但如果你们可以让我知道如何编辑它,或者知道任何解决方法,请告诉我。 (谢谢!你们是神!)
哦! PS - 为什么我的包版本似乎是唯一一个有这个问题的版本???为什么我的图表在默认情况下看起来不正确?
编辑:这不是这个包中唯一可疑的代码。我一直让事情中断并且无法正常工作,因为它似乎被记录在案(
Error in R[, nc] - coredata(Rf) : non-numeric argument to binary operator
似乎发生在所有其他函数调用中。)有人对这种类型的东西有更好的包有什么建议吗? 最佳答案
数据的子集化应该在传递给 charts.RollingRegression
函数之前完成。强大的 xts
提供了这个功能:charts.RollingRegression(Ret["2009-04::",], Rb["2009-04::",2, drop = FALSE], Rf = 0.001, na.pad = TRUE)
您可以通过 ?subset.xts
查看 R 中的帮助页面,了解有关如何使用 xts 进行子集化的更多信息。
关于r - PerformanceAnalytics charts.RollingRegression 绘制初始窗口值。我怎样才能让它不那样做?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/28968597/