如何从dose.p [R package 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

但是如何提取所有数字（剂量和SE）不仅可以计算EDx，还可以计算信心下限和上限？最后我想计算一下这样一张表：

  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

非常感谢您提前

基督教