Home / Questions / 106156
Accepted
Alejandro Carrera
Asked: 2020-09-30 14:32:33 +0800 CST 2020-09-30 14:32:33 +0800 CST

如何获取R中函数的源代码？（车包）

  • 772

R 中的“car”包有一个有趣的函数，叫做 outlierTest，我想分析它的源代码以了解它是如何工作的。我尝试了以下操作并收到错误：

getMethod("outlierTest")

    Error in getMethod("outlierTest") : 
     no generic function found for 'outlierTest'

这似乎让我得到了回应，但它将我重定向到一个没有显示源代码的页面：

methods(outlierTest)
    [1] outlierTest.lm*
    see '?methods' for accessing help and source code

最后，在控制台中仅输入“outlierTest”时，我没有得到任何响应：

outlierTest
    function (model, ...) 
    {
          UseMethod("outlierTest")
    }
<bytecode: 0x000000000cf4d3d0>
<environment: namespace:car>

我很感激这方面的任何指导。

r

2 Answers

  1. Best Answer

    您可以执行以下操作，至少使用此软件包可以正常工作：

    我们寻找方法：

    > methods("outlierTest")
[1] outlierTest.lm*
see '?methods' for accessing help and source code

    特别是根据文档outlierTest.lm具有的方法：*

    如果方法定义不是从定义方法的包命名空间导出的，则 S3 方法名称后跟星号 *

    要在这种情况下检索您的代码，我们可以尝试以下操作：

    getAnywhere('outlierTest.lm')

    如果没有，另一种方法是直接访问项目的 github，并且准时，您可以在此链接中检索此功能

    • 1

  2. 另一种方法是您可以使用三个“：”访问“未导出”函数（在您的命名空间中不可见）

    请记住，正如您所看到的，该函数outlierTest是通用的，并且要访问它，特别是它适用于类的对象，lm您必须添加一个点functiongenerica.clase：：

    car:::outlierTest.lm


function (model, cutoff = 0.05, n.max = 10, order = TRUE, labels = names(rstudent), 
    ...) 
{
    rstudent <- rstudent(model)
    if (length(rstudent) != length(labels)) 
        stop("Number of labels does not correspond to number of residuals.")
    else names(rstudent) <- labels
    df <- df.residual(model) - 1
    rstudent <- rstudent[!is.na(rstudent)]
    n <- length(rstudent)
    p <- if (class(model)[1] == "glm") 
        2 * (pnorm(abs(rstudent), lower.tail = FALSE))
    else 2 * (pt(abs(rstudent), df, lower.tail = FALSE))
    bp <- n * p
    ord <- if (order) 
        order(bp)
    else 1:n
    ord <- ord[bp[ord] <= cutoff]
    result <- if (length(ord) == 0) {
        which <- which.max(abs(rstudent))
        list(rstudent = rstudent[which], p = p[which], bonf.p = bp[which], 
            signif = FALSE, cutoff = cutoff)
    }
    else {
        if (length(ord) > n.max) 
            ord <- ord[1:n.max]
        result <- list(rstudent = rstudent[ord], p = p[ord], 
            bonf.p = bp[ord], signif = TRUE, cutoff = cutoff)
    }
    class(result) <- "outlierTest"
    result
}
<bytecode: 0x000000001e2258a0>
<environment: namespace:car>
    • 1