使用COEF和summary.lm與R中robcov（提取p值）

Question

我可以提取p值用於從醇我的斜率&截距對象這樣：

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)