분류 전체보기 850

R을 이용한 ARIMA(Autoregressive integrated moving average) 분석 예

ARIMA 모델 (Autoregressive integrated moving average) ARIMA (p,d,q) 에서 AR - p MA - q I - 차분 횟수 변태는 아니지만, 연도에 따른 여자들 치마길이 분석을 ARIMA(Autoregressive integrated moving average)를 해본다. > skirts sts1 sts2 acf(sts2, lag.max=20) arima 자동 > auto.arima(skirts) Series: skirts ARIMA(1,2,0) Coefficients: ar1 -0.2997 s.e. 0.1424 sigma^2 estimated as 388.7: log likelihood=-193.66 AIC=391.33 AICc=391.62 BIC=394.9 ..

R을 이용한 시계열분석 함수 이용 예

> kings kingsts kingsts Time Series: Start = 1 End = 42 Frequency = 1 [1] 60 43 67 50 56 42 50 65 68 43 65 34 47 34 49 41 13 35 53 56 16 43 69 59 48 59 86 55 68 51 33 49 67 77 81 67 71 81 68 70 [41] 77 56 > births birthts souvenir plot.ts(kingsts) lot() 함수를 취해서 plot할수도 있다. 평탄화 (TTR패키지 SMA 함수) > SMA(kingsts, n=3) 시계열 요소 죄다 decompose > result result plot()을 그리면 예쁘다.

R을 이용한 회귀분석 함수 이용 예

> set.seed(2) > x = runif(10,0,11) 0~11 까지 랜덤하게 10개를 뽑는다. > y = runif(10,11,20) > d = data.frame(x,y) > d x y 1 2.033705 15.97407 2 7.726114 13.15005 3 6.306590 17.84462 4 1.848571 12.62738 5 10.382233 14.64754 6 10.378225 18.68194 7 1.420749 19.78759 8 9.167937 13.03243 9 5.148204 15.00328 10 6.049821 11.67481 > result = lm(y~x, data=d)> result Call: lm(formula = y ~ x, data = d) Coefficients: ..