lm() 或 gls() 的正系数

【问题标题】：Positive coefficients for lm() or gls()lm() 或 gls() 的正系数
【发布时间】：2014-02-08 08:53:16
【问题描述】：

我想知道是否有人知道如何将lm 或gls 中的某些系数限制为正数？

例如，我希望lm(y ~ x1 + x2 + ... + xn) 或gls(y ~ x1 + x2 + ... + xn) 并且我希望x1 和x2 的系数为非负数。

【问题讨论】：

【解决方案1】：

一种选择是使用nls(...)，它允许指定参数的上限和下限。所以，类似：

# not tested...
fit <- nls(y~a*x1 + b*x2 +c*x3...,data=mydata, 
           start=c(a=1,b=1,c=1...), lower=c(a=0,b=0,c=NA,...), algorithm="port")

阅读文档here。

【讨论】：