r - 在 R 中评估 C 代码

Question

我想查看 R 中族背后的代码，例如:

make.link("logit")
make.link("identity")
make.link("probit")

现在在 R-repsones 中看到调用了用编程语言 C 编写的函数:

例子:

.Call(C_logit_link, mu)
.Call(C_logit_linkinv, eta)
.Call(C_logit_mu_eta, eta)

现在想访问此特定代码。有什么办法可以做到吗？

Answer 1

调用.Call的函数是编译代码中的调用入口点，所以如果你想完全理解这个函数，你必须查看编译代码的源代码

例如，您要查找的是 here

SEXP logit_link(SEXP mu)
{
    int i, n = LENGTH(mu);
    SEXP ans = PROTECT(duplicate(mu));
    double *rans = REAL(ans), *rmu=REAL(mu);

    if (!n || !isReal(mu))
    error(_("Argument %s must be a nonempty numeric vector"), "mu");
    for (i = 0; i < n; i++)
    rans[i] = log(x_d_omx(rmu[i]));
    UNPROTECT(1);
    return ans;
}

我建议您阅读这个优秀的答案 here关于如何获取任何 R 函数源代码。