r - rcpp 文件中有多个函数且没有匹配的函数
问题描述
我试图weights
在 R 中运行下面的函数。函数mN
和PsiN
下面的函数单独工作,我不需要将它们导出到 R 中(即它们的唯一目的是保持函数权重看起来更整洁)。
出于某种原因,仅mN
在函数内给我错误“没有匹配的函数调用” weights
。知道为什么吗?
#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
using namespace Rcpp;
using namespace arma;
arma::vec mN(double kappa, double ni, arma::vec m, arma::vec C) {
arma::vec mN;
double kappaN;
kappaN = kappa + ni;
mN = (kappa/kappaN) * m + (ni/kappaN) * C;
return mN;
}
arma::mat PsiN(double ni, double kappa, double kappaN,
arma::mat y, arma::vec indic, arma::vec C, arma::vec m,
arma::mat Psi) {
double r = indic.n_elem;
double p = y.n_cols;
double n = y.n_rows;
arma::mat S(p, p);
arma::mat PN(p, p);
if (r == 1) {
PN = Psi + kappa*ni/kappaN * ( (C - m ) * (C - m).t() );
} else {
for (int i = 0; i < n; i++) {
S += y.row(i).t() * y.row(i);
}
PN = Psi + S + kappa*ni/kappaN * ( (C - m ) * (C - m).t() ) ;
}
return PN;
}
// [[Rcpp::export]]
arma::vec weights(int alpha, int v, int kappa, int m,
arma::vec nj, arma::mat x, arma::mat Psi,
List C, List indic) {
int p = x.n_cols;
int n = x.n_rows;
int kappaN;
int vN;
int crp;
arma::vec Q;
arma::vec mu;
arma::mat Sigma;
for (int i = 0; i < n; i++) {
kappaN = kappa + nj[i];
vN = v + nj[i] - p + 1;
Sigma = PsiN(nj[i], kappa, kappaN, x, indic[i], C[i], m, Psi);
Sigma = Sigma*(kappaN + 1)/(kappaN*(vN - p + 1));
mu = mN(kappa, kappaN, nj[i], m, C[i]);
crp = log(nj[i]) - log(n + alpha - 1);
Q[i] = (crp - lgamma((vN + p)/2) - lgamma(vN/2) - p/2*(log(vN) + log(
datum::pi)) - 1/2*log_det(Sigma) - (vN + p)/2 * (log(1 + 1/v*(x - mu).t()*Sigma.i()*(x-mu)))) ;
}
return Q;
}
解决方案
您的函数签名mN
是
arma::vec mN(double, double, arma::vec, arma::vec)
但你称它为
mu = mN(kappa, kappaN, nj[i], m, C[i]);
其中kappa
is int
, kappaN
is int
, nj[i]
is double
, m
isint
并且C
是sList
的a SEXP
。
请注意,与 R 不同,您通常不能在 C++ 中随意混合双精度和整数。特别是,我还要指出像这样的声明
double r = indic.n_elem;
double p = y.n_cols;
double n = y.n_rows;
应该都是int
s或unsigned int
s。
推荐阅读
- mysql - 如何搜索关键字在两个表格之间的多个文本中出现的次数
- virtualbox - VirtualBox 一直要求我重新安装 Red Hat
- ios - 从循环播放列表中播放视频播放列表和单个视频剪辑,AVQueuePlayer iOS Swift
- java - 如何获取“主机”对象java
- json - 如何使用 JOLT 将整个 JSON 对象作为值映射到新的 JSON 键中?
- java - 如何从包含excel xls文件中的公式的单元格中获取单元格值
- amazon-web-services - Creating TWO way trusts between two AWS managed Active Directory
- android - 使用 gradle 模块作为具有特定风格的依赖项
- php - 图像移动到文件夹但不保存数据库中的行
- css - 在本机反应中将图像放在视图上