使用 r 模擬多變量時間序列數據

Question

我想模擬一個適合表中給出的年份的多元時間序列數據。 我怎樣才能做到這一點？

我使用以下代碼

Z <- rnorm(100, mean = 0, sd = 1.5)

# process simulation
X <- c()
for (i in 1:length(Z)) {

# my table
my_table <- 1982:2008 

Answer 1

您將需要序列的協方差或相關性。 如果沒有相關性，那將很容易。

# No correlation: 
my_table <- 1982:2008
num.series <- 4

sdata=matrix(rnorm(length(my_table)*num.series,0,sd=1.5),ncol=num.series)
sdf=data.frame(my_table,sdata)
summary(sdf)

# with covariance
library(MASS)
Sigma <- matrix(c(10,3,3,3,3,5,3,3,3,3,8,3,3,3,3,15)/10,ncol=4)
Sigma
smdata=mvrnorm(n = length(my_table), rep(0, num.series), Sigma)
scdf=data.frame(my_table,smdata)
summary(scdf)