是否可以从模型中调用nls参数？

Question

我使用“nls”来拟合某些药物 - 受体结合数据的逻辑方程式;

Y= 100/(1+B*Exp(-k*X))

这个过程似乎运作良好（见图表）

我有B和k的估计值（来自我的等式的常数），我想用这些来估计激动剂（药物）的Log EC50值。我可以通过重新排列等式轻松完成，Y=50，X=ln(1/B)/-k.我遇到的问题是将B和k的值放入重新排列的等式中 - 如何从我为我估算的B模型中调用参数（k和nls）？

Answer 1

使用coef()。

尝试

with(as.list(coef(model)),log(1/B)/(-k))

作为

的快捷方式

cc <- coef(model)
log(1/cc["B"])/(-cc["k"])