我刚开始学习R,因此陷入了琐碎的问题.我试图弄清楚如何在R中输出值.在C中我们只使用返回变量,但这似乎不是R的情况.假设我有一个带有4个输入参数的函数,将这些参数传递给ac执行所需计算的函数,现在如果我想在R中加载这个myfun并从c函数获取输出,我需要做什么？以下是我尝试使用的模板.
extern "C" {
SEXP myfun(SEXP S,SEXP A,SEXP B,SEXP C) { //will call this function from R.
SEXP rate,dir,list,list_names; //declare variables
PROTECT( rate = allocMatrix(REALSXP,10,2) ); //allocate 10x2 matrix of double type?
PROTECT( dir = allocVector(INTSXP,10) ); //allocated vector(10) of int type?
double* p_rate = REAL(rate); //why do I need pointers?
int* p_dir = INTEGER(dir);
//这里我称之为C函数,它计算向量< vector< double> > someVal和vector< int> someVal2.
现在我想将这些值传递给rate和dir.
for(int i =0; i < 10; i++){
rate[i][0] = someVal1[i][0];
rate[i][1] = someVal1[i][1];
dir[i] = someVal2[i];
}
但编译器不喜欢它.在这种情况下,我不确定如何使用p_rate和p_dir指针.
另外,我不太明白为什么下面的代码是这样的.它似乎将字符串,rate和dir传递给list_names,并将值传递给rate和dir.为什么我们不能简单地做cout并打印所需的值？
char *tag[2] = {"rate","dir"};
PROTECT(list_names = allocVector(STRSXP,2));
SET_STRING_ELT(list_names,mkChar(tag[0]));
SET_STRING_ELT(list_names,1,mkChar(tag[1]));
PROTECT(list = allocVector(VECSXP,2));
SET_VECTOR_ELT(list,rate);
SET_VECTOR_ELT(list,dir);
setAttrib(list,R_NamesSymbol,list_names);
UNPROTECT(8);
return (list);
}
我没有得到在R中返回变量的基本想法,这看起来很混乱.如果有人能指出适当的资源解释这些东西,我真的很感激.
感谢您的帮助！
