大家好,所以我的问题可能是统计数据或编程问题。我有两个xts时间序列,主要是重叠的时间段,我只是绘制他们的对数差异的回归:
logdiff <- merge.xts(diff(log(ts1)),diff(log(ts2)))
plot(logdiff[,1],logdiff[,2])
abline(lm(logdiff[,1]~logdiff[,2]),col=2)
给了我这个情节
所以只是在一个直观的层面上我宁愿回归线适合更广泛的数据点,即使这个结果它给我的技术上是最小二乘的正确数据点。是否有任何内置功能可以做到“更广泛的回归”,还是我不得不诉诸手工捏造?
答案 0 :(得分:7)
我认为你将y绘制为x的函数,但将x作为y的函数回归。
尝试abline(lm(logdiff[,2]~logdiff[,1]),col=2)
- 是的,使用列名而不是索引是一个好主意。