Plot RDW

#radial distance weight - plot
plot(st_geometry(data_merge), border="blue",col='gray', main="Radial Distance Weight dengan dmax=15")
plot(w.rdw.s, koord, add = TRUE, col = "red")

Plot IDW

#invers distance weight - plot
plot(st_geometry(data_merge), border="blue",col='gray', main="Invers Distance Weight dengan alpha=1")
plot(w.idw.s, koord, add = TRUE, col = "red")

Plot EDW

#exponential distance weight - plot
plot(st_geometry(data_merge), border="blue",col='gray', main="Exponential Distance Weight dengan alpha=1")
plot(w.edw.s, koord, add = TRUE, col = "red")

Plot K-Nearest Neighbor Distance

#k-nearest neighbor distance weight - plot
plot(st_geometry(data_merge), border="blue",col='gray', main="K-Nearest Neighbor Distance (k=5)")
plot(w.knn.s, koord, add = TRUE, col = "red")

Plot Rook Distance

plot(st_geometry(data_merge), border="blue",col='gray')
coords <- st_coordinates(st_centroid(st_geometry(desa.pegunungan.bint)))
plot(rook, koord, add = TRUE, col = "red")

Plot Queen Distance

plot(st_geometry(data_merge), border="blue",col='gray')
coords <- st_coordinates(st_centroid(st_geometry(desa.pegunungan.bint)))
plot(queen, coords, add = TRUE, col = "red")

Indeks Moran Global

W.rdw <- nb2listw(w.rdw, style='W')
W.rdw2 <- nb2listw(w.rdw, style='B')

moran(data_merge$KEK, W.rdw, n=length(W.rdw$neighbours), S0=Szero(W.rdw))

## $I
## [1] 0.2774202
## 
## $K
## [1] 18.9482

IM <- moran.test(data_merge$KEK, W.rdw,randomisation=T,
alternative="greater", zero.policy = FALSE)
IM

## 
##  Moran I test under randomisation
## 
## data:  data_merge$KEK  
## weights: W.rdw    
## 
## Moran I statistic standard deviate = 11.281, p-value < 2.2e-16
## alternative hypothesis: greater
## sample estimates:
## Moran I statistic       Expectation          Variance 
##      0.2774202493     -0.0036231884      0.0006206612

Indeks Moran Lokal

LM <- localmoran(data_merge$KEK, W.rdw, alternative = "greater")
head(LM)

##              Ii          E.Ii       Var.Ii        Z.Ii Pr(z > E(Ii))
## 1  0.2035971603 -7.376709e-04 0.0402428298  1.01858704     0.1541995
## 2  0.0934837596 -7.376709e-04 0.0285327391  0.55779929     0.2884907
## 3  0.2096692831 -8.612240e-03 0.1397526672  0.58389791     0.2796445
## 4 -0.0372393247 -5.062083e-05 0.0005608504 -1.57031767     0.9418294
## 5  0.1685610783 -7.376709e-04 0.0178872021  1.26585086     0.1027832
## 6 -0.0004164114 -3.425436e-03 0.0361045585  0.01583597     0.4936826

data_merge$lmI <- LM[, "Ii"] # local Moran's I
data_merge$lmZ <- LM[, "Z.Ii"] # z-scores
# p-values corresponding to alternative greater
data_merge$lmp <- LM[, "Pr(z > E(Ii))"]

p1 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p1

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

p2 <- tm_shape(data_merge) +
  tm_polygons(col = "lmp", title = "p-value",
              breaks = c(-Inf, 0.05, Inf)) +
  tm_layout(legend.outside = TRUE)
p2

Klasifikasi Warna :

• Less than 0.05 = Ada Autokorelasi Positif

• 0.05 or more = Tidak Ada Autokorelasi Positif

tmap_mode("view")

## tmap mode set to interactive viewing

p3 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p3

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ",
title = "Local Moran's I", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Negative SAC", "No SAC", "Positive SAC"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Rentang nilai untuk klasifikasi warna, dengan:

• -Inf hingga -1.96 mewakili Negative Spatial Autocorrelation (SAC).

• -1.96 hingga 1.96 mewakili No Spatial Autocorrelation (No SAC).

• 1.96 hingga Inf mewakili Positive Spatial Autocorrelation (Positive SAC)

# Autokorelasi spasial positif (Z-score > 1.96)
positive_autocorr <- data_merge[data_merge$lmZ > 1.96, ]
# Autokorelasi spasial negatif (Z-score < -1.96)
negative_autocorr <- data_merge[data_merge$lmZ < -1.96, ]
# Area tanpa autokorelasi spasial (Z-score antara -1.96 dan 1.96)
no_autocorr <- data_merge[data_merge$lmZ >= -1.96 & data_merge$lmZ <= 1.96, ]

moran.plot(data_merge$KEK, W.rdw)

Moran Plot memiliki kemiringan positif: ini menunjukkan adanya autokorelasi spasial positif, di mana lokasi dengan nilai tinggi cenderung berdekatan dengan lokasi bernilai tinggi, dan lokasi bernilai rendah cenderung berdekatan dengan lokasi bernilai rendah.

Getis-Ord Gi

Mendeteksi adanya wilayah yang hotspot/coldspot.

Local_M <- localG(data_merge$KEK, W.rdw, alternative = "greater")
Local_M

##   [1] -1.01858704 -0.55779929  0.58389791 -1.57031767 -1.26585086  0.01583597
##   [7] -0.29187721 -0.39699968 -1.28890004 -0.02435708 -0.69798031 -0.77770318
##  [13] -0.78609873  1.28923707 -1.35317274 -1.73344938  0.33085763 -1.01858704
##  [19] -0.18198005 -1.43174321 -0.78609873 -0.62135822 -1.34748412 -1.19784904
##  [25]  0.28639421  7.28303689  0.85162797  4.68939204  6.84520368 -0.90937568
##  [31] -0.79972206 -0.88373988 -0.63901149 -0.87540779 -1.73344938  0.89488138
##  [37] -0.78609873  0.26442350 -1.29561327 -0.77509302 -1.20968047  5.84660804
##  [43] -1.43712578 -1.37886105 -1.42167103 -1.92119908  0.52113001 -0.47605657
##  [49] -0.78651694 -0.36856416 -0.49958416 -0.64933258 -1.26638264  8.17571318
##  [55]  7.80623778  7.39264672 -0.55203869  5.25144564 -0.57123687  1.44240627
##  [61] -0.62124319  0.98681165 -0.14913046  0.33160457 -0.64362699 -1.01243401
##  [67] -1.56786372 -1.11069373 -0.55614890  6.93801841 -1.34748412 -1.34748412
##  [73] -0.31377048 -0.47689405  2.84064219 -0.90999658 -1.20238198 -1.01858704
##  [79] -1.86024912 -1.11787060 -0.12649199 -0.33221411  0.11947923 -1.20968047
##  [85] -1.17369127 -0.87716036  2.98092285 -0.88833621  6.96011843 -1.20968047
##  [91] -0.72268193  0.23311494 -0.55779929 -1.11787060 -1.73344938  1.07604490
##  [97] -0.82399472 -1.19776267  0.56450153 -1.01487395  0.69902352  0.58563925
## [103]  0.27296867  1.18523295 -1.20968047 -1.01858704  6.11230046 -0.18525077
## [109]  2.68899730 -0.38971324  7.03693285 -0.28437039  4.29245859 -0.69070753
## [115]  1.18523295  8.09201990 -0.78609873 -1.11787060 -0.77509302  3.90808405
## [121] -1.51372098  0.37895835 -1.92119908  7.77213613  3.10915718  0.84923421
## [127]  1.15701945 -1.92119908 -1.01858704  5.14794316 -1.66721494  0.99412244
## [133]  2.08441126 -0.52789255 -0.52789255 -0.90937568  3.03453971 -1.17078885
## [139] -1.26638264  0.09261845 -1.19030519 -0.88373988 -0.63260206 -0.98546246
## [145]  3.03453971  6.19149184 -0.12331651 -1.59877169 -0.96888758 -0.64067460
## [151] -1.08738024 -0.84490896  4.24127024  0.72723853 -1.42098051 -0.89050958
## [157]  0.74266089  0.49631440 -1.79772462 -1.07876005 -0.76646709 -1.57167816
## [163] -1.63699651  0.65240582 -1.03284927 -0.95639687 -0.82241181  0.82745842
## [169]  0.61955808  0.14563246 -1.28174329  0.83543909  0.09770559 -0.83944271
## [175] -0.88250534  1.13753378 -1.20968047 -0.42637502 -1.70804715 -0.88373988
## [181]  1.24606004 -0.94155138 -1.64096007 -0.69070753 -1.49476738  0.39513290
## [187]  0.23194836 -0.49745472  0.65240582 -1.09120711 -1.00388523 -1.20238198
## [193]  0.53025887 -0.59722757 -2.03895683 -0.78609873  0.83229286 -1.01858704
## [199] -0.62734701 -0.88373988 -0.50627720  5.14794316  0.36552307  0.09408189
## [205] -1.54919695 -1.42532709 -1.45397511  1.28923707 -1.98072503 -1.52611475
## [211]  5.09369832  3.34977143 -1.01858704 -0.88250534  0.99102489 -0.31884344
## [217]  0.31137561 -2.20694888 -0.78609873  0.65244978 -0.41009362 -0.07285891
## [223]  7.37869551  1.22842152 -0.31702708  0.99110659  0.16777281  1.24067261
## [229] -0.78609873 -1.42167103  6.44163321 -0.47216997  6.21896335  8.31934179
## [235] -1.11958522  0.70767470 -1.34748412  0.53506268 -1.01858704 -0.63183300
## [241] -1.01858704 -1.11958522  5.25703906  0.75237819 -0.30941180 -0.75986964
## [247] -1.42532709 -1.54876067 -0.62936479  2.55208567 -1.03187408 -0.41501916
## [253] -0.70354633 -2.09600737  7.67646100 -1.20968047 -1.01858704  5.18768135
## [259] -1.11787060 -1.00388523  0.21883148  2.07309136 -1.43124664 -0.82924147
## [265]  4.08227717 -1.00388523 -1.79772462  1.34713048  4.35430322  0.06585946
## [271] -0.63435466  0.99187727  5.51417690  3.19281325  3.11354752  0.19256523
## [277] -0.39534412
## attr(,"internals")
##                  Gi       E(Gi)        V(Gi)       Z(Gi) Pr(z > E(Gi))
##   [1,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
##   [2,] 0.0019524894 0.003623188 8.970993e-06 -0.55779929  7.115093e-01
##   [3,] 0.0047743056 0.003623188 3.886558e-06  0.58389791  2.796445e-01
##   [4,] 0.0010944185 0.003623188 2.593245e-06 -1.57031767  9.418294e-01
##   [5,] 0.0006212466 0.003623188 5.623924e-06 -1.26585086  8.972168e-01
##   [6,] 0.0036482335 0.003623188 2.501240e-06  0.01583597  4.936826e-01
##   [7,] 0.0022779043 0.003623188 2.124358e-05 -0.29187721  6.148097e-01
##   [8,] 0.0030633886 0.003623188 1.988315e-06 -0.39699968  6.543162e-01
##   [9,] 0.0012655024 0.003623188 3.346058e-06 -1.28890004  9.012836e-01
##  [10,] 0.0035755149 0.003623188 3.830931e-06 -0.02435708  5.097161e-01
##  [11,] 0.0024082569 0.003623188 3.029823e-06 -0.69798031  7.574052e-01
##  [12,] 0.0023769436 0.003623188 2.567905e-06 -0.77770318  7.816280e-01
##  [13,] 0.0000000000 0.003623188 2.124358e-05 -0.78609873  7.840952e-01
##  [14,] 0.0054157132 0.003623188 1.933147e-06  1.28923707  9.865784e-02
##  [15,] 0.0004141644 0.003623188 5.623924e-06 -1.35317274  9.119998e-01
##  [16,] 0.0000000000 0.003623188 4.368773e-06 -1.73344938  9.584921e-01
##  [17,] 0.0041701418 0.003623188 2.732866e-06  0.33085763  3.703760e-01
##  [18,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
##  [19,] 0.0032903063 0.003623188 3.346058e-06 -0.18198005  5.722008e-01
##  [20,] 0.0010790073 0.003623188 3.157663e-06 -1.43174321  9.238913e-01
##  [21,] 0.0000000000 0.003623188 2.124358e-05 -0.78609873  7.840952e-01
##  [22,] 0.0007593014 0.003623188 2.124358e-05 -0.62135822  7.328180e-01
##  [23,] 0.0005694761 0.003623188 5.135810e-06 -1.34748412  9.110878e-01
##  [24,] 0.0016566577 0.003623188 2.695239e-06 -1.19784904  8.845121e-01
##  [25,] 0.0041536864 0.003623188 3.431145e-06  0.28639421  3.872881e-01
##  [26,] 0.0138306138 0.003623188 1.964298e-06  7.28303689  1.631940e-13
##  [27,] 0.0054187192 0.003623188 4.445145e-06  0.85162797  1.972103e-01
##  [28,] 0.0124899275 0.003623188 3.575154e-06  4.68939204  1.370090e-06
##  [29,] 0.0127145582 0.003623188 1.763948e-06  6.84520368  3.818357e-12
##  [30,] 0.0000000000 0.003623188 1.587432e-05 -0.90937568  8.184241e-01
##  [31,] 0.0015186029 0.003623188 6.925562e-06 -0.79972206  7.880641e-01
##  [32,] 0.0009762447 0.003623188 8.970993e-06 -0.88373988  8.115817e-01
##  [33,] 0.0000000000 0.003623188 3.214878e-05 -0.63901149  7.385923e-01
##  [34,] 0.0009829620 0.003623188 9.096232e-06 -0.87540779  8.093240e-01
##  [35,] 0.0000000000 0.003623188 4.368773e-06 -1.73344938  9.584921e-01
##  [36,] 0.0048066876 0.003623188 1.749061e-06  0.89488138  1.854253e-01
##  [37,] 0.0000000000 0.003623188 2.124358e-05 -0.78609873  7.840952e-01
##  [38,] 0.0045662100 0.003623188 1.271871e-05  0.26442350  3.957268e-01
##  [39,] 0.0000000000 0.003623188 7.820438e-06 -1.29561327  9.024456e-01
##  [40,] 0.0013016596 0.003623188 8.970993e-06 -0.77509302  7.808577e-01
##  [41,] 0.0000000000 0.003623188 8.970993e-06 -1.20968047  8.867992e-01
##  [42,] 0.0125740319 0.003623188 2.343797e-06  5.84660804  2.508489e-09
##  [43,] 0.0011389522 0.003623188 2.988108e-06 -1.43712578  9.246589e-01
##  [44,] 0.0013016596 0.003623188 2.834700e-06 -1.37886105  9.160312e-01
##  [45,] 0.0008542141 0.003623188 3.793496e-06 -1.42167103  9.224391e-01
##  [46,] 0.0000000000 0.003623188 3.556617e-06 -1.92119908  9.726467e-01
##  [47,] 0.0047635850 0.003623188 4.788723e-06  0.52113001  3.011381e-01
##  [48,] 0.0026212320 0.003623188 4.429764e-06 -0.47605657  6.829830e-01
##  [49,] 0.0014236902 0.003623188 7.820438e-06 -0.78651694  7.842177e-01
##  [50,] 0.0029816514 0.003623188 3.029823e-06 -0.36856416  6.437737e-01
##  [51,] 0.0029681784 0.003623188 1.719011e-06 -0.49958416  6.913160e-01
##  [52,] 0.0015186029 0.003623188 1.050507e-05 -0.64933258  7.419383e-01
##  [53,] 0.0009762447 0.003623188 4.368773e-06 -1.26638264  8.973119e-01
##  [54,] 0.0162013730 0.003623188 2.366926e-06  8.17571318  1.470604e-16
##  [55,] 0.0156250000 0.003623188 2.363796e-06  7.80623778  2.946019e-15
##  [56,] 0.0176335981 0.003623188 3.591714e-06  7.39264672  7.196727e-14
##  [57,] 0.0027105922 0.003623188 2.732866e-06 -0.55203869  7.095391e-01
##  [58,] 0.0116628702 0.003623188 2.343797e-06  5.25144564  7.545504e-08
##  [59,] 0.0025458931 0.003623188 3.556617e-06 -0.57123687  7.160805e-01
##  [60,] 0.0058314351 0.003623188 2.343797e-06  1.44240627  7.459390e-02
##  [61,] 0.0011415525 0.003623188 1.595706e-05 -0.62124319  7.327802e-01
##  [62,] 0.0049233756 0.003623188 1.735974e-06  0.98681165  1.618675e-01
##  [63,] 0.0034168565 0.003623188 1.914256e-06 -0.14913046  5.592747e-01
##  [64,] 0.0042153439 0.003623188 3.188823e-06  0.33160457  3.700939e-01
##  [65,] 0.0022779043 0.003623188 4.368773e-06 -0.64362699  7.400913e-01
##  [66,] 0.0013287775 0.003623188 5.135810e-06 -1.01243401  8.443347e-01
##  [67,] 0.0005694761 0.003623188 3.793496e-06 -1.56786372  9.415435e-01
##  [68,] 0.0013016596 0.003623188 4.368773e-06 -1.11069373  8.666499e-01
##  [69,] 0.0027524677 0.003623188 2.451182e-06 -0.55614890  7.109455e-01
##  [70,] 0.0150134603 0.003623188 2.695239e-06  6.93801841  1.988190e-12
##  [71,] 0.0005694761 0.003623188 5.135810e-06 -1.34748412  9.110878e-01
##  [72,] 0.0005694761 0.003623188 5.135810e-06 -1.34748412  9.110878e-01
##  [73,] 0.0031890661 0.003623188 1.914256e-06 -0.31377048  6.231523e-01
##  [74,] 0.0028202625 0.003623188 2.834700e-06 -0.47689405  6.832812e-01
##  [75,] 0.0100607441 0.003623188 5.135810e-06  2.84064219  2.251140e-03
##  [76,] 0.0018507973 0.003623188 3.793496e-06 -0.90999658  8.185878e-01
##  [77,] 0.0012813212 0.003623188 3.793496e-06 -1.20238198  8.853922e-01
##  [78,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
##  [79,] 0.0000000000 0.003623188 3.793496e-06 -1.86024912  9.685749e-01
##  [80,] 0.0000000000 0.003623188 1.050507e-05 -1.11787060  8.681889e-01
##  [81,] 0.0032903063 0.003623188 6.925562e-06 -0.12649199  5.503288e-01
##  [82,] 0.0026212320 0.003623188 9.096232e-06 -0.33221411  6.301362e-01
##  [83,] 0.0045871560 0.003623188 6.509385e-05  0.11947923  4.524478e-01
##  [84,] 0.0000000000 0.003623188 8.970993e-06 -1.20968047  8.867992e-01
##  [85,] 0.0018525053 0.003623188 2.276009e-06 -1.17369127  8.797406e-01
##  [86,] 0.0017897820 0.003623188 4.368773e-06 -0.87716036  8.098002e-01
##  [87,] 0.0089740480 0.003623188 3.222149e-06  2.98092285  1.436906e-03
##  [88,] 0.0011389522 0.003623188 7.820438e-06 -0.88833621  8.128200e-01
##  [89,] 0.0187489049 0.003623188 4.722790e-06  6.96011843  1.699934e-12
##  [90,] 0.0000000000 0.003623188 8.970993e-06 -1.20968047  8.867992e-01
##  [91,] 0.0018223235 0.003623188 6.209661e-06 -0.72268193  7.650623e-01
##  [92,] 0.0040496077 0.003623188 3.346058e-06  0.23311494  4.078361e-01
##  [93,] 0.0019524894 0.003623188 8.970993e-06 -0.55779929  7.115093e-01
##  [94,] 0.0000000000 0.003623188 1.050507e-05 -1.11787060  8.681889e-01
##  [95,] 0.0000000000 0.003623188 4.368773e-06 -1.73344938  9.584921e-01
##  [96,] 0.0051223242 0.003623188 1.940980e-06  1.07604490  1.409536e-01
##  [97,] 0.0022935780 0.003623188 2.603754e-06 -0.82399472  7.950287e-01
##  [98,] 0.0012842466 0.003623188 3.813268e-06 -1.19776267  8.844953e-01
##  [99,] 0.0043323140 0.003623188 1.578033e-06  0.56450153  2.862064e-01
## [100,] 0.0021026809 0.003623188 2.244672e-06 -1.01487395  8.449170e-01
## [101,] 0.0045558087 0.003623188 1.780025e-06  0.69902352  2.422687e-01
## [102,] 0.0043988270 0.003623188 1.754117e-06  0.58563925  2.790590e-01
## [103,] 0.0040606121 0.003623188 2.567905e-06  0.27296867  3.924387e-01
## [104,] 0.0052410128 0.003623188 1.863182e-06  1.18523295  1.179627e-01
## [105,] 0.0000000000 0.003623188 8.970993e-06 -1.20968047  8.867992e-01
## [106,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
## [107,] 0.0127900603 0.003623188 2.249225e-06  6.11230046  4.910253e-10
## [108,] 0.0033029613 0.003623188 2.988108e-06 -0.18525077  5.734838e-01
## [109,] 0.0097701149 0.003623188 5.225590e-06  2.68899730  3.583350e-03
## [110,] 0.0028604119 0.003623188 3.830931e-06 -0.38971324  6.516257e-01
## [111,] 0.0143963554 0.003623188 2.343797e-06  7.03693285  9.825875e-13
## [112,] 0.0032059394 0.003623188 2.152890e-06 -0.28437039  6.119367e-01
## [113,] 0.0097056735 0.003623188 2.007936e-06  4.29245859  8.835276e-06
## [114,] 0.0022779043 0.003623188 3.793496e-06 -0.69070753  7.551253e-01
## [115,] 0.0052410128 0.003623188 1.863182e-06  1.18523295  1.179627e-01
## [116,] 0.0148741419 0.003623188 1.933147e-06  8.09201990  2.934166e-16
## [117,] 0.0000000000 0.003623188 2.124358e-05 -0.78609873  7.840952e-01
## [118,] 0.0000000000 0.003623188 1.050507e-05 -1.11787060  8.681889e-01
## [119,] 0.0013016596 0.003623188 8.970993e-06 -0.77509302  7.808577e-01
## [120,] 0.0110126242 0.003623188 3.575154e-06  3.90808405  4.651547e-05
## [121,] 0.0011911852 0.003623188 2.581289e-06 -1.51372098  9.349517e-01
## [122,] 0.0051369863 0.003623188 1.595706e-05  0.37895835  3.523594e-01
## [123,] 0.0000000000 0.003623188 3.556617e-06 -1.92119908  9.726467e-01
## [124,] 0.0146958305 0.003623188 2.029653e-06  7.77213613  3.858672e-15
## [125,] 0.0089977221 0.003623188 2.988108e-06  3.10915718  9.381094e-04
## [126,] 0.0047479402 0.003623188 1.754117e-06  0.84923421  1.978755e-01
## [127,] 0.0054472477 0.003623188 2.485402e-06  1.15701945  1.236322e-01
## [128,] 0.0000000000 0.003623188 3.556617e-06 -1.92119908  9.726467e-01
## [129,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
## [130,] 0.0130644343 0.003623188 3.363498e-06  5.14794316  1.316791e-07
## [131,] 0.0000000000 0.003623188 4.722790e-06 -1.66721494  9.522642e-01
## [132,] 0.0054416603 0.003623188 3.346058e-06  0.99412244  1.600816e-01
## [133,] 0.0067461013 0.003623188 2.244672e-06  2.08441126  1.856139e-02
## [134,] 0.0026575550 0.003623188 3.346058e-06 -0.52789255  7.012130e-01
## [135,] 0.0026575550 0.003623188 3.346058e-06 -0.52789255  7.012130e-01
## [136,] 0.0000000000 0.003623188 1.587432e-05 -0.90937568  8.184241e-01
## [137,] 0.0177880897 0.003623188 2.178921e-05  3.03453971  1.204515e-03
## [138,] 0.0016382700 0.003623188 2.874274e-06 -1.17078885  8.791582e-01
## [139,] 0.0009762447 0.003623188 4.368773e-06 -1.26638264  8.973119e-01
## [140,] 0.0037844037 0.003623188 3.029823e-06  0.09261845  4.631033e-01
## [141,] 0.0018398458 0.003623188 2.244672e-06 -1.19030519  8.830368e-01
## [142,] 0.0009762447 0.003623188 8.970993e-06 -0.88373988  8.115817e-01
## [143,] 0.0000000000 0.003623188 3.280353e-05 -0.63260206  7.365032e-01
## [144,] 0.0010170353 0.003623188 6.993904e-06 -0.98546246  8.378015e-01
## [145,] 0.0177880897 0.003623188 2.178921e-05  3.03453971  1.204515e-03
## [146,] 0.0131116390 0.003623188 2.348551e-06  6.19149184  2.979871e-10
## [147,] 0.0032978184 0.003623188 6.961658e-06 -0.12331651  5.490718e-01
## [148,] 0.0000000000 0.003623188 5.135810e-06 -1.59877169  9.450643e-01
## [149,] 0.0016704632 0.003623188 4.061959e-06 -0.96888758  8.336994e-01
## [150,] 0.0000000000 0.003623188 3.198209e-05 -0.64067460  7.391329e-01
## [151,] 0.0010377750 0.003623188 5.653236e-06 -1.08738024  8.615656e-01
## [152,] 0.0017084282 0.003623188 5.135810e-06 -0.84490896  8.009192e-01
## [153,] 0.0094912680 0.003623188 1.914256e-06  4.24127024  1.111292e-05
## [154,] 0.0049154778 0.003623188 3.157663e-06  0.72723853  2.335399e-01
## [155,] 0.0007593014 0.003623188 4.061959e-06 -1.42098051  9.223388e-01
## [156,] 0.0004555809 0.003623188 1.265277e-05 -0.89050958  8.134038e-01
## [157,] 0.0045871560 0.003623188 1.684782e-06  0.74266089  2.288435e-01
## [158,] 0.0042694497 0.003623188 1.695519e-06  0.49631440  3.098363e-01
## [159,] 0.0000000000 0.003623188 4.061959e-06 -1.79772462  9.638897e-01
## [160,] 0.0007645260 0.003623188 7.022246e-06 -1.07876005  8.596526e-01
## [161,] 0.0011389522 0.003623188 1.050507e-05 -0.76646709  7.783008e-01
## [162,] 0.0004555809 0.003623188 4.061959e-06 -1.57167816  9.419874e-01
## [163,] 0.0005359775 0.003623188 3.556617e-06 -1.63699651  9.491844e-01
## [164,] 0.0044560944 0.003623188 1.629882e-06  0.65240582  2.570697e-01
## [165,] 0.0014643671 0.003623188 4.368773e-06 -1.03284927  8.491628e-01
## [166,] 0.0021902926 0.003623188 2.244672e-06 -0.95639687  8.305641e-01
## [167,] 0.0024406118 0.003623188 2.067664e-06 -0.82241181  7.945787e-01
## [168,] 0.0046898030 0.003623188 1.661586e-06  0.82745842  2.039886e-01
## [169,] 0.0044218143 0.003623188 1.661586e-06  0.61955808  2.677744e-01
## [170,] 0.0038210008 0.003623188 1.844976e-06  0.14563246  4.421058e-01
## [171,] 0.0012059494 0.003623188 3.556617e-06 -1.28174329  9.000336e-01
## [172,] 0.0047261607 0.003623188 1.743009e-06  0.83543909  2.017353e-01
## [173,] 0.0037421535 0.003623188 1.482520e-06  0.09770559  4.610830e-01
## [174,] 0.0023655160 0.003623188 2.244672e-06 -0.83944271  7.993895e-01
## [175,] 0.0021743632 0.003623188 2.695239e-06 -0.88250534  8.112482e-01
## [176,] 0.0051790471 0.003623188 1.870732e-06  1.13753378  1.276576e-01
## [177,] 0.0000000000 0.003623188 8.970993e-06 -1.20968047  8.867992e-01
## [178,] 0.0030100879 0.003623188 2.067664e-06 -0.42637502  6.650827e-01
## [179,] 0.0004019831 0.003623188 3.556617e-06 -1.70804715  9.561862e-01
## [180,] 0.0009762447 0.003623188 8.970993e-06 -0.88373988  8.115817e-01
## [181,] 0.0062541145 0.003623188 4.457993e-06  1.24606004  1.063712e-01
## [182,] 0.0015770107 0.003623188 4.722790e-06 -0.94155138  8.267888e-01
## [183,] 0.0004271071 0.003623188 3.793496e-06 -1.64096007  9.495971e-01
## [184,] 0.0022779043 0.003623188 3.793496e-06 -0.69070753  7.551253e-01
## [185,] 0.0007118451 0.003623188 3.793496e-06 -1.49476738  9.325124e-01
## [186,] 0.0041042974 0.003623188 1.482520e-06  0.39513290  3.463724e-01
## [187,] 0.0039863326 0.003623188 2.451182e-06  0.23194836  4.082891e-01
## [188,] 0.0028882093 0.003623188 2.182945e-06 -0.49745472  6.905658e-01
## [189,] 0.0044560944 0.003623188 1.629882e-06  0.65240582  2.570697e-01
## [190,] 0.0010354111 0.003623188 5.623924e-06 -1.09120711  8.624091e-01
## [191,] 0.0012424933 0.003623188 5.623924e-06 -1.00388523  8.422830e-01
## [192,] 0.0012813212 0.003623188 3.793496e-06 -1.20238198  8.853922e-01
## [193,] 0.0042954767 0.003623188 1.607442e-06  0.53025887  2.979662e-01
## [194,] 0.0024481567 0.003623188 3.870967e-06 -0.59722757  7.248223e-01
## [195,] 0.0000000000 0.003623188 3.157663e-06 -2.03895683  9.792728e-01
## [196,] 0.0000000000 0.003623188 2.124358e-05 -0.78609873  7.840952e-01
## [197,] 0.0054319257 0.003623188 4.722790e-06  0.83229286  2.026218e-01
## [198,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
## [199,] 0.0013761468 0.003623188 1.282941e-05 -0.62734701  7.347841e-01
## [200,] 0.0009762447 0.003623188 8.970993e-06 -0.88373988  8.115817e-01
## [201,] 0.0018223235 0.003623188 1.265277e-05 -0.50627720  6.936690e-01
## [202,] 0.0130644343 0.003623188 3.363498e-06  5.14794316  1.316791e-07
## [203,] 0.0043173233 0.003623188 3.606269e-06  0.36552307  3.573605e-01
## [204,] 0.0039049788 0.003623188 8.970993e-06  0.09408189  4.625221e-01
## [205,] 0.0006761325 0.003623188 3.618791e-06 -1.54919695  9.393328e-01
## [206,] 0.0005256702 0.003623188 4.722790e-06 -1.42532709  9.229687e-01
## [207,] 0.0000000000 0.003623188 6.209661e-06 -1.45397511  9.270234e-01
## [208,] 0.0054157132 0.003623188 1.933147e-06  1.28923707  9.865784e-02
## [209,] 0.0000000000 0.003623188 3.346058e-06 -1.98072503  9.761889e-01
## [210,] 0.0012338648 0.003623188 2.451182e-06 -1.52611475  9.365094e-01
## [211,] 0.0117856789 0.003623188 2.567905e-06  5.09369832  1.755727e-07
## [212,] 0.0097666159 0.003623188 3.363498e-06  3.34977143  4.043914e-04
## [213,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
## [214,] 0.0021743632 0.003623188 2.695239e-06 -0.88250534  8.112482e-01
## [215,] 0.0049568966 0.003623188 1.811142e-06  0.99102489  1.608367e-01
## [216,] 0.0027841053 0.003623188 6.925562e-06 -0.31884344  6.250774e-01
## [217,] 0.0045558087 0.003623188 8.970993e-06  0.31137561  3.777575e-01
## [218,] 0.0000000000 0.003623188 2.695239e-06 -2.20694888  9.863412e-01
## [219,] 0.0000000000 0.003623188 2.124358e-05 -0.78609873  7.840952e-01
## [220,] 0.0044727568 0.003623188 1.695519e-06  0.65244978  2.570555e-01
## [221,] 0.0027660267 0.003623188 4.368773e-06 -0.41009362  6.591314e-01
## [222,] 0.0034403670 0.003623188 6.296351e-06 -0.07285891  5.290408e-01
## [223,] 0.0145903479 0.003623188 2.209170e-06  7.37869551  7.992397e-14
## [224,] 0.0061162080 0.003623188 4.118666e-06  1.22842152  1.096444e-01
## [225,] 0.0030751708 0.003623188 2.988108e-06 -0.31702708  6.243885e-01
## [226,] 0.0056947608 0.003623188 4.368773e-06  0.99110659  1.608168e-01
## [227,] 0.0045766590 0.003623188 3.229769e-05  0.16777281  4.333810e-01
## [228,] 0.0053516820 0.003623188 1.940980e-06  1.24067261  1.073634e-01
## [229,] 0.0000000000 0.003623188 2.124358e-05 -0.78609873  7.840952e-01
## [230,] 0.0008542141 0.003623188 3.793496e-06 -1.42167103  9.224391e-01
## [231,] 0.0136130536 0.003623188 2.405068e-06  6.44163321  5.909732e-11
## [232,] 0.0022935780 0.003623188 7.929615e-06 -0.47216997  6.815973e-01
## [233,] 0.0123446630 0.003623188 1.966726e-06  6.21896335  2.502251e-10
## [234,] 0.0171277997 0.003623188 2.635034e-06  8.31934179  4.422657e-17
## [235,] 0.0013667426 0.003623188 4.061959e-06 -1.11958522  8.685547e-01
## [236,] 0.0049577918 0.003623188 3.556617e-06  0.70767470  2.395736e-01
## [237,] 0.0005694761 0.003623188 5.135810e-06 -1.34748412  9.110878e-01
## [238,] 0.0044608960 0.003623188 2.451182e-06  0.53506268  2.963033e-01
## [239,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
## [240,] 0.0013698630 0.003623188 1.271871e-05 -0.63183300  7.362520e-01
## [241,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
## [242,] 0.0013667426 0.003623188 4.061959e-06 -1.11958522  8.685547e-01
## [243,] 0.0124742380 0.003623188 2.834700e-06  5.25703906  7.319662e-08
## [244,] 0.0046269932 0.003623188 1.780025e-06  0.75237819  2.259118e-01
## [245,] 0.0028032620 0.003623188 7.022246e-06 -0.30941180  6.214959e-01
## [246,] 0.0009132420 0.003623188 1.271871e-05 -0.75986964  7.763337e-01
## [247,] 0.0005256702 0.003623188 4.722790e-06 -1.42532709  9.229687e-01
## [248,] 0.0005760369 0.003623188 3.870967e-06 -1.54876067  9.392804e-01
## [249,] 0.0026997385 0.003623188 2.152890e-06 -0.62936479  7.354449e-01
## [250,] 0.0139211137 0.003623188 1.628205e-05  2.55208567  5.354009e-03
## [251,] 0.0015290520 0.003623188 4.118666e-06 -1.03187408  8.489344e-01
## [252,] 0.0025310048 0.003623188 6.925562e-06 -0.41501916  6.609361e-01
## [253,] 0.0017717034 0.003623188 6.925562e-06 -0.70354633  7.591423e-01
## [254,] 0.0000000000 0.003623188 2.988108e-06 -2.09600737  9.819592e-01
## [255,] 0.0172641170 0.003623188 3.157663e-06  7.67646100  8.177191e-15
## [256,] 0.0000000000 0.003623188 8.970993e-06 -1.20968047  8.867992e-01
## [257,] 0.0000000000 0.003623188 1.265277e-05 -1.01858704  8.458005e-01
## [258,] 0.0130249867 0.003623188 3.284545e-06  5.18768135  1.064643e-07
## [259,] 0.0000000000 0.003623188 1.050507e-05 -1.11787060  8.681889e-01
## [260,] 0.0012424933 0.003623188 5.623924e-06 -1.00388523  8.422830e-01
## [261,] 0.0039151481 0.003623188 1.780025e-06  0.21883148  4.133907e-01
## [262,] 0.0066649793 0.003623188 2.152890e-06  2.07309136  1.908189e-02
## [263,] 0.0003796507 0.003623188 5.135810e-06 -1.43124664  9.238202e-01
## [264,] 0.0016566577 0.003623188 5.623924e-06 -0.82924147  7.965161e-01
## [265,] 0.0118721461 0.003623188 4.083130e-06  4.08227717  2.229828e-05
## [266,] 0.0012424933 0.003623188 5.623924e-06 -1.00388523  8.422830e-01
## [267,] 0.0000000000 0.003623188 4.061959e-06 -1.79772462  9.638897e-01
## [268,] 0.0061637411 0.003623188 3.556617e-06  1.34713048  8.896911e-02
## [269,] 0.0104403948 0.003623188 2.451182e-06  4.35430322  6.674542e-06
## [270,] 0.0037965072 0.003623188 6.925562e-06  0.06585946  4.737449e-01
## [271,] 0.0013667426 0.003623188 1.265277e-05 -0.63435466  7.370753e-01
## [272,] 0.0054937693 0.003623188 3.556617e-06  0.99187727  1.606287e-01
## [273,] 0.0115522291 0.003623188 2.067664e-06  5.51417690  1.752080e-08
## [274,] 0.0185185185 0.003623188 2.176473e-05  3.19281325  7.044702e-04
## [275,] 0.0097031963 0.003623188 3.813268e-06  3.11354752  9.242637e-04
## [276,] 0.0039179954 0.003623188 2.343797e-06  0.19256523  4.236497e-01
## [277,] 0.0030861930 0.003623188 1.844976e-06 -0.39534412  6.537055e-01
## attr(,"cluster")
##   [1] Low  Low  High High Low  High Low  Low  Low  High High Low  Low  High Low 
##  [16] Low  High Low  Low  Low  Low  Low  Low  Low  High High High Low  High Low 
##  [31] Low  Low  Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [46] Low  High High Low  High Low  Low  Low  High High High High Low  Low  Low 
##  [61] Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [76] Low  Low  Low  Low  Low  Low  High High Low  High Low  High Low  Low  Low 
##  [91] Low  Low  Low  Low  Low  High High Low  High Low  Low  High Low  High Low 
## [106] Low  High Low  High High Low  Low  High Low  High High Low  Low  Low  Low 
## [121] Low  Low  Low  High Low  High High Low  Low  Low  Low  Low  Low  Low  Low 
## [136] Low  High High Low  High Low  Low  High High High High Low  Low  Low  Low 
## [151] Low  Low  Low  Low  Low  Low  High High Low  High Low  Low  Low  High Low 
## [166] Low  Low  Low  Low  Low  Low  High High Low  Low  High Low  Low  Low  Low 
## [181] High Low  Low  Low  Low  High Low  High High Low  Low  Low  Low  High Low 
## [196] Low  Low  Low  High Low  Low  Low  High Low  High Low  Low  High Low  Low 
## [211] Low  Low  Low  Low  High Low  Low  Low  Low  High Low  High High High Low 
## [226] Low  High High Low  Low  High High High High Low  Low  Low  Low  Low  Low 
## [241] Low  Low  Low  Low  High Low  Low  High Low  High High Low  Low  Low  Low 
## [256] Low  Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
## [271] Low  Low  Low  High Low  Low  Low 
## Levels: Low High
## attr(,"gstari")
## [1] FALSE
## attr(,"call")
## localG(x = data_merge$KEK, listw = W.rdw, alternative = "greater")
## attr(,"class")
## [1] "localG"

