Qgamma 함수는 Rcpp를 통해 호출 될 때 0을 반환합니다.

Question

RedHat에서 R 2.15.2 및 Rcpp 10.4 (며칠 전에 업그레이드 됨)를 실행 중입니다.

make를 통해 컴파일 된 내 C++ 프로그램에서 qgamma 함수를 호출하면 0이 반환되지만 동일한 프로그램의 qnorm은 올바른 값 (뒤따라 오는 소스)을 반환합니다. cppFunction 컴파일을 통해 호출하면 올바른 대답.

C++ 소스 :

#include <iostream> 
using namespace std; 

#include <Rcpp.h> 

int main() { 
    cout << R::qnorm(0.3, 1.0, 1.0, 1, 0) << endl; 
    cout << R::qgamma(0.3, 5.0, 5.0, 1, 0) << endl; 
    return 0; 
} 

와 연관된 컴파일 메시지 및 실행 : 되 물론

[jbowma1@smartrepl-app00 src]$ make test 
g++ -I/usr/include/R -I/usr/lib64/R/library/Rcpp/include -L/usr/lib64/R/lib -lR -L/usr/lib64/R/library/Rcpp/lib -lRcpp -Wl,-rpath,/usr/lib64/R/library/Rcpp/lib -o test test.cpp 
[jbowma1@smartrepl-app00 src]$ ./test 
0.475599 
0 
[jbowma1@smartrepl-app00 src]$ 

문제 "0"것을 상기 출력의 두 번째 행이다. cppFunction 버전 이제

: 단지 R 프롬프트에서 qgamma를 호출하는 경우

가

> foosrc <- "double foo() {return R::qgamma(0.3, 5.0, 5.0, 1, 0);}" 
> bar <- cppFunction(foosrc) 
> bar() 
[1] 18.16805 

(. 같은 컴퓨터가, 자연이 하나가 RStudio에서 실행)

주, 세 번째 매개 변수는 필요가 같은 대답을 얻기 위해 거꾸로 (이 예에서는 0.2).

qnorm이 실제로 작동하지 않는다면 RcppExamples 디렉토리에서 실제로 복사되는 makefile을 엉망으로 만들 수 있습니다. dgamma도 작동하지만 pgamma은 그렇지 않습니다. 다른 분포, 예를 들어, 음 이항 또한 작용한다.

의견을 보내 주시면 대단히 감사하겠습니다.

Answer 1

첫 번째 예는 "우스운 냄새"입니다. main()과 Rcpp.h으로 빌드 할 수 없습니다. Rcpp.h은 개의 애드온을 의미합니다. 여기에는 메인을 제공하는 R에로드됩니다.

지금, 당신은 당신이, R을 포함 할 것이다

1. RInside을 사용할 수 있으며 Rcpp.h를 사용하여 얻을
2. Rcpp
의 나머지 부분 또는 당신은 또한 R에서 제공하는 독립 수학 라이브러리 Rmath을 사용할 수 있습니다 (내 패키지를 통해 Debian/Ubuntu의 경우 r-mathlib 패키지로 제공됩니다).

여기에서 올바른 컨텍스트를 사용하고 있기 때문에 두 번째 예제가 작동합니다. 첫 번째는 그렇지 않습니다. 그 의미에서 당신의 비교는 꺼져 있습니다.

편집 다음은 qbeta 않는 내가 주위를 어슬렁했다 샘플 프로그램입니다. 여기에 어떤 Rcpp, 쓰기 R 확장에 설명 된대로 R의의 MathLib 단지 외부 사용 :

는

// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; 
//   compile-command: "gcc -s -Wall -O3 -I/usr/share/R/include 
//        -o rmath_qbeta rmath_qbeta.c -lRmath -lm" -*- 

#include <stdio.h> 

#define MATHLIB_STANDALONE 1 
#include <Rmath.h> 

int main(void) { 

    double x = 0.25, a = 0.8, b = 2.0; 

    printf("qbeta: %f %f %f %f\n", 
      qbeta(x, a, b, 0, 0), 
      qbeta(log(x), a, b, 0, 1), 
      qbeta(x, a, b, 1, 0), 
      qbeta(log(x), a, b, 1, 1)); 

    return 0; 
} 

은/링크 경로가 데비안/우분투 작업을 포함하지 않습니다; 주석 달린 첫 번째 세 줄은 한 줄에서 들여 쓰여졌습니다 (이맥스 사용을위한 것입니다).