1

将 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. 1649当我在专有程序中拟合相同的数据时，我得到了一个类似的图，但是 我尝试了 log10 和计算的（正确的） EC50，但无论如何它都关闭了45.5log2ED50

解决方案

所以我想通了，模型是正确的，但显然在使用时ED()可以选择type="absolute"，type="relative"其中“相对”是默认值；在我的情况下，鉴于这些响应值是绝对的，改变这个参数给了我权利ED50s

4

0 回答 0