계수에 선형 제한이 적용되는 R의 피팅 모델


16

계수를 바인딩하는 하나 이상의 정확한 선형 제한이있는 경우 R로 모델 공식을 어떻게 정의해야합니까? 예를 들어, 간단한 선형 회귀 모델에서 b1 = 2 * b0을 알고 있다고 가정하십시오.

감사합니다!

답변:


16

모델이

와이()=β0+β1엑스1()+β2엑스2()+ε()

예를 들어 다음과 같이 계수를 제한 할 계획입니다.

β1=2β2

제한을 삽입하고 얻을 원래 회귀 모델을 다시 작성

와이()=β0+2β2엑스1()+β2엑스2()+ε()

와이()=β0+β2(2엑스1()+엑스2())+ε()

()=2엑스1()+엑스2()

와이()=β0+β2()+ε()

이 방법으로 등호 의 수가 알 수없는 매개 변수의 수를 같은 수만큼 줄 이므로 정확한 제한을 처리 할 수 ​​있습니다 .

R 공식으로 연주하면 I () 함수로 직접 할 수 있습니다

lm(formula = Y ~ I(1 + 2*X1) + X2 + X3 - 1, data = <your data>) 
lm(formula = Y ~ I(2*X1 + X2) + X3, data = <your data>)

이것은 분명하지만 질문은 b0과 b1 사이의 제한을 암시하는 것이 었습니다. 또한 새로운 변수 Z = 2X + 1을 작성하고 절편없이 모형을 적합시켜야합니까?
George Dontas

2
나는 데는 보통 생각 나는 대신 수식 평가, 즉 Y ~ I (1 개 + 2 * X1) + X2 + X3-1의 사용
mpiktas

@ gd047 : 코드 조각으로 업데이트되었습니다. 그렇습니다. @mpiktas :이 변경됩니다, 그래 그것은 짧습니다;)
Dmitrij Celov

4
이것은 일반적인 이론적 접근 방식에 대한 좋은 해답이지만 R에서 이러한 가설을 실제로 쉽게 구현할 수있는 방법은 여러 모델을 추정 할 필요가 없다는 장점 linearHypothesis()이 있습니다 ( car패키지 참조) .
Jake Westfall
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.