我正在使用 drc 包来计算我使用的抑制剂的 ED50。我得到了这个值:
concentration:
[1] 1.0e-04 1.0e-04 1.0e-04 1.0e-01 1.0e-01 1.0e-01 2.5e-01 2.5e-01 2.5e-01 5.0e-01 5.0e-01
[12] 5.0e-01 1.0e+03
str(concentration)
num [1:13] 0.0001 0.0001 0.0001 0.1 0.1 0.1 0.25 0.25 0.25 0.5 ...
viability:
[1] 2576959.0000 3241331.0000 3095945.0000 2476539.0000 2667830.0000 2768508.0000 1175715.0000
[8] 1511373.0000 1603243.0000 1025978.0000 1200592.0000 1202150.0000 0.0001
str(viability)
num [1:13] 2576959 3241331 3095945 2476539 2667830 ...
model_logistic<-drm(viability~concentration,fct = LL.4(names = c("Slope","Lower Limit","Upper Limit","IC50")))
ED(model_logistic, c(50), interval = "delta")
Estimate Std. Error Lower Upper:
e:1:50 0.28926 0.066065 0.139813 0.438714
我想检查 Estimation +/- 1.96 SE 是否是 CI 的上限/下限。但这不起作用,因为它宁愿更小。我检查了系数,计算出的 CI 为 2.26 SE。这里发生了什么?为什么我不能继续使用我知道的公式来计算 CI?