빅분기 실기연습(2)

R 프로그램_회귀 분석

프로그램 시작 전에 자신의 디렉토리 확인 및 설정부터...

---

getwd()

setwd("C:\\myRproject")

오늘 할 것은 근무 연수에 따른 연봉 그래프 그리고 1차 함수로 FITTING. 이걸 1차 회귀분석이라고 함

두 변수를 하나로 묶기 위해서는 data.frame(x변수, y변수)

그리고 간단한 통계 지수 보려면 summary()

---

year=c(26,16,20,7,22,15,29)

annual_salary=c(1267,887,1022,511,1193,795,1713)

Data=data.frame(year,annual_salary)

summary(Data)

> summary(Data)

      year       annual_salary 

 Min.   : 7.00   Min.   : 511  

 1st Qu.:15.50   1st Qu.: 841  

 Median :20.00   Median :1022  

 Mean   :19.29   Mean   :1055  

 3rd Qu.:24.00   3rd Qu.:1230  

 Max.   :29.00   Max.   :1713  

1차식 만들려면 ls=lm(y변수~x변수, data = ...)

그리고 plot과 abline

기울기와 y 절편 알려면 summary(ls)

---

plot(year,annual_salary)

abline(ls)

summary(ls)

> summary(ls)

Call:

lm(formula = annual_salary ~ year, data = Data)

Residuals:

        1         2         3         4         5 

-127.1104   -2.6907  -69.4586   75.2871    0.6575 

        6         7 

 -44.2487  167.5637 

Coefficients:

            Estimate Std. Error t value Pr(>|t|)    

(Intercept)   82.619    120.415   0.686 0.523151    

year          50.442      5.885   8.572 0.000356 ***

---

Signif. codes:  

0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 106.5 on 5 degrees of freedom

Multiple R-squared:  0.9363,Adjusted R-squared:  0.9235 

F-statistic: 73.48 on 1 and 5 DF,  p-value: 0.0003561

여기서 y 절편은 82.619, 기울기는 50.442

코드 정리

---

getwd()

setwd("C:\\myRproject")

year=c(26,16,20,7,22,15,29)

annual_salary=c(1267,887,1022,511,1193,795,1713)

Data=data.frame(year,annual_salary)

summary(Data)

ls=lm(annual_salary~year,data=Data)

plot(year,annual_salary)

abline(ls)

summary(ls)