data_merge$lmZ2 <- Local_M
tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ2",
title = "Local Gi", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Coldspot", "No SAC", "Hotspot"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Indeks Moran Global

W.idw <- mat2listw(w.idw, style='W')
moran(data_merge$KEK, W.idw, n=length(W.idw$neighbours), S0=Szero(W.idw))

## $I
## [1] 0.1130384
## 
## $K
## [1] 18.9482

IM <- moran.test(data_merge$KEK, W.idw,randomisation=T,
alternative="greater", zero.policy = FALSE)
IM

## 
##  Moran I test under randomisation
## 
## data:  data_merge$KEK  
## weights: W.idw    
## 
## Moran I statistic standard deviate = 19.51, p-value < 2.2e-16
## alternative hypothesis: greater
## sample estimates:
## Moran I statistic       Expectation          Variance 
##      1.130384e-01     -3.623188e-03      3.575541e-05

Indeks Moran Lokal

LM <- localmoran(data_merge$KEK, W.idw, alternative = "greater")
head(LM)

##            Ii          E.Ii       Var.Ii       Z.Ii Pr(z > E(Ii))
## 1  0.01947620 -7.376709e-04 2.496782e-04  1.2792607     0.1004026
## 2  0.01712903 -7.376709e-04 4.332487e-04  0.8583723     0.1953435
## 3  0.03770718 -8.612240e-03 8.644418e-03  0.4981903     0.3091750
## 4 -0.01523860 -5.062083e-05 4.568421e-05 -2.2470710     0.9876823
## 5  0.01424062 -7.376709e-04 1.565552e-03  0.3785546     0.3525093
## 6 -0.05410789 -3.425436e-03 2.881584e-03 -0.9441523     0.8274541

data_merge$lmI <- LM[, "Ii"] # local Moran's I
data_merge$lmZ <- LM[, "Z.Ii"] # z-scores
# p-values corresponding to alternative greater
data_merge$lmp <- LM[, "Pr(z > E(Ii))"]

p1 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p1

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

p2 <- tm_shape(data_merge) +
  tm_polygons(col = "lmp", title = "p-value",
              breaks = c(-Inf, 0.05, Inf)) +
  tm_layout(legend.outside = TRUE)
p2

Klasifikasi Warna :

• Less than 0.05 = Ada Autokorelasi Positif

• 0.05 or more = Tidak Ada Autokorelasi Positif

tmap_mode("view")

## tmap mode set to interactive viewing

p3 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p3

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ",
title = "Local Moran's I", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Negative SAC", "No SAC", "Positive SAC"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Rentang nilai untuk klasifikasi warna, dengan:

• -Inf hingga -1.96 mewakili Negative Spatial Autocorrelation (SAC).

• -1.96 hingga 1.96 mewakili No Spatial Autocorrelation (No SAC).

• 1.96 hingga Inf mewakili Positive Spatial Autocorrelation (Positive SAC)

# Autokorelasi spasial positif (Z-score > 1.96)
positive_autocorr <- data_merge[data_merge$lmZ > 1.96, ]
# Autokorelasi spasial negatif (Z-score < -1.96)
negative_autocorr <- data_merge[data_merge$lmZ < -1.96, ]
# Area tanpa autokorelasi spasial (Z-score antara -1.96 dan 1.96)
no_autocorr <- data_merge[data_merge$lmZ >= -1.96 & data_merge$lmZ <= 1.96, ]

moran.plot(data_merge$KEK, W.idw)

Moran Plot memiliki kemiringan positif: ini menunjukkan adanya autokorelasi spasial positif, di mana lokasi dengan nilai tinggi cenderung berdekatan dengan lokasi bernilai tinggi, dan lokasi bernilai rendah cenderung berdekatan dengan lokasi bernilai rendah.

Getis-Ord Gi

Mendeteksi adanya wilayah yang hotspot/coldspot.

Local_M <- localG(data_merge$KEK, W.idw, alternative = "greater")
Local_M

