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 설명서를 읽으십시오.