时间序列分析-R语言-随机游走以及回归画图

###随机游走以及回归
library(TSA)
win.graph(width=4.8, height=3,pointsize=10)
data(rwalk)
plot(rwalk,type='o',ylab='数值',main='随机游走')
model1=lm(rwalk~time(rwalk))
summary(model1)
abline(model1)
Call:
lm(formula = rwalk ~ time(rwalk))

Residuals:
     Min       1Q   Median       3Q      Max 
-2.70045 -0.79782  0.06391  0.63064  2.22128 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) -1.007888   0.297245  -3.391  0.00126 ** 
time(rwalk)  0.134087   0.008475  15.822  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 1.137 on 58 degrees of freedom
Multiple R-squared:  0.8119,	Adjusted R-squared:  0.8086 
F-statistic: 250.3 on 1 and 58 DF,  p-value: < 2.2e-16