Rcpp:Rcpp 中是否有实现 fisher.test() [英] Rcpp: Is there an implementation fisher.test() in Rcpp
问题描述
是否有 fisher 的实现.test 用于 Rcpp?
Is there an implementation of fisher.test for Rcpp?
推荐答案
目前在 Rcpp 中没有 fisher.test
函数的实现.特定的 组件由于计算量大,R 函数中的 test 是用 C 编写的.非常欢迎您在 Rcpp 中重新实现该测试.
There is no current implementation of the fisher.test
function within Rcpp. A specific component of the test in the R function is written in C due to the computational intensity of it. You are more than welcome to reimplement the test in Rcpp.
尽管有一些注意事项,SEXP
对象中没有 factor
的表示.因此,factor
不是 Rcpp 支持的东西.因此,在将对象传递到 C++ 之前,您必须将对象转换为 integer
或 character
类型.
A few notes on that though, there is no representation of a factor
in a SEXP
object. Thus, factor
not something Rcpp that supports. As a result, you would have to convert the object into an integer
or character
type before passing it into C++.
为了将它用于 Rcpp,您必须从 C++ 调用 fisher.test
R 函数.也就是说,您必须将数据传输回 R.
In order to use it for Rcpp, you would have to call the fisher.test
R function from C++. That is, you would have to transfer the data back into R.
例如
#include <Rcpp.h>
//' @title Accessing R's fisher.test function from Rcpp
// [[Rcpp::export]]
Rcpp::List fisher_test_cpp(const Rcpp::NumericMatrix& x, double conf_level = 0.95){
// Obtain environment containing function
Rcpp::Environment base("package:stats");
// Make function callable from C++
Rcpp::Function fisher_test = base["fisher.test"];
// Call the function and receive its list output
Rcpp::List test_out = fisher_test(Rcpp::_["x"] = x,
Rcpp::_["conf.level"] = conf_level);
// Return test object in list structure
return test_out;
}
/***R
Job = matrix(c(1,2,1,0, 3,3,6,1, 10,10,14,9, 6,7,12,11), 4, 4,
dimnames = list(income = c("< 15k", "15-25k", "25-40k", "> 40k"),
satisfaction = c("VeryD", "LittleD", "ModerateS", "VeryS")))
fisher.test(Job)
fisher_test_cpp(Job)
*/
注意 cpp 函数以列表形式返回对象:
Note the cpp function returns objects in list form under:
List of 7
$ p.value : num 0.783
$ alternative: chr "two.sided"
$ method : chr "Fisher's Exact Test for Count Data"
$ data.name1 : chr "structure(c(1, 2, 1, 0, 3, 3, 6, 1, 10, 10, 14, 9, 6, 7, 12, "
$ data.name2 : chr "11), .Dim = c(4L, 4L), .Dimnames = structure(list(income = c(\"< 15k\", "
$ data.name3 : chr "\"15-25k\", \"25-40k\", \"> 40k\"), satisfaction = c(\"VeryD\", \"LittleD\", "
$ data.name4 : chr "\"ModerateS\", \"VeryS\")), .Names = c(\"income\", \"satisfaction\")))"
- attr(*, "class")= chr "htest"
可以使用以下方法访问这些值:
These values can be accessed using:
double p_value = test_out[0];
std::string alternative = test_out[1];
std::string method = test_out[2];
等等……
这篇关于Rcpp:Rcpp 中是否有实现 fisher.test()的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持IT屋!