能否从模型中调用nls参数？

Question

我使用“nls”将逻辑方程拟合到某些药物受体结合数据上；

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

这个过程似乎很顺利（见图表）

我有从方程中获得的B和k（常数）的估计值，并希望使用它们来估计激动剂（药物）的Log EC50值。当Y=50，X=ln(1/B)/-k.时，我可以轻松地通过重新排列方程来做到这一点。问题是如何从为我估计它们的nls模型中调用参数（B和k），以将这些值放入重新排列的方程中？

- G Boyd

---

如果您提供了一个包含样本输入数据的可重现示例，那么帮助您会更容易。通常，您可以使用predict()函数来完成这些操作。 - MrFlick

注意 - 如果您想为此数量获取置信区间，您可能需要研究Delta方法。 - Dason

Answer 1

使用coef()函数。

尝试一下。

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

作为...的快捷方式

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

- Ben Bolker

---

抱歉，本，我真的在尝试！但它一直告诉我我不能再接受另一个答案三分钟。我现在已经接受了答案。再次感谢你的帮助。 - G Boyd