##   [1] -1.279260750 -0.858372256  0.498190301 -2.247071035 -0.378554578
##   [6] -0.944152310 -0.445399843 -1.335026152 -1.342082350  0.744299352
##  [11]  1.519684712 -1.128685757 -1.965131204  0.037175336 -1.874462497
##  [16] -1.538558714  0.069012474 -2.334062306 -1.337771804 -0.966991992
##  [21] -1.742218298  0.706155517 -0.921037316 -1.399989894 -0.380411424
##  [26]  9.205054744  2.840505309  2.967987317  8.501967106  1.426663219
##  [31] -0.863944428 -1.079396956  0.765370579 -1.036157166 -2.147189882
##  [36]  0.823439153 -2.567232844  1.058102689 -1.821734140  3.001472729
##  [41] -2.718289190  6.491453116 -2.208040269 -1.782312906 -1.061195999
##  [46] -2.019429092  0.952501441  0.130633393 -0.438156149  0.298404916
##  [51] -1.264165546 -0.077940618 -1.013052073  7.563328764  7.376048210
##  [56]  4.340330375  0.908895467  4.218090193 -1.612076347  2.146202888
##  [61] -0.046290344  0.777349745 -1.572477291 -0.935125932 -1.455709849
##  [66] -2.717576217 -2.421762524 -0.227920935 -1.564158399  7.503615227
##  [71] -0.883901260 -2.794685626 -1.697679432 -1.424297190  4.223952204
##  [76] -0.948733146 -0.927777475 -2.415343483 -1.637852982 -2.446478769
##  [81]  0.062525874  0.402563819  0.875134685 -1.723207590 -1.617785706
##  [86] -1.460493985  2.644620600 -0.431214470  4.353153538 -2.611782713
##  [91]  1.690457552 -0.572639009 -1.583779210 -2.326321381 -2.122455136
##  [96] -0.232324341  1.093174033 -1.887317209  0.394542616 -1.026138031
## [101] -0.358439653 -0.032565620 -1.065765762  0.386152326 -2.825179785
## [106]  1.994676022  7.308647338  0.923058127  5.194217116  1.505947645
## [111]  7.316738091 -1.164392532  3.265044507  0.274872961  0.248745542
## [116]  8.346874793 -0.716339422 -0.583566787 -0.279381031  2.883450210
## [121] -1.727569901  0.723633731 -1.860799589  5.377814359  3.423276668
## [126]  0.798427592 -0.521400766 -1.573348281 -2.524270275  2.844982985
## [131] -0.935547517  1.178701686  2.201451944  0.747524521  1.045744453
## [136] -2.041285849  2.279086671 -1.090559227 -1.083463940 -1.236553605
## [141] -1.245802448 -1.006401495 -0.040369109 -1.235646624  2.477926949
## [146]  5.830820571  0.333748091 -2.721119460 -1.901408034 -1.414094216
## [151] -0.963806457 -2.280089104  2.510614483  2.560679019 -2.100017175
## [156] -1.959980297  0.424675359  0.721463859 -1.354323673 -1.346788995
## [161] -1.729479575 -1.274292375 -2.253575001  0.466993983 -1.391607981
## [166] -1.443639958 -1.778009201  0.282096124  0.711009506 -1.209717290
## [171] -1.189612718  0.213381750  0.001225376 -1.423699812 -1.331842300
## [176]  0.873461759 -2.024547054 -1.369017508 -1.248919071 -1.013682217
## [181]  0.774546014 -0.727015670 -1.250525301 -1.512589045 -1.860207162
## [186]  0.302798270 -1.300504668 -1.362469531  0.362145699 -1.615166780
## [191] -1.183334504 -1.420208092 -0.277278866  0.019728424 -1.878159660
## [196] -2.622440323  2.589770513 -2.424945352  1.166575793  1.719192732
## [201]  1.026557434  2.981694299 -0.728563683 -0.282120734 -1.710431603
## [206] -0.829470762 -1.717934391  0.903327863 -1.922285721 -1.737278534
## [211]  4.433874661  3.860209942 -1.253096027 -1.364812074  0.134248818
## [216]  2.887265801 -0.437747491 -1.787219688 -2.437041889  0.976196432
## [221]  1.845241352  0.809227596  9.471827139  0.042493479  0.794417351
## [226] -0.264011037  0.790693397  0.384532678  0.881116028 -0.945670818
## [231]  8.267264481 -0.255133205  5.699069175  8.860900586 -1.010803911
## [236] -0.305863704 -0.657598383  1.315675988 -2.596691945 -0.141439232
## [241] -2.439236336 -1.506585133  5.363868209 -0.125247414  0.206311776
## [246] -0.192328123 -0.969068531 -1.701622822 -1.431436298  1.710955099
## [251] -1.569744480 -0.133112133 -0.324663931 -1.776417488  6.932767785
## [256] -2.620057522 -2.405378751  4.826071805 -2.887894769 -1.084841748
## [261] -0.952133658  1.779765160 -2.213717410 -0.731161252  3.470935484
## [266] -0.680944973 -1.317017057 -0.479413326  3.914262736  0.515654321
## [271] -1.504825386  3.363270356  4.695874383  2.613428193  3.254599682
## [276] -1.102239827 -1.637584630
## attr(,"internals")
##                 Gi       E(Gi)        V(Gi)        Z(Gi) Pr(z > E(Gi))
##   [1,] 0.003264764 0.003623188 7.850146e-08 -1.279260750  8.995974e-01
##   [2,] 0.003306383 0.003623188 1.362179e-07 -0.858372256  8.046565e-01
##   [3,] 0.003867456 0.003623188 2.404035e-07  0.498190301  3.091750e-01
##   [4,] 0.002590431 0.003623188 2.112334e-07 -2.247071035  9.876823e-01
##   [5,] 0.003357599 0.003623188 4.922259e-07 -0.378554578  6.474907e-01
##   [6,] 0.003201342 0.003623188 1.996294e-07 -0.944152310  8.274541e-01
##   [7,] 0.003502011 0.003623188 7.401870e-08 -0.445399843  6.719846e-01
##   [8,] 0.002819573 0.003623188 3.623409e-07 -1.335026152  9.090661e-01
##   [9,] 0.002857655 0.003623188 3.253645e-07 -1.342082350  9.102154e-01
##  [10,] 0.003973958 0.003623188 2.221002e-07  0.744299352  2.283477e-01
##  [11,] 0.004372810 0.003623188 2.433198e-07  1.519684712  6.429512e-02
##  [12,] 0.003079691 0.003623188 2.318725e-07 -1.128685757  8.704848e-01
##  [13,] 0.002916060 0.003623188 1.294831e-07 -1.965131204  9.753005e-01
##  [14,] 0.003644947 0.003623188 3.425871e-07  0.037175336  4.851726e-01
##  [15,] 0.002851907 0.003623188 1.693058e-07 -1.874462497  9.695666e-01
##  [16,] 0.003039477 0.003623188 1.439356e-07 -1.538558714  9.380440e-01
##  [17,] 0.003656450 0.003623188 2.322919e-07  0.069012474  4.724898e-01
##  [18,] 0.002519735 0.003623188 2.235028e-07 -2.334062306  9.902038e-01
##  [19,] 0.003019504 0.003623188 2.036368e-07 -1.337771804  9.095146e-01
##  [20,] 0.002959359 0.003623188 4.712679e-07 -0.966991992  8.332260e-01
##  [21,] 0.003035170 0.003623188 1.139140e-07 -1.742218298  9.592649e-01
##  [22,] 0.003919132 0.003623188 1.756373e-07  0.706155517  2.400457e-01
##  [23,] 0.002574841 0.003623188 1.295556e-06 -0.921037316  8.214845e-01
##  [24,] 0.002563427 0.003623188 5.730152e-07 -1.399989894  9.192418e-01
##  [25,] 0.003411023 0.003623188 3.110591e-07 -0.380411424  6.481800e-01
##  [26,] 0.010441002 0.003623188 5.485768e-07  9.205054744  1.707486e-20
##  [27,] 0.005418063 0.003623188 3.992795e-07  2.840505309  2.252106e-03
##  [28,] 0.006093493 0.003623188 6.927506e-07  2.967987317  1.498783e-03
##  [29,] 0.009240104 0.003623188 4.364726e-07  8.501967106  9.320232e-18
##  [30,] 0.004180027 0.003623188 1.523406e-07  1.426663219  7.683850e-02
##  [31,] 0.002824102 0.003623188 8.554923e-07 -0.863944428  8.061908e-01
##  [32,] 0.002647619 0.003623188 8.168715e-07 -1.079396956  8.597946e-01
##  [33,] 0.003903557 0.003623188 1.341888e-07  0.765370579  2.220255e-01
##  [34,] 0.003001162 0.003623188 3.603843e-07 -1.036157166  8.499356e-01
##  [35,] 0.002871545 0.003623188 1.225414e-07 -2.147189882  9.841109e-01
##  [36,] 0.004239391 0.003623188 5.599944e-07  0.823439153  2.051292e-01
##  [37,] 0.002599038 0.003623188 1.591465e-07 -2.567232844  9.948743e-01
##  [38,] 0.004044242 0.003623188 1.583506e-07  1.058102689  1.450043e-01
##  [39,] 0.002987172 0.003623188 1.218895e-07 -1.821734140  9.657523e-01
##  [40,] 0.004509197 0.003623188 8.713795e-08  3.001472729  1.343386e-03
##  [41,] 0.002461073 0.003623188 1.827711e-07 -2.718289190  9.967190e-01
##  [42,] 0.008278671 0.003623188 5.143343e-07  6.491453116  4.250618e-11
##  [43,] 0.002598469 0.003623188 2.153753e-07 -2.208040269  9.863793e-01
##  [44,] 0.002597499 0.003623188 3.311796e-07 -1.782312906  9.626509e-01
##  [45,] 0.002678792 0.003623188 7.919852e-07 -1.061195999  8.556996e-01
##  [46,] 0.002797767 0.003623188 1.670686e-07 -2.019429092  9.782787e-01
##  [47,] 0.003994882 0.003623188 1.522790e-07  0.952501441  1.704214e-01
##  [48,] 0.003675902 0.003623188 1.628328e-07  0.130633393  4.480327e-01
##  [49,] 0.003417009 0.003623188 2.214280e-07 -0.438156149  6.693635e-01
##  [50,] 0.003822122 0.003623188 4.444294e-07  0.298404916  3.826971e-01
##  [51,] 0.002841633 0.003623188 3.822182e-07 -1.264165546  8.969147e-01
##  [52,] 0.003576076 0.003623188 3.653761e-07 -0.077940618  5.310624e-01
##  [53,] 0.002658496 0.003623188 9.068053e-07 -1.013052073  8.444824e-01
##  [54,] 0.008688315 0.003623188 4.484919e-07  7.563328764  1.964415e-14
##  [55,] 0.010138411 0.003623188 7.802091e-07  7.376048210  8.152850e-14
##  [56,] 0.006161218 0.003623188 3.419384e-07  4.340330375  7.113433e-06
##  [57,] 0.004149751 0.003623188 3.356389e-07  0.908895467  1.817027e-01
##  [58,] 0.005783310 0.003623188 2.622555e-07  4.218090193  1.231901e-05
##  [59,] 0.002881154 0.003623188 2.118736e-07 -1.612076347  9.465273e-01
##  [60,] 0.004752518 0.003623188 2.768855e-07  2.146202888  1.592840e-02
##  [61,] 0.003586045 0.003623188 6.438632e-07 -0.046290344  5.184606e-01
##  [62,] 0.004131960 0.003623188 4.283633e-07  0.777349745  2.184762e-01
##  [63,] 0.002890784 0.003623188 2.169367e-07 -1.572477291  9.420800e-01
##  [64,] 0.003245409 0.003623188 1.632064e-07 -0.935125932  8.251383e-01
##  [65,] 0.003032949 0.003623188 1.644020e-07 -1.455709849  9.272636e-01
##  [66,] 0.002564804 0.003623188 1.516782e-07 -2.717576217  9.967119e-01
##  [67,] 0.002560725 0.003623188 1.924709e-07 -2.421762524  9.922773e-01
##  [68,] 0.003544521 0.003623188 1.191305e-07 -0.227920935  5.901461e-01
##  [69,] 0.002832068 0.003623188 2.558134e-07 -1.564158399  9.411098e-01
##  [70,] 0.008386100 0.003623188 4.029062e-07  7.503615227  3.104068e-14
##  [71,] 0.002659536 0.003623188 1.188593e-06 -0.883901260  8.116252e-01
##  [72,] 0.002489559 0.003623188 1.645418e-07 -2.794685626  9.974025e-01
##  [73,] 0.002858738 0.003623188 2.027622e-07 -1.697679432  9.552159e-01
##  [74,] 0.003127503 0.003623188 1.211189e-07 -1.424297190  9.228198e-01
##  [75,] 0.006367503 0.003623188 4.221142e-07  4.223952204  1.200275e-05
##  [76,] 0.002835674 0.003623188 6.890160e-07 -0.948733146  8.286218e-01
##  [77,] 0.002559629 0.003623188 1.314122e-06 -0.927777475  8.232385e-01
##  [78,] 0.002480736 0.003623188 2.237269e-07 -2.415343483  9.921398e-01
##  [79,] 0.002729858 0.003623188 2.974913e-07 -1.637852982  9.492738e-01
##  [80,] 0.002459097 0.003623188 2.264076e-07 -2.446478769  9.927870e-01
##  [81,] 0.003657731 0.003623188 3.051974e-07  0.062525874  4.750720e-01
##  [82,] 0.003889857 0.003623188 4.388085e-07  0.402563819  3.436346e-01
##  [83,] 0.003886244 0.003623188 9.035323e-08  0.875134685  1.907503e-01
##  [84,] 0.003077487 0.003623188 1.002848e-07 -1.723207590  9.575745e-01
##  [85,] 0.002492390 0.003623188 4.885722e-07 -1.617785706  9.471456e-01
##  [86,] 0.002917777 0.003623188 2.332845e-07 -1.460493985  9.279228e-01
##  [87,] 0.005891550 0.003623188 7.356948e-07  2.644620600  4.089128e-03
##  [88,] 0.003434088 0.003623188 1.923082e-07 -0.431214470  6.668438e-01
##  [89,] 0.005852043 0.003623188 2.621535e-07  4.353153538  6.709655e-06
##  [90,] 0.002482404 0.003623188 1.907800e-07 -2.611782713  9.954964e-01
##  [91,] 0.004371851 0.003623188 1.961390e-07  1.690457552  4.547023e-02
##  [92,] 0.003417722 0.003623188 1.287424e-07 -0.572639009  7.165554e-01
##  [93,] 0.003067253 0.003623188 1.232137e-07 -1.583779210  9.433780e-01
##  [94,] 0.002818507 0.003623188 1.196488e-07 -2.326321381  9.899993e-01
##  [95,] 0.002899419 0.003623188 1.162848e-07 -2.122455136  9.831002e-01
##  [96,] 0.003500267 0.003623188 2.799382e-07 -0.232324341  5.918569e-01
##  [97,] 0.004259077 0.003623188 3.383640e-07  1.093174033  1.371587e-01
##  [98,] 0.002635378 0.003623188 2.739411e-07 -1.887317209  9.704412e-01
##  [99,] 0.003904238 0.003623188 5.074329e-07  0.394542616  3.465902e-01
## [100,] 0.003007799 0.003623188 3.596569e-07 -1.026138031  8.475867e-01
## [101,] 0.003446852 0.003623188 2.420191e-07 -0.358439653  6.399928e-01
## [102,] 0.003600062 0.003623188 5.043105e-07 -0.032565620  5.129895e-01
## [103,] 0.003181574 0.003623188 1.716970e-07 -1.065765762  8.567352e-01
## [104,] 0.003924696 0.003623188 6.096470e-07  0.386152326  3.496919e-01
## [105,] 0.002465215 0.003623188 1.679984e-07 -2.825179785  9.976373e-01
## [106,] 0.004445547 0.003623188 1.699723e-07  1.994676022  2.303911e-02
## [107,] 0.010810854 0.003623188 9.671676e-07  7.308647338  1.349225e-13
## [108,] 0.004242354 0.003623188 4.499414e-07  0.923058127  1.779885e-01
## [109,] 0.006693784 0.003623188 3.494662e-07  5.194217116  1.027916e-07
## [110,] 0.004283269 0.003623188 1.921210e-07  1.505947645  6.604031e-02
## [111,] 0.008971956 0.003623188 5.344074e-07  7.316738091  1.270355e-13
## [112,] 0.002946040 0.003623188 3.381966e-07 -1.164392532  8.778675e-01
## [113,] 0.005698876 0.003623188 4.041530e-07  3.265044507  5.472344e-04
## [114,] 0.003736806 0.003623188 1.708547e-07  0.274872961  3.917069e-01
## [115,] 0.003809600 0.003623188 5.616132e-07  0.248745542  4.017788e-01
## [116,] 0.008526883 0.003623188 3.451429e-07  8.346874793  3.504648e-17
## [117,] 0.003423280 0.003623188 7.787960e-08 -0.716339422  7.631091e-01
## [118,] 0.003433748 0.003623188 1.053817e-07 -0.583566787  7.202441e-01
## [119,] 0.003487160 0.003623188 2.370640e-07 -0.279381031  6.100238e-01
## [120,] 0.006047363 0.003623188 7.068103e-07  2.883450210  1.966724e-03
## [121,] 0.002588193 0.003623188 3.589266e-07 -1.727569901  9.579673e-01
## [122,] 0.003902975 0.003623188 1.494915e-07  0.723633731  2.346453e-01
## [123,] 0.002976344 0.003623188 1.208371e-07 -1.860799589  9.686138e-01
## [124,] 0.006885111 0.003623188 3.679050e-07  5.377814359  3.769775e-08
## [125,] 0.006626706 0.003623188 7.697972e-07  3.423276668  3.093553e-04
## [126,] 0.004180602 0.003623188 4.873989e-07  0.798427592  2.123112e-01
## [127,] 0.003380609 0.003623188 2.164535e-07 -0.521400766  6.989562e-01
## [128,] 0.002782685 0.003623188 2.853835e-07 -1.573348281  9.421809e-01
## [129,] 0.002543583 0.003623188 1.829187e-07 -2.524270275  9.942031e-01
## [130,] 0.006142012 0.003623188 7.838559e-07  2.844982985  2.220690e-03
## [131,] 0.003256724 0.003623188 1.534373e-07 -0.935547517  8.252469e-01
## [132,] 0.004278536 0.003623188 3.091264e-07  1.178701686  1.192585e-01
## [133,] 0.004630709 0.003623188 2.094544e-07  2.201451944  1.385202e-02
## [134,] 0.003991811 0.003623188 2.431715e-07  0.747524521  2.273735e-01
## [135,] 0.004054003 0.003623188 1.697191e-07  1.045744453  1.478395e-01
## [136,] 0.002929487 0.003623188 1.154882e-07 -2.041285849  9.793888e-01
## [137,] 0.005085656 0.003623188 4.117661e-07  2.279086671  1.133096e-02
## [138,] 0.003280608 0.003623188 9.867947e-08 -1.090559227  8.622666e-01
## [139,] 0.002738829 0.003623188 6.662367e-07 -1.083463940  8.606987e-01
## [140,] 0.003069013 0.003623188 2.008484e-07 -1.236553605  8.918736e-01
## [141,] 0.002789801 0.003623188 4.475027e-07 -1.245802448  8.935815e-01
## [142,] 0.002621928 0.003623188 9.898090e-07 -1.006401495  8.428888e-01
## [143,] 0.003610984 0.003623188 9.140093e-08 -0.040369109  5.161006e-01
## [144,] 0.002633268 0.003623188 6.418184e-07 -1.235646624  8.917050e-01
## [145,] 0.005372467 0.003623188 4.983576e-07  2.477926949  6.607410e-03
## [146,] 0.008764351 0.003623188 7.774335e-07  5.830820571  2.757773e-09
## [147,] 0.003851769 0.003623188 4.690725e-07  0.333748091  3.692848e-01
## [148,] 0.002478660 0.003623188 1.769123e-07 -2.721119460  9.967469e-01
## [149,] 0.002772346 0.003623188 2.002387e-07 -1.901408034  9.713757e-01
## [150,] 0.003240201 0.003623188 7.335217e-08 -1.414094216  9.213329e-01
## [151,] 0.003119865 0.003623188 2.727183e-07 -0.963806457  8.324285e-01
## [152,] 0.002650207 0.003623188 1.820980e-07 -2.280089104  9.886988e-01
## [153,] 0.005188523 0.003623188 3.887355e-07  2.510614483  6.026062e-03
## [154,] 0.004819421 0.003623188 2.182332e-07  2.560679019  5.223391e-03
## [155,] 0.002499229 0.003623188 2.864543e-07 -2.100017175  9.821363e-01
## [156,] 0.003029095 0.003623188 9.187681e-08 -1.959980297  9.750010e-01
## [157,] 0.004032546 0.003623188 9.291630e-07  0.424675359  3.355367e-01
## [158,] 0.004141735 0.003623188 5.165909e-07  0.721463859  2.353121e-01
## [159,] 0.002935436 0.003623188 2.578808e-07 -1.354323673  9.121834e-01
## [160,] 0.003000349 0.003623188 2.138713e-07 -1.346788995  9.109759e-01
## [161,] 0.003080059 0.003623188 9.862255e-08 -1.729479575  9.581384e-01
## [162,] 0.002760392 0.003623188 4.584359e-07 -1.274292375  8.987201e-01
## [163,] 0.002552872 0.003623188 2.255694e-07 -2.253575001  9.878885e-01
## [164,] 0.004027349 0.003623188 7.490064e-07  0.466993983  3.202521e-01
## [165,] 0.002820124 0.003623188 3.330170e-07 -1.391607981  9.179794e-01
## [166,] 0.002868198 0.003623188 2.735052e-07 -1.443639958  9.255799e-01
## [167,] 0.002722616 0.003623188 2.565484e-07 -1.778009201  9.622988e-01
## [168,] 0.003830418 0.003623188 5.396481e-07  0.282096124  3.889349e-01
## [169,] 0.004088944 0.003623188 4.291071e-07  0.711009506  2.385392e-01
## [170,] 0.002980895 0.003623188 2.819027e-07 -1.209717290  8.868063e-01
## [171,] 0.002879976 0.003623188 3.903145e-07 -1.189612718  8.829007e-01
## [172,] 0.003812083 0.003623188 7.836555e-07  0.213381750  4.155146e-01
## [173,] 0.003623863 0.003623188 3.034490e-07  0.001225376  4.995111e-01
## [174,] 0.002685041 0.003623188 4.342152e-07 -1.423699812  9.227333e-01
## [175,] 0.002678632 0.003623188 5.029796e-07 -1.331842300  9.085440e-01
## [176,] 0.004321061 0.003623188 6.383588e-07  0.873461759  1.912057e-01
## [177,] 0.002928879 0.003623188 1.176116e-07 -2.024547054  9.785431e-01
## [178,] 0.002784719 0.003623188 3.751085e-07 -1.369017508  9.145031e-01
## [179,] 0.002735997 0.003623188 5.046221e-07 -1.248919071  8.941527e-01
## [180,] 0.002551568 0.003623188 1.117578e-06 -1.013682217  8.446328e-01
## [181,] 0.004004266 0.003623188 2.420653e-07  0.774546014  2.193040e-01
## [182,] 0.003104951 0.003623188 5.081244e-07 -0.727015670  7.663918e-01
## [183,] 0.002827487 0.003623188 4.048695e-07 -1.250525301  8.944461e-01
## [184,] 0.002878501 0.003623188 2.423852e-07 -1.512589045  9.348080e-01
## [185,] 0.002879933 0.003623188 1.596443e-07 -1.860207162  9.685719e-01
## [186,] 0.003822004 0.003623188 4.311165e-07  0.302798270  3.810218e-01
## [187,] 0.003122572 0.003623188 1.481790e-07 -1.300504668  9.032860e-01
## [188,] 0.002953844 0.003623188 2.413494e-07 -1.362469531  9.134751e-01
## [189,] 0.003906455 0.003623188 6.118196e-07  0.362145699  3.586216e-01
## [190,] 0.002839555 0.003623188 2.353918e-07 -1.615166780  9.468627e-01
## [191,] 0.003036936 0.003623188 2.454447e-07 -1.183334504  8.816617e-01
## [192,] 0.002708653 0.003623188 4.146650e-07 -1.420208092  9.222264e-01
## [193,] 0.003449945 0.003623188 3.903731e-07 -0.277278866  6.092170e-01
## [194,] 0.003633139 0.003623188 2.543853e-07  0.019728424  4.921300e-01
## [195,] 0.002795802 0.003623188 1.940670e-07 -1.878159660  9.698203e-01
## [196,] 0.002534022 0.003623188 1.724953e-07 -2.622440323  9.956349e-01
## [197,] 0.004871319 0.003623188 2.322724e-07  2.589770513  4.801997e-03
## [198,] 0.002471014 0.003623188 2.257525e-07 -2.424945352  9.923447e-01
## [199,] 0.004035460 0.003623188 1.248936e-07  1.166575793  1.216909e-01
## [200,] 0.004346188 0.003623188 1.768590e-07  1.719192732  4.278964e-02
## [201,] 0.004016112 0.003623188 1.465038e-07  1.026557434  1.523145e-01
## [202,] 0.006142914 0.003623188 7.141351e-07  2.981694299  1.433290e-03
## [203,] 0.003306609 0.003623188 1.888124e-07 -0.728563683  7.668657e-01
## [204,] 0.003470225 0.003623188 2.939731e-07 -0.282120734  6.110745e-01
## [205,] 0.002441968 0.003623188 4.769257e-07 -1.710431603  9.564070e-01
## [206,] 0.002865108 0.003623188 8.352727e-07 -0.829470762  7.965810e-01
## [207,] 0.003058239 0.003623188 1.081448e-07 -1.717934391  9.570957e-01
## [208,] 0.004314858 0.003623188 5.862823e-07  0.903327863  1.831760e-01
## [209,] 0.002863741 0.003623188 1.560844e-07 -1.922285721  9.727151e-01
## [210,] 0.002631928 0.003623188 3.255640e-07 -1.737278534  9.588310e-01
## [211,] 0.006318310 0.003623188 3.694791e-07  4.433874661  4.627722e-06
## [212,] 0.006342688 0.003623188 4.963135e-07  3.860209942  5.664482e-05
## [213,] 0.003275535 0.003623188 7.697029e-08 -1.253096027  8.949146e-01
## [214,] 0.002549829 0.003623188 6.185069e-07 -1.364812074  9.138439e-01
## [215,] 0.003713394 0.003623188 4.514894e-07  0.134248818  4.466029e-01
## [216,] 0.004986399 0.003623188 2.229218e-07  2.887265801  1.943029e-03
## [217,] 0.003415077 0.003623188 2.260179e-07 -0.437747491  6.692153e-01
## [218,] 0.002708191 0.003623188 2.621102e-07 -1.787219688  9.630490e-01
## [219,] 0.002869224 0.003623188 9.571408e-08 -2.437041889  9.925960e-01
## [220,] 0.004342868 0.003623188 5.435049e-07  0.976196432  1.644836e-01
## [221,] 0.004440111 0.003623188 1.959999e-07  1.845241352  3.250122e-02
## [222,] 0.003936871 0.003623188 1.502584e-07  0.809227596  2.091921e-01
## [223,] 0.012910371 0.003623188 9.613919e-07  9.471827139  1.374954e-21
## [224,] 0.003641590 0.003623188 1.875334e-07  0.042493479  4.830527e-01
## [225,] 0.004120772 0.003623188 3.923146e-07  0.794417351  2.134763e-01
## [226,] 0.003506915 0.003623188 1.939607e-07 -0.264011037  6.041143e-01
## [227,] 0.003871336 0.003623188 9.849267e-08  0.790693397  2.145615e-01
## [228,] 0.003894381 0.003623188 4.973795e-07  0.384532678  3.502918e-01
## [229,] 0.003981396 0.003623188 1.652740e-07  0.881116028  1.891275e-01
## [230,] 0.002524974 0.003623188 1.348633e-06 -0.945670818  8.278417e-01
## [231,] 0.011832109 0.003623188 9.859354e-07  8.267264481  6.853383e-17
## [232,] 0.003518406 0.003623188 1.686703e-07 -0.255133205  6.006899e-01
## [233,] 0.006760110 0.003623188 3.029699e-07  5.699069175  6.023167e-09
## [234,] 0.010347322 0.003623188 5.758601e-07  8.860900586  3.968531e-19
## [235,] 0.002834449 0.003623188 6.088824e-07 -1.010803911  8.439449e-01
## [236,] 0.003492777 0.003623188 1.817930e-07 -0.305863704  6.201458e-01
## [237,] 0.003080286 0.003623188 6.815884e-07 -0.657598383  7.446019e-01
## [238,] 0.004364730 0.003623188 3.176675e-07  1.315675988  9.414141e-02
## [239,] 0.002709057 0.003623188 1.239299e-07 -2.596691945  9.952937e-01
## [240,] 0.003564105 0.003623188 1.744985e-07 -0.141439232  5.562385e-01
## [241,] 0.002493004 0.003623188 2.146797e-07 -2.439236336  9.926408e-01
## [242,] 0.002716132 0.003623188 3.624778e-07 -1.506585133  9.340415e-01
## [243,] 0.006801488 0.003623188 3.511021e-07  5.363868209  4.072917e-08
## [244,] 0.003538600 0.003623188 4.561305e-07 -0.125247414  5.498362e-01
## [245,] 0.003781286 0.003623188 5.872236e-07  0.206311776  4.182737e-01
## [246,] 0.003459695 0.003623188 7.226324e-07 -0.192328123  5.762574e-01
## [247,] 0.002592740 0.003623188 1.130691e-06 -0.969068531  8.337445e-01
## [248,] 0.002425577 0.003623188 4.953421e-07 -1.701622822  9.555870e-01
## [249,] 0.002652663 0.003623188 4.596949e-07 -1.431436298  9.238474e-01
## [250,] 0.004413045 0.003623188 2.131175e-07  1.710955099  4.354470e-02
## [251,] 0.002358037 0.003623188 6.495716e-07 -1.569744480  9.417627e-01
## [252,] 0.003562193 0.003623188 2.099740e-07 -0.133112133  5.529476e-01
## [253,] 0.003481578 0.003623188 1.902499e-07 -0.324663931  6.272823e-01
## [254,] 0.002727134 0.003623188 2.544360e-07 -1.776417488  9.621679e-01
## [255,] 0.007509962 0.003623188 3.143150e-07  6.932767785  2.063427e-12
## [256,] 0.002429681 0.003623188 2.075051e-07 -2.620057522  9.956043e-01
## [257,] 0.002576202 0.003623188 1.894590e-07 -2.405378751  9.919221e-01
## [258,] 0.006314745 0.003623188 3.110424e-07  4.826071805  6.962616e-07
## [259,] 0.002440193 0.003623188 1.678043e-07 -2.887894769  9.980609e-01
## [260,] 0.003143552 0.003623188 1.954752e-07 -1.084841748  8.610041e-01
## [261,] 0.003102651 0.003623188 2.988873e-07 -0.952133658  8.294854e-01
## [262,] 0.004656808 0.003623188 3.372843e-07  1.779765160  3.755720e-02
## [263,] 0.002464975 0.003623188 2.737365e-07 -2.213717410  9.865759e-01
## [264,] 0.003380033 0.003623188 1.105964e-07 -0.731161252  7.676597e-01
## [265,] 0.006101388 0.003623188 5.097760e-07  3.470935484  2.593243e-04
## [266,] 0.003302020 0.003623188 2.224542e-07 -0.680944973  7.520468e-01
## [267,] 0.003108669 0.003623188 1.526235e-07 -1.317017057  9.060835e-01
## [268,] 0.003428315 0.003623188 1.652292e-07 -0.479413326  6.841777e-01
## [269,] 0.005647746 0.003623188 2.675225e-07  3.914262736  4.534043e-05
## [270,] 0.003839124 0.003623188 1.753595e-07  0.515654321  3.030479e-01
## [271,] 0.003156709 0.003623188 9.609337e-08 -1.504825386  9.338155e-01
## [272,] 0.005303071 0.003623188 2.494793e-07  3.363270356  3.851245e-04
## [273,] 0.006092536 0.003623188 2.765231e-07  4.695874383  1.327343e-06
## [274,] 0.005146711 0.003623188 3.398419e-07  2.613428193  4.481946e-03
## [275,] 0.005862888 0.003623188 4.735702e-07  3.254599682  5.677615e-04
## [276,] 0.003027790 0.003623188 2.917855e-07 -1.102239827  8.648213e-01
## [277,] 0.002872439 0.003623188 2.101757e-07 -1.637584630  9.492458e-01
## attr(,"cluster")
##   [1] Low  Low  High High Low  High Low  Low  Low  High High Low  Low  High Low 
##  [16] Low  High Low  Low  Low  Low  Low  Low  Low  High High High Low  High Low 
##  [31] Low  Low  Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [46] Low  High High Low  High Low  Low  Low  High High High High Low  Low  Low 
##  [61] Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [76] Low  Low  Low  Low  Low  Low  High High Low  High Low  High Low  Low  Low 
##  [91] Low  Low  Low  Low  Low  High High Low  High Low  Low  High Low  High Low 
## [106] Low  High Low  High High Low  Low  High Low  High High Low  Low  Low  Low 
## [121] Low  Low  Low  High Low  High High Low  Low  Low  Low  Low  Low  Low  Low 
## [136] Low  High High Low  High Low  Low  High High High High Low  Low  Low  Low 
## [151] Low  Low  Low  Low  Low  Low  High High Low  High Low  Low  Low  High Low 
## [166] Low  Low  Low  Low  Low  Low  High High Low  Low  High Low  Low  Low  Low 
## [181] High Low  Low  Low  Low  High Low  High High Low  Low  Low  Low  High Low 
## [196] Low  Low  Low  High Low  Low  Low  High Low  High Low  Low  High Low  Low 
## [211] Low  Low  Low  Low  High Low  Low  Low  Low  High Low  High High High Low 
## [226] Low  High High Low  Low  High High High High Low  Low  Low  Low  Low  Low 
## [241] Low  Low  Low  Low  High Low  Low  High Low  High High Low  Low  Low  Low 
## [256] Low  Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
## [271] Low  Low  Low  High Low  Low  Low 
## Levels: Low High
## attr(,"gstari")
## [1] FALSE
## attr(,"call")
## localG(x = data_merge$KEK, listw = W.idw, alternative = "greater")
## attr(,"class")
## [1] "localG"

