假设我有一个向量:
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
会发生什么。我能得到我想要的东西吗?
答案 0 :(得分:0)
没有!请注意,rmvnorm
会在10x50矩阵的行中返回10个随机向量。因此,您添加Q + vars
,R将逐列回收Q
以获得10x50矩阵。但你可以用这个得到你想要的东西:
t(t(vars)+Q)
其行是随mean==Q
多变量法线密度和cov.mat
给出的方差 - 协方差矩阵的随机观察值。