> #membaca data peta

> setwd("D:/Idris")

> peta=readOGR(dsn="SULAWESI", layer="sul-sel")

OGR data source with driver: ESRI Shapefile

Source: "D:\Idris\SULAWESI", layer: "sul-sel"

with 24 features

It has 11 fields

Integer64 fields read as strings: Y

> plot(peta)

> text(peta,peta@data$ID,cex=0.7)

> #matriks pembobot

> queen<-poly2nb(peta,row.names(peta@data))

> plot(queen,coordinates(peta),col='red', add=TRUE)

> wqueen<-nb2mat(queen,style="W")

> wqueen

[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]

0 0.0000000 0.2000000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

1 0.1428571 0.0000000 0.00 0.0000000 0.1428571 0.0000000 0.0000000 0.0000000 0.0000000

2 0.0000000 0.0000000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

3 0.0000000 0.0000000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000

4 0.0000000 0.1428571 0.00 0.0000000 0.0000000 0.1428571 0.0000000 0.0000000 0.0000000

5 0.0000000 0.0000000 0.00 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000 0.0000000

6 0.0000000 0.0000000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 0.0000000

7 0.0000000 0.0000000 0.00 0.0000000 0.0000000 0.0000000 0.3333333 0.0000000 0.3333333

8 0.0000000 0.0000000 0.00 0.1428571 0.0000000 0.0000000 0.0000000 0.1428571 0.0000000

9 0.0000000 0.0000000 0.00 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000 0.0000000

10 0.0000000 0.2500000 0.00 0.0000000 0.2500000 0.0000000 0.0000000 0.0000000 0.0000000

11 0.0000000 0.0000000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.5000000

12 0.3333333 0.3333333 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

13 0.3333333 0.0000000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

14 0.0000000 0.0000000 0.00 0.2500000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

15 0.0000000 0.0000000 1.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

16 0.1428571 0.0000000 0.00 0.1428571 0.0000000 0.0000000 0.0000000 0.0000000 0.1428571

17 0.0000000 0.2500000 0.25 0.0000000 0.2500000 0.0000000 0.0000000 0.0000000 0.0000000

18 0.2500000 0.2500000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

19 0.0000000 0.0000000 0.00 0.0000000 0.3333333 0.3333333 0.0000000 0.0000000 0.0000000

20 0.0000000 0.0000000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000 0.2500000

21 0.0000000 0.0000000 0.00 0.2500000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000

22 0.0000000 0.2500000 0.00 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000

23 0.0000000 0.0000000 0.25 0.0000000 0.2500000 0.2500000 0.0000000 0.0000000 0.0000000

[,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18]

0 0.0000000 0.0000000 0.0000000 0.2000000 0.2000000 0.0000000 0.0000000 0.2000000

0.0000000

1 0.0000000 0.1428571 0.0000000 0.1428571 0.0000000 0.0000000 0.0000000 0.0000000

0.1428571
2 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.3333333 0.0000000
0.3333333

3 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000 0.0000000 0.2500000

0.0000000

4 0.1428571 0.1428571 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.1428571

5 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

6 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

7 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

8 0.0000000 0.0000000 0.1428571 0.0000000 0.0000000 0.0000000 0.0000000 0.1428571

0.0000000

9 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

10 0.2500000 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

11 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

12 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

13 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.3333333 0.0000000 0.3333333

0.0000000

14 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000 0.2500000

0.0000000

15 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

16 0.0000000 0.0000000 0.0000000 0.0000000 0.1428571 0.1428571 0.0000000 0.0000000

0.0000000

17 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

18 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000

0.0000000

19 0.3333333 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000

20 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000
21 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000
0.0000000

22 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000

0.0000000

23 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.2500000

[,19] [,20] [,21] [,22] [,23] [,24]

0 0.2000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

1 0.1428571 0.0000000 0.0000000 0.0000000 0.1428571 0.0000000

2 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.3333333

3 0.0000000 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000

4 0.0000000 0.1428571 0.0000000 0.0000000 0.0000000 0.1428571

5 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000 0.3333333

6 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

7 0.0000000 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000

