- univariate_eQTL <- function(mask2, x, y, numRegulators, threads = 1L) {
- .Call('GOAL_univariate_eQTL', PACKAGE = 'GOAL', mask2, x, y, numRegulators, threads)
- }
复制代码
这里用的.Call函数,是说明引用了内置函数GOAL_univariate_eQTL.
在src文件里找关于C或者cpp的文件。然后在RcppExports.cpp发现了这个函数:
- // univariate_eQTL
- arma::vec univariate_eQTL(arma::umat mask2, arma::mat x, arma::mat y, arma::vec numRegulators, int threads = 1);
- RcppExport SEXP GOAL_univariate_eQTL(SEXP mask2SEXP, SEXP xSEXP, SEXP ySEXP, SEXP numRegulatorsSEXP, SEXP threadsSEXP) {
- BEGIN_RCPP
- SEXP __sexp_result;
- {
- Rcpp::RNGScope __rngScope;
- Rcpp::traits::input_parameter< arma::umat >::type mask2(mask2SEXP );
- Rcpp::traits::input_parameter< arma::mat >::type x(xSEXP );
- Rcpp::traits::input_parameter< arma::mat >::type y(ySEXP );
- Rcpp::traits::input_parameter< arma::vec >::type numRegulators(numRegulatorsSEXP );
- Rcpp::traits::input_parameter< int >::type threads(threadsSEXP );
- arma::vec __result = univariate_eQTL(mask2, x, y, numRegulators, threads);
- PROTECT(__sexp_result = Rcpp::wrap(__result));
- }
- UNPROTECT(1);
- return __sexp_result;
- END_RCPP
- }
复制代码