Data = read.table(file.choose(), header = TRUE, sep = "\t")
Data
##               Provinsi Formal
## 1                 ACEH  36.09
## 2       SUMATERA_UTARA  42.74
## 3       SUMATERA_BARAT  35.67
## 4                 RIAU  48.00
## 5                JAMBI  42.73
## 6     SUMATERA_SELATAN  36.87
## 7             BENGKULU  34.83
## 8              LAMPUNG  35.28
## 9  KEP_BANGKA_BELITUNG  48.71
## 10            KEP_RIAU  65.75
## 11         DKI_JAKARTA  63.37
## 12          JAWA_BARAT  45.05
## 13         JAWA_TENGAH  41.21
## 14       DI_YOGYAKARTA  46.03
## 15          JAWA_TIMUR  38.27
## 16              BANTEN  54.88
## 17                BALI  49.53
## 18 NUSA_TENGGARA_BARAT  31.58
## 19 NUSA_TENGGARA_TIMUR  31.56
## 20    KALIMANTAN_BARAT  43.61
## 21   KALIMANTAN_TENGAH  49.85
## 22  KALIMANTAN_SELATAN  46.52
## 23    KALIMANTAN_TIMUR  56.78
## 24    KALIMANTAN_UTARA  51.56
## 25      SULAWESI_UTARA  45.44
## 26     SULAWESI_TENGAH  35.43
## 27    SULAWESI_SELATAN  39.55
## 28   SULAWESI_TENGGARA  38.87
## 29           GORONTALO  38.86
## 30      SULAWESI_BARAT  29.37
## 31              MALUKU  33.94
## 32        MALUKU_UTARA  35.49
## 33         PAPUA_BARAT  36.81
## 34    PAPUA_BARAT_DAYA  44.34
## 35               PAPUA  47.12
## 36       PAPUA_SELATAN  34.12
## 37        PAPUA_TENGAH  15.28
## 38    PAPUA_PEGUNUNGAN   4.70

#VISUALISASI DATA UNTUK SELURUH PROVINSI DI INDONESIA

#Histogram

library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.5.3
ggplot(Data, aes(x = Formal)) +
  geom_histogram(
    binwidth = 10,
    boundary = 0,
    fill = "lightblue",
    color = "black"
  ) +
  scale_x_continuous(
    breaks = seq(0, 70, by = 10),
    limits = c(0, 70)
  ) +
  labs(
    title = "Persentase Tenaga Kerja Formal di Indonesia",
    x = "Tenaga Kerja Formal (%)",
    y = "Freq"
  ) +
  theme_classic() +
  theme(
    axis.text.x = element_text(angle = 0, hjust = 0.5),
    plot.title = element_text(
      hjust = 0.5,
      face = "bold"
    )
  )

#Boxplot

ggplot(Data, aes(x = "", y = Formal)) +
  geom_boxplot(fill = "tomato") +
  labs(
    title = "Box Plot Tenaga Kerja Formal Indonesia 2025",
    x = "",
    y = "Tenaga Kerja Formal (%)"
  ) +
  theme_minimal() +
  theme(
    plot.title = element_text(hjust = 0.5, face = "bold")
  )

#VISUALISASI DATA UNTUK PULAU KALIMANTAN

library(dplyr)
## Warning: package 'dplyr' was built under R version 4.5.3
## 
## 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
library(tidyverse)
## Warning: package 'tidyverse' was built under R version 4.5.3
## Warning: package 'tibble' was built under R version 4.5.3
## Warning: package 'tidyr' was built under R version 4.5.3
## Warning: package 'readr' was built under R version 4.5.3
## Warning: package 'purrr' was built under R version 4.5.3
## Warning: package 'stringr' was built under R version 4.5.2
## Warning: package 'forcats' was built under R version 4.5.3
## Warning: package 'lubridate' was built under R version 4.5.3
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats   1.0.1     ✔ stringr   1.6.0
## ✔ lubridate 1.9.5     ✔ tibble    3.3.1
## ✔ purrr     1.2.2     ✔ tidyr     1.3.2
## ✔ readr     2.2.0
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(ggplot2)

Data_Kalimantan <- Data %>%
  filter(str_detect(Provinsi, "KALIMANTAN"))

Data_Kalimantan
##             Provinsi Formal
## 1   KALIMANTAN_BARAT  43.61
## 2  KALIMANTAN_TENGAH  49.85
## 3 KALIMANTAN_SELATAN  46.52
## 4   KALIMANTAN_TIMUR  56.78
## 5   KALIMANTAN_UTARA  51.56
ggplot(Data_Kalimantan, aes(x = Provinsi, y = Formal, fill = Provinsi)) +
  geom_col() +
  labs(
    title = "Tenaga Kerja Formal Pulau Kalimantan Tahun 2025",
    x = "Provinsi",
    y = "Tenaga Kerja Formal (%)"
  ) +
  theme_minimal() +
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1),
    plot.title = element_text(hjust = 0.5, face = "bold"),
    legend.position = "none"
  )