8 0.0000000 0.0000000 0.1428571 0.1428571 0.1428571 0.0000000

9 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000 0.0000000

10 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

11 0.0000000 0.0000000 0.5000000 0.0000000 0.0000000 0.0000000

12 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

13 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

14 0.0000000 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000

15 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

16 0.1428571 0.0000000 0.0000000 0.0000000 0.1428571 0.0000000

17 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000

18 0.0000000 0.0000000 0.0000000 0.0000000 0.2500000 0.0000000

19 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

20 0.0000000 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000

21 0.0000000 0.0000000 0.2500000 0.0000000 0.0000000 0.0000000

22 0.2500000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

23 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

attr(,"call")
nb2mat(neighbours = queen, style = "W")

> pembobotspasial<-mat2listw(wqueen,style = "W")

> #memanggil data dalam bentuk csv

> Data <- read.csv("DATAFIX.csv", sep=";")

> Data

Kab.Kota Tahun APM JMMTP KP PP AL IPM

1 Bantaeng 2017 56.46 16916 469 856592 52.68 67.27

2 Bantaeng 2018 51.74 2846 471 867837 76.39 67.76

3 Bantaeng 2019 52.46 11200 474 908574 91.31 68.30

4 Barru 2017 60.26 3984 147 777788 68.09 69.56

5 Barru 2018 62.06 3821 148 890120 61.68 70.05

6 Barru 2019 62.70 3888 148 845883 78.39 70.60

7 Bone 2017 49.42 11391 165 649605 72.02 64.16

8 Bone 2018 50.96 19068 166 685530 77.77 65.04

9 Bone 2019 51.14 19263 166 742907 67.76 65.67

10 Bulukumba 2017 64.50 8008 360 752592 71.13 67.08

11 Bulukumba 2018 66.66 9167 362 858083 68.51 67.70

12 Bulukumba 2019 65.88 3777 364 835098 74.52 68.28

13 Enrekang 2017 69.37 6984 114 673795 57.60 71.44

14 Enrekang 2018 67.90 7145 115 767873 34.79 72.15

15 Enrekang 2019 68.85 7084 116 724410 58.92 72.66

16 Gowa 2017 59.70 4026 397 864763 82.39 68.33

17 Gowa 2018 59.04 16425 404 1008146 63.64 68.87

18 Gowa 2019 59.92 16551 410 935435 72.09 69.66

19 Jeneponto 2017 53.48 9067 398 667407 93.50 62.67

20 Jeneponto 2018 51.71 6754 401 671299 96.51 63.33

21 Jeneponto 2019 50.85 7811 403 698594 84.94 64.00

22 Luwu 2017 64.95 12512 119 698582 74.22 69.02

23 Luwu 2018 67.46 10204 120 878710 48.91 69.60

24 Luwu 2019 67.11 12037 121 852041 57.16 70.39

25 Luwu Timur 2017 64.47 9983 41 959369 82.79 71.46

26 Luwu Timur 2018 60.16 9960 42 1071700 53.43 72.16

27 Luwu Timur 2019 61.11 7818 43 1195872 30.47 72.80

28 Luwu Utara 2017 58.41 9783 41 727129 70.95 68.35

29 Luwu Utara 2018 60.66 8263 41 927567 60.29 68.79

30 Luwu Utara 2019 61.05 9342 42 908519 49.47 69.46

31 Makassar 2017 62.47 38788 8471 1536017 99.72 81.13

32 Makassar 2018 59.36 32507 8580 1576103 93.73 81.73

33 Makassar 2019 59.24 38478 8686 1513094 99.04 82.25

34 Maros 2017 62.77 7920 214 994814 80.13 68.42

35 Maros 2018 66.57 8144 216 1023702 64.41 68.94

36 Maros 2019 65.67 9828 218 1102522 62.58 69.50

37 Palopo 2017 66.42 5970 715 1100514 94.56 76.71

38 Palopo 2018 67.99 5733 730 1285357 78.07 77.30

39 Palopo 2019 67.40 5700 746 1188645 84.52 77.98

40 Pangkep 2017 58.62 5248 296 860812 78.37 67.25

41 Pangkep 2018 58.67 7944 299 1001377 77.71 67.71

