angka-perceraian-dan-pernikahan-di-jateng-2023.R

library(readxl)

## Warning: package 'readxl' was built under R version 4.5.3

# Membaca data Excel
data <- read_excel(file.choose(), skip = 4)

# Memberi nama kolom
names(data) <- c("Kabupaten_Kota", "Pernikahan", "Perceraian")

# Melihat data
View(data)

# Mengeluarkan Provinsi Jawa Tengah
data_35 <- data[data$Kabupaten_Kota != "3300 PROVINSI JAWA TENGAH", ]

# Mengubah variabel menjadi numerik
pernikahan <- as.numeric(data_35$Pernikahan)
perceraian <- as.numeric(data_35$Perceraian)

# Mengecek jumlah data
length(pernikahan)

## [1] 35

length(perceraian)

## [1] 35

# Q1, Q2, Q3
Q1_pernikahan <- quantile(pernikahan, 0.25)
Q2_pernikahan <- quantile(pernikahan, 0.50)
Q3_pernikahan <- quantile(pernikahan, 0.75)

# IQR
IQR_pernikahan <- Q3_pernikahan - Q1_pernikahan

# Batas bawah dan batas atas
BB_pernikahan <- Q1_pernikahan - 1.5 * IQR_pernikahan
BA_pernikahan <- Q3_pernikahan + 1.5 * IQR_pernikahan

# Menampilkan hasil
Q1_pernikahan

##  25% 
## 5645

Q2_pernikahan

##  50% 
## 6771

Q3_pernikahan

##    75% 
## 9122.5

IQR_pernikahan

##    75% 
## 3477.5

BB_pernikahan

##    25% 
## 428.75

BA_pernikahan

##      75% 
## 14338.75

pencilan_pernikahan <- pernikahan[
  pernikahan < BB_pernikahan |
    pernikahan > BA_pernikahan
]

pencilan_pernikahan

## [1] 15156 15505

boxplot(pernikahan,
        main = "Boxplot Jumlah Pernikahan",
        ylab = "Jumlah Pernikahan")

# Q1, Q2, Q3
Q1_perceraian <- quantile(perceraian, 0.25)
Q2_perceraian <- quantile(perceraian, 0.50)
Q3_perceraian <- quantile(perceraian, 0.75)

# IQR
IQR_perceraian <- Q3_perceraian - Q1_perceraian

# Batas bawah dan batas atas
BB_perceraian <- Q1_perceraian - 1.5 * IQR_perceraian
BA_perceraian <- Q3_perceraian + 1.5 * IQR_perceraian

# Menampilkan hasil
Q1_perceraian

##    25% 
## 1356.5

Q2_perceraian

##  50% 
## 2060

Q3_perceraian

##    75% 
## 2619.5

IQR_perceraian

##  75% 
## 1263

BB_perceraian

##  25% 
## -538

BA_perceraian

##  75% 
## 4514

pencilan_perceraian <- perceraian[
  perceraian < BB_perceraian |
    perceraian > BA_perceraian
]

pencilan_perceraian

## [1] 5922 4579 5041

boxplot(perceraian,
        main = "Boxplot Jumlah Perceraian",
        ylab = "Jumlah Perceraian")