我目前正在使用nlme
执行混合效果回归。
我想通过为调用中的参数提供上限和下限来执行约束优化nlme
。
这可能吗?
如果您将上限定义为与传递给 的参数upper_bounds
的向量长度相同的向量,则可以通过以下方式设置上限:start
nlme
nlmeControl()
nlme(..., control = nlmeControl(opt = "nlminb", upper = upper_bounds))
提出问题时这可能是不可能的,但它在nlme
版本 3.1.137 中有效。
这里有两种简单的方法,不会弄乱 nlme 参数:1) 在边界上拟合一组模型并选择最适合的模型,以及 2) 使用参数的转换版本,将实数映射到所需的区间。