data_merge$lmZ2 <- Local_M
tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ2",
title = "Local Gi", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Coldspot", "No SAC", "Hotspot"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Indeks Moran Global

W.edw <- mat2listw(w.edw, style='W')

moran(data_merge$KEK, W.edw, n=length(W.edw$neighbours), S0=Szero(W.edw))

## $I
## [1] 0.009386819
## 
## $K
## [1] 18.9482

IM <- moran.test(data_merge$KEK, W.edw,randomisation=T,
alternative="greater", zero.policy = FALSE)
IM

## 
##  Moran I test under randomisation
## 
## data:  data_merge$KEK  
## weights: W.edw    
## 
## Moran I statistic standard deviate = 10.967, p-value < 2.2e-16
## alternative hypothesis: greater
## sample estimates:
## Moran I statistic       Expectation          Variance 
##      9.386819e-03     -3.623188e-03      1.407236e-06

Indeks Moran Lokal

LM <- localmoran(data_merge$KEK, W.edw, alternative = "greater")
head(LM)

##             Ii          E.Ii       Var.Ii        Z.Ii Pr(z > E(Ii))
## 1 -0.001034241 -7.376709e-04 2.454606e-05 -0.05985991    0.52386640
## 2  0.007641443 -7.376709e-04 3.520212e-05  1.41225688    0.07893717
## 3 -0.040977015 -8.612240e-03 3.279923e-04 -1.78706688    0.96303666
## 4 -0.006417541 -5.062083e-05 3.953247e-06 -3.20222902    0.99931816
## 5 -0.010114682 -7.376709e-04 2.936083e-05 -1.73053432    0.95823257
## 6 -0.037849995 -3.425436e-03 1.771435e-04 -2.58646031    0.99515163

data_merge$lmI <- LM[, "Ii"] # local Moran's I
data_merge$lmZ <- LM[, "Z.Ii"] # z-scores
# p-values corresponding to alternative greater
data_merge$lmp <- LM[, "Pr(z > E(Ii))"]

p1 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p1

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

p2 <- tm_shape(data_merge) +
  tm_polygons(col = "lmp", title = "p-value",
              breaks = c(-Inf, 0.05, Inf)) +
  tm_layout(legend.outside = TRUE)
p2

Klasifikasi Warna :

• Less than 0.05 = Ada Autokorelasi Positif

• 0.05 or more = Tidak Ada Autokorelasi Positif

tmap_mode("view")

## tmap mode set to interactive viewing

p3 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p3

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ",
title = "Local Moran's I", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Negative SAC", "No SAC", "Positive SAC"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Rentang nilai untuk klasifikasi warna, dengan:

• -Inf hingga -1.96 mewakili Negative Spatial Autocorrelation (SAC).

• -1.96 hingga 1.96 mewakili No Spatial Autocorrelation (No SAC).

• 1.96 hingga Inf mewakili Positive Spatial Autocorrelation (Positive SAC)

# Autokorelasi spasial positif (Z-score > 1.96)
positive_autocorr <- data_merge[data_merge$lmZ > 1.96, ]
# Autokorelasi spasial negatif (Z-score < -1.96)
negative_autocorr <- data_merge[data_merge$lmZ < -1.96, ]
# Area tanpa autokorelasi spasial (Z-score antara -1.96 dan 1.96)
no_autocorr <- data_merge[data_merge$lmZ >= -1.96 & data_merge$lmZ <= 1.96, ]

moran.plot(data_merge$KEK, W.edw)

Moran Plot memiliki kemiringan positif: ini menunjukkan adanya autokorelasi spasial positif, di mana lokasi dengan nilai tinggi cenderung berdekatan dengan lokasi bernilai tinggi, dan lokasi bernilai rendah cenderung berdekatan dengan lokasi bernilai rendah.

Getis-Ord Gi

Mendeteksi adanya wilayah yang hotspot/coldspot.

Local_M <- localG(data_merge$KEK, W.edw, alternative = "greater")
Local_M

##   [1]  0.059859910 -1.412256882 -1.787066878 -3.202229023  1.730534316
##   [6] -2.586460310  0.980751097 -2.796257615 -0.844877511  2.177682673
##  [11]  3.079466274 -2.604289935 -2.401632818 -2.606632497 -2.268903207
##  [16] -0.110190810 -2.240551217 -3.346579052 -1.122310943 -3.142132277
##  [21] -1.852290278  2.736859633  0.838983983 -3.007012114 -2.269731104
##  [26]  4.709400326  4.217485031  4.729960152  4.690150951  3.271893044
##  [31] -2.110744110 -2.208729697  2.432382222 -1.360245088 -1.933054415
##  [36] -2.621324425 -3.521086622  0.840734772 -1.255117513  3.001201339
##  [41] -3.506632311  4.433426818 -3.217615639 -3.126183630  0.065898893
##  [46] -2.020214500  2.249542006  1.427445864  1.047492680  2.562498908
##  [51] -2.732050239  1.153010945  0.511789997  4.633960537  4.425780603
##  [56]  4.392438001  2.672585267  4.171461118 -2.728567149  3.500157166
##  [61]  0.998518838 -2.616910922 -2.978632684 -2.577971361 -2.824696188
##  [66] -3.348788243 -3.317467204  1.484241743 -2.944385984  4.953249970
##  [71]  0.993446311 -3.403337877 -3.056548119 -2.148202764  4.774699562
##  [76] -3.032349006 -0.178940734 -3.499667112 -1.629642040 -3.416244219
##  [81]  1.789051765  1.968967821  2.266291673 -0.701157480 -2.978781780
##  [86] -2.415039014  4.412685963  1.177469361  4.209233053 -3.443929199
##  [91]  3.352347318 -1.570636179 -2.081316129 -2.755509644 -1.328623007
##  [96] -2.603263679  2.939626254 -3.198240983 -2.773643603 -2.780740755
## [101] -2.484364195 -2.722999513 -2.427445957 -2.550887506 -3.445691917
## [106]  3.611886448  4.377251839  2.677190459  4.631525051  2.776797261
## [111]  4.631687411 -2.639546666  3.935419295  1.646024681 -2.568732905
## [116]  4.930214042  0.341190696  0.999069008  1.009477072  4.642132274
## [121] -3.115294627  2.035352027 -2.221205723  4.429704634  4.673374571
## [126] -2.695808699 -2.726096968 -1.595560602 -3.444956718  4.734229503
## [131]  0.767342064 -1.986133982  3.297942807  2.360223726  2.570435688
## [136] -2.322475560  0.267285596 -0.291072389  0.296932167 -2.499292677
## [141] -2.949871629 -2.356042657 -0.509926571 -2.134622956  0.376654633
## [146]  4.153106256  1.934624254 -3.407621678 -2.977385763 -1.203937612
## [151] -1.569963794 -3.324794720  3.767336800  3.631240086 -3.240553553
## [156] -1.816894269 -2.754681208 -2.715299656  0.002372012 -1.561750204
## [161] -2.639888685 -0.605144344 -3.250570615 -2.747392982 -3.199829321
## [166] -2.813089760 -3.096432181 -2.498130844 -2.602564650 -2.940726071
## [171] -3.050840273 -2.770376532 -2.831947227 -3.070123000 -2.879694032
## [176] -2.573518445 -1.751439517 -2.966041797 -0.421808035 -2.301783360
## [181] -1.731984849 -3.139508871 -0.158154343 -2.406110527 -2.326877818
## [186] -2.797773281 -2.454817050 -2.765177718 -2.715747572 -2.926911473
## [191] -1.831990617 -3.022299500 -2.753917328  1.802109049 -1.168467983
## [196] -3.511475040  3.918740252 -3.514984868  2.691295875  3.385706558
## [201]  2.819849742  4.705631534 -1.207521606 -0.942373515 -3.311973588
## [206]  1.185879077 -0.453199207 -2.542740911 -0.903323842 -3.095037635
## [211]  4.141270545  4.460503582 -0.489372446 -2.946426851 -2.611401093
## [216]  4.025454523 -1.043311858 -1.732648685 -2.758489966 -2.697899363
## [221]  3.307982804  2.320664292  4.633400858 -2.004967657  2.434602355
## [226] -1.320092680  2.047973952 -2.572929166  2.964624020 -0.088843272
## [231]  4.453245526  1.394420775  4.280947315  4.636890280 -0.374817274
## [236] -2.093047177  1.437045066  3.058612304 -3.097683157  1.288743263
## [241] -3.439505031 -2.984901869  4.324182662 -2.514820351  1.903562331
## [246]  1.029770441  0.674596745 -3.302819450 -3.022793958  0.386106277
## [251] -3.160192677  1.686862729  1.420105560 -1.509251829  4.731218406
## [256] -3.501977308 -3.273098645  4.224606622 -3.495215922 -0.568736764
## [261] -2.820992747  3.473789512 -3.300691952  0.667340416  4.728376573
## [266] -1.228539452  0.310524618 -2.089664026  4.193898656  2.047278240
## [271] -1.425895386  4.096328730  4.136514714  0.476949048  4.475762836
## [276] -2.565941350 -2.929000324
## attr(,"internals")
##                 Gi       E(Gi)        V(Gi)        Z(Gi) Pr(z > E(Gi))
##   [1,] 0.003628447 0.003623188 7.717538e-09  0.059859910  4.761336e-01
##   [2,] 0.003474613 0.003623188 1.106792e-08 -1.412256882  9.210628e-01
##   [3,] 0.003452511 0.003623188 9.121551e-09 -1.787066878  9.630367e-01
##   [4,] 0.003190248 0.003623188 1.827892e-08 -3.202229023  9.993182e-01
##   [5,] 0.003789458 0.003623188 9.231354e-09  1.730534316  4.176743e-02
##   [6,] 0.003336662 0.003623188 1.227209e-08 -2.586460310  9.951516e-01
##   [7,] 0.003715867 0.003623188 8.929754e-09  0.980751097  1.633577e-01
##   [8,] 0.003288784 0.003623188 1.430174e-08 -2.796257615  9.974151e-01
##   [9,] 0.003553856 0.003623188 6.734213e-09 -0.844877511  8.009104e-01
##  [10,] 0.003856089 0.003623188 1.143809e-08  2.177682673  1.471483e-02
##  [11,] 0.004015364 0.003623188 1.621851e-08  3.079466274  1.036859e-03
##  [12,] 0.003315618 0.003623188 1.394792e-08 -2.604289935  9.953968e-01
##  [13,] 0.003335828 0.003623188 1.431667e-08 -2.401632818  9.918390e-01
##  [14,] 0.003325749 0.003623188 1.302082e-08 -2.606632497  9.954281e-01
##  [15,] 0.003429303 0.003623188 7.302245e-09 -2.268903207  9.883629e-01
##  [16,] 0.003614502 0.003623188 6.214288e-09 -0.110190810  5.438710e-01
##  [17,] 0.003388145 0.003623188 1.100491e-08 -2.240551217  9.874724e-01
##  [18,] 0.003122964 0.003623188 2.234225e-08 -3.346579052  9.995909e-01
##  [19,] 0.003529577 0.003623188 6.957178e-09 -1.122310943  8.691349e-01
##  [20,] 0.003165719 0.003623188 2.119705e-08 -3.142132277  9.991614e-01
##  [21,] 0.003406168 0.003623188 1.372724e-08 -1.852290278  9.680079e-01
##  [22,] 0.004062651 0.003623188 2.578336e-08  2.736859633  3.101438e-03
##  [23,] 0.003704004 0.003623188 9.278682e-09  0.838983983  2.007391e-01
##  [24,] 0.003228939 0.003623188 1.718980e-08 -3.007012114  9.986809e-01
##  [25,] 0.003398103 0.003623188 9.834396e-09 -2.269731104  9.883880e-01
##  [26,] 0.004278964 0.003623188 1.939009e-08  4.709400326  1.242233e-06
##  [27,] 0.004290194 0.003623188 2.501217e-08  4.217485031  1.235211e-05
##  [28,] 0.004333481 0.003623188 2.255066e-08  4.729960152  1.122820e-06
##  [29,] 0.004291943 0.003623188 2.033106e-08  4.690150951  1.365018e-06
##  [30,] 0.004142379 0.003623188 2.517997e-08  3.271893044  5.341499e-04
##  [31,] 0.003406139 0.003623188 1.057419e-08 -2.110744110  9.826028e-01
##  [32,] 0.003395152 0.003623188 1.065918e-08 -2.208729697  9.864033e-01
##  [33,] 0.004003701 0.003623188 2.447233e-08  2.432382222  7.499934e-03
##  [34,] 0.003486929 0.003623188 1.003460e-08 -1.360245088  9.131238e-01
##  [35,] 0.003467163 0.003623188 6.514819e-09 -1.933054415  9.733853e-01
##  [36,] 0.003315584 0.003623188 1.377027e-08 -2.621324425  9.956206e-01
##  [37,] 0.003032051 0.003623188 2.818536e-08 -3.521086622  9.997851e-01
##  [38,] 0.003726736 0.003623188 1.516911e-08  0.840734772  2.002483e-01
##  [39,] 0.003520954 0.003623188 6.634787e-09 -1.255117513  8.952819e-01
##  [40,] 0.003923539 0.003623188 1.001539e-08  3.001201339  1.344583e-03
##  [41,] 0.003056947 0.003623188 2.607489e-08 -3.506632311  9.997731e-01
##  [42,] 0.004188912 0.003623188 1.628284e-08  4.433426818  4.637350e-06
##  [43,] 0.003163261 0.003623188 2.043194e-08 -3.217615639  9.993537e-01
##  [44,] 0.003200525 0.003623188 1.827934e-08 -3.126183630  9.991145e-01
##  [45,] 0.003629634 0.003623188 9.566206e-09  0.065898893  4.737292e-01
##  [46,] 0.003455984 0.003623188 6.850132e-09 -2.020214500  9.783194e-01
##  [47,] 0.003884110 0.003623188 1.345338e-08  2.249542006  1.223902e-02
##  [48,] 0.003757591 0.003623188 8.865308e-09  1.427445864  7.672571e-02
##  [49,] 0.003752723 0.003623188 1.529226e-08  1.047492680  1.474362e-01
##  [50,] 0.003905287 0.003623188 1.211925e-08  2.562498908  5.196096e-03
##  [51,] 0.003304288 0.003623188 1.362489e-08 -2.732050239  9.968529e-01
##  [52,] 0.003776722 0.003623188 1.773133e-08  1.153010945  1.244529e-01
##  [53,] 0.003672677 0.003623188 9.350171e-09  0.511789997  3.043990e-01
##  [54,] 0.004211322 0.003623188 1.610823e-08  4.633960537  1.793677e-06
##  [55,] 0.004200388 0.003623188 1.700876e-08  4.425780603  4.804708e-06
##  [56,] 0.004162933 0.003623188 1.509961e-08  4.392438001  5.604330e-06
##  [57,] 0.003935368 0.003623188 1.364410e-08  2.672585267  3.763462e-03
##  [58,] 0.004202741 0.003623188 1.930234e-08  4.171461118  1.513264e-05
##  [59,] 0.003289251 0.003623188 1.497825e-08 -2.728567149  9.968195e-01
##  [60,] 0.004064269 0.003623188 1.588038e-08  3.500157166  2.324920e-04
##  [61,] 0.003753889 0.003623188 1.713323e-08  0.998518838  1.590139e-01
##  [62,] 0.003300542 0.003623188 1.520108e-08 -2.616910922  9.955635e-01
##  [63,] 0.003236481 0.003623188 1.685512e-08 -2.978632684  9.985523e-01
##  [64,] 0.003363009 0.003623188 1.018570e-08 -2.577971361  9.950309e-01
##  [65,] 0.003270358 0.003623188 1.560234e-08 -2.824696188  9.976337e-01
##  [66,] 0.003136263 0.003623188 2.114216e-08 -3.348788243  9.995942e-01
##  [67,] 0.003128685 0.003623188 2.221904e-08 -3.317467204  9.995458e-01
##  [68,] 0.003752701 0.003623188 7.614081e-09  1.484241743  6.887240e-02
##  [69,] 0.003223381 0.003623188 1.843791e-08 -2.944385984  9.983820e-01
##  [70,] 0.004322135 0.003623188 1.991164e-08  4.953249970  3.649206e-07
##  [71,] 0.003718811 0.003623188 9.264644e-09  0.993446311  1.602463e-01
##  [72,] 0.003118037 0.003623188 2.203096e-08 -3.403337877  9.996672e-01
##  [73,] 0.003203496 0.003623188 1.885380e-08 -3.056548119  9.988805e-01
##  [74,] 0.003426517 0.003623188 8.381719e-09 -2.148202764  9.841512e-01
##  [75,] 0.004308463 0.003623188 2.059857e-08  4.774699562  8.998787e-07
##  [76,] 0.003221203 0.003623188 1.757365e-08 -3.032349006  9.987867e-01
##  [77,] 0.003605530 0.003623188 9.738050e-09 -0.178940734  5.710079e-01
##  [78,] 0.003048395 0.003623188 2.697551e-08 -3.499667112  9.997671e-01
##  [79,] 0.003490273 0.003623188 6.652252e-09 -1.629642040  9.484114e-01
##  [80,] 0.003101782 0.003623188 2.329455e-08 -3.416244219  9.996825e-01
##  [81,] 0.003853786 0.003623188 1.661366e-08  1.789051765  3.680324e-02
##  [82,] 0.003898657 0.003623188 1.957346e-08  1.968967821  2.447839e-02
##  [83,] 0.003991095 0.003623188 2.635380e-08  2.266291673  1.171677e-02
##  [84,] 0.003566659 0.003623188 6.499990e-09 -0.701157480  7.583976e-01
##  [85,] 0.003247476 0.003623188 1.590870e-08 -2.978781780  9.985530e-01
##  [86,] 0.003404224 0.003623188 8.220547e-09 -2.415039014  9.921332e-01
##  [87,] 0.004307777 0.003623188 2.406874e-08  4.412685963  5.104802e-06
##  [88,] 0.003767284 0.003623188 1.497620e-08  1.177469361  1.195041e-01
##  [89,] 0.004116280 0.003623188 1.372303e-08  4.209233053  1.281195e-05
##  [90,] 0.003081641 0.003623188 2.472656e-08 -3.443929199  9.997133e-01
##  [91,] 0.004103525 0.003623188 2.053020e-08  3.352347318  4.006471e-04
##  [92,] 0.003484141 0.003623188 7.837449e-09 -1.570636179  9.418664e-01
##  [93,] 0.003387529 0.003623188 1.282020e-08 -2.081316129  9.812975e-01
##  [94,] 0.003272571 0.003623188 1.619058e-08 -2.755509644  9.970700e-01
##  [95,] 0.003520420 0.003623188 5.983002e-09 -1.328623007  9.080138e-01
##  [96,] 0.003340673 0.003623188 1.177740e-08 -2.603263679  9.953830e-01
##  [97,] 0.003975526 0.003623188 1.436596e-08  2.939626254  1.643042e-03
##  [98,] 0.003183902 0.003623188 1.886568e-08 -3.198240983  9.993087e-01
##  [99,] 0.003266112 0.003623188 1.657377e-08 -2.773643603  9.972284e-01
## [100,] 0.003284031 0.003623188 1.487580e-08 -2.780740755  9.972882e-01
## [101,] 0.003362690 0.003623188 1.099465e-08 -2.484364195  9.935108e-01
## [102,] 0.003304377 0.003623188 1.370792e-08 -2.722999513  9.967654e-01
## [103,] 0.003371804 0.003623188 1.072453e-08 -2.427445957  9.923972e-01
## [104,] 0.003331202 0.003623188 1.310213e-08 -2.550887506  9.946275e-01
## [105,] 0.003099887 0.003623188 2.306489e-08 -3.445691917  9.997152e-01
## [106,] 0.004194761 0.003623188 2.504234e-08  3.611886448  1.519888e-04
## [107,] 0.004204025 0.003623188 1.760778e-08  4.377251839  6.009254e-06
## [108,] 0.003921447 0.003623188 1.241157e-08  2.677190459  3.712121e-03
## [109,] 0.004269445 0.003623188 1.946983e-08  4.631525051  1.814910e-06
## [110,] 0.003970336 0.003623188 1.562936e-08  2.776797261  2.744870e-03
## [111,] 0.004211702 0.003623188 1.614489e-08  4.631687411  1.813487e-06
## [112,] 0.003324514 0.003623188 1.280380e-08 -2.639546666  9.958492e-01
## [113,] 0.004113690 0.003623188 1.553454e-08  3.935419295  4.152579e-05
## [114,] 0.003773874 0.003623188 8.380519e-09  1.646024681  4.987934e-02
## [115,] 0.003329233 0.003623188 1.309555e-08 -2.568732905  9.948964e-01
## [116,] 0.004342047 0.003623188 2.125963e-08  4.930214042  4.106978e-07
## [117,] 0.003657379 0.003623188 1.004192e-08  0.341190696  3.664800e-01
## [118,] 0.003731435 0.003623188 1.173923e-08  0.999069008  1.588806e-01
## [119,] 0.003750570 0.003623188 1.592294e-08  1.009477072  1.563729e-01
## [120,] 0.004337232 0.003623188 2.365998e-08  4.642132274  1.724160e-06
## [121,] 0.003206621 0.003623188 1.788019e-08 -3.115294627  9.990812e-01
## [122,] 0.003918994 0.003623188 2.112200e-08  2.035352027  2.090773e-02
## [123,] 0.003426692 0.003623188 7.825828e-09 -2.221205723  9.868315e-01
## [124,] 0.004248029 0.003623188 1.989708e-08  4.429704634  4.718112e-06
## [125,] 0.004339344 0.003623188 2.348301e-08  4.673374571  1.481453e-06
## [126,] 0.003292561 0.003623188 1.504182e-08 -2.695808699  9.964891e-01
## [127,] 0.003302568 0.003623188 1.383249e-08 -2.726096968  9.967956e-01
## [128,] 0.003491829 0.003623188 6.777949e-09 -1.595560602  9.447065e-01
## [129,] 0.003074638 0.003623188 2.535513e-08 -3.444956718  9.997144e-01
## [130,] 0.004338675 0.003623188 2.284041e-08  4.734229503  1.099444e-06
## [131,] 0.003686788 0.003623188 6.869515e-09  0.767342064  2.214391e-01
## [132,] 0.003425201 0.003623188 9.937101e-09 -1.986133982  9.764908e-01
## [133,] 0.003990965 0.003623188 1.243602e-08  3.297942807  4.869799e-04
## [134,] 0.003861375 0.003623188 1.018422e-08  2.360223726  9.131958e-03
## [135,] 0.003880464 0.003623188 1.001807e-08  2.570435688  5.078534e-03
## [136,] 0.003337761 0.003623188 1.510390e-08 -2.322475560  9.898963e-01
## [137,] 0.003656533 0.003623188 1.556338e-08  0.267285596  3.946246e-01
## [138,] 0.003597649 0.003623188 7.698764e-09 -0.291072389  6.145020e-01
## [139,] 0.003652006 0.003623188 9.418737e-09  0.296932167  3.832591e-01
## [140,] 0.003338112 0.003623188 1.301033e-08 -2.499292677  9.937779e-01
## [141,] 0.003245628 0.003623188 1.638199e-08 -2.949871629  9.984105e-01
## [142,] 0.003376448 0.003623188 1.096771e-08 -2.356042657  9.907646e-01
## [143,] 0.003567802 0.003623188 1.179735e-08 -0.509926571  6.949485e-01
## [144,] 0.003401816 0.003623188 1.075484e-08 -2.134622956  9.836041e-01
## [145,] 0.003670579 0.003623188 1.583084e-08  0.376654633  3.532151e-01
## [146,] 0.004171088 0.003623188 1.740427e-08  4.153106256  1.639962e-05
## [147,] 0.003880897 0.003623188 1.774447e-08  1.934624254  2.651821e-02
## [148,] 0.003107170 0.003623188 2.293123e-08 -3.407621678  9.996723e-01
## [149,] 0.003233301 0.003623188 1.714776e-08 -2.977385763  9.985464e-01
## [150,] 0.003488791 0.003623188 1.246167e-08 -1.203937612  8.856932e-01
## [151,] 0.003466366 0.003623188 9.977888e-09 -1.569963794  9.417882e-01
## [152,] 0.003115574 0.003623188 2.330974e-08 -3.324794720  9.995576e-01
## [153,] 0.004086281 0.003623188 1.511006e-08  3.767336800  8.249917e-05
## [154,] 0.004112364 0.003623188 1.814762e-08  3.631240086  1.410313e-04
## [155,] 0.003165751 0.003623188 1.992618e-08 -3.240553553  9.994035e-01
## [156,] 0.003464056 0.003623188 7.671143e-09 -1.816894269  9.653833e-01
## [157,] 0.003269800 0.003623188 1.645742e-08 -2.754681208  9.970625e-01
## [158,] 0.003296479 0.003623188 1.447734e-08 -2.715299656  9.966892e-01
## [159,] 0.003623379 0.003623188 6.483450e-09  0.002372012  4.990537e-01
## [160,] 0.003456764 0.003623188 1.135559e-08 -1.561750204  9.408266e-01
## [161,] 0.003332821 0.003623188 1.209830e-08 -2.639888685  9.958533e-01
## [162,] 0.003573933 0.003623188 6.625106e-09 -0.605144344  7.274585e-01
## [163,] 0.003167623 0.003623188 1.964183e-08 -3.250570615  9.994241e-01
## [164,] 0.003271004 0.003623188 1.643231e-08 -2.747392982  9.969964e-01
## [165,] 0.003146037 0.003623188 2.223615e-08 -3.199829321  9.993125e-01
## [166,] 0.003283664 0.003623188 1.456714e-08 -2.813089760  9.975466e-01
## [167,] 0.003185693 0.003623188 1.996292e-08 -3.096432181  9.990207e-01
## [168,] 0.003344579 0.003623188 1.243834e-08 -2.498130844  9.937575e-01
## [169,] 0.003316867 0.003623188 1.385324e-08 -2.602564650  9.953735e-01
## [170,] 0.003222839 0.003623188 1.853398e-08 -2.940726071  9.983628e-01
## [171,] 0.003190126 0.003623188 2.014940e-08 -3.050840273  9.988590e-01
## [172,] 0.003267789 0.003623188 1.645719e-08 -2.770376532  9.972004e-01
## [173,] 0.003261435 0.003623188 1.631755e-08 -2.831947227  9.976867e-01
## [174,] 0.003185452 0.003623188 2.032891e-08 -3.070123000  9.989301e-01
## [175,] 0.003259410 0.003623188 1.595806e-08 -2.879694032  9.980097e-01
## [176,] 0.003320624 0.003623188 1.382228e-08 -2.573518445  9.949665e-01
## [177,] 0.003479928 0.003623188 6.690533e-09 -1.751439517  9.600649e-01
## [178,] 0.003214816 0.003623188 1.895655e-08 -2.966041797  9.984917e-01
## [179,] 0.003589017 0.003623188 6.562839e-09 -0.421808035  6.634174e-01
## [180,] 0.003382630 0.003623188 1.092221e-08 -2.301783360  9.893263e-01
## [181,] 0.003453436 0.003623188 9.606055e-09 -1.731984849  9.583619e-01
## [182,] 0.003161722 0.003623188 2.160516e-08 -3.139508871  9.991538e-01
## [183,] 0.003610388 0.003623188 6.550990e-09 -0.158154343  5.628324e-01
## [184,] 0.003405935 0.003623188 8.152676e-09 -2.406110527  9.919383e-01
## [185,] 0.003419426 0.003623188 7.668367e-09 -2.326877818  9.900141e-01
## [186,] 0.003263471 0.003623188 1.653091e-08 -2.797773281  9.974272e-01
## [187,] 0.003388675 0.003623188 9.126367e-09 -2.454817050  9.929522e-01
## [188,] 0.003297645 0.003623188 1.386030e-08 -2.765177718  9.971554e-01
## [189,] 0.003281137 0.003623188 1.586372e-08 -2.715747572  9.966937e-01
## [190,] 0.003239719 0.003623188 1.716490e-08 -2.926911473  9.982883e-01
## [191,] 0.003442196 0.003623188 9.760537e-09 -1.831990617  9.665236e-01
## [192,] 0.003202477 0.003623188 1.937730e-08 -3.022299500  9.987457e-01
## [193,] 0.003266507 0.003623188 1.677482e-08 -2.753917328  9.970557e-01
## [194,] 0.003802264 0.003623188 9.874421e-09  1.802109049  3.576413e-02
## [195,] 0.003529969 0.003623188 6.364677e-09 -1.168467983  8.786910e-01
## [196,] 0.003045309 0.003623188 2.708289e-08 -3.511475040  9.997772e-01
## [197,] 0.004224976 0.003623188 2.358268e-08  3.918740252  4.450649e-05
## [198,] 0.003039313 0.003623188 2.759267e-08 -3.514984868  9.997801e-01
## [199,] 0.004012219 0.003623188 2.089516e-08  2.691295875  3.558752e-03
## [200,] 0.004130141 0.003623188 2.242004e-08  3.385706558  3.549763e-04
## [201,] 0.004058932 0.003623188 2.387866e-08  2.819849742  2.402307e-03
## [202,] 0.004340183 0.003623188 2.321644e-08  4.705631534  1.265406e-06
## [203,] 0.003505611 0.003623188 9.481072e-09 -1.207521606  8.863843e-01
## [204,] 0.003529759 0.003623188 9.829192e-09 -0.942373515  8.269993e-01
## [205,] 0.003132290 0.003623188 2.196899e-08 -3.311973588  9.995368e-01
## [206,] 0.003736951 0.003623188 9.202824e-09  1.185879077  1.178350e-01
## [207,] 0.003587870 0.003623188 6.073373e-09 -0.453199207  6.747973e-01
## [208,] 0.003323409 0.003623188 1.389955e-08 -2.542740911  9.945007e-01
## [209,] 0.003553070 0.003623188 6.025342e-09 -0.903323842  8.168230e-01
## [210,] 0.003210821 0.003623188 1.775161e-08 -3.095037635  9.990161e-01
## [211,] 0.004126374 0.003623188 1.476347e-08  4.141270545  1.726936e-05
## [212,] 0.004311297 0.003623188 2.379831e-08  4.460503582  4.088365e-06
## [213,] 0.003578851 0.003623188 8.208373e-09 -0.489372446  6.877110e-01
## [214,] 0.003243853 0.003623188 1.657507e-08 -2.946426851  9.983927e-01
## [215,] 0.003328243 0.003623188 1.275662e-08 -2.611401093  9.954914e-01
## [216,] 0.004257420 0.003623188 2.482368e-08  4.025454523  2.843268e-05
## [217,] 0.003515716 0.003623188 1.061116e-08 -1.043311858  8.515981e-01
## [218,] 0.003484813 0.003623188 6.378181e-09 -1.732648685  9.584209e-01
## [219,] 0.003261476 0.003623188 1.719427e-08 -2.758489966  9.970965e-01
## [220,] 0.003295104 0.003623188 1.478837e-08 -2.697899363  9.965111e-01
## [221,] 0.004075367 0.003623188 1.868503e-08  3.307982804  4.698529e-04
## [222,] 0.003920007 0.003623188 1.635908e-08  2.320664292  1.015248e-02
## [223,] 0.004247453 0.003623188 1.815257e-08  4.633400858  1.798535e-06
## [224,] 0.003414723 0.003623188 1.081065e-08 -2.004967657  9.775167e-01
## [225,] 0.003881566 0.003623188 1.126298e-08  2.434602355  7.454081e-03
## [226,] 0.003491355 0.003623188 9.973414e-09 -1.320092680  9.065980e-01
## [227,] 0.003939310 0.003623188 2.382651e-08  2.047973952  2.028128e-02
## [228,] 0.003327616 0.003623188 1.319691e-08 -2.572929166  9.949579e-01
## [229,] 0.004102117 0.003623188 2.609765e-08  2.964624020  1.515265e-03
## [230,] 0.003614451 0.003623188 9.672816e-09 -0.088843272  5.353968e-01
## [231,] 0.004219565 0.003623188 1.793447e-08  4.453245526  4.229097e-06
## [232,] 0.003810547 0.003623188 1.805337e-08  1.394420775  8.159529e-02
## [233,] 0.004250492 0.003623188 2.147220e-08  4.280947315  9.304971e-06
## [234,] 0.004241498 0.003623188 1.778108e-08  4.636890280  1.768450e-06
## [235,] 0.003586178 0.003623188 9.750318e-09 -0.374817274  6.461018e-01
## [236,] 0.003419819 0.003623188 9.440895e-09 -2.093047177  9.818275e-01
## [237,] 0.003760597 0.003623188 9.142918e-09  1.437045066  7.535259e-02
## [238,] 0.003973835 0.003623188 1.314285e-08  3.058612304  1.111824e-03
## [239,] 0.003195968 0.003623188 1.902083e-08 -3.097683157  9.990248e-01
## [240,] 0.003795351 0.003623188 1.784614e-08  1.288743263  9.874368e-02
## [241,] 0.003087776 0.003623188 2.423174e-08 -3.439505031  9.997086e-01
## [242,] 0.003217228 0.003623188 1.849726e-08 -2.984901869  9.985817e-01
## [243,] 0.004145204 0.003623188 1.457331e-08  4.324182662  7.654920e-06
## [244,] 0.003346967 0.003623188 1.206422e-08 -2.514820351  9.940453e-01
## [245,] 0.003884057 0.003623188 1.878063e-08  1.903562331  2.848360e-02
## [246,] 0.003759370 0.003623188 1.748870e-08  1.029770441  1.515589e-01
## [247,] 0.003688249 0.003623188 9.301445e-09  0.674596745  2.499660e-01
## [248,] 0.003138878 0.003623188 2.150193e-08 -3.302819450  9.995214e-01
## [249,] 0.003199197 0.003623188 1.967417e-08 -3.022793958  9.987477e-01
## [250,] 0.003670889 0.003623188 1.526273e-08  0.386106277  3.497090e-01
## [251,] 0.003199035 0.003623188 1.801437e-08 -3.160192677  9.992117e-01
## [252,] 0.003827606 0.003623188 1.468514e-08  1.686862729  4.581488e-02
## [253,] 0.003796321 0.003623188 1.486331e-08  1.420105560  7.778848e-02
## [254,] 0.003504579 0.003623188 6.176076e-09 -1.509251829  9.343828e-01
## [255,] 0.004228869 0.003623188 1.638860e-08  4.731218406  1.115881e-06
## [256,] 0.003055630 0.003623188 2.626599e-08 -3.501977308  9.997691e-01
## [257,] 0.003148380 0.003623188 2.104357e-08 -3.273098645  9.994681e-01
## [258,] 0.004182843 0.003623188 1.754959e-08  4.224606622  1.196793e-05
## [259,] 0.003072040 0.003623188 2.486502e-08 -3.495215922  9.997632e-01
## [260,] 0.003565682 0.003623188 1.022361e-08 -0.568736764  7.152326e-01
## [261,] 0.003254265 0.003623188 1.710282e-08 -2.820992747  9.976062e-01
## [262,] 0.004042835 0.003623188 1.459356e-08  3.473789512  2.565817e-04
## [263,] 0.003146881 0.003623188 2.082402e-08 -3.300691952  9.995178e-01
## [264,] 0.003678907 0.003623188 6.971093e-09  0.667340416  2.522774e-01
## [265,] 0.004321649 0.003623188 2.182026e-08  4.728376573  1.131611e-06
## [266,] 0.003496240 0.003623188 1.067763e-08 -1.228539452  8.903777e-01
## [267,] 0.003648217 0.003623188 6.496373e-09  0.310524618  3.780810e-01
## [268,] 0.003396812 0.003623188 1.173571e-08 -2.089664026  9.816760e-01
## [269,] 0.004227433 0.003623188 2.075817e-08  4.193898656  1.371003e-05
## [270,] 0.003874655 0.003623188 1.508712e-08  2.047278240  2.031539e-02
## [271,] 0.003466390 0.003623188 1.209230e-08 -1.425895386  9.230507e-01
## [272,] 0.004233498 0.003623188 2.219785e-08  4.096328730  2.098769e-05
## [273,] 0.004169463 0.003623188 1.744028e-08  4.136514714  1.763104e-05
## [274,] 0.003682848 0.003623188 1.564624e-08  0.476949048  3.166992e-01
## [275,] 0.004323750 0.003623188 2.449961e-08  4.475762836  3.806948e-06
## [276,] 0.003341728 0.003623188 1.203210e-08 -2.565941350  9.948552e-01
## [277,] 0.003253967 0.003623188 1.589044e-08 -2.929000324  9.982997e-01
## attr(,"cluster")
##   [1] Low  Low  High High Low  High Low  Low  Low  High High Low  Low  High Low 
##  [16] Low  High Low  Low  Low  Low  Low  Low  Low  High High High Low  High Low 
##  [31] Low  Low  Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [46] Low  High High Low  High Low  Low  Low  High High High High Low  Low  Low 
##  [61] Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [76] Low  Low  Low  Low  Low  Low  High High Low  High Low  High Low  Low  Low 
##  [91] Low  Low  Low  Low  Low  High High Low  High Low  Low  High Low  High Low 
## [106] Low  High Low  High High Low  Low  High Low  High High Low  Low  Low  Low 
## [121] Low  Low  Low  High Low  High High Low  Low  Low  Low  Low  Low  Low  Low 
## [136] Low  High High Low  High Low  Low  High High High High Low  Low  Low  Low 
## [151] Low  Low  Low  Low  Low  Low  High High Low  High Low  Low  Low  High Low 
## [166] Low  Low  Low  Low  Low  Low  High High Low  Low  High Low  Low  Low  Low 
## [181] High Low  Low  Low  Low  High Low  High High Low  Low  Low  Low  High Low 
## [196] Low  Low  Low  High Low  Low  Low  High Low  High Low  Low  High Low  Low 
## [211] Low  Low  Low  Low  High Low  Low  Low  Low  High Low  High High High Low 
## [226] Low  High High Low  Low  High High High High Low  Low  Low  Low  Low  Low 
## [241] Low  Low  Low  Low  High Low  Low  High Low  High High Low  Low  Low  Low 
## [256] Low  Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
## [271] Low  Low  Low  High Low  Low  Low 
## Levels: Low High
## attr(,"gstari")
## [1] FALSE
## attr(,"call")
## localG(x = data_merge$KEK, listw = W.edw, alternative = "greater")
## attr(,"class")
## [1] "localG"

