1

この方法で、olsオブジェクトから勾配と切片のp値を抽出できます。

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

しかし、robcovオブジェクトで同じことを試してみると、summary.lmは、robcovモデルではなく、元のモデル（m1）からのp値を提供します。

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

これは、robcovヘルプページの警告に関連している必要があると思います。

警告

調整されたolsフィットには、print.olsで印刷された修正された標準誤差がありません。これを取得するには、sqrt（diag（adjfit $ var））を使用します。ここで、adjfitはrobcovの結果です。

でもどうすればいいのかわかりません。

robcovオブジェクトからp値を抽出する方法はありますか？（私は本当にスロープ用のものだけに興味があります、それが違いを生むなら...）

4

1 に答える 1

0

ハッキングし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)

于 2012-02-15T20:29:34.630 に答える