- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我在r中有一段代码运行得很慢,因此我希望在rcpp中重写该函数,但是,这是我第一次尝试使用Rcpp,而且我无法获取代码来进行编译或运行。
在R中,我要重写的部分是
drawsamples<-lapply(1:numstudies,function(v){
temp<-t(mapply(rmvn,n=1,mu=finalmu_b2_l[[v]],sigma=finalcov_b2_l[[v]]))
rownames(temp)<-id[[v]]
temp
})
numstudies
,并且每个元素的长度应为
n[v]
列的
ql
行的
matrix。这些矩阵的每一行都应来自多元正态分布,但是每一行的绘制将由不同的均值 vector
finalmu_b2_l
和不同的协方差矩阵
finalcov_b2_l
控制。
finalmu_b2_l
和
finalcov_b2_l
都是嵌套列表,其长度为
numstudies
的顶层,其每个元素都是一个长度为
n[v]
的列表。
#include <RcppDist.h>
#include <RcppArmadillo.h>
using namespace Rcpp;
// [[Rcpp::depends(RcppDist,RcppArmadillo)]]
// [[Rcpp::export]]
arma::mat draw_b(const int n,
const int ql,
const Rcpp::List mu,
const Rcpp::List cov) {
arma::mat draws = arma::zeros(n,ql);
for (int iter = 1; iter <= n; iter++) {
draws.row(iter) = rmvnorm(1, mu[iter], cov[iter]);
}
return draws;
}
sourceCpp
函数通过获取代码来保存代码所在的.cpp文件),都会收到错误消息:
invalid initialization of reference of type 'const::mat&'
draws
矩阵还是填充它的问题?
#include <RcppDist.h>
#include <RcppArmadillo.h>
using namespace Rcpp;
// [[Rcpp::depends(RcppDist,RcppArmadillo)]]
//' @keywords internal
// [[Rcpp::export]]
arma::mat draw_b(const int & n,
const int & ql,
const Rcpp::List & mu,
const Rcpp::List & cov) {
arma::mat draws = arma::zeros(n,ql);
for (int iter = 0; iter <= n; iter++ ) {
draws.row(iter) = rmvnorm(1, mu_temp, cov_temp);
}
return draws;
}
error: Mat::init(): requested size is not compatible with column vector layout
Error in draw_b(n = n, ql = ql, mu = mu_example, cov = cov_example) :
Mat::init(): requested size is not compatible with column vector layout
list(n = 10, ql = 2, mu_example = list(c(0.342909965003207, -0.788070875792469
), c(-0.00499810116271365, 0.0114865660452949), c(-0.149753928200309,
0.344162379034614), c(0.335176829763227, -0.770298692761465),
c(0.254206123984596, -0.584212951520601), c(0.379893097582703,
-0.873064992779416), c(0.137231089484867, -0.315382566602526
), c(0.405123380985852, -0.931048876501857), c(-0.00505917031396947,
0.0116269143128456), c(-0.0743318653279181, 0.170828451158346
)), cov_example = list(structure(c(0.199912910315971, -0.459437048770092,
-0.459437048770092, 4.49223135519527), .Dim = c(2L, 2L)), structure(c(0.199912910315971,
-0.459437048770092, -0.459437048770092, 4.49223135519527), .Dim = c(2L,
2L)), structure(c(0.199912910315971, -0.459437048770092, -0.459437048770092,
4.49223135519527), .Dim = c(2L, 2L)), structure(c(0.199912910315971,
-0.459437048770092, -0.459437048770092, 4.49223135519527), .Dim = c(2L,
2L)), structure(c(0.199912910315971, -0.459437048770092, -0.459437048770092,
4.49223135519527), .Dim = c(2L, 2L)), structure(c(0.199912910315971,
-0.459437048770092, -0.459437048770092, 4.49223135519527), .Dim = c(2L,
2L)), structure(c(0.199912910315971, -0.459437048770092, -0.459437048770092,
4.49223135519527), .Dim = c(2L, 2L)), structure(c(0.199912910315971,
-0.459437048770092, -0.459437048770092, 4.49223135519527), .Dim = c(2L,
2L)), structure(c(0.199912910315971, -0.459437048770092, -0.459437048770092,
4.49223135519527), .Dim = c(2L, 2L)), structure(c(0.199912910315971,
-0.459437048770092, -0.459437048770092, 4.49223135519527), .Dim = c(2L,
2L))))
drawsampletest<-draw_b(n=n,
ql=ql,
mu=mu_example,
cov=cov_example)
#include <RcppDist.h>
#include <RcppArmadillo.h>
using namespace Rcpp;
// [[Rcpp::depends(RcppDist,RcppArmadillo)]]
//' @keywords internal
// [[Rcpp::export]]
arma::mat draw_b(const int & n,
const int & ql,
const Rcpp::List & mu,
const Rcpp::List & cov) {
arma::mat draws = arma::zeros(n,ql);
for (int iter = 0; iter <= n; iter++ ) {
arma::rowvec mu_temp = mu[iter];
arma::mat cov_temp = cov[iter];
draws.row(iter) = rmvnorm(1, mu_temp, cov_temp);
}
return draws;
}
最佳答案
这是一个基本设置,应执行您想要的操作:
#include <RcppDist.h>
#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppDist,RcppArmadillo)]]
// [[Rcpp::export]]
arma::mat draw_b(const int ql,
const Rcpp::List& mu,
const Rcpp::List& cov) {
int n = mu.size();
arma::mat draws = arma::zeros(n, ql);
for ( arma::uword iter = 0; iter < n; iter++ ) {
draws.row(iter) = rmvnorm(1, mu[iter], cov[iter]);
}
return draws;
}
// [[Rcpp::export]]
Rcpp::List get_list_of_draws(Rcpp::List mu, Rcpp::List Sigma, int ql) {
int numstudies = mu.size();
Rcpp::List res(numstudies);
for ( int iter = 0; iter < numstudies; ++iter ) {
Rcpp::List mu_temp = mu[iter];
Rcpp::List cov_temp = Sigma[iter];
res[iter] = draw_b(ql, mu_temp, cov_temp);
}
return res;
}
res <- get_list_of_draws(mu_example, cov_example, ql)
str(res)
# List of 1
# $ : num [1:10, 1:2] -0.0156 -0.4717 -0.8134 0.5489 0.1215 ...
mu_example
和
cov_example
时,我将它们包装在
list()
中,正如您所说的,它们应该是列表列表。)
关于c++ - 带有列表参数的RCPP函数返回矩阵,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/60986805/
我想将函数参数中的默认值设置为 Rcpp::Function 参数。 只是简单的赋值,Rcpp::Function func = mean , 不可能。它返回错误:no viable conversi
我正在处理需要逐元素矩阵乘法的代码。我试图在 Rcpp 中实现这一点,因为代码需要一些昂贵的循环。我对 Rcpp 还很陌生,可能会遗漏一些东西,但我无法使逐元素矩阵乘法工作。 // [[Rcpp::e
在 C++ 中,我们可以声明一个变量作为引用。 int a = 10; int& b = a; 如果我们设置 b=15 , a 也会改变。 我想在 Rcpp 中做类似的事情。 List X = obj
我正在阅读很棒的 Rcpp vignette关于使用 Rcpp 模块公开 C++ 类和函数。在这种情况下,是否可以创建一个 Rcpp 函数,该函数具有一个类型为 Uniform 的类作为参数之一,并且
我在 R 中有一个命名列表: l = list(a=1, b=2) 我想在 Rcpp 中使用这个列表,并迭代值和名称。理想情况下,它可能类似于(为简洁起见使用 C++11 格式): void prin
这个问题在这里已经有了答案: Rcpp - sourceCpp - undefined symbol (2 个答案) 关闭 4 年前。 我现有的 C 代码由三个文件组成:头文件(“.h”文件)、库文
我目前正在为类作业编写模拟退火算法(“解决”背包问题),并想在 Rcpp 中完成(我必须使用 R,而 Rcpp 更快)。 Rcpp 一直给我以下错误 invalid static_cast from
根据我的理解,在 Rcpp 和 C++ 之间转换 vector 会创建新 vector ,如下所示。我的理解对吗? 将 Rcpp vector 转换为 C++ vector 时,我们使用 Rcpp::
我想将参数的默认值设置为 NULL在Rcpp如果参数不是NULL,则函数并根据参数进行一些计算.这种代码的一个例子是 #include using namespace Rcpp; // [[Rcpp
任何人都可以解释以下行为吗? 当声明一个新的NumericMatrix时,y,作为原始矩阵,x,乘以一个标量,c,标量/矩阵乘法的顺序很重要。如果我将左侧的标量与右侧的矩阵相乘(例如 NumericM
有一种方法可以使用 NA 值初始化数值向量,例如。 NumericVector x(10,NumericVector::get_na()) 有没有类似的方法可以将矩阵初始化为 NA 值? 最佳答案 这
这可能是一个非常简单的问题,但我不知道哪里出了问题。 我有一个传递给 Rcpp 函数的列表,该列表的第一个元素是一个 data.frame。 我如何获取该 data.frame? bar = list
我正在尝试开发一个使用 Sundials 的 R 包用于求解微分方程的 C 库。为了不让用户安装库,我将库的源代码放在我的包中。 我已将库中的所有头文件放入 /inst/include/sundial
我正在研究一个同时使用 Rcpp::IntegerVector (行/列指针)和模板化 std::vector 的 Rcpp 稀疏矩阵类。基本原理是,在极大的稀疏矩阵中深度复制整数指针 vector
我想将一个R函数翻译成Rcpp,一个简单的测试代码如下,但我不知道如何处理默认设置为NULL的参数。 test t=R_NilValue, Rcpp
我想将一个R函数翻译成Rcpp,一个简单的测试代码如下,但我不知道如何处理默认设置为NULL的参数。 test t=R_NilValue, Rcpp
我想公开一个 C++ 类和一个将该类的对象作为 R 参数的函数。我必须遵循简化的示例。我使用创建了一个包 Rscript -e 'Rcpp::Rcpp.package.skeleton("soq")'
我想用 Rcpp 编写一个 C++ 函数,它使用 hypred 包中的 C 函数,它在 CRAN here 上. 我读了using C function from other package in R
[我在别处将其草拟为评论,但决定创建一个适当的问题...] 在 Rcpp 中使用数据帧时,就代码结构而言,目前被认为是“最佳实践”的是什么?从 R 到 C++ 代码的输入数据帧“传输”非常容易,但是如
我正在尝试使用 Rcpp::CharacterMatrix 并将每一行转换为 Rcpp::List 中它自己的元素。 但是,我为此编写的函数有一个奇怪的行为,即列表的每个条目都对应于矩阵的最后一行。为
我是一名优秀的程序员,十分优秀!