42 Pangkep 2019 57.88 8053 302 927619 86.58 68.29

43 Parepare 2017 61.52 3952 1431 952831 96.99 76.68

44 Parepare 2018 62.67 3941 1447 1055320 72.74 77.19

45 Parepare 2019 62.35 3940 1462 1186848 69.82 77.62

46 Pinrang 2017 60.07 6517 190 855502 85.28 69.90

47 Pinrang 2018 60.78 6111 191 904833 54.69 70.62

48 Pinrang 2019 60.36 6250 192 903867 53.45 71.12

49 Selayar 2017 50.11 9812 147 941869 73.13 65.39

50 Selayar 2018 50.52 3265 149 936760 77.84 66.04

51 Selayar 2019 49.60 3430 150 1003365 85.45 66.91

52 Sidrap 2017 55.42 8240 157 840741 94.69 69.84

53 Sidrap 2018 57.96 6875 159 1006166 58.35 70.60

54 Sidrap 2019 58.45 7844 160 1009325 72.46 71.05

55 Sinjai 2017 72.62 19256 294 737954 76.38 65.80

56 Sinjai 2018 70.29 6991 296 773423 76.49 66.24

57 Sinjai 2019 69.22 7517 298 835167 85.12 67.05

58 Soppeng 2017 65.97 7640 167 607629 73.56 66.67

59 Soppeng 2018 68.88 4581 167 817740 57.70 67.60

60 Soppeng 2019 69.41 4994 167 843584 83.89 68.26

61 Takalar 2017 55.50 7439 517 849474 98.18 65.48

62 Takalar 2018 58.24 5301 522 976008 84.15 66.07

63 Takalar 2019 59.13 8806 527 903084 65.53 66.94

64 Tana Toraja 2017 48.85 7099 113 730696 72.28 66.82

65 Tana Toraja 2018 50.83 7455 113 800153 30.15 67.66

66 Tana Toraja 2019 51.58 7380 114 907188 28.54 68.25

67 Toraja Utara 2017 63.72 7791 198 670879 69.74 67.90

68 Toraja Utara 2018 68.13 7409 200 819426 42.23 68.49

69 Toraja Utara 2019 68.86 9975 201 825524 55.67 69.23

70 Wajo 2017 55.37 3320 158 942581 88.04 68.18

71 Wajo 2018 54.62 7888 158 966218 70.50 68.57

72 Wajo 2019 55.96 7818 159 1074893 69.42 69.05

> Idriss=pdata.frame(Data,c("Kab.Kota","Tahun"))

> Idris=Idriss[order(Idriss$Tahun),]

> #Analisis Regresi

> OLS <- plm(APM ~ JMMTP+KP+PP+AL+IPM, model = "within", effect="individual",

+ data = Idris, index = c("Tahun","Kab.Kota"))

> summary(OLS)

Oneway (individual) effect Within Model

Call:

plm(formula = APM ~ JMMTP + KP + PP + AL + IPM, data = Idris,

effect = "individual", model = "within", index = c("Tahun",

"Kab.Kota"))

Balanced Panel: n = 24, T = 3, N = 72

Residuals:

Min. 1st Qu. Median 3rd Qu. Max.

-2.24038 -0.67003 0.11849 0.50961 3.60971

Coefficients:

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

JMMTP 1.4010e-04 6.2586e-05 2.2386 0.03041 *

KP -1.5273e-02 8.9910e-03 -1.6987 0.09661 .

PP 6.2015e-06 4.1572e-06 1.4917 0.14307

AL -1.4063e-02 1.8174e-02 -0.7738 0.44328

IPM 9.7793e-02 4.5211e-01 0.2163 0.82977

---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Total Sum of Squares: 106.78

Residual Sum of Squares: 73.635

R-Squared: 0.31043

Adj. R-Squared: -0.13859

F-statistic: 3.87155 on 5 and 43 DF, p-value: 0.0055164

> summary(fixef(OLS))

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

Bantaeng 48.268 27.213 1.7737 0.08320 .

Barru 52.310 29.079 1.7989 0.07906 .

Bone 41.087 27.093 1.5165 0.13670

