如何在R包中用R与testthat对Rcpp函数做单元测试

核心问题拆解 #

你碰到的问题根源有两个：一是Rcpp函数没正确导出到包的命名空间，导致测试时找不到函数；二是测试环境未正确加载包，无法调用包内的Rcpp函数。以下是纯R侧的解决方案，无需用C++端测试工具：

步骤1：确保Rcpp函数正确导出 #

1. 在你的C++源文件（比如src/timesTwo.cpp）里，给函数加上// [[Rcpp::export]]标签：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector timesTwo(NumericVector x) {
  return x * 2;
}

2. 切换到包的根目录，在R控制台运行compileAttributes()，这会自动生成R/RcppExports.R文件——里面包含将Rcpp函数导出到包命名空间的代码，从此R侧就能直接调用timesTwo()了。

步骤2：配置testthat测试框架 #

1. 修改包的DESCRIPTION文件，添加testthat依赖：

Suggests: testthat (>= 3.0.0)
Config/testthat/edition: 3

2. 在包根目录创建tests/testthat.R文件，内容如下：

library(testthat)
library(你的包名) # 替换成你的包的实际名称

test_check("你的包名")

3. 创建tests/testthat文件夹，用来存放测试脚本。

步骤3：编写对比测试脚本 #

在tests/testthat/test-timesTwo.R中编写测试代码，用R实现做参考，对比Rcpp函数的输出：

# 定义R语言参考实现
timesTwo_r <- function(x) {
  x * 2
}

test_that("timesTwo函数输出与R参考实现一致", {
  # 生成测试输入
  test_input <- c(1, 2, 3, 4.5, -2)
  
  # 分别调用两个函数
  r_result <- timesTwo_r(test_input)
  rcpp_result <- timesTwo(test_input)
  
  # 对比结果（允许微小浮点误差）
  expect_equal(rcpp_result, r_result, tolerance = 1e-6)
  
  # 额外测试边界情况
  expect_equal(timesTwo(numeric(0)), timesTwo_r(numeric(0)))
  expect_equal(timesTwo(NA_real_), timesTwo_r(NA_real_))
})

步骤4：运行测试验证 #

在R控制台执行devtools::test()，或者在终端运行R CMD check，此时测试脚本会正确加载包并调用timesTwo()函数，完成结果对比验证。

内容的提问来源于stack exchange，提问作者titfortat