寻求评论良好的mgcv源代码版本

Question

为交叉发布道歉(在R-help上),但是SO会获得更多观点,我希望知道的人可以在这里找到问题.

我正在寻找包含mgcv的各种功能的评论很多的版本,以便我可以为我正在进行的项目修改它的一部分.特别是我正在寻找

• testStat
• summary.gam
• 刘备
• SIMF

我知道我可以通过键入mgcv ::: what来找到它们.但是有很多嵌套if语句,很难遵循.代码中的注释描述了每个步骤的确切内容,这将使我的生活变得更加轻松.

我在哪里可以找到更详细的代码版本？这样的事情存在吗？

编辑:回应删除的评论:我不知道它可以在github上找到,之前没有使用它.FFR,就在这里:https: //github.com/cran/mgcv

Answer 1

您是否在GitHub上搜索过任何源代码？那里有评论代码.我testStat在mgcv.r文件中找到了.

以下是一些帮助我理解复杂代码的建议:

(1)将源写入文本文件,以便将其分隔,移动,等等.

> sink("testStat.txt")
> mgcv:::testStat
> sink()

(2)将其body分成一个部分列表

> x <- as.list(body(mgcv:::testStat))
> x
[[1]]
`{`

[[2]]
qrx <- qr(X, tol = 0)

 ...

[[19]]
if (nu > 0) {
    if (k1 == 1) 
        rank1 <- val <- 1
    else {
        val <- rep(1, k1)
        rp <- nu + 1
        val[k] <- (rp + sqrt(rp * (2 - rp)))/2
        val[k1] <- (rp - val[k])
    }
    if (res.df <= 0) 
        pval <- liu2(d, val)
    else pval <- simf(d, val, res.df)
} else {
    pval <- 2
}

 ...

[[21]]
list(stat = d, pval = min(1, pval), rank = rank)