使用Rcpp中的rWishart生成Wishart分布？

Question

我正在尝试让R的rWishart的统计数据在Rcpp中工作。 这是一个测试示例：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector rWishart_cpp(int p, int df, NumericMatrix Sigma) {
    return rWishart(p, df, Sigma);
}

我收到的错误是：“ test.cpp：12：12：错误：使用未声明的标识符'rWishart'返回rWishart（p，df，Sigma）;”

从其中一个小插曲来看，我认为d / r函数是由Rmath.h提供的Rcpp实现的？

无论如何，有什么方法可以在Rcpp中生成wishart分布吗？

Answer 1

这里的问题似乎是您假设Wishart Rmath.h ：

edd@max:~$ grep -i wish /usr/include/Rmath.h         # no Wishart
edd@max:~$ grep -i weib /usr/include/Rmath.h         # but Weibull has an example
#define dweibull        Rf_dweibull
#define pweibull        Rf_pweibull
#define qweibull        Rf_qweibull
#define rweibull        Rf_rweibull
        /* Weibull Distribution */
double  dweibull(double, double, double, int);
double  pweibull(double, double, double, int, int);
double  qweibull(double, double, double, int, int);
double  rweibull(double, double);
edd@max:~$ 

最好的选择是从某个地方获取另一个实现，然后通过Rcpp进行编码。

编辑：或者根据在rdocumentation.org上两分钟的研究，您可能可以遵循stats::rWishart的文档，并将其编码为多元法线的函数。 我将基于RcppArmadillo ....