Bulukumba 59.559 27.534 2.1631 0.03614 *

Enrekang 58.655 30.347 1.9328 0.05986 .

Gowa 52.464 27.739 1.8914 0.06533 .

Jeneponto 47.913 25.846 1.8538 0.07064 .

Luwu 55.726 29.044 1.9187 0.06168 .

Luwu Timur 48.314 29.762 1.6234 0.11182

Luwu Utara 48.206 28.809 1.6733 0.10153

Makassar 170.074 72.659 2.3407 0.02395 *

Maros 54.869 27.936 1.9641 0.05600 .

Palopo 63.865 30.309 2.1071 0.04098 *

Pangkep 50.709 27.451 1.8473 0.07159 .

Parepare 70.694 30.198 2.3410 0.02394 *

Pinrang 50.939 29.056 1.7532 0.08670 .

Selayar 40.262 27.032 1.4894 0.14368

Sidrap 46.886 28.980 1.6179 0.11300

Sinjai 63.429 27.180 2.3337 0.02436 *

Soppeng 59.549 28.071 2.1214 0.03969 *

Takalar 53.641 26.344 2.0361 0.04793 *

Tana Toraja 40.092 28.110 1.4263 0.16101

Toraja Utara 58.073 28.404 2.0446 0.04705 *

Wajo 45.039 28.037 1.6064 0.11551

---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> RR=APM ~ JMMTP+KP+PP+AL+IPM

> slmtest(RR, data=Idriss, listw=pembobotspasial,model = "within",

+ test=c("lml"))

LM test for spatial lag dependence

data: formula (within transformation)

LM = 5.3433, df = 1, p-value = 0.0208

alternative hypothesis: spatial lag dependence

> slmtest(RR, data=Idriss, listw=pembobotspasial,model = "within",

+ test=c("lme"))

LM test for spatial error dependence

data: formula (within transformation)

LM = 4.7433, df = 1, p-value = 0.02941

alternative hypothesis: spatial error dependence

> R=APM ~ JMMTP+KP+PP+AL+IPM

> SARF <- spml(R, data = Idris, listw=pembobotspasial, effect = "individual",

+ model="within", spatial.error="none", lag=TRUE)

There were 38 warnings (use warnings() to see them)

> summary(SARF)

Spatial panel fixed effects lag model

Call:

spml(formula = R, data = Idris, listw = pembobotspasial, model = "within",

effect = "individual", lag = TRUE, spatial.error = "none")

Residuals:

Min. 1st Qu. Median 3rd Qu. Max.

-1.743150 -0.573702 0.072248 0.519633 2.962942

Spatial autoregressive coefficient:

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

lambda -0.50861 0.13445 -3.7829 0.000155 ***

Coefficients:

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

JMMTP 1.3473e-04 4.4005e-05 3.0617 0.002201 **

KP -1.4350e-02 6.3215e-03 -2.2700 0.023208 *

PP 4.1134e-06 2.9320e-06 1.4029 0.160643

AL -2.3042e-02 1.2863e-02 -1.7913 0.073239 .

IPM 3.9940e-01 3.2727e-01 1.2204 0.222310

---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> effects.splm(SARF)

Intercept:

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

(Intercept) 69.122 19.431 3.5573 0.0003747 ***

Spatial fixed effects:

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

1 -7.84075 19.13279 -0.4098 0.68195

2 -3.37198 20.44413 -0.1649 0.86899

3 -16.26400 19.04781 -0.8539 0.39319

4 0.62603 19.35820 0.0323 0.97420

5 0.36813 21.33561 0.0173 0.98623

6 -2.58363 19.50229 -0.1325 0.89461

7 -4.57162 18.17142 -0.2516 0.80136

8 -3.14265 20.41963 -0.1539 0.87769

9 -9.13242 20.92459 -0.4364 0.66251

10 -5.92670 20.25448 -0.2926 0.76982

11 105.44896 51.08361 2.0642 0.03899 *

12 -2.01179 19.64050 -0.1024 0.91841

13 3.80251 21.30911 0.1784 0.85837

14 -8.30558 19.29954 -0.4304 0.66694

