#Uji Independensi Chi Square 
Wonosobo<-matrix(c(83, 99, 48, 67), byrow=T, nrow=2)
dimnames(Wonosobo)<-
  list(Kecamatan=c("Kepil", "Kalibawang"),
       JumlahPelangganTeleponPerTahun =c("2019", "2020"))
addmargins(Wonosobo)
##             JumlahPelangganTeleponPerTahun
## Kecamatan    2019 2020 Sum
##   Kepil        83   99 182
##   Kalibawang   48   67 115
##   Sum         131  166 297
chisq.test(Wonosobo)
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  Wonosobo
## X-squared = 0.28468, df = 1, p-value = 0.5937
#Koefisien Kontingensi
N<-sum(Wonosobo)
xsquared<-0.28468
KoefisienKontingensi<-sqrt((xsquared)/(N+xsquared))
KoefisienKontingensi
## [1] 0.03094512