0

gnls我正在尝试在nlme包中使用 R 中的逻辑增长模型。

我之前已经成功拟合了一个模型：

mod1 <- gnls(Weight ~ I(A/(1+exp(b + v0*Age + v1*Sum.T))),
                        data = df,
                        start = c(A= 13.157132, b= 3, v0= 0.16, v1= -0.0059),
                        na.action=na.omit)

但是，我现在希望限制b它不被模型拟合，因此尝试拟合第二个模型：

mod2 <- gnls(Weight ~ I(A/(1+exp(log((A/1.022)-1) + v0*Age + v1*Sum.T))),
                       data = df,
                       start = c(A= 13.157132, v0= 0.16, v1= -0.0059),
                       na.action=na.omit)

此模型返回错误：

Error in gnls(Weight ~ A/(1 + exp(log((A/1.022) - 1) + v0 * Age +  : 
approximate covariance matrix for parameter estimates not of full rank

Warning messages:
1: In log((A/1.022) - 1) : NaNs produced

搜索错误表明问题是由模型中的对称性引起的，具体问题的解决方案涉及使用不同参数调整公式。不幸的是，我的统计知识不足以 a) 完全理解问题或 b) 自己调整公式。

至于警告消息（总共有 15 条，都是一样的），我不明白它们为什么会出现，因为模型的这一部分单独工作（带有示例编号）。

任何这些查询的任何帮助将不胜感激。

4

1 回答 1

0

用户知道我最终用一个简单的解决方案（在朋友的帮助下）解决了这个问题可能会提供有用的信息。

由于 exp(a+b) = exp(a)*exp(b)，方程可以改写为：

Weight ~ I(A/(1+((A/1.022)-1) * exp(v0*Age + v1*Sum.T))

适合没有任何问题。一般来说，以不同的形式写出方程似乎是答案。

于 2015-01-27T14:30:17.000 回答