我想让我的gecpLU函数像arma::LU一样工作,输入L、U和P,arma::LU函数修改输入的L、U和P矩阵,而不是返回L、U和P。arma::colvec example(arma::colvec X) { return X;我知道当暴露在R中时,这不会改变输入,因为arma对象是R对象的副本::mat L, arma::mat U, arma::mat P, arma::mat Q, arm
我想使用Rcpp提高一些R代码的速度。然而,我对C++的了解很少。因此,我查看了Rcpp提供的文档,以及Dirk Eddelbuttel网站上提供的其他文档。读完所有的东西后,我试着执行一个我用R写的简单循环。不幸的是,我不能这么做。下面是R函数: p = ncol(beta) H = matrix(