没有用于调用sqrt的匹配函数

Question

在Rcpp中编写逐列计算样本偏斜度的函数我在使用sqrt((-函数时遇到了麻烦。我知道sqrt(x(适用于NumericVector类型(在一个单独的文件中测试(，但在我的代码中(我试图传递doubles(，它不起作用。

这是我的代码：

#include <Rcpp.h>
#include <cmath>
using namespace Rcpp;
// [[Rcpp::export]]
NumericVector colSkew(NumericMatrix x) {
int nc = x.ncol();
int nr = x.nrow();
NumericVector colS(nc);
for(int i = 0; i < nc; i++){
double cMean = mean(x( _ ,i));
double xSq = 0;
double cSt = 0;
for(int j = 0; j < nr; j++){
xSq += std::pow(x(j,i), 2.0);
cSt += std::pow(x(j,i) - cMean, 3.0);
}
double colMsq = nr * std::pow(cMean, 2.0);
double cTT = sqrt((xSq - colMsq)) / (nr - 1);
double colVar = cTT / (nr - 1);
double colNew = nr * std::pow(colVar, 3);
colS[i] = cSt / colNew;
}
return(colS);
}

我尝试了std::sqrt()；对sqrt的调用不明确"；错误对于那个，Getting：error C2668:'；sqrt'；：对重载函数的调用不明确，也http://dirk.eddelbuettel.com/code/rcpp.examples.html的

inline static double sqrt_double( double x ){ return ::sqrt( x ); }

帮助。(后者本身存在"无匹配函数"或"模糊"问题(。代码本身可以编译，但不能正常工作(我可以调用函数，但结果不好(。

可复制示例：

# First load the cpp function, however you want to;
Create example data:
A = matrix(rchisq(1000, 5), nrow = 100)
library(timeDate)
skew.1 = apply(A, 2, skewness)
skew.2 = colSkew(A)
# A custom R-function which should mimic the cpp function
colSkew.r = function(x){

nc = ncol(x)
nr = nrow(x)
colS = numeric(nc)
cMean = colMeans(x)
xSq = colSums(x^2)
cSt = 0
for(i in 1:nc){
cSt = sum((x[,i]-cMean[i])^3)
colMsq = nr * cMean[i]^2
cTT = sqrt((xSq[i] - colMsq) / (nr - 1))
colNew = nr * cTT^3
colS[i] = cSt / colNew
}
return(colS)
}
skew.3 = colSkew.r(A)

Answer 1

简单地说，我们可以避免一些反复出现的常见初学者错误：

• 不包括不需要的头(大部分是无害的，但Rcpp已经带来了数学头(

• 不要使用全局扁平化的命名空间Rcpp，尤其是在符号可见性上出现编译错误之后

有了这些以及所需的最小更改，它就构建并运行了(并返回了一些毫无意义的东西，但我将其留给您处理：(

还有其他方法可以解决这一问题，因为std::符号和Rcpp符号通常与不同的签名非常愉快地共存。

代码

#include <Rcpp.h>
// [[Rcpp::export]]
Rcpp::NumericVector colSkew(Rcpp::NumericMatrix x) {
int nc = x.ncol();
int nr = x.nrow();
Rcpp::NumericVector colS(nc);
for(int i = 0; i < nc; i++){
double cMean = Rcpp::mean(x( Rcpp::_ ,i));
double xSq = 0;
double cSt = 0;
for(int j = 0; j < nr; j++){
xSq += std::pow(x(j,i), 2.0);
cSt += std::pow(x(j,i) - cMean, 3.0);
}
double colMsq = nr * std::pow(cMean, 2.0);
double cTT = std::sqrt((xSq - colMsq)) / (nr - 1);
double colVar = cTT / (nr - 1);
double colNew = nr * std::pow(colVar, 3);
colS[i] = cSt / colNew;
}
return(colS);
}
/*** R
set.seed(42)
colSkew( matrix(rnorm(100), 100, 1) )
*/

输出

> sourceCpp("answer.cpp")
> set.seed(42)
> colSkew( matrix(rnorm(100), 100, 1) )
[1] -448250877
>