drc 封装中的反向 ED95 和 ED5

Question

我的目标是确定治疗注射部位与治疗目标之间的距离，成功概率为 0.95。

结果变量是一个二进制变量（成功：1/失败：0）

我使用 Dixon 上下方法测试了六个距离：0、2、4、6、8 和 10 毫米。这是我的数据：

第 1 列：使用的距离

第 2 列：成功次数

第 3 列：患者总数

data <- data.frame(1:6,1:6,1:6)
data[,1] <- c(0, 2, 4, 6, 8, 10)
data[,2] <- c(2, 12, 3, 2, 1, 0)
data[,3] <- c(2, 12, 15, 8, 4, 1)
names(data) <- c("Distance", "Success", "Total")

我在 Windows Vista Os 上使用 DRC 包 2.3-96 和 R ver 3.1.2 构建了一个模型：

library(drc)
model <- drm(Success/Total~Distance, weights=Total, 
             data=data, fct=LL.2(), type="binomial")

summary(model)
plot(model, bp=.5, legend=FALSE
     , xlab=paste("Distance"), ylab="Probability of success", lwd=2,
     cex=1.2, cex.axis=1.2, cex.lab=1.2, log = "")

似乎一切正常

但是在估计 ED 95（有效剂量 95：成功概率为 0.95 所需的距离）时，我认为该 ED95 与 ED5（有效剂量 5：成功概率为 0.05 所需的距离）相反) :

ED(model, 95, interval="delta")
ED(model, 5, interval="delta")

ED95: 8.0780 SE: 2.0723 CI 95 % (4.0165 ; 12.139)

ED5 : 1.58440 SE: 0.46413 CI 95 % (0.67472 ; 2.4941)

Answer 1

drc 包中的 ED 值默认是相对于控制水平计算的。在我们的例子中，我们正在寻找相对于上限计算的 ED 值。

所以我们必须将 reference 值从 "control"（默认值）更改为 "upper" :

ED(model, 95, interval="delta", reference = "upper")

非常感谢 Christian Ritz

drc 封装中的反向 ED95 和 ED5

Inverted ED95 and ED5 in drc package

model

r