TM KOMSTAT 2

data <- read.table(file.choose(), header = TRUE, sep = "\t")

data <- data[order(data$IPM), ]

par(mar = c(5, 12, 4, 2))

barplot(
  data$IPM,
  names.arg = data$Provinsi,
  horiz = TRUE,
  las = 1,
  cex.names = 0.65,
  space = 0.8,
  col = "steelblue",
  border = NA,
  main = "Indeks Pembangunan Manusia Menurut Provinsi Tahun 2025",
  xlab = "IPM",
  xlim = c(0, 90)
)

par(mar = c(5, 4, 4, 2))

hist(
  data$IPM,
  main = "Distribusi Indeks Pembangunan Manusia 38 Provinsi Tahun 2025",
  xlab = "IPM",
  ylab = "Frekuensi",
  breaks = 8,
  col = "lightblue",
  border = "white"
)

kalimantan <- data[data$Provinsi %in% c(
  "Kalimantan Barat",
  "Kalimantan Tengah",
  "Kalimantan Selatan",
  "Kalimantan Timur",
  "Kalimantan Utara"
), ]

kalimantan <- kalimantan[order(kalimantan$IPM), ]

barplot(
  kalimantan$IPM,
  names.arg = kalimantan$Provinsi,
  las = 2,
  cex.names = 0.8,
  space = 0.7,
  col = "steelblue",
  border = NA,
  main = "IPM Provinsi di Pulau Kalimantan Tahun 2025",
  xlab = "Provinsi",
  ylab = "IPM",
  ylim = c(0, 90)
)