R от C — самый простой Helloworld
Какова самая простая функция C для запуска интерпретатора R, передачи в маленьком выражении (например, 2+2), и уход результат? Я пытаюсь скомпилировать с MingW в Windows.
3 ответа
Вы хотите вызвать R из C?
Посмотрите раздел 8.1 в руководстве Writing R Extensions. Вам также следует заглянуть в каталог "tests" (скачайте исходный пакет, извлеките его, и вы получите каталог tests). Похожий вопрос был задан ранее на R-Help и здесь был пример:
#include <Rinternals.h>
#include <Rembedded.h>
SEXP hello() {
return mkString("Hello, world!\n");
}
int main(int argc, char **argv) {
SEXP x;
Rf_initEmbeddedR(argc, argv);
x = hello();
return x == NULL; /* i.e. 0 on success */
}
Простой пример из руководства R выглядит так:
#include <Rembedded.h>
int main(int ac, char **av)
{
/* do some setup */
Rf_initEmbeddedR(argc, argv);
/* do some more setup */
/* submit some code to R, which is done interactively via
run_Rmainloop();
A possible substitute for a pseudo-console is
R_ReplDLLinit();
while(R_ReplDLLdo1() > 0) {
add user actions here if desired
}
*/
Rf_endEmbeddedR(0);
/* final tidying up after R is shutdown */
return 0;
}
Кстати, вы можете рассмотреть возможность использования Rinside вместо этого: Дирк предоставляет хороший пример "hello world" на домашней странице проекта.
Если вас интересует вызов языка C из R, вот мой первоначальный ответ:
Это не совсем "hello world", но вот несколько хороших ресурсов:
- Джей Эмерсон недавно выступил с докладом о разработке пакетов R в Нью-Йоркской группе useR, и он предоставил несколько очень хороших примеров использования языка C из R. Посмотрите статью из этого обсуждения на его сайте, начиная со страницы 9. Весь соответствующий исходный код находится здесь: http://www.stat.yale.edu/~jay/Rmeetup/MyToolkitWithC/.
- Курс, прочитанный в Гарварде Гопи Госвами в 2005 году: C-C++-R (в статистике). Он включает обширные примеры и исходный код.
Я думаю, вы не можете сделать ничего лучше, чем встроенный пакет (который поддерживает C, C ++ и Fortran):
library(inline)
fun <- cfunction(signature(x="ANY"),
body='printf("Hello, world\\n"); return R_NilValue;')
res <- fun(NULL)
который напечатает вам "Hello, World". И вы даже не знаете, где / как / когда вызываются компилятор и компоновщик. [R_NilValue - это NULL-версия SEXP для R, а используемая здесь подпись .Call () требует, чтобы вы возвращали SEXP - см. Руководство «Написание R-расширений», которого вы не можете здесь избежать. ]
Затем вы возьмете такой код и обернете его в пакет. Мы добились большого успеха при использовании inline для модульных тестов Rcpp (более 200, и сейчас их количество растет) и некоторых примеров.
Да, и этот встроенный пример будет работать в любой ОС. Даже Windoze при условии, что у вас установлена цепочка инструментов для сборки пакетов R, в PATH и т. Д.
Edit: Я неправильно прочитал вопрос. По сути, вам нужно то, что делает внешний интерфейс littler (с использованием чистого C) и то, что классы RInside вынесены за скобки для C ++.
Джефф и я никогда не беспокоились о переносе littler в Windoze, но RInside действительно работал там в самой последней версии.Таким образом, вы должны иметь возможность копаться в рецептах сборки и создать вариант RInside только для C, чтобы вы могли передавать выражение во встроенный процесс R. Я подозреваю, что вам все еще нужно что-то вроде Rcpp в качестве подсказки, поскольку в противном случае это будет утомительно.
Edit 2: И, как упоминает Шейн, действительно есть несколько примеров в источниках R в tests / Embedding / вместе с Makefile.win. Возможно, это самый простой старт, если вы хотите узнать о внутреннем устройстве R.
Вот, пожалуйста. Это основная функция, но вы должны быть в состоянии адаптировать ее к более общим целям. Этот пример строит выражение R из вызовов C, а также из строки C. Для компиляции на windows вы справитесь сами, но я привел шаги компиляции на linux:
/* simple.c */
#include <Rinternals.h>
#include <Rembedded.h>
#include <R_ext/Parse.h>
int
main(int argc, char *argv[])
{
char *localArgs[] = {"R", "--no-save","--silent"};
SEXP e, tmp, ret;
ParseStatus status;
int i;
Rf_initEmbeddedR(3, localArgs);
/* EXAMPLE #1 */
/* Create the R expressions "rnorm(10)" with the R API.*/
PROTECT(e = allocVector(LANGSXP, 2));
tmp = findFun(install("rnorm"), R_GlobalEnv);
SETCAR(e, tmp);
SETCADR(e, ScalarInteger(10));
/* Call it, and store the result in ret */
PROTECT(ret = R_tryEval(e, R_GlobalEnv, NULL));
/* Print out ret */
printf("EXAMPLE #1 Output: ");
for (i=0; i<length(ret); i++){
printf("%f ",REAL(ret)[i]);
}
printf("\n");
UNPROTECT(2);
/* EXAMPLE 2*/
/* Parse and eval the R expression "rnorm(10)" from a string */
PROTECT(tmp = mkString("rnorm(10)"));
PROTECT(e = R_ParseVector(tmp, -1, &status, R_NilValue));
PROTECT(ret = R_tryEval(VECTOR_ELT(e,0), R_GlobalEnv, NULL));
/* And print. */
printf("EXAMPLE #2 Output: ");
for (i=0; i<length(ret); i++){
printf("%f ",REAL(ret)[i]);
}
printf("\n");
UNPROTECT(3);
Rf_endEmbeddedR(0);
return(0);
}
Шаги компиляции:
$ gcc -I/usr/share/R/include/ -c -ggdb simple.c
$ gcc -o simple simple.o -L/usr/lib/R/lib -lR
$ LD_LIBRARY_PATH=/usr/lib/R/lib R_HOME=/usr/lib/R ./simple
EXAMPLE #1 Output: 0.164351 -0.052308 -1.102335 -0.924609 -0.649887 0.605908 0.130604 0.243198 -2.489826 1.353731
EXAMPLE #2 Output: -1.532387 -1.126142 -0.330926 0.672688 -1.150783 -0.848974 1.617413 -0.086969 -1.334659 -0.313699