Question

我正在使用caper包运行系统发育分析，其中回归函数（使用系统发育独立对比）是crunch。 crunch函数使用名为caper的{{1}}包内部的对象。

该模型通过以下方式启动：

caic

当我运行crunchMod <- crunch(y ~ f(x), data = comparison)时，我的格式与summary(crunchMod)函数生成的摘要相同。

但是，在尝试通过输入lm()开始检查模型假设时，我收到以下错误：

rstandard(crunchMod)

阅读第19-20页的http://cran.r-project.org/web/packages/caper/vignettes/caper.pdf，我发现Error in UseMethod("rstandard") : no applicable method for 'rstandard' applied to an object of class "caic"使用包装器来检查回归假设。但是，这些是图形检查：

plot(crunchMod)

有没有人知道如何使用我自己的包装器访问标准化残差，还是允许我访问p值而不是图形图像？

Answer 1

这个非常简单，但只有在if() crunch().语句中查看caic的第二个> summary.caic function (object, ...) { summary(object$mod, ...) } <environment: namespace:caper>语句后，在圆圈中转了一圈，它只是一个子集整个摘要/模型

names(summary(crunchMod))
# [1] "call"          "terms"         "residuals"     "coefficients" 
# [5] "aliased"       "sigma"         "df"            "r.squared"    
# [9] "adj.r.squared" "fstatistic"    "cov.unscaled"

您可以看到整个模型的名称显示其余部分中有用的统计信息。

mod

请注意，只有lm()部分继承自crunchMod您可以使用

查看

> sapply(crunchMod, is)
$contrast.data
[1] "list"   "vector"

$mod
[1] "lm"       "oldClass"

$data
[1] "comparative.data"

个对象的所有继承性

{{1}}

R：基本函数不能使用包Caper中的对象

1 个答案: