Question

包名为pareto。这是.c目录中的src文件：

#include <R.h>
#include <math.h>
#include "Rinternals.h"
#include "R_ext/Rdynload.h"

static void dpareto(double *x, double *a, double *b, int *n, int *islog,
                           double *den){
    int length = n[0];
    int i;
    int isLog = islog[0];
    for (i = 0; i < length; i++){
        if (a[i] > 0 && b[i] > 0) {
            if (x[i] > a[i])
                den[i] = log(b[i]) + b[i] * log(a[i]) - (b[i] + 1) * log(x[i]);
            else
                den[i] =  R_NegInf;
            if (!isLog)
                den[i] = exp(den[i]);
        }
        else {
            den[i] = R_NaN;
        }
    }
}

static R_CMethodDef DotCEntries[] = {
    {"dpareto", (DL_FUNC) dpareto, 6},
    {NULL}
};

void R_init_pareto(DllInfo *info)
{
    R_registerRoutines(info, DotCEntries, NULL, NULL, NULL);
}

在R目录中，相应的.R文件为：

#' @useDynLib pareto
#' @export
dpareto <- function(x, a, b, log = FALSE) {
    nx <- length(x)
    na <- length(a)
    nb <- length(b)
    n <- max(nx, na, nb)
    if (nx < n) x <- rep(x, length.out = n)
    if (na < n) a <- rep(a, length.out = n)
    if (nb < n) b <- rep(b, length.out = n)
    rt <- .C("dpareto", as.double(x), as.double(a), as.double(b), as.integer(n),
             as.integer(log), den = double(n), PACKAGE="pareto")
    rt$den

}

roxygen记录后，NAMESPACE有：

export(dpareto)
useDynLib(pareto)

但是包无法通过检查，R会生成错误消息：

 "dpareto" not available for .C() for package "pareto"
Calls: dpareto -> .C

我真的无法弄清楚我犯了哪一步。

Answer 1

您已将static关键字添加到dpareto函数定义中。这意味着该函数不会被导出，所以R不会看到它。删除static然后重试。

Answer 2

琐碎的错误。就在void R_init_pareto我使用了错误的包名。这么傻的错误。

如何在R包中调用C函数

2 个答案: