lm(数据〜时间)和tslm(数据〜趋势)之间有什么区别

时间:2018-01-01 23:29:26

标签: r regression time-series forecasting

我观察到两种方法的结果都不同。为什么是这样?我知道lm上发生了什么,但无法弄清楚tslm下发生了什么。

> library(forecast)
> set.seed(2)
> tts <- ts(100*runif(1200)+seq(1:1200)*0.1, frequency=12, start=c(2000,1))
> lm(tts~time(tts))

Call:
lm(formula = tts ~ time(tts))

Coefficients:
(Intercept)    time(tts)  
  -2400.365        1.225  

> tslm(tts~trend)

Call:
tslm(formula = tts ~ trend)

Coefficients:
(Intercept)        trend  
    48.9350       0.1021  

1 个答案:

答案 0 :(得分:1)

运行以下三个命令:

predict(lm(tts~time(tts)))
predict(tslm(tts~time(tts)))
all.equal(predict(lm(tts~time(tts))), predict(tslm(tts~trend)))

你会说服自己,他们是完全相同的。如果输出相同,那么lm回归的X变量,即

time(tts) 

必须是

的线性转换
trend

最简单的猜测:

tmp <- time(tts)*12
lm(tts~tmp)

具有与tslm系数相同的系数。所以趋势只是

12*time(tts)

即。趋势是从0年开始经过的时间的(整数)计数,以月为单位。 time(tts)是从0年开始经过的时间计数,以年为单位。