LM 계수 값을 강제로 입력하십시오.

Question

나사에 대한 무게에 관한 프로젝트를 진행하고 있습니다.

다음과 같이 여러 가지 조합을 반복하여 찾아 냈습니다.

head style  measurement intercept  slope 
a    1   x    1.0 
a    2   x    2.2 
a    3   x    4.1 
b    1   x    1.2 
b    2   x    2.0 
b    3   x    4.1 

동일한 측정 값을 가진 것들이 동일한 기울기를 가져야한다는 사실을 알고 있다면. 그래서 a-1은 b-1과 같은 기울기를 가져야합니다. 슬로프를 세트 번호로 설정 한 다른 작품을 실행할 수있는 방법이 있습니까?

Answer 1

nls 기능을 사용하여이를 수행 할 수 있습니다. 슬로프의 범위 제한을 슬로프를 원하는대로 정확하게 설정하면됩니다.

data(mtcars) 
df <- mtcars 

fitlm <- lm(hp~wt,data=df) 
summary(fitlm) 

fitnls <- nls(hp~a+b*wt,algorithm="port",data=df, 
       start=c(a=0,b=50),lower=c(a=-Inf,b=50),upper=c(a=Inf,b=50)) 
summary(fitnls) 

Answer 2

그것은 아직도 당신이 요구하는지 완전히 명확하지 않다,하지만의 원래 데이터를 이런 식으로 뭔가 보이는 가정하자 : 당신이 지금까지했던 그래서

head_style measurement x y 
a   1   . . 
a   1   . . 
a   1   . . 
a   2   . . 
... 
b   1   . . 

이 적합을 여기 예제 measurement과 head_style의 각 조합에 별도의 모델 (y~x 포함)을 사용하십시오. 이제 경사가 measurement에만 의존하고 가로 채기가 measurement과 head_slope에 의존하기를 원한다고 가정 해보십시오.

먼저 확인 measurement을하는 요소이다

my_dat <- transform(my_dat,measurement=factor(measurement)) 

이제 하나의 lm 모델에 맞게 :

lm(y~head_style:measurement + measurement:x, data=my_data) 

이것은 head_style 및 measurement하고 다른의 각 조합에 대해 당신에게 별도의 절편을 제공한다을 기울기는 각 값 measurement에 대해