data_merge$lmZ2 <- Local_M
tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ2",
title = "Local Gi", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Coldspot", "No SAC", "Hotspot"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Indeks Moran Global

W.knn <- nb2listw(w.knn, style='W')

moran(data_merge$KEK, W.knn, n=length(W.knn$neighbours), S0=Szero(W.knn))

## $I
## [1] 0.4490796
## 
## $K
## [1] 18.9482

IM <- moran.test(data_merge$KEK, W.knn,randomisation=T,
alternative="greater", zero.policy = FALSE)
IM

## 
##  Moran I test under randomisation
## 
## data:  data_merge$KEK  
## weights: W.knn    
## 
## Moran I statistic standard deviate = 13.196, p-value < 2.2e-16
## alternative hypothesis: greater
## sample estimates:
## Moran I statistic       Expectation          Variance 
##       0.449079634      -0.003623188       0.001176860

Indeks Moran Lokal

LM <- localmoran(data_merge$KEK, W.knn, alternative = "greater")
head(LM)

##            Ii          E.Ii      Var.Ii       Z.Ii Pr(z > E(Ii))
## 1  0.20359716 -7.376709e-04 0.040242830  1.0185870    0.15419953
## 2  0.04943840 -7.376709e-04 0.040242830  0.2501223    0.40124639
## 3  0.97234254 -8.612240e-03 0.466128896  1.4367989    0.07538758
## 4 -0.04660356 -5.062083e-05 0.002763463 -0.8855647    0.81207393
## 5  0.20359716 -7.376709e-04 0.040242830  1.0185870    0.15419953
## 6 -0.10653468 -3.425436e-03 0.186368293 -0.2388425    0.59438614

data_merge$lmI <- LM[, "Ii"] # local Moran's I
data_merge$lmZ <- LM[, "Z.Ii"] # z-scores
# p-values corresponding to alternative greater
data_merge$lmp <- LM[, "Pr(z > E(Ii))"]

p1 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p1

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

p2 <- tm_shape(data_merge) +
  tm_polygons(col = "lmp", title = "p-value",
              breaks = c(-Inf, 0.05, Inf)) +
  tm_layout(legend.outside = TRUE)
p2

Klasifikasi Warna :

• Less than 0.05 = Ada Autokorelasi Positif

• 0.05 or more = Tidak Ada Autokorelasi Positif

tmap_mode("view")

## tmap mode set to interactive viewing

p3 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p3

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ",
title = "Local Moran's I", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Negative SAC", "No SAC", "Positive SAC"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Rentang nilai untuk klasifikasi warna, dengan:

• -Inf hingga -1.96 mewakili Negative Spatial Autocorrelation (SAC).

• -1.96 hingga 1.96 mewakili No Spatial Autocorrelation (No SAC).

• 1.96 hingga Inf mewakili Positive Spatial Autocorrelation (Positive SAC)

# Autokorelasi spasial positif (Z-score > 1.96)
positive_autocorr <- data_merge[data_merge$lmZ > 1.96, ]
# Autokorelasi spasial negatif (Z-score < -1.96)
negative_autocorr <- data_merge[data_merge$lmZ < -1.96, ]
# Area tanpa autokorelasi spasial (Z-score antara -1.96 dan 1.96)
no_autocorr <- data_merge[data_merge$lmZ >= -1.96 & data_merge$lmZ <= 1.96, ]

moran.plot(data_merge$KEK, W.knn)

Moran Plot memiliki kemiringan positif: ini menunjukkan adanya autokorelasi spasial positif, di mana lokasi dengan nilai tinggi cenderung berdekatan dengan lokasi bernilai tinggi, dan lokasi bernilai rendah cenderung berdekatan dengan lokasi bernilai rendah.

Getis-Ord Gi

Mendeteksi adanya wilayah yang hotspot/coldspot.

Local_M <- localG(data_merge$KEK, W.knn, alternative = "greater")
Local_M

