我想知道是否有人知道如何限制 lm
中的某些系数或 gls
要积极?
例如,我想要 lm(y ~ x1 + x2 + ... + xn)
或 gls(y ~ x1 + x2 + ... + xn)
我想要,比如说 x1
的系数和 x2
为非负数。
最佳答案
一种选择是使用 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 .
关于r - lm() 或 gls() 的正系数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/21212782/