在 R 中将 coef 和 summary.lm 与 robcov 一起使用（提取 p 值）

Question

我可以提取斜率和斜率的 p 值以这种方式从 ols 对象拦截：

library(rms)
m1 <- ols(wt ~ cyl, data= mtcars, x= TRUE, y= TRUE)
coef(summary.lm(m1))

但是当我尝试使用 robcov 对象执行相同的操作时，summary.lm 为我提供原始模型 (m1) 的 p 值，而不是 robcov 模型：

m2 <- robcov(m1)
m2
coef(summary.lm(m2))

我认为这必须与来自 robcov 帮助页面的警告，

警告

调整后的 ols 拟合没有打印更正的标准误差 与 print.ols。使用 sqrt(diag(adjfit$var)) 来获取此值，其中 adjfit 是 robcov 的结果。

但我不知道怎么做。

有没有办法从 robcov 对象中提取 p 值？ （我真的只对斜坡感兴趣，如果这有影响的话......）

Answer 1

通过 print.ols 和 prModFit 进行黑客攻击，我想出了这个。

errordf <- m2$df.residual
beta <- m2$coefficients
se <- sqrt(diag(m2$var))
Z <- beta/se
P <- 2 * (1 - pt(abs(Z), errordf))

将 m2 更改为另一个 robcov 模型。

通过比较 P 与 print(m2) 的结果来亲自尝试一下