##   [1] -1.018587043 -0.250122285  1.436798857 -0.885564710 -1.018587043
##   [6] -0.238842496 -0.378199745 -1.018587043 -1.018587043  0.394769811
##  [11]  0.012992394 -0.250122285 -0.634354664  0.650836715 -1.018587043
##  [16] -1.018587043  1.421739075 -1.018587043 -1.018587043 -0.506277205
##  [21] -1.018587043 -0.506277205 -1.018587043 -0.890509584  0.940073008
##  [26]  8.180644436  0.015321089 -0.503796354  7.404348110 -1.018587043
##  [31] -0.762432124 -0.762432124 -0.247723069 -0.883482767 -1.018587043
##  [36]  1.168581248 -1.018587043  0.264423501 -1.018587043 -1.018587043
##  [41] -1.018587043  2.823736747 -1.018587043 -1.018587043 -1.018587043
##  [46] -1.018587043  0.269128154 -0.115075486 -0.506277205 -0.627347007
##  [51] -0.634354664 -0.634354664 -1.018587043  4.491840277  5.844470515
##  [56]  3.211505756  0.012992394  1.286807231 -0.634354664 -0.250122285
##  [61] -0.759869639  1.034937071 -0.634354664 -0.117363998 -0.634354664
##  [66] -0.506277205 -1.018587043 -1.018587043 -1.018587043  5.001053562
##  [71] -1.018587043 -1.018587043 -0.634354664 -1.018587043 -0.122044826
##  [76] -0.378199745 -1.018587043 -1.018587043 -1.018587043 -1.018587043
##  [81] -0.506277205 -0.371211247 -0.499279127 -1.018587043 -1.011550647
##  [86] -0.762432124 -0.110592440 -0.506277205  3.464124045 -1.018587043
##  [91] -0.634354664  1.286807231 -0.634354664 -1.018587043 -1.018587043
##  [96]  0.397196034 -0.243143367 -0.503796354  0.781399674 -0.250122285
## [101]  0.006032634  0.145907672 -0.378199745  0.522803263 -1.018587043
## [106] -1.018587043  7.120955420  0.134110094  4.500221706  0.522803263
## [111]  4.104511344 -0.634354664  1.291003975  0.006032634  0.522803263
## [116]  5.772174798 -1.018587043 -1.018587043 -0.634354664  0.008350216
## [121] -0.887906281  0.520496786 -1.018587043  2.366061489  0.134110094
## [126]  1.043658065  0.525263914 -1.018587043 -1.018587043 -0.503796354
## [131] -1.018587043  2.311426908 -0.762432124  0.006032634 -0.634354664
## [136] -1.018587043  3.116200157 -1.011550647 -1.018587043  0.012992394
## [141] -0.634354664 -0.762432124  1.312845745 -1.013598162  3.116200157
## [146]  4.052336521 -0.247723069 -1.018587043 -0.634354664 -0.634354664
## [151] -0.375759712 -0.506277205 -0.762432124 -0.634354664 -1.018587043
## [156] -0.890509584  0.781399674  1.171908121 -1.018587043 -1.011550647
## [161] -0.634354664 -1.018587043 -0.762432124  0.781399674  0.134110094
## [166] -0.634354664 -1.018587043  0.646419932  1.414884690 -0.634354664
## [171] -0.378199745  0.397196034  1.677874835 -0.506277205 -0.634354664
## [176]  1.549806955 -1.018587043 -1.018587043 -1.018587043 -0.762432124
## [181]  2.069658514  0.134110094 -1.018587043 -0.762432124 -1.018587043
## [186]  1.165603314 -0.762432124  0.012992394  0.397196034 -1.018587043
## [191] -0.250122285 -0.378199745  0.518342473 -0.238842496 -1.018587043
## [196] -1.018587043 -1.018587043 -1.018587043 -0.627347007 -1.018587043
## [201] -0.506277205  0.008350216  0.525263914  0.518342473 -0.369098964
## [206] -1.018587043 -1.018587043  1.547070879 -1.018587043 -0.762432124
## [211]  1.542962150  0.392460143 -1.018587043 -0.634354664 -0.112818929
## [216] -0.378199745 -0.506277205 -1.018587043 -1.018587043  1.428408233
## [221] -0.378199745 -0.243143367  8.180644436  1.421739075  0.262187553
## [226]  0.774497392  0.010669454  0.653331794 -0.890509584 -1.018587043
## [231]  4.946076411 -0.883482767  1.404495747  7.145582668 -1.018587043
## [236]  0.774497392 -1.018587043 -0.250122285 -1.018587043 -0.631832997
## [241] -1.018587043 -1.018587043  1.542962150  0.646419932 -0.371211247
## [246] -0.759869639 -1.018587043 -0.495342608 -1.018587043  3.116200157
## [251] -0.883482767 -0.506277205 -0.634354664 -1.018587043  3.207969126
## [256] -1.018587043 -1.018587043  2.637592870 -1.018587043 -1.018587043
## [261] -0.250122285 -0.378199745 -1.018587043  0.006032634 -0.247723069
## [266] -0.250122285 -1.018587043  0.134110094  1.542962150 -0.122044826
## [271] -0.634354664 -0.890509584  1.542962150  3.236990420  0.392460143
## [276] -0.378199745 -1.018587043
## attr(,"internals")
##                  Gi       E(Gi)        V(Gi)        Z(Gi) Pr(z > E(Gi))
##   [1,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##   [2,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
##   [3,] 0.0087962963 0.003623188 1.296317e-05  1.436798857  7.538758e-02
##   [4,] 0.0004576659 0.003623188 1.277763e-05 -0.885564710  8.120739e-01
##   [5,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##   [6,] 0.0027649770 0.003623188 1.291116e-05 -0.238842496  5.943861e-01
##   [7,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
##   [8,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##   [9,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [10,] 0.0050343249 0.003623188 1.277763e-05  0.394769811  3.465064e-01
##  [11,] 0.0036697248 0.003623188 1.282941e-05  0.012992394  4.948169e-01
##  [12,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
##  [13,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [14,] 0.0059496568 0.003623188 1.277763e-05  0.650836715  2.575759e-01
##  [15,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [16,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [17,] 0.0087155963 0.003623188 1.282941e-05  1.421739075  7.755101e-02
##  [18,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [19,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [20,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [21,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [22,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [23,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [24,] 0.0004555809 0.003623188 1.265277e-05 -0.890509584  8.134038e-01
##  [25,] 0.0070093458 0.003623188 1.297451e-05  0.940073008  1.735901e-01
##  [26,] 0.0331002331 0.003623188 1.298353e-05  8.180644436  1.411649e-16
##  [27,] 0.0036781609 0.003623188 1.287395e-05  0.015321089  4.938880e-01
##  [28,] 0.0018264840 0.003623188 1.271871e-05 -0.503796354  6.927978e-01
##  [29,] 0.0303030303 0.003623188 1.298353e-05  7.404348110  6.589814e-14
##  [30,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [31,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
##  [32,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
##  [33,] 0.0027397260 0.003623188 1.271871e-05 -0.247723069  5.978257e-01
##  [34,] 0.0004587156 0.003623188 1.282941e-05 -0.883482767  8.115123e-01
##  [35,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [36,] 0.0078160920 0.003623188 1.287395e-05  1.168581248  1.212862e-01
##  [37,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [38,] 0.0045662100 0.003623188 1.271871e-05  0.264423501  3.957268e-01
##  [39,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [40,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [41,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [42,] 0.0136674260 0.003623188 1.265277e-05  2.823736747  2.373368e-03
##  [43,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [44,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [45,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [46,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [47,] 0.0045871560 0.003623188 1.282941e-05  0.269128154  3.939155e-01
##  [48,] 0.0032110092 0.003623188 1.282941e-05 -0.115075486  5.458074e-01
##  [49,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [50,] 0.0013761468 0.003623188 1.282941e-05 -0.627347007  7.347841e-01
##  [51,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [52,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [53,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [54,] 0.0196796339 0.003623188 1.277763e-05  4.491840277  3.530519e-06
##  [55,] 0.0240384615 0.003623188 1.220169e-05  5.844470515  2.540910e-09
##  [56,] 0.0151029748 0.003623188 1.277763e-05  3.211505756  6.602066e-04
##  [57,] 0.0036697248 0.003623188 1.282941e-05  0.012992394  4.948169e-01
##  [58,] 0.0082004556 0.003623188 1.265277e-05  1.286807231  9.908074e-02
##  [59,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [60,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
##  [61,] 0.0009132420 0.003623188 1.271871e-05 -0.759869639  7.763337e-01
##  [62,] 0.0073226545 0.003623188 1.277763e-05  1.034937071  1.503492e-01
##  [63,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [64,] 0.0032036613 0.003623188 1.277763e-05 -0.117363998  5.467142e-01
##  [65,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [66,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [67,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [68,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [69,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [70,] 0.0214123007 0.003623188 1.265277e-05  5.001053562  2.850893e-07
##  [71,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [72,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [73,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [74,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [75,] 0.0031890661 0.003623188 1.265277e-05 -0.122044826  5.485682e-01
##  [76,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
##  [77,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [78,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [79,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [80,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [81,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [82,] 0.0022935780 0.003623188 1.282941e-05 -0.371211247  6.447599e-01
##  [83,] 0.0018348624 0.003623188 1.282941e-05 -0.499279127  6.912086e-01
##  [84,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [85,] 0.0000000000 0.003623188 1.282941e-05 -1.011550647  8.441235e-01
##  [86,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
##  [87,] 0.0032258065 0.003623188 1.291116e-05 -0.110592440  5.440302e-01
##  [88,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [89,] 0.0159453303 0.003623188 1.265277e-05  3.464124045  2.659806e-04
##  [90,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [91,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [92,] 0.0082004556 0.003623188 1.265277e-05  1.286807231  9.908074e-02
##  [93,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [94,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [95,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [96,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
##  [97,] 0.0027522936 0.003623188 1.282941e-05 -0.243143367  5.960528e-01
##  [98,] 0.0018264840 0.003623188 1.271871e-05 -0.503796354  6.927978e-01
##  [99,] 0.0064220183 0.003623188 1.282941e-05  0.781399674  2.172837e-01
## [100,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
## [101,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [102,] 0.0041474654 0.003623188 1.291116e-05  0.145907672  4.419971e-01
## [103,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
## [104,] 0.0054919908 0.003623188 1.277763e-05  0.522803263  3.005556e-01
## [105,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [106,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [107,] 0.0289786223 0.003623188 1.267843e-05  7.120955420  5.359073e-13
## [108,] 0.0041002278 0.003623188 1.265277e-05  0.134110094  4.466578e-01
## [109,] 0.0197701149 0.003623188 1.287395e-05  4.500221706  3.394131e-06
## [110,] 0.0054919908 0.003623188 1.277763e-05  0.522803263  3.005556e-01
## [111,] 0.0182232346 0.003623188 1.265277e-05  4.104511344  2.025852e-05
## [112,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [113,] 0.0082379863 0.003623188 1.277763e-05  1.291003975  9.835115e-02
## [114,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [115,] 0.0054919908 0.003623188 1.277763e-05  0.522803263  3.005556e-01
## [116,] 0.0242562929 0.003623188 1.277763e-05  5.772174798  3.912744e-09
## [117,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [118,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [119,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [120,] 0.0036529680 0.003623188 1.271871e-05  0.008350216  4.966688e-01
## [121,] 0.0004566210 0.003623188 1.271871e-05 -0.887906281  8.127044e-01
## [122,] 0.0054794521 0.003623188 1.271871e-05  0.520496786  3.013587e-01
## [123,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [124,] 0.0119617225 0.003623188 1.242016e-05  2.366061489  8.989229e-03
## [125,] 0.0041002278 0.003623188 1.265277e-05  0.134110094  4.466578e-01
## [126,] 0.0073732719 0.003623188 1.291116e-05  1.043658065  1.483218e-01
## [127,] 0.0055045872 0.003623188 1.282941e-05  0.525263914  2.996999e-01
## [128,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [129,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [130,] 0.0018264840 0.003623188 1.271871e-05 -0.503796354  6.927978e-01
## [131,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [132,] 0.0118451025 0.003623188 1.265277e-05  2.311426908  1.040464e-02
## [133,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [134,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [135,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [136,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [137,] 0.0148491879 0.003623188 1.297775e-05  3.116200157  9.159894e-04
## [138,] 0.0000000000 0.003623188 1.282941e-05 -1.011550647  8.441235e-01
## [139,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [140,] 0.0036697248 0.003623188 1.282941e-05  0.012992394  4.948169e-01
## [141,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [142,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [143,] 0.0083526682 0.003623188 1.297775e-05  1.312845745  9.461747e-02
## [144,] 0.0000000000 0.003623188 1.277763e-05 -1.013598162  8.446127e-01
## [145,] 0.0148491879 0.003623188 1.297775e-05  3.116200157  9.159894e-04
## [146,] 0.0180522565 0.003623188 1.267843e-05  4.052336521  2.535433e-05
## [147,] 0.0027397260 0.003623188 1.271871e-05 -0.247723069  5.978257e-01
## [148,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [149,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [150,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [151,] 0.0022831050 0.003623188 1.271871e-05 -0.375759712  6.464522e-01
## [152,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [153,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [154,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [155,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [156,] 0.0004555809 0.003623188 1.265277e-05 -0.890509584  8.134038e-01
## [157,] 0.0064220183 0.003623188 1.282941e-05  0.781399674  2.172837e-01
## [158,] 0.0078341014 0.003623188 1.291116e-05  1.171908121  1.206170e-01
## [159,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [160,] 0.0000000000 0.003623188 1.282941e-05 -1.011550647  8.441235e-01
## [161,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [162,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [163,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [164,] 0.0064220183 0.003623188 1.282941e-05  0.781399674  2.172837e-01
## [165,] 0.0041002278 0.003623188 1.265277e-05  0.134110094  4.466578e-01
## [166,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [167,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [168,] 0.0059225513 0.003623188 1.265277e-05  0.646419932  2.590037e-01
## [169,] 0.0086560364 0.003623188 1.265277e-05  1.414884690  7.855115e-02
## [170,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [171,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
## [172,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [173,] 0.0096330275 0.003623188 1.282941e-05  1.677874835  4.668577e-02
## [174,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [175,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [176,] 0.0091743119 0.003623188 1.282941e-05  1.549806955  6.059393e-02
## [177,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [178,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [179,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [180,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [181,] 0.0110599078 0.003623188 1.291116e-05  2.069658514  1.924217e-02
## [182,] 0.0041002278 0.003623188 1.265277e-05  0.134110094  4.466578e-01
## [183,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [184,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [185,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [186,] 0.0077981651 0.003623188 1.282941e-05  1.165603314  1.218874e-01
## [187,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [188,] 0.0036697248 0.003623188 1.282941e-05  0.012992394  4.948169e-01
## [189,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [190,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [191,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
## [192,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
## [193,] 0.0054669704 0.003623188 1.265277e-05  0.518342473  3.021097e-01
## [194,] 0.0027649770 0.003623188 1.291116e-05 -0.238842496  5.943861e-01
## [195,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [196,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [197,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [198,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [199,] 0.0013761468 0.003623188 1.282941e-05 -0.627347007  7.347841e-01
## [200,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [201,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [202,] 0.0036529680 0.003623188 1.271871e-05  0.008350216  4.966688e-01
## [203,] 0.0055045872 0.003623188 1.282941e-05  0.525263914  2.996999e-01
## [204,] 0.0054669704 0.003623188 1.265277e-05  0.518342473  3.021097e-01
## [205,] 0.0022988506 0.003623188 1.287395e-05 -0.369098964  6.439730e-01
## [206,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [207,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [208,] 0.0091533181 0.003623188 1.277763e-05  1.547070879  6.092308e-02
## [209,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [210,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [211,] 0.0091116173 0.003623188 1.265277e-05  1.542962150  6.141998e-02
## [212,] 0.0050228311 0.003623188 1.271871e-05  0.392460143  3.473591e-01
## [213,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [214,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [215,] 0.0032183908 0.003623188 1.287395e-05 -0.112818929  5.449129e-01
## [216,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
## [217,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [218,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [219,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [220,] 0.0087557604 0.003623188 1.291116e-05  1.428408233  7.658720e-02
## [221,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
## [222,] 0.0027522936 0.003623188 1.282941e-05 -0.243143367  5.960528e-01
## [223,] 0.0331002331 0.003623188 1.298353e-05  8.180644436  1.411649e-16
## [224,] 0.0087155963 0.003623188 1.282941e-05  1.421739075  7.755101e-02
## [225,] 0.0045558087 0.003623188 1.265277e-05  0.262187553  3.965884e-01
## [226,] 0.0063781321 0.003623188 1.265277e-05  0.774497392  2.193184e-01
## [227,] 0.0036613272 0.003623188 1.277763e-05  0.010669454  4.957436e-01
## [228,] 0.0059633028 0.003623188 1.282941e-05  0.653331794  2.567712e-01
## [229,] 0.0004555809 0.003623188 1.265277e-05 -0.890509584  8.134038e-01
## [230,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [231,] 0.0214452214 0.003623188 1.298353e-05  4.946076411  3.786212e-07
## [232,] 0.0004587156 0.003623188 1.282941e-05 -0.883482767  8.115123e-01
## [233,] 0.0085918854 0.003623188 1.251538e-05  1.404495747  8.008564e-02
## [234,] 0.0293706294 0.003623188 1.298353e-05  7.145582668  4.480747e-13
## [235,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [236,] 0.0063781321 0.003623188 1.265277e-05  0.774497392  2.193184e-01
## [237,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [238,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
## [239,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [240,] 0.0013698630 0.003623188 1.271871e-05 -0.631832997  7.362520e-01
## [241,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [242,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [243,] 0.0091116173 0.003623188 1.265277e-05  1.542962150  6.141998e-02
## [244,] 0.0059225513 0.003623188 1.265277e-05  0.646419932  2.590037e-01
## [245,] 0.0022935780 0.003623188 1.282941e-05 -0.371211247  6.447599e-01
## [246,] 0.0009132420 0.003623188 1.271871e-05 -0.759869639  7.763337e-01
## [247,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [248,] 0.0018433180 0.003623188 1.291116e-05 -0.495342608  6.898209e-01
## [249,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [250,] 0.0148491879 0.003623188 1.297775e-05  3.116200157  9.159894e-04
## [251,] 0.0004587156 0.003623188 1.282941e-05 -0.883482767  8.115123e-01
## [252,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [253,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [254,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [255,] 0.0150341686 0.003623188 1.265277e-05  3.207969126  6.683792e-04
## [256,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [257,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [258,] 0.0129186603 0.003623188 1.242016e-05  2.637592870  4.174838e-03
## [259,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [260,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [261,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
## [262,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
## [263,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [264,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [265,] 0.0027397260 0.003623188 1.271871e-05 -0.247723069  5.978257e-01
## [266,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
## [267,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [268,] 0.0041002278 0.003623188 1.265277e-05  0.134110094  4.466578e-01
## [269,] 0.0091116173 0.003623188 1.265277e-05  1.542962150  6.141998e-02
## [270,] 0.0031890661 0.003623188 1.265277e-05 -0.122044826  5.485682e-01
## [271,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [272,] 0.0004555809 0.003623188 1.265277e-05 -0.890509584  8.134038e-01
## [273,] 0.0091116173 0.003623188 1.265277e-05  1.542962150  6.141998e-02
## [274,] 0.0152777778 0.003623188 1.296317e-05  3.236990420  6.039873e-04
## [275,] 0.0050228311 0.003623188 1.271871e-05  0.392460143  3.473591e-01
## [276,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
## [277,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## attr(,"cluster")
##   [1] Low  Low  High High Low  High Low  Low  Low  High High Low  Low  High Low 
##  [16] Low  High Low  Low  Low  Low  Low  Low  Low  High High High Low  High Low 
##  [31] Low  Low  Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [46] Low  High High Low  High Low  Low  Low  High High High High Low  Low  Low 
##  [61] Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [76] Low  Low  Low  Low  Low  Low  High High Low  High Low  High Low  Low  Low 
##  [91] Low  Low  Low  Low  Low  High High Low  High Low  Low  High Low  High Low 
## [106] Low  High Low  High High Low  Low  High Low  High High Low  Low  Low  Low 
## [121] Low  Low  Low  High Low  High High Low  Low  Low  Low  Low  Low  Low  Low 
## [136] Low  High High Low  High Low  Low  High High High High Low  Low  Low  Low 
## [151] Low  Low  Low  Low  Low  Low  High High Low  High Low  Low  Low  High Low 
## [166] Low  Low  Low  Low  Low  Low  High High Low  Low  High Low  Low  Low  Low 
## [181] High Low  Low  Low  Low  High Low  High High Low  Low  Low  Low  High Low 
## [196] Low  Low  Low  High Low  Low  Low  High Low  High Low  Low  High Low  Low 
## [211] Low  Low  Low  Low  High Low  Low  Low  Low  High Low  High High High Low 
## [226] Low  High High Low  Low  High High High High Low  Low  Low  Low  Low  Low 
## [241] Low  Low  Low  Low  High Low  Low  High Low  High High Low  Low  Low  Low 
## [256] Low  Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
## [271] Low  Low  Low  High Low  Low  Low 
## Levels: Low High
## attr(,"gstari")
## [1] FALSE
## attr(,"call")
## localG(x = data_merge$KEK, listw = W.knn, alternative = "greater")
## attr(,"class")
## [1] "localG"

data_merge$lmZ2 <- Local_M
tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ2",
title = "Local Gi", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Coldspot", "No SAC", "Hotspot"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Indeks Moran Global

W.rook <- nb2listw(rook, style='W')

moran(data_merge$KEK, W.rook, n=length(W.rook$neighbours), S0=Szero(W.rook))

## $I
## [1] 0.3487822
## 
## $K
## [1] 18.9482

IM <- moran.test(data_merge$KEK, W.rook,randomisation=T,
alternative="greater", zero.policy = FALSE)
IM

## 
##  Moran I test under randomisation
## 
## data:  data_merge$KEK  
## weights: W.rook    
## 
## Moran I statistic standard deviate = 9.6447, p-value < 2.2e-16
## alternative hypothesis: greater
## sample estimates:
## Moran I statistic       Expectation          Variance 
##       0.348782225      -0.003623188       0.001335086

Indeks Moran Lokal

LM <- localmoran(data_merge$KEK, W.rook, alternative = "greater")
head(LM)

##            Ii          E.Ii      Var.Ii       Z.Ii Pr(z > E(Ii))
## 1  0.02695691 -7.376709e-04 0.024873336  0.1756013    0.43030357
## 2 -0.07902723 -7.376709e-04 0.040242830 -0.3902650    0.65182970
## 3  1.38934374 -8.612240e-03 0.584811161  1.8280399    0.03377178
## 4 -0.03650774 -5.062083e-05 0.001356238 -0.9899526    0.83890135
## 5  0.07513153 -7.376709e-04 0.067566375  0.2918772    0.38519026
## 6 -0.08280637 -3.425436e-03 0.132137772 -0.2183749    0.58643149

data_merge$lmI <- LM[, "Ii"] # local Moran's I
data_merge$lmZ <- LM[, "Z.Ii"] # z-scores
# p-values corresponding to alternative greater
data_merge$lmp <- LM[, "Pr(z > E(Ii))"]

p1 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p1

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

p2 <- tm_shape(data_merge) +
  tm_polygons(col = "lmp", title = "p-value",
              breaks = c(-Inf, 0.05, Inf)) +
  tm_layout(legend.outside = TRUE)
p2

Klasifikasi Warna :

• Less than 0.05 = Ada Autokorelasi Positif

• 0.05 or more = Tidak Ada Autokorelasi Positif

tmap_mode("view")

## tmap mode set to interactive viewing

p3 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p3

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ",
title = "Local Moran's I", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Negative SAC", "No SAC", "Positive SAC"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Rentang nilai untuk klasifikasi warna, dengan:

• -Inf hingga -1.96 mewakili Negative Spatial Autocorrelation (SAC).

• -1.96 hingga 1.96 mewakili No Spatial Autocorrelation (No SAC).

• 1.96 hingga Inf mewakili Positive Spatial Autocorrelation (Positive SAC)

# Autokorelasi spasial positif (Z-score > 1.96)
positive_autocorr <- data_merge[data_merge$lmZ > 1.96, ]
# Autokorelasi spasial negatif (Z-score < -1.96)
negative_autocorr <- data_merge[data_merge$lmZ < -1.96, ]
# Area tanpa autokorelasi spasial (Z-score antara -1.96 dan 1.96)
no_autocorr <- data_merge[data_merge$lmZ >= -1.96 & data_merge$lmZ <= 1.96, ]

moran.plot(data_merge$KEK, W.rook)

Moran Plot memiliki kemiringan positif: ini menunjukkan adanya autokorelasi spasial positif, di mana lokasi dengan nilai tinggi cenderung berdekatan dengan lokasi bernilai tinggi, dan lokasi bernilai rendah cenderung berdekatan dengan lokasi bernilai rendah.

Getis-Ord Gi

Mendeteksi adanya wilayah yang hotspot/coldspot.

Local_M <- localG(data_merge$KEK, W.rook, alternative = "greater")
Local_M

