Question

我正在尝试运行下面的函数weights。函数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;
}

Answer 1

mN的功能签名是

arma::vec mN(double, double, arma::vec, arma::vec)

但你称之为

mu     = mN(kappa, kappaN, nj[i], m, C[i]);

其中kappa为int，kappaN为int，nj[i]为double，m为int }和C是List的{{1}}。

请注意，与R不同，您通常不能将C ++中的双精度和整数混合使用。特别是，我还要指出像

这样的声明

SEXP

应该都是double r = indic.n_elem; double p = y.n_cols; double n = y.n_rows;或int s。

Answer 2

定义中的函数签名：

arma::vec mN(double kappa, double ni, arma::vec m, arma::vec C)

与其用法不符：

mu     = mN(kappa, kappaN, nj[i], m, C[i]);

您必须决定该函数是否需要4个或5个参数。

Rcpp在文件中有多个函数，没有匹配函数

2 个答案: