Моделирование данных многомерного временного ряда с использованием 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)