r - lm() 或 gls() 的正系数

Question

我想知道是否有人知道如何限制 lm 中的某些系数或 gls要积极？

例如，我想要 lm(y ~ x1 + x2 + ... + xn)或 gls(y ~ x1 + x2 + ... + xn)我想要，比如说 x1 的系数和 x2为非负数。

Answer 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 .