将drm模型拟合到我的剂量反应数据时,获得的ED50
关闭了
log(ED50)
,log2(ED50)
,log10(ED50)
# here is some actual data
test_response <- c(0,1.130912987,-2.0159852,6.158574058,5.679161469,12.97,20.46711451,35.55009039,66.31837697)
test_dose <- c(0,0.006103516,0.024414063,0.09765625,0.39,1.5625,6.25,25,100)
test_df <- data.frame(dose=test_dose, response=test_response)
# fitting a Hill formula
m1 <- drm(response ~ dose, data=test_df, fct=LL.4())
# optaining ED50
ED(m1, 50)
我估计有ED50
的{{1}}。当我在专有程序中拟合相同的数据时,我得到了相似的图,但1649
的{{1}}(正确)
我尝试了log10以及计算出的EC50
中的45.5
,但是无论如何它都是关闭的
解决方案
因此我发现模型是正确的,但是显然在使用log2
时,可以在ED50
和ED()
之间进行选择,其中“相对”是默认值;在我的情况下,鉴于这些响应值是绝对的,因此更改此参数给了我正确的type="absolute"