15 8.47997 21.23105 0.3994 0.68959

16 -11.48137 20.42790 -0.5620 0.57409

17 -14.27374 19.00514 -0.7510 0.45262

18 -10.90654 20.37497 -0.5353 0.59245

19 5.94967 19.10900 0.3114 0.75553

20 4.08955 19.73560 0.2072 0.83584

21 -2.23371 18.52169 -0.1206 0.90401

22 -15.93082 19.76311 -0.8061 0.42019

23 1.32584 19.96962 0.0664 0.94707

24 -12.11337 19.71209 -0.6145 0.53888

---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> SEMF <- spml(R, data = Idris, listw = pembobotspasial,effect = "individual",

+ model="within", spatial.error="b",lag = FALSE)

There were 34 warnings (use warnings() to see them)

> summary(SEMF)

Spatial panel fixed effects error model

Call:

spml(formula = R, data = Idris, listw = pembobotspasial, model = "within",

effect = "individual", lag = FALSE, spatial.error = "b")

Residuals:

Min. 1st Qu. Median 3rd Qu. Max.

-2.38262 -0.71992 0.10047 0.64073 3.52436

Spatial error parameter:

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

rho -0.54335 0.14564 -3.7307 0.0001909 ***

Coefficients:

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

JMMTP 1.2270e-04 4.0253e-05 3.0483 0.002302 **

KP -1.1367e-02 6.0834e-03 -1.8685 0.061693 .

PP 3.8178e-06 2.9342e-06 1.3011 0.193208

AL -1.9360e-02 1.1571e-02 -1.6731 0.094310 .

IPM 1.4342e-01 2.6982e-01 0.5315 0.595044

---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> effects.splm(SEMF)

Intercept:

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

(Intercept) 54.544 19.441 2.8057 0.005021 **

Spatial fixed effects:

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

1 -8.54888 19.14218 -0.4466 0.6552

2 -3.57590 20.45416 -0.1748 0.8612

3 -14.74489 19.05715 -0.7737 0.4391

4 2.95547 19.36770 0.1526 0.8787

5 2.48361 21.34608 0.1163 0.9074

6 -3.97156 19.51186 -0.2035 0.8387

7 -8.84482 18.18034 -0.4865 0.6266

8 -0.01523 20.42965 -0.0007 0.9994

9 -6.66626 20.93485 -0.3184 0.7502

10 -7.12731 20.26442 -0.3517 0.7251

11 83.12266 51.10867 1.6264 0.1039

12 -0.66926 19.65013 -0.0341 0.9728

13 6.33535 21.31957 0.2972 0.7663

14 -5.32562 19.30901 -0.2758 0.7827

15 10.00895 21.24147 0.4712 0.6375

16 -5.00211 20.43792 -0.2447 0.8067

17 -15.07684 19.01447 -0.7929 0.4278

18 -8.69341 20.38497 -0.4265 0.6698

19 7.18109 19.11838 0.3756 0.7072

20 3.55518 19.74529 0.1801 0.8571

21 -3.23072 18.53077 -0.1743 0.8616

22 -15.68249 19.77280 -0.7931 0.4277

23 1.90363 19.97942 0.0953 0.9241

24 -10.37064 19.72176 -0.5258 0.5990

> #R-Square

> summary(SARF)$rsqr

[1] 0.9779748

> summary(SEMF)$rsqr

[1] 0.972926

> #sigma

> summary(SARF)$sigma2

[1] 0.8464543

> summary(SEMF)$sigma2

[1] 0.8472853

> #uji asumsi

> residual=SARF$residuals

> #normalitas

> normtest::jb.norm.test(residual)

Jarque-Bera test for normality

data: residual

JB = 0.72395, p-value = 0.6405

> nortest::lillie.test(residual)

Lilliefors (Kolmogorov-Smirnov) normality test

data: residual

D = 0.05956, p-value = 0.7619

> #autokorelasi

> Box.test(residual)

Box-Pierce test

data: residual

X-squared = 0.31862, df = 1, p-value = 0.5724

> #heteroskedastisitas

> Yduga=SARF$fitted.values

> plot(Yduga,residual)