如何提高产生multi-norm random variable的速度

swimming

我把程序改了，原来10分钟现在只要1分多钟，我很开心，大体上就是把循环改成矩阵操作。但是下面一段程序大家有什么办法改吗？

我有一千个矩阵存在list中，还有一千个向量，用他们产生一千个随机变量

<br />
for (i in seq(1, 1000, 1)) {<br />
    beta.hat[, i] <- rmvnorm(1, mean = temp.beta.hat[, index[i]],<br />
        sigma = Sigma.square[i] * solve(temp.temp1[[index[i]]]))<br />
}<br />

</p>

怎样不写成循环形式？

temp.temp1是list

还有我的

砸不起作用。。。。。
</p>

rtist

It is the inversion that is slow, not the generation of random numbers.

swimming

回复第2楼的 Rtist：谢谢有什么办法没，还有我的

砸不起作用？
</p>

rtist

Like finding cholesky decomposition, or sth similar.

The backward slash should be forward slash.

swimming

回复第4楼的 Rtist：有解决办法吗？还有可以不写成循环形式吗？ you must be in USA. Where, haha

rtist

无语了[s:12]

swimming

回复第6楼的 Rtist：是不是英文我没看懂，解释下吧，大哥

foison

回复第1楼的 swimming：用mapply

回复第6楼的 Rtist：想起我们这里以前一个故事，一个老头子突然失忆了，被发现时没有任何身份证明，仅能说流利但是有口音的英语，请来语言专家来分析结果说来自荷兰的，当时我就无语了。后来老头子醒悟了，说来自苏格兰...

Ihavenothing

回复第1楼的 swimming：

code标签的写法是[code][/code]而不是[code][\code]。

swimming

回复第8楼的 Feng Li：多谢，我试试，

swimming

回复第8楼的 Feng Li：但是我前面所有mean存在一个矩阵中，也可以用mapply逐行做吗

idoget

回复第11楼的 swimming：

<br />
mu=list(rnorm(4),rnorm(4))<br />
sig=list(cor(matrix(rnorm(40),10,4)),cor(matrix(rnorm(40),10,4)))<br />
mapply(rmvnorm,mean=mu,sigma=sig,MoreArgs=list(n=1),SIMPLIFY=F)<br />

</p>

longoR

回复第8楼的 Feng Li：lost...