我有50个矩阵(50,000行乘10列),我正在尝试为蒙特卡洛模拟填充。我在一个循环中创建了所有的矩阵,它们被称为mCMatrix1,mCMatrix2等。我想在一个循环中填充矩阵,如下所示:for (j in 50000){
num <- mu + tR %*% rnorm(10) # returns a 10 rowmCMatrixC"i"[]= num[,1] # basically rotates
我有一个漂亮的三角n x n矩阵,比如L (对于下三角),我想要解一个类似的系统其中B和X是n x k矩阵(即:我想要解一个多个右手边的三角形线性系统)。虽然有许多功能可以解决类似的问题::采用密集格式的三角形矩阵,并对多个右手边进行求解。
我真正需要的是两者的结合。我希望有一个函数接受打包的三角格式,如stpsv(),也可以接受多个右手边,如strsm()。但是,似乎没有这样的功能随时