使用DRM的剂量反应模型错误的ED50

时间:2019-08-20 07:45:46

标签: r drc

将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时,可以在ED50ED()之间进行选择,其中“相对”是默认值;在我的情况下,鉴于这些响应值是绝对的,因此更改此参数给了我正确的type="absolute"

0 个答案:

没有答案
相关问题