Rcpp：如何将R函数和Rcpp函数组合在一起做一个包

Question

假设我在名为 test.cpp

的文件中有以下 C++ 代码

#include <Rcpp.h>

//[[Rcpp::export]]
Rcpp::NumericMatrix MyAbar (const Rcpp::NumericMatrix & x, int T){
    unsigned int outrows = x.nrow(), i = 0, j = 0;
    double d;
    Rcpp::NumericMatrix out(outrows,outrows);
    // Rcpp::LogicalVector comp;

    for (i = 0; i < outrows - 1; i++){
        Rcpp::NumericVector v1 = x.row(i);
        Rcpp::NumericVector ans(outrows);
        for (j = i + 1; j < outrows ; j ++){
            d = mean(Rcpp::runif( T ) < x(i,j));
            out(j,i)=d;
            out(i,j)=d;
        }
    }
    return out;
}

我知道用下面的命令，我可以拥有自己的包

Rcpp.package.skeleton("test",cpp_files = "~/Desktop/test.cpp")

但是，如果我想将调用 Rcpp 函数的以下 R 函数组合到包中怎么办

random = function(A, T){
    if (!is.matrix(A)){
        A = Reduce("+",A)/T
    } 
    # global constant and threshold
    n = nrow(A)
    B_0 = 3
    w = min(sqrt(n),sqrt(T * log(n)))
    q = B_0 * log(n) / (sqrt(n) * w)

    A2 = MyAbar(A) 
    diag(A2) <- NA 
    K = A2 <= rowQuantiles(A2, probs=q, na.rm =TRUE)
    diag(K) = FALSE 
    P = K %*% A * ( 1/(rowSums(K) + 1e-10))

    return( (P + t(P))*0.5 )
}

我怎样才能做到？

Answer 1

所以你问的是如何制作R包？有很多好的教程。

第一个近似值：

将您的文件复制到文件 R/random.R
处理函数的帮助文件，手动编写 man/random.Rd 或学习包 roxygen2
确保您知道 NAMESPACE 的用途并且 DESCRIPTION 是正确的

Rcpp：如何将R函数和Rcpp函数组合在一起做一个包

Rcpp: how to combine the R function and Rcpp function together to make a package

r

rcpp

rcpp11

r-package