library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.5.3
data <- data.frame(
  provinsi = c("Aceh","Sumatera Utara","Sumatera Barat","Riau","Jambi",
               "Sumatera Selatan","Bengkulu","Lampung","Kep. Bangka Belitung",
               "Kepulauan Riau","DKI Jakarta","Jawa Barat","Jawa Tengah",
               "DI Yogyakarta","Jawa Timur","Banten","Bali","Nusa Tenggara Barat",
               "Nusa Tenggara Timur","Kalimantan Barat","Kalimantan Tengah",
               "Kalimantan Selatan","Kalimantan Timur","Kalimantan Utara",
               "Sulawesi Utara","Sulawesi Tengah","Sulawesi Selatan",
               "Sulawesi Tenggara","Gorontalo","Sulawesi Barat","Maluku",
               "Maluku Utara","Papua Barat","Papua Barat Daya","Papua",
               "Papua Selatan","Papua Tengah","Papua Pegunungan"),
  rls = c(10.42, 10.29, 10.03, 9.80, 9.32, 9.28, 9.71, 9.00, 9.09, 10.78, 11.58,
          9.46, 8.62, 10.58, 8.81, 9.82, 10.07, 8.87, 8.88, 8.54, 9.36, 9.17, 
          10.30, 9.69, 10.10, 9.38, 9.38, 9.90, 8.78, 8.91, 10.63, 9.78, 9.90,
          10.60, 10.65, 8.87, 6.27, 4.76 )
)

# Cek data
str(data)
## 'data.frame':    38 obs. of  2 variables:
##  $ provinsi: chr  "Aceh" "Sumatera Utara" "Sumatera Barat" "Riau" ...
##  $ rls     : num  10.42 10.29 10.03 9.8 9.32 ...
nrow(data)              
## [1] 38
summary(data$rls)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   4.760   8.932   9.575   9.457  10.092  11.580
# 4. Visualisasi 1: Histogram
ggplot(data, aes(x = rls)) +
  geom_histogram(bins = 8, fill = "steelblue", color = "white") +
  labs(title = "Distribusi Rata-rata Lama Sekolah 38 Provinsi",
       x = "Rata-rata Lama Sekolah (tahun)", y = "Frekuensi")

# Boxplot
ggplot(data, aes(y = rls)) +
  geom_boxplot(fill = "orange") +
  labs(title = "Boxplot Rata-rata Lama Sekolah 38 Provinsi",
       y = "Rata-rata Lama Sekolah (tahun)")

# Visualisasi khusus Kalimantan
kalimantan <- subset(data, provinsi %in% c("Kalimantan Barat",
                                           "Kalimantan Tengah","Kalimantan Selatan","Kalimantan Timur",
                                           "Kalimantan Utara"))

ggplot(kalimantan, aes(x = reorder(provinsi, rls), y = rls)) +
  geom_col(fill = "darkgreen") +
  geom_text(aes(label = rls), vjust = -0.3) +
  labs(title = "Rata-rata Lama Sekolah Provinsi di Pulau Kalimantan",
       x = "Provinsi", y = "Rata-rata Lama Sekolah (tahun)")