##   [1] -0.175601343  0.390265013  1.828039900 -0.989952578 -0.291877205
##   [6] -0.218374907 -0.480581224 -0.786098730 -0.766467085  0.533996866
##  [11] -0.173142404 -0.291877205 -0.786098730  0.523902035 -1.117870599
##  [16] -1.018587043  0.763863545 -0.909375678 -0.634354664 -0.127136697
##  [21] -0.786098730 -0.621358222 -0.291877205 -0.621358222  0.550882783
##  [26]  9.077185239 -0.449674614  2.379361417  5.167857453 -0.634354664
##  [31] -0.456617713 -0.786098730 -0.061094178 -0.760172919 -1.295613270
##  [36]  1.809281336 -0.640674601  1.160679998 -1.117870599 -1.179740170
##  [41] -1.018587043  2.396164542 -0.649332580 -0.623512709 -0.291877205
##  [46] -1.295613270 -0.117074734  0.526114228 -0.766444193 -0.173142404
##  [51] -0.766467085 -0.456617713 -0.291877205  3.683673298  3.989061729
##  [56]  0.809666782  0.269128154  1.224819495 -0.415063571 -0.250122285
##  [61] -0.454682106  1.275834530 -1.018587043 -0.576777591 -0.532198076
##  [66] -1.091974737 -0.506277205 -1.209680475 -1.088328523  4.504585626
##  [71] -0.291877205 -1.117870599 -0.649332580 -1.209680475  0.037603812
##  [76] -0.378199745 -0.291877205 -0.786098730 -0.909375678 -0.909375678
##  [81] -0.297929066  0.526114228 -0.032104298 -1.018587043 -0.903093714
##  [86] -0.623512709  2.248500327 -0.506277205 -0.337649739 -1.018587043
##  [91] -0.532198076  1.341953999 -0.684697673 -1.117870599 -1.117870599
##  [96]  0.545946531 -0.223575372 -0.049699455  0.240272639 -0.043344556
## [101] -0.415063571  0.434384744 -0.340505568  0.523902035 -1.117870599
## [106] -0.909375678  3.558294459  0.390265013  2.681923891  0.316777802
## [111]  3.136194059  0.006032634 -1.013598162 -0.194718255  0.381019662
## [116]  7.552579583  0.234076200 -0.342333802 -0.775093022  2.056936495
## [121] -0.880703586  0.407294521 -1.117870599  2.100022453  2.949774409
## [126]  0.817259344 -0.115075486 -1.117870599 -0.909375678  3.003766040
## [131] -0.909375678  1.414884690 -0.992386748 -0.180794562 -1.018587043
## [136] -1.209680475  2.608199220 -1.318391544 -0.291877205 -0.166735022
## [141] -0.337649739 -0.786098730  0.948707013 -0.904921697  3.034539714
## [146]  2.069378039  0.378958346 -1.295613270 -0.883739885 -0.786098730
## [151] -0.631832997 -0.337649739 -0.623512709 -0.634354664 -1.117870599
## [156] -0.873805045  0.383193434  1.389755737 -0.883739885 -0.903093714
## [161] -0.180794562 -0.766467085 -0.992386748  0.397196034  0.164911686
## [166] -0.532198076 -0.775093022  0.377007684 -0.798640360 -1.018587043
## [171]  0.377007684 -0.045568949  0.042972570 -0.506277205 -0.250122285
## [176]  1.549806955 -0.786098730 -1.117870599 -0.480581224 -0.786098730
## [181]  1.776111383  0.696565845 -0.480581224 -0.762432124 -1.018587043
## [186]  0.240272639 -0.883601590 -0.371211247  0.397196034 -0.786098730
## [191] -0.480581224 -0.786098730 -0.340505568 -0.327168183 -1.117870599
## [196] -0.909375678 -0.909375678 -1.018587043 -0.120216783 -0.766467085
## [201] -0.506277205  2.665133284  0.012992394  0.662870654 -1.165960889
## [206] -0.291877205 -1.117870599  1.162970523 -1.295613270 -0.623512709
## [211] -0.883601590  3.236677019  0.006032634 -0.449152431  0.752001323
## [216]  1.810492018  0.518342473 -0.909375678 -0.786098730  1.428408233
## [221] -0.532198076 -0.114936635  5.851755458  0.397196034  0.519939169
## [226]  0.170608952  0.041170976  0.653331794 -0.786098730 -0.291877205
## [231]  4.946076411 -0.370358323  0.540729809  3.781631922 -0.051786770
## [236]  0.390265013 -0.291877205 -0.415063571 -0.909375678 -0.646580053
## [241] -1.117870599 -0.786098730 -0.762432124 -0.194718255  0.653331794
## [246] -0.619370113 -0.291877205 -0.118342671 -0.909375678  1.874751242
## [251] -0.875407794 -0.909375678 -0.634354664 -1.117870599  5.236678165
## [256] -1.209680475 -0.909375678  4.191974305 -1.018587043 -0.480581224
## [261] -0.762432124 -0.766467085 -0.909375678 -0.318843435 -0.619370113
## [266] -0.480581224 -1.018587043 -0.063660057  2.396164542 -0.480581224
## [271] -0.634354664  1.224819495  1.508495334  3.192813251  0.521844279
## [276] -0.194718255 -0.992386748
## attr(,"internals")
##                  Gi       E(Gi)        V(Gi)        Z(Gi) Pr(z > E(Gi))
##   [1,] 0.0031321185 0.003623188 7.820438e-06 -0.175601343  5.696964e-01
##   [2,] 0.0050113895 0.003623188 1.265277e-05  0.390265013  3.481703e-01
##   [3,] 0.0109953704 0.003623188 1.626375e-05  1.828039900  3.377178e-02
##   [4,] 0.0011441648 0.003623188 6.270939e-06 -0.989952578  8.389014e-01
##   [5,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##   [6,] 0.0029624753 0.003623188 9.154199e-06 -0.218374907  5.864315e-01
##   [7,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
##   [8,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##   [9,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
##  [10,] 0.0052304675 0.003623188 9.059519e-06  0.533996866  2.966719e-01
##  [11,] 0.0030581040 0.003623188 1.065172e-05 -0.173142404  5.687303e-01
##  [12,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [13,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [14,] 0.0057208238 0.003623188 1.603097e-05  0.523902035  3.001733e-01
##  [15,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [16,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [17,] 0.0061162080 0.003623188 1.065172e-05  0.763863545  2.224743e-01
##  [18,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
##  [19,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [20,] 0.0030372058 0.003623188 2.124358e-05 -0.127136697  5.505839e-01
##  [21,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [22,] 0.0007593014 0.003623188 2.124358e-05 -0.621358222  7.328180e-01
##  [23,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [24,] 0.0007593014 0.003623188 2.124358e-05 -0.621358222  7.328180e-01
##  [25,] 0.0056074766 0.003623188 1.297451e-05  0.550882783  2.908570e-01
##  [26,] 0.0254291163 0.003623188 5.770943e-06  9.077185239  5.570944e-20
##  [27,] 0.0015325670 0.003623188 2.161494e-05 -0.449674614  6.735275e-01
##  [28,] 0.0131278539 0.003623188 1.595706e-05  2.379361417  8.671331e-03
##  [29,] 0.0205905206 0.003623188 1.077969e-05  5.167857453  1.183964e-07
##  [30,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [31,] 0.0015186029 0.003623188 2.124358e-05 -0.456617713  6.760271e-01
##  [32,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [33,] 0.0034246575 0.003623188 1.055982e-05 -0.061094178  5.243579e-01
##  [34,] 0.0005733945 0.003623188 1.609593e-05 -0.760172919  7.764244e-01
##  [35,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
##  [36,] 0.0101149425 0.003623188 1.287395e-05  1.809281336  3.520365e-02
##  [37,] 0.0000000000 0.003623188 3.198209e-05 -0.640674601  7.391329e-01
##  [38,] 0.0077625571 0.003623188 1.271871e-05  1.160679998  1.228860e-01
##  [39,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [40,] 0.0006833713 0.003623188 6.209661e-06 -1.179740170  8.809482e-01
##  [41,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [42,] 0.0113895216 0.003623188 1.050507e-05  2.396164542  8.283826e-03
##  [43,] 0.0015186029 0.003623188 1.050507e-05 -0.649332580  7.419383e-01
##  [44,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
##  [45,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [46,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
##  [47,] 0.0033129460 0.003623188 7.022246e-06 -0.117074734  5.465996e-01
##  [48,] 0.0057339450 0.003623188 1.609593e-05  0.526114228  2.994044e-01
##  [49,] 0.0005694761 0.003623188 1.587432e-05 -0.766444193  7.782940e-01
##  [50,] 0.0030581040 0.003623188 1.065172e-05 -0.173142404  5.687303e-01
##  [51,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
##  [52,] 0.0015186029 0.003623188 2.124358e-05 -0.456617713  6.760271e-01
##  [53,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [54,] 0.0147106898 0.003623188 9.059519e-06  3.683673298  1.149484e-04
##  [55,] 0.0192307692 0.003623188 1.530840e-05  3.989061729  3.316758e-05
##  [56,] 0.0068649886 0.003623188 1.603097e-05  0.809666782  2.090659e-01
##  [57,] 0.0045871560 0.003623188 1.282941e-05  0.269128154  3.939155e-01
##  [58,] 0.0075930144 0.003623188 1.050507e-05  1.224819495  1.103216e-01
##  [59,] 0.0022779043 0.003623188 1.050507e-05 -0.415063571  6.609523e-01
##  [60,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
##  [61,] 0.0015220700 0.003623188 2.135430e-05 -0.454682106  6.753310e-01
##  [62,] 0.0059014814 0.003623188 3.188823e-06  1.275834530  1.010070e-01
##  [63,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [64,] 0.0020022883 0.003623188 7.897611e-06 -0.576777591  7.179551e-01
##  [65,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
##  [66,] 0.0005694761 0.003623188 7.820438e-06 -1.091974737  8.625779e-01
##  [67,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [68,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
##  [69,] 0.0009111617 0.003623188 6.209661e-06 -1.088328523  8.617749e-01
##  [70,] 0.0182232346 0.003623188 1.050507e-05  4.504585626  3.325129e-06
##  [71,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [72,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [73,] 0.0015186029 0.003623188 1.050507e-05 -0.649332580  7.419383e-01
##  [74,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
##  [75,] 0.0037965072 0.003623188 2.124358e-05  0.037603812  4.850018e-01
##  [76,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
##  [77,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [78,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [79,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
##  [80,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
##  [81,] 0.0026575550 0.003623188 1.050507e-05 -0.297929066  6.171213e-01
##  [82,] 0.0057339450 0.003623188 1.609593e-05  0.526114228  2.994044e-01
##  [83,] 0.0034403670 0.003623188 3.242857e-05 -0.032104298  5.128056e-01
##  [84,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [85,] 0.0000000000 0.003623188 1.609593e-05 -0.903093714  8.167619e-01
##  [86,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
##  [87,] 0.0126728111 0.003623188 1.619851e-05  2.248500327  1.227215e-02
##  [88,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [89,] 0.0022779043 0.003623188 1.587432e-05 -0.337649739  6.321864e-01
##  [90,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [91,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
##  [92,] 0.0079726651 0.003623188 1.050507e-05  1.341953999  8.980546e-02
##  [93,] 0.0017084282 0.003623188 7.820438e-06 -0.684697673  7.532326e-01
##  [94,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [95,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [96,] 0.0051605505 0.003623188 7.929615e-06  0.545946531  2.925513e-01
##  [97,] 0.0029488860 0.003623188 9.096232e-06 -0.223575372  5.884561e-01
##  [98,] 0.0034246575 0.003623188 1.595706e-05 -0.049699455  5.198191e-01
##  [99,] 0.0045871560 0.003623188 1.609593e-05  0.240272639  4.050595e-01
## [100,] 0.0035203976 0.003623188 5.623924e-06 -0.043344556  5.172866e-01
## [101,] 0.0022779043 0.003623188 1.050507e-05 -0.415063571  6.609523e-01
## [102,] 0.0049374589 0.003623188 9.154199e-06  0.434384744  3.320045e-01
## [103,] 0.0026033192 0.003623188 8.970993e-06 -0.340505568  6.332621e-01
## [104,] 0.0057208238 0.003623188 1.603097e-05  0.523902035  3.001733e-01
## [105,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [106,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [107,] 0.0178147268 0.003623188 1.590652e-05  3.558294459  1.866354e-04
## [108,] 0.0050113895 0.003623188 1.265277e-05  0.390265013  3.481703e-01
## [109,] 0.0160919540 0.003623188 2.161494e-05  2.681923891  3.660006e-03
## [110,] 0.0045766590 0.003623188 9.059519e-06  0.316777802  3.757061e-01
## [111,] 0.0130165962 0.003623188 8.970993e-06  3.136194059  8.557795e-04
## [112,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [113,] 0.0000000000 0.003623188 1.277763e-05 -1.013598162  8.446127e-01
## [114,] 0.0028473804 0.003623188 1.587432e-05 -0.194718255  5.771932e-01
## [115,] 0.0051487414 0.003623188 1.603097e-05  0.381019662  3.515943e-01
## [116,] 0.0282227307 0.003623188 1.060873e-05  7.552579583  2.133604e-14
## [117,] 0.0045558087 0.003623188 1.587432e-05  0.234076200  4.074629e-01
## [118,] 0.0028473804 0.003623188 5.135810e-06 -0.342333802  6.339501e-01
## [119,] 0.0013016596 0.003623188 8.970993e-06 -0.775093022  7.808577e-01
## [120,] 0.0109589041 0.003623188 1.271871e-05  2.056936495  1.984617e-02
## [121,] 0.0009784736 0.003623188 9.017750e-06 -0.880703586  8.107609e-01
## [122,] 0.0049467275 0.003623188 1.055982e-05  0.407294521  3.418958e-01
## [123,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [124,] 0.0103668262 0.003623188 1.031195e-05  2.100022453  1.786343e-02
## [125,] 0.0153758542 0.003623188 1.587432e-05  2.949774409  1.590030e-03
## [126,] 0.0069124424 0.003623188 1.619851e-05  0.817259344  2.068901e-01
## [127,] 0.0032110092 0.003623188 1.282941e-05 -0.115075486  5.458074e-01
## [128,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [129,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [130,] 0.0175038052 0.003623188 2.135430e-05  3.003766040  1.333301e-03
## [131,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [132,] 0.0086560364 0.003623188 1.265277e-05  1.414884690  7.855115e-02
## [133,] 0.0006508298 0.003623188 8.970993e-06 -0.992386748  8.394955e-01
## [134,] 0.0030372058 0.003623188 1.050507e-05 -0.180794562  5.717356e-01
## [135,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [136,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
## [137,] 0.0185614849 0.003623188 3.280353e-05  2.608199220  4.550998e-03
## [138,] 0.0014958117 0.003623188 2.603754e-06 -1.318391544  9.063137e-01
## [139,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [140,] 0.0031536697 0.003623188 7.929615e-06 -0.166735022  5.662107e-01
## [141,] 0.0022779043 0.003623188 1.587432e-05 -0.337649739  6.321864e-01
## [142,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [143,] 0.0057112261 0.003623188 4.844093e-06  0.948707013  1.713848e-01
## [144,] 0.0000000000 0.003623188 1.603097e-05 -0.904921697  8.172466e-01
## [145,] 0.0177880897 0.003623188 2.178921e-05  3.034539714  1.204515e-03
## [146,] 0.0118764846 0.003623188 1.590652e-05  2.069378039  1.925531e-02
## [147,] 0.0051369863 0.003623188 1.595706e-05  0.378958346  3.523594e-01
## [148,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
## [149,] 0.0009762447 0.003623188 8.970993e-06 -0.883739885  8.115817e-01
## [150,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [151,] 0.0013698630 0.003623188 1.271871e-05 -0.631832997  7.362520e-01
## [152,] 0.0022779043 0.003623188 1.587432e-05 -0.337649739  6.321864e-01
## [153,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
## [154,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [155,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [156,] 0.0020248038 0.003623188 3.346058e-06 -0.873805045  8.088878e-01
## [157,] 0.0051605505 0.003623188 1.609593e-05  0.383193434  3.507882e-01
## [158,] 0.0092165899 0.003623188 1.619851e-05  1.389755737  8.230153e-02
## [159,] 0.0009762447 0.003623188 8.970993e-06 -0.883739885  8.115817e-01
## [160,] 0.0000000000 0.003623188 1.609593e-05 -0.903093714  8.167619e-01
## [161,] 0.0030372058 0.003623188 1.050507e-05 -0.180794562  5.717356e-01
## [162,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
## [163,] 0.0006508298 0.003623188 8.970993e-06 -0.992386748  8.394955e-01
## [164,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [165,] 0.0045558087 0.003623188 3.198209e-05  0.164911686  4.345067e-01
## [166,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
## [167,] 0.0013016596 0.003623188 8.970993e-06 -0.775093022  7.808577e-01
## [168,] 0.0051252847 0.003623188 1.587432e-05  0.377007684  3.530839e-01
## [169,] 0.0023728170 0.003623188 2.451182e-06 -0.798640360  7.877505e-01
## [170,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [171,] 0.0051252847 0.003623188 1.587432e-05  0.377007684  3.530839e-01
## [172,] 0.0034403670 0.003623188 1.609593e-05 -0.045568949  5.181731e-01
## [173,] 0.0038226300 0.003623188 2.154015e-05  0.042972570  4.828617e-01
## [174,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [175,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
## [176,] 0.0091743119 0.003623188 1.282941e-05  1.549806955  6.059393e-02
## [177,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [178,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [179,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [180,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [181,] 0.0086405530 0.003623188 7.980148e-06  1.776111383  3.785728e-02
## [182,] 0.0068337130 0.003623188 2.124358e-05  0.696565845  2.430373e-01
## [183,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [184,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [185,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [186,] 0.0045871560 0.003623188 1.609593e-05  0.240272639  4.050595e-01
## [187,] 0.0007593014 0.003623188 1.050507e-05 -0.883601590  8.115443e-01
## [188,] 0.0022935780 0.003623188 1.282941e-05 -0.371211247  6.447599e-01
## [189,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [190,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [191,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [192,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [193,] 0.0026033192 0.003623188 8.970993e-06 -0.340505568  6.332621e-01
## [194,] 0.0026333114 0.003623188 9.154199e-06 -0.327168183  6.282297e-01
## [195,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [196,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [197,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [198,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [199,] 0.0033361134 0.003623188 5.702437e-06 -0.120216783  5.478443e-01
## [200,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
## [201,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [202,] 0.0142694064 0.003623188 1.595706e-05  2.665133284  3.847891e-03
## [203,] 0.0036697248 0.003623188 1.282941e-05  0.012992394  4.948169e-01
## [204,] 0.0062642369 0.003623188 1.587432e-05  0.662870654  2.537067e-01
## [205,] 0.0009578544 0.003623188 5.225590e-06 -1.165960889  8.781849e-01
## [206,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [207,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [208,] 0.0077803204 0.003623188 1.277763e-05  1.162970523  1.224207e-01
## [209,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
## [210,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
## [211,] 0.0007593014 0.003623188 1.050507e-05 -0.883601590  8.115443e-01
## [212,] 0.0165525114 0.003623188 1.595706e-05  3.236677019  6.046510e-04
## [213,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [214,] 0.0022779043 0.003623188 8.970993e-06 -0.449152431  6.733392e-01
## [215,] 0.0056194125 0.003623188 7.046629e-06  0.752001323  2.260251e-01
## [216,] 0.0094912680 0.003623188 1.050507e-05  1.810492018  3.510976e-02
## [217,] 0.0054669704 0.003623188 1.265277e-05  0.518342473  3.021097e-01
## [218,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [219,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [220,] 0.0087557604 0.003623188 1.291116e-05  1.428408233  7.658720e-02
## [221,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
## [222,] 0.0032765400 0.003623188 9.096232e-06 -0.114936635  5.457523e-01
## [223,] 0.0247086247 0.003623188 1.298353e-05  5.851755458  2.432057e-09
## [224,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [225,] 0.0056947608 0.003623188 1.587432e-05  0.519939169  3.015530e-01
## [226,] 0.0041761579 0.003623188 1.050507e-05  0.170608952  4.322656e-01
## [227,] 0.0038138825 0.003623188 2.145321e-05  0.041170976  4.835798e-01
## [228,] 0.0059633028 0.003623188 1.282941e-05  0.653331794  2.567712e-01
## [229,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [230,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [231,] 0.0214452214 0.003623188 1.298353e-05  4.946076411  3.786212e-07
## [232,] 0.0025802752 0.003623188 7.929615e-06 -0.370358323  6.444422e-01
## [233,] 0.0050384513 0.003623188 6.850364e-06  0.540729809  2.943469e-01
## [234,] 0.0172494172 0.003623188 1.298353e-05  3.781631922  7.790181e-05
## [235,] 0.0034168565 0.003623188 1.587432e-05 -0.051786770  5.206507e-01
## [236,] 0.0050113895 0.003623188 1.265277e-05  0.390265013  3.481703e-01
## [237,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [238,] 0.0022779043 0.003623188 1.050507e-05 -0.415063571  6.609523e-01
## [239,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [240,] 0.0015220700 0.003623188 1.055982e-05 -0.646580053  7.410481e-01
## [241,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [242,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [243,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [244,] 0.0028473804 0.003623188 1.587432e-05 -0.194718255  5.771932e-01
## [245,] 0.0059633028 0.003623188 1.282941e-05  0.653331794  2.567712e-01
## [246,] 0.0007610350 0.003623188 2.135430e-05 -0.619370113  7.321637e-01
## [247,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [248,] 0.0030721966 0.003623188 2.167742e-05 -0.118342671  5.471019e-01
## [249,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [250,] 0.0123743233 0.003623188 2.178921e-05  1.874751242  3.041348e-02
## [251,] 0.0009829620 0.003623188 9.096232e-06 -0.875407794  8.093240e-01
## [252,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [253,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [254,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [255,] 0.0244874715 0.003623188 1.587432e-05  5.236678165  8.174621e-08
## [256,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
## [257,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [258,] 0.0160628845 0.003623188 8.806074e-06  4.191974305  1.382687e-05
## [259,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [260,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [261,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [262,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
## [263,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [264,] 0.0027841053 0.003623188 6.925562e-06 -0.318843435  6.250774e-01
## [265,] 0.0007610350 0.003623188 2.135430e-05 -0.619370113  7.321637e-01
## [266,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [267,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [268,] 0.0034168565 0.003623188 1.050507e-05 -0.063660057  5.253795e-01
## [269,] 0.0113895216 0.003623188 1.050507e-05  2.396164542  8.283826e-03
## [270,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [271,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [272,] 0.0075930144 0.003623188 1.050507e-05  1.224819495  1.103216e-01
## [273,] 0.0075930144 0.003623188 6.925562e-06  1.508495334  6.571390e-02
## [274,] 0.0185185185 0.003623188 2.176473e-05  3.192813251  7.044702e-04
## [275,] 0.0057077626 0.003623188 1.595706e-05  0.521844279  3.008894e-01
## [276,] 0.0028473804 0.003623188 1.587432e-05 -0.194718255  5.771932e-01
## [277,] 0.0006508298 0.003623188 8.970993e-06 -0.992386748  8.394955e-01
## attr(,"cluster")
##   [1] Low  Low  High High Low  High Low  Low  Low  High High Low  Low  High Low 
##  [16] Low  High Low  Low  Low  Low  Low  Low  Low  High High High Low  High Low 
##  [31] Low  Low  Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [46] Low  High High Low  High Low  Low  Low  High High High High Low  Low  Low 
##  [61] Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [76] Low  Low  Low  Low  Low  Low  High High Low  High Low  High Low  Low  Low 
##  [91] Low  Low  Low  Low  Low  High High Low  High Low  Low  High Low  High Low 
## [106] Low  High Low  High High Low  Low  High Low  High High Low  Low  Low  Low 
## [121] Low  Low  Low  High Low  High High Low  Low  Low  Low  Low  Low  Low  Low 
## [136] Low  High High Low  High Low  Low  High High High High Low  Low  Low  Low 
## [151] Low  Low  Low  Low  Low  Low  High High Low  High Low  Low  Low  High Low 
## [166] Low  Low  Low  Low  Low  Low  High High Low  Low  High Low  Low  Low  Low 
## [181] High Low  Low  Low  Low  High Low  High High Low  Low  Low  Low  High Low 
## [196] Low  Low  Low  High Low  Low  Low  High Low  High Low  Low  High Low  Low 
## [211] Low  Low  Low  Low  High Low  Low  Low  Low  High Low  High High High Low 
## [226] Low  High High Low  Low  High High High High Low  Low  Low  Low  Low  Low 
## [241] Low  Low  Low  Low  High Low  Low  High Low  High High Low  Low  Low  Low 
## [256] Low  Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
## [271] Low  Low  Low  High Low  Low  Low 
## Levels: Low High
## attr(,"gstari")
## [1] FALSE
## attr(,"call")
## localG(x = data_merge$KEK, listw = W.rook, alternative = "greater")
## attr(,"class")
## [1] "localG"

data_merge$lmZ2 <- Local_M
tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ2",
title = "Local Gi", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Coldspot", "No SAC", "Hotspot"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Indeks Moran Global

W.queen <- nb2listw(queen, style='W')
W.queen2 <- nb2listw(queen, style='B')

moran(data_merge$KEK, W.queen, n=length(W.queen$neighbours), S0=Szero(W.queen))

## $I
## [1] 0.3491721
## 
## $K
## [1] 18.9482

IM <- moran.test(data_merge$KEK, W.queen,randomisation=T,
alternative="greater", zero.policy = FALSE)
IM

## 
##  Moran I test under randomisation
## 
## data:  data_merge$KEK  
## weights: W.queen    
## 
## Moran I statistic standard deviate = 9.6888, p-value < 2.2e-16
## alternative hypothesis: greater
## sample estimates:
## Moran I statistic       Expectation          Variance 
##       0.349172064      -0.003623188       0.001325887

Indeks Moran Lokal

LM <- localmoran(data_merge$KEK, W.queen, alternative = "greater")
head(LM)

##            Ii          E.Ii      Var.Ii       Z.Ii Pr(z > E(Ii))
## 1  0.02695691 -7.376709e-04 0.024873336  0.1756013    0.43030357
## 2 -0.07902723 -7.376709e-04 0.040242830 -0.3902650    0.65182970
## 3  1.38934374 -8.612240e-03 0.584811161  1.8280399    0.03377178
## 4 -0.03803741 -5.062083e-05 0.001228309 -1.0838741    0.86078969
## 5  0.07513153 -7.376709e-04 0.067566375  0.2918772    0.38519026
## 6 -0.08280637 -3.425436e-03 0.132137772 -0.2183749    0.58643149

data_merge$lmI <- LM[, "Ii"] # local Moran's I
data_merge$lmZ <- LM[, "Z.Ii"] # z-scores
# p-values corresponding to alternative greater
data_merge$lmp <- LM[, "Pr(z > E(Ii))"]

p1 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p1

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

p2 <- tm_shape(data_merge) +
  tm_polygons(col = "lmp", title = "p-value",
              breaks = c(-Inf, 0.05, Inf)) +
  tm_layout(legend.outside = TRUE)
p2

Klasifikasi Warna :

• Less than 0.05 = Ada Autokorelasi Positif

• 0.05 or more = Tidak Ada Autokorelasi Positif

tmap_mode("view")

## tmap mode set to interactive viewing

p3 <- tm_shape(data_merge) +
  tm_polygons(col = "lmI", title = "Local Moran's I",
              style = "quantile") +
  tm_layout(legend.outside = TRUE)
p3

## Variable(s) "lmI" contains positive and negative values, so midpoint is set to 0. Set midpoint = NA to show the full spectrum of the color palette.

tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ",
title = "Local Moran's I", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Negative SAC", "No SAC", "Positive SAC"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)

Rentang nilai untuk klasifikasi warna, dengan:

• -Inf hingga -1.96 mewakili Negative Spatial Autocorrelation (SAC).

• -1.96 hingga 1.96 mewakili No Spatial Autocorrelation (No SAC).

• 1.96 hingga Inf mewakili Positive Spatial Autocorrelation (Positive SAC)

# Autokorelasi spasial positif (Z-score > 1.96)
positive_autocorr <- data_merge[data_merge$lmZ > 1.96, ]
# Autokorelasi spasial negatif (Z-score < -1.96)
negative_autocorr <- data_merge[data_merge$lmZ < -1.96, ]
# Area tanpa autokorelasi spasial (Z-score antara -1.96 dan 1.96)
no_autocorr <- data_merge[data_merge$lmZ >= -1.96 & data_merge$lmZ <= 1.96, ]

moran.plot(data_merge$KEK, W.queen)

Moran Plot memiliki kemiringan positif: ini menunjukkan adanya autokorelasi spasial positif, di mana lokasi dengan nilai tinggi cenderung berdekatan dengan lokasi bernilai tinggi, dan lokasi bernilai rendah cenderung berdekatan dengan lokasi bernilai rendah.

Getis-Ord Gi

Mendeteksi adanya wilayah yang hotspot/coldspot.

Local_M <- localG(data_merge$KEK, W.queen, alternative = "greater")
Local_M

##   [1] -0.175601343  0.390265013  1.828039900 -1.083874118 -0.291877205
##   [6] -0.218374907 -0.480581224 -0.786098730 -0.766467085  0.533996866
##  [11] -0.173142404 -0.291877205 -0.786098730  0.523902035 -1.117870599
##  [16] -1.018587043  0.763863545 -0.909375678 -0.634354664 -0.127136697
##  [21] -0.786098730 -0.621358222 -0.291877205 -0.621358222  0.550882783
##  [26]  9.077185239  0.099479164  2.379361417  5.167857453 -0.634354664
##  [31] -0.456617713 -0.786098730 -0.061094178 -0.760172919 -1.295613270
##  [36]  1.809281336 -0.640674601  1.160679998 -1.117870599 -1.179740170
##  [41] -1.018587043  2.396164542 -0.649332580 -0.623512709 -0.291877205
##  [46] -1.295613270 -0.117074734  0.526114228 -0.766444193 -0.173142404
##  [51] -0.766467085 -0.456617713 -0.291877205  3.683673298  3.989061729
##  [56]  0.809666782  0.269128154  1.224819495 -0.557799295 -0.250122285
##  [61] -0.454682106  1.275834530 -1.018587043 -0.576777591 -0.532198076
##  [66] -0.895897783 -0.506277205 -1.209680475 -1.088328523  4.504585626
##  [71] -0.291877205 -1.117870599 -0.649332580 -1.209680475  0.037603812
##  [76] -0.378199745 -0.291877205 -0.786098730 -0.909375678 -0.909375678
##  [81] -0.297929066  0.526114228 -0.032104298 -1.018587043 -0.903093714
##  [86] -0.623512709  2.326158626 -0.506277205 -0.337649739 -1.018587043
##  [91] -0.532198076  1.341953999 -0.799722058 -1.117870599 -1.117870599
##  [96]  0.545946531 -0.223575372 -0.049699455  0.240272639 -0.043344556
## [101] -0.415063571  0.434384744 -0.340505568  0.523902035 -1.117870599
## [106] -0.909375678  3.558294459  0.390265013  2.681923891  0.316777802
## [111]  3.136194059  0.006032634 -1.013598162 -0.194718255  0.381019662
## [116]  7.552579583  0.234076200 -0.342333802 -0.775093022  2.056936495
## [121] -0.880703586  0.407294521 -1.117870599  2.100022453  2.949774409
## [126]  0.817259344 -0.115075486 -1.117870599 -0.909375678  3.003766040
## [131] -0.909375678  1.414884690 -0.992386748 -0.180794562 -1.018587043
## [136] -1.209680475  2.608199220 -1.318391544 -0.291877205 -0.309411796
## [141] -0.337649739 -0.786098730  0.948707013 -0.904921697  3.034539714
## [146]  2.069378039  0.378958346 -1.295613270 -0.607370609 -0.786098730
## [151] -0.631832997 -0.337649739 -0.623512709 -0.634354664 -1.117870599
## [156] -0.873805045  0.383193434  1.389755737 -0.883739885 -0.903093714
## [161] -0.180794562 -0.766467085 -0.992386748  0.397196034  0.164911686
## [166] -0.532198076 -0.775093022  0.377007684 -0.798640360 -1.018587043
## [171]  0.377007684 -0.045568949  0.042972570 -0.506277205 -0.250122285
## [176]  1.549806955 -0.786098730 -1.117870599 -0.480581224 -0.786098730
## [181]  1.776111383  0.696565845 -0.480581224 -0.762432124 -1.018587043
## [186]  0.240272639 -0.883601590 -0.371211247  0.397196034 -0.786098730
## [191] -0.480581224 -0.786098730 -0.340505568 -0.327168183 -1.117870599
## [196] -0.909375678 -0.909375678 -1.018587043 -0.120216783 -0.766467085
## [201] -0.506277205  2.665133284  0.012992394  0.662870654 -1.165960889
## [206] -0.291877205 -1.117870599  1.162970523 -1.295613270 -0.623512709
## [211] -0.883601590  2.825156350  0.006032634 -0.449152431  0.752001323
## [216]  1.810492018  0.518342473 -0.909375678 -0.786098730  1.428408233
## [221] -0.532198076 -0.114936635  5.851755458  0.397196034  0.519939169
## [226]  0.170608952  0.041170976  0.653331794 -0.786098730 -0.291877205
## [231]  4.946076411 -0.370358323  0.540729809  3.781631922 -0.051786770
## [236]  0.390265013 -0.291877205 -0.415063571 -0.909375678 -0.646580053
## [241] -1.117870599 -0.786098730 -0.762432124 -0.194718255  0.653331794
## [246] -0.619370113 -0.291877205 -0.118342671 -0.909375678  1.874751242
## [251] -0.981228225 -0.909375678 -0.634354664 -1.117870599  5.236678165
## [256] -1.209680475 -0.909375678  4.191974305 -1.018587043 -0.480581224
## [261] -0.762432124 -0.766467085 -0.909375678 -0.318843435 -0.619370113
## [266] -0.480581224 -1.018587043 -0.063660057  2.396164542 -0.480581224
## [271] -0.634354664  1.224819495  1.508495334  3.192813251  0.392460143
## [276] -0.194718255 -0.992386748
## attr(,"internals")
##                  Gi       E(Gi)        V(Gi)        Z(Gi) Pr(z > E(Gi))
##   [1,] 0.0031321185 0.003623188 7.820438e-06 -0.175601343  5.696964e-01
##   [2,] 0.0050113895 0.003623188 1.265277e-05  0.390265013  3.481703e-01
##   [3,] 0.0109953704 0.003623188 1.626375e-05  1.828039900  3.377178e-02
##   [4,] 0.0010401498 0.003623188 5.679422e-06 -1.083874118  8.607897e-01
##   [5,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##   [6,] 0.0029624753 0.003623188 9.154199e-06 -0.218374907  5.864315e-01
##   [7,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
##   [8,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##   [9,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
##  [10,] 0.0052304675 0.003623188 9.059519e-06  0.533996866  2.966719e-01
##  [11,] 0.0030581040 0.003623188 1.065172e-05 -0.173142404  5.687303e-01
##  [12,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [13,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [14,] 0.0057208238 0.003623188 1.603097e-05  0.523902035  3.001733e-01
##  [15,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [16,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [17,] 0.0061162080 0.003623188 1.065172e-05  0.763863545  2.224743e-01
##  [18,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
##  [19,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [20,] 0.0030372058 0.003623188 2.124358e-05 -0.127136697  5.505839e-01
##  [21,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [22,] 0.0007593014 0.003623188 2.124358e-05 -0.621358222  7.328180e-01
##  [23,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [24,] 0.0007593014 0.003623188 2.124358e-05 -0.621358222  7.328180e-01
##  [25,] 0.0056074766 0.003623188 1.297451e-05  0.550882783  2.908570e-01
##  [26,] 0.0254291163 0.003623188 5.770943e-06  9.077185239  5.570944e-20
##  [27,] 0.0040229885 0.003623188 1.615182e-05  0.099479164  4.603789e-01
##  [28,] 0.0131278539 0.003623188 1.595706e-05  2.379361417  8.671331e-03
##  [29,] 0.0205905206 0.003623188 1.077969e-05  5.167857453  1.183964e-07
##  [30,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
##  [31,] 0.0015186029 0.003623188 2.124358e-05 -0.456617713  6.760271e-01
##  [32,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [33,] 0.0034246575 0.003623188 1.055982e-05 -0.061094178  5.243579e-01
##  [34,] 0.0005733945 0.003623188 1.609593e-05 -0.760172919  7.764244e-01
##  [35,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
##  [36,] 0.0101149425 0.003623188 1.287395e-05  1.809281336  3.520365e-02
##  [37,] 0.0000000000 0.003623188 3.198209e-05 -0.640674601  7.391329e-01
##  [38,] 0.0077625571 0.003623188 1.271871e-05  1.160679998  1.228860e-01
##  [39,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [40,] 0.0006833713 0.003623188 6.209661e-06 -1.179740170  8.809482e-01
##  [41,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [42,] 0.0113895216 0.003623188 1.050507e-05  2.396164542  8.283826e-03
##  [43,] 0.0015186029 0.003623188 1.050507e-05 -0.649332580  7.419383e-01
##  [44,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
##  [45,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [46,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
##  [47,] 0.0033129460 0.003623188 7.022246e-06 -0.117074734  5.465996e-01
##  [48,] 0.0057339450 0.003623188 1.609593e-05  0.526114228  2.994044e-01
##  [49,] 0.0005694761 0.003623188 1.587432e-05 -0.766444193  7.782940e-01
##  [50,] 0.0030581040 0.003623188 1.065172e-05 -0.173142404  5.687303e-01
##  [51,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
##  [52,] 0.0015186029 0.003623188 2.124358e-05 -0.456617713  6.760271e-01
##  [53,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [54,] 0.0147106898 0.003623188 9.059519e-06  3.683673298  1.149484e-04
##  [55,] 0.0192307692 0.003623188 1.530840e-05  3.989061729  3.316758e-05
##  [56,] 0.0068649886 0.003623188 1.603097e-05  0.809666782  2.090659e-01
##  [57,] 0.0045871560 0.003623188 1.282941e-05  0.269128154  3.939155e-01
##  [58,] 0.0075930144 0.003623188 1.050507e-05  1.224819495  1.103216e-01
##  [59,] 0.0019524894 0.003623188 8.970993e-06 -0.557799295  7.115093e-01
##  [60,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
##  [61,] 0.0015220700 0.003623188 2.135430e-05 -0.454682106  6.753310e-01
##  [62,] 0.0059014814 0.003623188 3.188823e-06  1.275834530  1.010070e-01
##  [63,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [64,] 0.0020022883 0.003623188 7.897611e-06 -0.576777591  7.179551e-01
##  [65,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
##  [66,] 0.0012655024 0.003623188 6.925562e-06 -0.895897783  8.148463e-01
##  [67,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [68,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
##  [69,] 0.0009111617 0.003623188 6.209661e-06 -1.088328523  8.617749e-01
##  [70,] 0.0182232346 0.003623188 1.050507e-05  4.504585626  3.325129e-06
##  [71,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [72,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [73,] 0.0015186029 0.003623188 1.050507e-05 -0.649332580  7.419383e-01
##  [74,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
##  [75,] 0.0037965072 0.003623188 2.124358e-05  0.037603812  4.850018e-01
##  [76,] 0.0022779043 0.003623188 1.265277e-05 -0.378199745  6.473589e-01
##  [77,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
##  [78,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
##  [79,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
##  [80,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
##  [81,] 0.0026575550 0.003623188 1.050507e-05 -0.297929066  6.171213e-01
##  [82,] 0.0057339450 0.003623188 1.609593e-05  0.526114228  2.994044e-01
##  [83,] 0.0034403670 0.003623188 3.242857e-05 -0.032104298  5.128056e-01
##  [84,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [85,] 0.0000000000 0.003623188 1.609593e-05 -0.903093714  8.167619e-01
##  [86,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
##  [87,] 0.0119815668 0.003623188 1.291116e-05  2.326158626  1.000504e-02
##  [88,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
##  [89,] 0.0022779043 0.003623188 1.587432e-05 -0.337649739  6.321864e-01
##  [90,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
##  [91,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
##  [92,] 0.0079726651 0.003623188 1.050507e-05  1.341953999  8.980546e-02
##  [93,] 0.0015186029 0.003623188 6.925562e-06 -0.799722058  7.880641e-01
##  [94,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [95,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
##  [96,] 0.0051605505 0.003623188 7.929615e-06  0.545946531  2.925513e-01
##  [97,] 0.0029488860 0.003623188 9.096232e-06 -0.223575372  5.884561e-01
##  [98,] 0.0034246575 0.003623188 1.595706e-05 -0.049699455  5.198191e-01
##  [99,] 0.0045871560 0.003623188 1.609593e-05  0.240272639  4.050595e-01
## [100,] 0.0035203976 0.003623188 5.623924e-06 -0.043344556  5.172866e-01
## [101,] 0.0022779043 0.003623188 1.050507e-05 -0.415063571  6.609523e-01
## [102,] 0.0049374589 0.003623188 9.154199e-06  0.434384744  3.320045e-01
## [103,] 0.0026033192 0.003623188 8.970993e-06 -0.340505568  6.332621e-01
## [104,] 0.0057208238 0.003623188 1.603097e-05  0.523902035  3.001733e-01
## [105,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [106,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [107,] 0.0178147268 0.003623188 1.590652e-05  3.558294459  1.866354e-04
## [108,] 0.0050113895 0.003623188 1.265277e-05  0.390265013  3.481703e-01
## [109,] 0.0160919540 0.003623188 2.161494e-05  2.681923891  3.660006e-03
## [110,] 0.0045766590 0.003623188 9.059519e-06  0.316777802  3.757061e-01
## [111,] 0.0130165962 0.003623188 8.970993e-06  3.136194059  8.557795e-04
## [112,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [113,] 0.0000000000 0.003623188 1.277763e-05 -1.013598162  8.446127e-01
## [114,] 0.0028473804 0.003623188 1.587432e-05 -0.194718255  5.771932e-01
## [115,] 0.0051487414 0.003623188 1.603097e-05  0.381019662  3.515943e-01
## [116,] 0.0282227307 0.003623188 1.060873e-05  7.552579583  2.133604e-14
## [117,] 0.0045558087 0.003623188 1.587432e-05  0.234076200  4.074629e-01
## [118,] 0.0028473804 0.003623188 5.135810e-06 -0.342333802  6.339501e-01
## [119,] 0.0013016596 0.003623188 8.970993e-06 -0.775093022  7.808577e-01
## [120,] 0.0109589041 0.003623188 1.271871e-05  2.056936495  1.984617e-02
## [121,] 0.0009784736 0.003623188 9.017750e-06 -0.880703586  8.107609e-01
## [122,] 0.0049467275 0.003623188 1.055982e-05  0.407294521  3.418958e-01
## [123,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [124,] 0.0103668262 0.003623188 1.031195e-05  2.100022453  1.786343e-02
## [125,] 0.0153758542 0.003623188 1.587432e-05  2.949774409  1.590030e-03
## [126,] 0.0069124424 0.003623188 1.619851e-05  0.817259344  2.068901e-01
## [127,] 0.0032110092 0.003623188 1.282941e-05 -0.115075486  5.458074e-01
## [128,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [129,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [130,] 0.0175038052 0.003623188 2.135430e-05  3.003766040  1.333301e-03
## [131,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [132,] 0.0086560364 0.003623188 1.265277e-05  1.414884690  7.855115e-02
## [133,] 0.0006508298 0.003623188 8.970993e-06 -0.992386748  8.394955e-01
## [134,] 0.0030372058 0.003623188 1.050507e-05 -0.180794562  5.717356e-01
## [135,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [136,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
## [137,] 0.0185614849 0.003623188 3.280353e-05  2.608199220  4.550998e-03
## [138,] 0.0014958117 0.003623188 2.603754e-06 -1.318391544  9.063137e-01
## [139,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [140,] 0.0028032620 0.003623188 7.022246e-06 -0.309411796  6.214959e-01
## [141,] 0.0022779043 0.003623188 1.587432e-05 -0.337649739  6.321864e-01
## [142,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [143,] 0.0057112261 0.003623188 4.844093e-06  0.948707013  1.713848e-01
## [144,] 0.0000000000 0.003623188 1.603097e-05 -0.904921697  8.172466e-01
## [145,] 0.0177880897 0.003623188 2.178921e-05  3.034539714  1.204515e-03
## [146,] 0.0118764846 0.003623188 1.590652e-05  2.069378039  1.925531e-02
## [147,] 0.0051369863 0.003623188 1.595706e-05  0.378958346  3.523594e-01
## [148,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
## [149,] 0.0020248038 0.003623188 6.925562e-06 -0.607370609  7.281975e-01
## [150,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [151,] 0.0013698630 0.003623188 1.271871e-05 -0.631832997  7.362520e-01
## [152,] 0.0022779043 0.003623188 1.587432e-05 -0.337649739  6.321864e-01
## [153,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
## [154,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [155,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [156,] 0.0020248038 0.003623188 3.346058e-06 -0.873805045  8.088878e-01
## [157,] 0.0051605505 0.003623188 1.609593e-05  0.383193434  3.507882e-01
## [158,] 0.0092165899 0.003623188 1.619851e-05  1.389755737  8.230153e-02
## [159,] 0.0009762447 0.003623188 8.970993e-06 -0.883739885  8.115817e-01
## [160,] 0.0000000000 0.003623188 1.609593e-05 -0.903093714  8.167619e-01
## [161,] 0.0030372058 0.003623188 1.050507e-05 -0.180794562  5.717356e-01
## [162,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
## [163,] 0.0006508298 0.003623188 8.970993e-06 -0.992386748  8.394955e-01
## [164,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [165,] 0.0045558087 0.003623188 3.198209e-05  0.164911686  4.345067e-01
## [166,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
## [167,] 0.0013016596 0.003623188 8.970993e-06 -0.775093022  7.808577e-01
## [168,] 0.0051252847 0.003623188 1.587432e-05  0.377007684  3.530839e-01
## [169,] 0.0023728170 0.003623188 2.451182e-06 -0.798640360  7.877505e-01
## [170,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [171,] 0.0051252847 0.003623188 1.587432e-05  0.377007684  3.530839e-01
## [172,] 0.0034403670 0.003623188 1.609593e-05 -0.045568949  5.181731e-01
## [173,] 0.0038226300 0.003623188 2.154015e-05  0.042972570  4.828617e-01
## [174,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [175,] 0.0027334852 0.003623188 1.265277e-05 -0.250122285  5.987536e-01
## [176,] 0.0091743119 0.003623188 1.282941e-05  1.549806955  6.059393e-02
## [177,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [178,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [179,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [180,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [181,] 0.0086405530 0.003623188 7.980148e-06  1.776111383  3.785728e-02
## [182,] 0.0068337130 0.003623188 2.124358e-05  0.696565845  2.430373e-01
## [183,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [184,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [185,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [186,] 0.0045871560 0.003623188 1.609593e-05  0.240272639  4.050595e-01
## [187,] 0.0007593014 0.003623188 1.050507e-05 -0.883601590  8.115443e-01
## [188,] 0.0022935780 0.003623188 1.282941e-05 -0.371211247  6.447599e-01
## [189,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [190,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [191,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [192,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [193,] 0.0026033192 0.003623188 8.970993e-06 -0.340505568  6.332621e-01
## [194,] 0.0026333114 0.003623188 9.154199e-06 -0.327168183  6.282297e-01
## [195,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [196,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [197,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [198,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [199,] 0.0033361134 0.003623188 5.702437e-06 -0.120216783  5.478443e-01
## [200,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
## [201,] 0.0018223235 0.003623188 1.265277e-05 -0.506277205  6.936690e-01
## [202,] 0.0142694064 0.003623188 1.595706e-05  2.665133284  3.847891e-03
## [203,] 0.0036697248 0.003623188 1.282941e-05  0.012992394  4.948169e-01
## [204,] 0.0062642369 0.003623188 1.587432e-05  0.662870654  2.537067e-01
## [205,] 0.0009578544 0.003623188 5.225590e-06 -1.165960889  8.781849e-01
## [206,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [207,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [208,] 0.0077803204 0.003623188 1.277763e-05  1.162970523  1.224207e-01
## [209,] 0.0000000000 0.003623188 7.820438e-06 -1.295613270  9.024456e-01
## [210,] 0.0011389522 0.003623188 1.587432e-05 -0.623512709  7.335262e-01
## [211,] 0.0007593014 0.003623188 1.050507e-05 -0.883601590  8.115443e-01
## [212,] 0.0136986301 0.003623188 1.271871e-05  2.825156350  2.362878e-03
## [213,] 0.0036446469 0.003623188 1.265277e-05  0.006032634  4.975933e-01
## [214,] 0.0022779043 0.003623188 8.970993e-06 -0.449152431  6.733392e-01
## [215,] 0.0056194125 0.003623188 7.046629e-06  0.752001323  2.260251e-01
## [216,] 0.0094912680 0.003623188 1.050507e-05  1.810492018  3.510976e-02
## [217,] 0.0054669704 0.003623188 1.265277e-05  0.518342473  3.021097e-01
## [218,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [219,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [220,] 0.0087557604 0.003623188 1.291116e-05  1.428408233  7.658720e-02
## [221,] 0.0018982536 0.003623188 1.050507e-05 -0.532198076  7.027056e-01
## [222,] 0.0032765400 0.003623188 9.096232e-06 -0.114936635  5.457523e-01
## [223,] 0.0247086247 0.003623188 1.298353e-05  5.851755458  2.432057e-09
## [224,] 0.0050458716 0.003623188 1.282941e-05  0.397196034  3.456115e-01
## [225,] 0.0056947608 0.003623188 1.587432e-05  0.519939169  3.015530e-01
## [226,] 0.0041761579 0.003623188 1.050507e-05  0.170608952  4.322656e-01
## [227,] 0.0038138825 0.003623188 2.145321e-05  0.041170976  4.835798e-01
## [228,] 0.0059633028 0.003623188 1.282941e-05  0.653331794  2.567712e-01
## [229,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [230,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [231,] 0.0214452214 0.003623188 1.298353e-05  4.946076411  3.786212e-07
## [232,] 0.0025802752 0.003623188 7.929615e-06 -0.370358323  6.444422e-01
## [233,] 0.0050384513 0.003623188 6.850364e-06  0.540729809  2.943469e-01
## [234,] 0.0172494172 0.003623188 1.298353e-05  3.781631922  7.790181e-05
## [235,] 0.0034168565 0.003623188 1.587432e-05 -0.051786770  5.206507e-01
## [236,] 0.0050113895 0.003623188 1.265277e-05  0.390265013  3.481703e-01
## [237,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [238,] 0.0022779043 0.003623188 1.050507e-05 -0.415063571  6.609523e-01
## [239,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [240,] 0.0015220700 0.003623188 1.055982e-05 -0.646580053  7.410481e-01
## [241,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [242,] 0.0000000000 0.003623188 2.124358e-05 -0.786098730  7.840952e-01
## [243,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [244,] 0.0028473804 0.003623188 1.587432e-05 -0.194718255  5.771932e-01
## [245,] 0.0059633028 0.003623188 1.282941e-05  0.653331794  2.567712e-01
## [246,] 0.0007610350 0.003623188 2.135430e-05 -0.619370113  7.321637e-01
## [247,] 0.0022779043 0.003623188 2.124358e-05 -0.291877205  6.148097e-01
## [248,] 0.0030721966 0.003623188 2.167742e-05 -0.118342671  5.471019e-01
## [249,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [250,] 0.0123743233 0.003623188 2.178921e-05  1.874751242  3.041348e-02
## [251,] 0.0008600917 0.003623188 7.929615e-06 -0.981228225  8.367599e-01
## [252,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [253,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [254,] 0.0000000000 0.003623188 1.050507e-05 -1.117870599  8.681889e-01
## [255,] 0.0244874715 0.003623188 1.587432e-05  5.236678165  8.174621e-08
## [256,] 0.0000000000 0.003623188 8.970993e-06 -1.209680475  8.867992e-01
## [257,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [258,] 0.0160628845 0.003623188 8.806074e-06  4.191974305  1.382687e-05
## [259,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [260,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [261,] 0.0009111617 0.003623188 1.265277e-05 -0.762432124  7.770989e-01
## [262,] 0.0011389522 0.003623188 1.050507e-05 -0.766467085  7.783008e-01
## [263,] 0.0000000000 0.003623188 1.587432e-05 -0.909375678  8.184241e-01
## [264,] 0.0027841053 0.003623188 6.925562e-06 -0.318843435  6.250774e-01
## [265,] 0.0007610350 0.003623188 2.135430e-05 -0.619370113  7.321637e-01
## [266,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [267,] 0.0000000000 0.003623188 1.265277e-05 -1.018587043  8.458005e-01
## [268,] 0.0034168565 0.003623188 1.050507e-05 -0.063660057  5.253795e-01
## [269,] 0.0113895216 0.003623188 1.050507e-05  2.396164542  8.283826e-03
## [270,] 0.0017084282 0.003623188 1.587432e-05 -0.480581224  6.845929e-01
## [271,] 0.0013667426 0.003623188 1.265277e-05 -0.634354664  7.370753e-01
## [272,] 0.0075930144 0.003623188 1.050507e-05  1.224819495  1.103216e-01
## [273,] 0.0075930144 0.003623188 6.925562e-06  1.508495334  6.571390e-02
## [274,] 0.0185185185 0.003623188 2.176473e-05  3.192813251  7.044702e-04
## [275,] 0.0050228311 0.003623188 1.271871e-05  0.392460143  3.473591e-01
## [276,] 0.0028473804 0.003623188 1.587432e-05 -0.194718255  5.771932e-01
## [277,] 0.0006508298 0.003623188 8.970993e-06 -0.992386748  8.394955e-01
## attr(,"cluster")
##   [1] Low  Low  High High Low  High Low  Low  Low  High High Low  Low  High Low 
##  [16] Low  High Low  Low  Low  Low  Low  Low  Low  High High High Low  High Low 
##  [31] Low  Low  Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [46] Low  High High Low  High Low  Low  Low  High High High High Low  Low  Low 
##  [61] Low  High Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
##  [76] Low  Low  Low  Low  Low  Low  High High Low  High Low  High Low  Low  Low 
##  [91] Low  Low  Low  Low  Low  High High Low  High Low  Low  High Low  High Low 
## [106] Low  High Low  High High Low  Low  High Low  High High Low  Low  Low  Low 
## [121] Low  Low  Low  High Low  High High Low  Low  Low  Low  Low  Low  Low  Low 
## [136] Low  High High Low  High Low  Low  High High High High Low  Low  Low  Low 
## [151] Low  Low  Low  Low  Low  Low  High High Low  High Low  Low  Low  High Low 
## [166] Low  Low  Low  Low  Low  Low  High High Low  Low  High Low  Low  Low  Low 
## [181] High Low  Low  Low  Low  High Low  High High Low  Low  Low  Low  High Low 
## [196] Low  Low  Low  High Low  Low  Low  High Low  High Low  Low  High Low  Low 
## [211] Low  Low  Low  Low  High Low  Low  Low  Low  High Low  High High High Low 
## [226] Low  High High Low  Low  High High High High Low  Low  Low  Low  Low  Low 
## [241] Low  Low  Low  Low  High Low  Low  High Low  High High Low  Low  Low  Low 
## [256] Low  Low  High Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low  Low 
## [271] Low  Low  Low  High Low  Low  Low 
## Levels: Low High
## attr(,"gstari")
## [1] FALSE
## attr(,"call")
## localG(x = data_merge$KEK, listw = W.queen, alternative = "greater")
## attr(,"class")
## [1] "localG"

data_merge$lmZ2 <- Local_M
tmap_mode("plot")

## tmap mode set to plotting

tm_shape(data_merge) + tm_polygons(col = "lmZ2",
title = "Local Gi", style = "fixed",
breaks = c(-Inf, -1.96, 1.96, Inf),
labels = c("Coldspot", "No SAC", "Hotspot"),
palette =  c("blue", "white", "red")) +
tm_layout(legend.outside = TRUE)