TM 2 Komputasi Statistika

Data=read.table(file.choose(),header = T) 
Data

##               Provinsi JKN_2025
## 1                 ACEH    98.18
## 2       SUMATERA_UTARA    68.69
## 3       SUMATERA_BARAT    78.19
## 4                 RIAU    73.30
## 5                JAMBI    61.04
## 6     SUMATERA_SELATAN    73.33
## 7             BENGKULU    76.20
## 8              LAMPUNG    71.90
## 9  KEP_BANGKA_BELITUNG    84.11
## 10            KEP_RIAU    80.00
## 11         DKI_JAKARTA    91.69
## 12          JAWA_BARAT    70.62
## 13         JAWA_TENGAH    72.96
## 14       DI_YOGYAKARTA    89.27
## 15          JAWA_TIMUR    73.86
## 16              BANTEN    74.49
## 17                BALI    89.15
## 18 NUSA_TENGGARA_BARAT    71.73
## 19 NUSA_TENGGARA_TIMUR    77.45
## 20    KALIMANTAN_BARAT    68.92
## 21   KALIMANTAN_TENGAH    67.89
## 22  KALIMANTAN_SELATAN    74.32
## 23    KALIMANTAN_TIMUR    83.00
## 24    KALIMANTAN_UTARA    88.40
## 25      SULAWESI_UTARA    83.46
## 26     SULAWESI_TENGAH    80.86
## 27    SULAWESI_SELATAN    86.56
## 28   SULAWESI_TENGGARA    84.19
## 29           GORONTALO    87.21
## 30      SULAWESI_BARAT    89.32
## 31              MALUKU    65.53
## 32        MALUKU_UTARA    73.41
## 33         PAPUA_BARAT    76.74
## 34    PAPUA_BARAT_DAYA    84.13
## 35               PAPUA    79.90
## 36       PAPUA_SELATAN    92.19
## 37        PAPUA_TENGAH    70.83
## 38    PAPUA_PEGUNUNGAN    96.31

library(ggplot2)

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

ggplot(Data, aes(x = JKN_2025)) +
  geom_histogram(
    binwidth = 5,
    boundary = 60,
    fill = "pink",
    color = "black"
  ) +
  scale_x_continuous(
    breaks = seq(60, 100, by = 5),
    limits = c(60, 100)
  ) +
  labs(
    title = "Persentase Penduduk yang Memiliki JKN di Indonesia Tahun 2025",
    x = "Persentase Penduduk Memiliki JKN (%)",
    y = "Freq"
  ) +
  theme_classic() +
  theme(
    axis.text.x = element_text(angle = 0, hjust = 0.5),
    plot.title = element_text(
      hjust = 0.5,
      face = "bold"
    )
  )

ggplot(Data, aes(x = "", y = JKN_2025)) +
  geom_boxplot(fill = "red") +
  labs(
    title = "Box Plot Persentase Penduduk yang Memiliki JKN Tahun 2025",
    x = "",
    y = "JKN (%)"
  ) +
  theme_minimal() +
  theme(
    plot.title = element_text(
      hjust = 0.5,
      face = "bold"
    )
  )

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)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats   1.0.1     ✔ stringr   1.5.2
## ✔ lubridate 1.9.4     ✔ tibble    3.3.0
## ✔ purrr     1.1.0     ✔ tidyr     1.3.1
## ✔ readr     2.1.5

## ── 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_starts(Provinsi, "KALIMANTAN"))

Data_Kalimantan

##             Provinsi JKN_2025
## 1   KALIMANTAN_BARAT    68.92
## 2  KALIMANTAN_TENGAH    67.89
## 3 KALIMANTAN_SELATAN    74.32
## 4   KALIMANTAN_TIMUR    83.00
## 5   KALIMANTAN_UTARA    88.40

ggplot(
  Data_Kalimantan,
  aes(x = Provinsi, y = JKN_2025, fill = Provinsi)
) +
  geom_col() +
  labs(
    title = "Persentase Penduduk yang Memiliki JKN di Pulau Kalimantan Tahun 2025",
    x = "Provinsi",
    y = "JKN (%)"
  ) +
  theme_minimal() +
  theme(
    axis.text.x = element_text(
      angle = 45,
      hjust = 1
    ),
    plot.title = element_text(
      hjust = 0.5,
      face = "bold"
    ),
    legend.position = "none"
  )