r - dose.p からデータを抽出する方法 [R パッケージ MASS] 原子ベクトルの結果

Question

dose.p を使用して、glm モデル (プロビット モデル) から EDx (ED50 など) データを計算しました。

glm.logit <- glm(cbind(nok,nges-nok) ~ log(dosis), family=binomial(lk), data=edx.data)
r <- dose.p(glm.logit,p=seq(0.1,0.9,0.2))

r は次のとおりです。

             Dose         SE
p = 0.1: 0.866650 0.10578072
p = 0.3: 1.301613 0.06405342
p = 0.5: 1.574622 0.05168480
p = 0.7: 1.847632 0.05971840
p = 0.9: 2.282595 0.09898567

exp(r) は正しい EDx を与えます:

> exp(r)
             Dose         SE
p = 0.1: 2.378928 0.10578072
p = 0.3: 3.675219 0.06405342
p = 0.5: 4.828918 0.05168480
p = 0.7: 6.344778 0.05971840
p = 0.9: 9.802082 0.09898567

しかし、EDx だけでなく信頼限界の下限と上限も計算できるように、すべての数値 (用量と SE) を抽出するにはどうすればよいでしょうか? 最後に、そのようなテーブルを計算したいと思います:

  p          Dose        Lower limit       Upper limit
0.1    exp(2.379)   exp(2.379-0.106)  exp(2.379+0.106)
0.2 ...

これらのデータを抽出しようとしましたが、エラー メッセージが表示されました。

> se <- r$SE
    Error in r$SE : $ operator is invalid for atomic vectors

事前にどうもありがとうございました

キリスト教徒