#⽣成两个⻓度为1000的时间序列
set.seed(20140623)
N <- 1000
x <- cumsum(rnorm(N))
gamma <- 0.7
#固定随机数种⼦
#定义模拟的⻓度
#模拟⼀个正态随机游⾛
#设置初始的参数值
y <- gamma * x + rnorm(N)
#模拟协整序列
plot(x, type='l')
lines(y,col="red")

生活、成长记录
#⽣成两个⻓度为1000的时间序列
set.seed(20140623)
N <- 1000
x <- cumsum(rnorm(N))
gamma <- 0.7
#固定随机数种⼦
#定义模拟的⻓度
#模拟⼀个正态随机游⾛
#设置初始的参数值
y <- gamma * x + rnorm(N)
#模拟协整序列
plot(x, type='l')
lines(y,col="red")