从现有变量创建相关错误项

时间:2013-09-11 13:06:48

标签: r

假设我有一个向量:

Q<-rnorm(50,mean=0,sd=30)

我想创建10个由Q+error term组成的变量,这样10个变量中的每个变量都有相关的误差项,导致10个变量之间的相关性,比如说.8。

通常我会使用以下方法创建相关变量:

vars<-rmvnorm(n=10, mean=c(rep(0,50)),cov.mat, method="chol")

如果我添加Q+vars会发生什么。我能得到我想要的东西吗?

1 个答案:

答案 0 :(得分:0)

没有!请注意,rmvnorm会在10x50矩阵的中返回10个随机向量。因此,您添加Q + vars,R将逐列回收Q以获得10x50矩阵。但你可以用这个得到你想要的东西:

t(t(vars)+Q)

其行是随mean==Q多变量法线密度和cov.mat给出的方差 - 协方差矩阵的随机观察值。