TM Komstat 2

library(readxl)
library(ggplot2)

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

library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

data_air <- read_excel(file.choose())

## New names:
## • `` -> `...2`

colnames(data_air) <- c("Provinsi", "Persentase")

data_air <- data_air %>%
  filter(!is.na(Provinsi)) %>%
  filter(!is.na(Persentase)) %>%
  filter(Provinsi != "INDONESIA")

data_air$Persentase <- as.numeric(data_air$Persentase)

ggplot(data_air, aes(x = reorder(Provinsi, Persentase), y = Persentase)) +
  geom_bar(stat = "identity", fill = "deeppink") +
  coord_flip() +
  labs(
    title = "Persentase Rumah Tangga Memiliki Akses Air Minum Layak 2025",
    x = "Provinsi",
    y = "Persentase (%)"
  ) +
  theme_minimal()

ggplot(data_air, aes(y = Persentase)) +
  geom_boxplot(fill = "blue4", color = "darkorchid4", outlier.color = "deeppink3") +
  labs(
    title = "Boxplot Distribusi Akses Air Minum Layak 38 Provinsi",
    y = "Persentase (%)"
  ) +
  theme_minimal()

data_kalimantan <- data_air %>%
  filter(grepl("KALIMANTAN", Provinsi))

ggplot(data_kalimantan, aes(x = reorder(Provinsi, -Persentase), y = Persentase)) +
  geom_bar(stat = "identity", fill = "green") +
  geom_text(aes(label = paste0(Persentase, "%")), vjust = -0.5, fontface = "bold") +
  labs(
    title = "Persentase Akses Air Minum Layak di Pulau Kalimantan 2025",
    x = "Provinsi",
    y = "Persentase (%)"
  ) +
  ylim(0, 110) +
  theme_minimal()