r - 是否可以从模型中调用 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"])