我正在使用 R 中的 nls 程序来拟合逻辑增长模型。在他们的 SSlogis 函数中,José Pinheiro 和 Douglas Bates 选择了公式
Asym / (1 + exp((xmid-input) / scal))
为他们的模型。由于我对此类模型的数值属性相当缺乏经验,我想知道:
有人可以解释为什么作者选择这个公式而不是可能的替代方案吗?特别是,生态学家似乎更喜欢具有初始人口、承载能力和增长率的模型。上面的公式是否具有有利的数值特性?
似乎当模型指定错误并且数据实际上与时间相当线性时,该公式通常无法收敛。这样的问题可以避免吗?
参数正交性是这里的关键问题还是模型的其他方面更重要?
扩展此模型以允许灵活拦截是否微不足道?以下模型会提供合理的数值属性吗?
Intercept + (Asym - Intercept ) / (1 + exp((xmid-input) / scal))当然,我对替代方案持开放态度,只要它允许在截距、已实现 50% 增长的位置和渐近线方面具有一定的灵活性。