특정 범위 내의 여러 열에 대한 선형 회귀를 플롯

Question

여러 개의 열이있는 목록이 있습니다. 가장 왼쪽의 열은 x 축 데이터이고 나머지는 다양한 y 데이터 집합입니다.

각 열에 대해 선형 회귀를 생성하고 싶지만 x의 특정 범위에 대해서만 선형 회귀를 0,0으로 강제합니다. 그리고 마지막으로 상기 선형 회귀의 선 기울기를 추출한다. 아래는 내 데이터의 예입니다. 여기

x  y1  y2  y3  y4 
1 2.495 -1e-04 -1e-03 -1e-04 0e+00 
2 2.995 -2e-04 -7e-04 -2e-04 0e+00 
3 3.495 2e-04 -2e-04 1e-04 2e-04 
4 3.995 0e+00 -3e-04 0e+00 1e-04 
5 4.495 0e+00 -3e-04 0e+00 0e+00 
6 4.995 0e+00 -3e-04 0e+00 1e-04 
7 5.495 1e-04 -2e-04 0e+00 1e-04 
8 5.995 -1e-04 -1e-04 0e+00 3e-04 
9 6.495 1e-04 1e-04 2e-04 3e-04 
10 6.995 1e-04 0e+00 0e+00 3e-04 

예컨대 I 0,0 교차하도록 회귀 강요, 예를 들어 행 (X)의 3 내지 10의 범위에 대한 선형 회귀를 구하는 것이다. 마지막으로 y1, y2 등에 대한 선형 회귀 식의 기울기를 구한다.

나는 abind 및 lm 함수를 알고 있지만 0,0 교차를 강제하지 않고 단일 회귀를 플로팅 할 수만 있습니다. 감사합니다

Answer 1

회귀 모델은 절편을 사용하지 않고 회귀선을 원점을지나 가게합니다. 이것은 R에 -1을 모델에 추가하여 얻을 수 있습니다. 당신이 원하는 x의 범위를 얻으려면, 단순히 그 범위에 대한 관찰을 부분 집합 및 다른 독립 변수를 통해 루프 lapply를 사용

dat2 <- dat[3:9,] 
fits <- lapply(dat2[,-1], function(x) lm(x ~ dat2[,1] - 1)) 

ggplot2와 데이터 및 회귀 선을 플로팅. reshape2 패키지의 melt을 사용하면 데이터를 길게 만들 수 있습니다. 이 카테고리의 데이터를 파괴 ggplot2하는 데 도움이 :

library(ggplot2); library(reshape2) 
mdat <- melt(dat2,id="x") 
ggplot(data = mdat) + 
    geom_point(aes(y = value, x = x, color=variable)) + 
    geom_smooth(data=mdat,aes(x=x,y=value,color = variable),formula = y ~ x - 1, 
       method = "lm", se = F)