TM 2 Prak Komstat

Data=read.table(file.choose(),header = T, sep = "\t")
Data

##                     Provinsi Rumah_Sakit
## 1                       Aceh          77
## 2             Sumatera Utara         185
## 3             Sumatera Barat          53
## 4                       Riau          67
## 5                      Jambi          42
## 6           Sumatera Selatan          72
## 7                   Bengkulu          25
## 8                    Lampung          64
## 9  Kepulauan Bangka Belitung          25
## 10            Kepulauan Riau          34
## 11               DKI Jakarta         142
## 12                Jawa Barat         382
## 13               Jawa Tengah         317
## 14             DI Yogyakarta          64
## 15                Jawa Timur         358
## 16                    Banten         111
## 17                      Bali          66
## 18       Nusa Tenggara Barat          39
## 19       Nusa Tenggara Timur          64
## 20          Kalimantan Barat          51
## 21         Kalimantan Tengah          31
## 22        Kalimantan Selatan          44
## 23          Kalimantan Timur          54
## 24          Kalimantan Utara          17
## 25            Sulawesi Utara          52
## 26           Sulawesi Tengah          37
## 27          Sulawesi Selatan         100
## 28         Sulawesi Tenggara          39
## 29                 Gorontalo          21
## 30            Sulawesi Barat          14
## 31                    Maluku          28
## 32              Maluku Utara          21
## 33               Papua Barat          14
## 34          Papua Barat Daya          12
## 35                     Papua          18
## 36             Papua Selatan           8
## 37              Papua Tengah          14
## 38          Papua Pegunungan           8

length(Data$Rumah_Sakit)

## [1] 38

mean(Data$Rumah_Sakit)

## [1] 72.89474

median(Data$Rumah_Sakit)

## [1] 43

range(Data$Rumah_Sakit)

## [1]   8 382

max(Data$Rumah_Sakit)-min(Data$Rumah_Sakit)

## [1] 374

max(Data$Rumah_Sakit)

## [1] 382

min(Data$Rumah_Sakit)

## [1] 8

var(Data$Rumah_Sakit)

## [1] 8304.367

sd(Data$Rumah_Sakit)

## [1] 91.1283

quantile(Data$Rumah_Sakit)

##     0%    25%    50%    75%   100% 
##   8.00  22.00  43.00  66.75 382.00

library(ggplot2)
ggplot(Data, aes(x = Rumah_Sakit)) +
  geom_histogram(
    binwidth = 50,
    boundary = 0,
    fill = "pink",
    color = "antiquewhite1"
  ) +
  scale_x_continuous(
    breaks = seq(0, 400, by = 50),
    limits = c(0, 400)
  ) +
  labs(
    title = "Jumlah Rumah Sakit di Indonesia Tahun 2025",
    x = "Jumlah Rumah Sakit",
    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 = Rumah_Sakit)) +
  geom_boxplot(fill = "red") +
  labs(title = "Box Plot Jumlah Rumah Sakit Indonesia 2025", x = "", y = "Jumlah Rumah Sakit") +
  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)

## 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_starts(Provinsi, "Kalimantan"))
Data_Kalimantan

##             Provinsi Rumah_Sakit
## 1   Kalimantan Barat          51
## 2  Kalimantan Tengah          31
## 3 Kalimantan Selatan          44
## 4   Kalimantan Timur          54
## 5   Kalimantan Utara          17

ggplot(Data_Kalimantan, aes(x = Provinsi, y = Rumah_Sakit, fill = Provinsi)) +
  geom_col() +
  labs(
    title = "Jumlah Rumah Sakit Pulau Kalimantan Tahun 2025",
    x = "Provinsi",
    y = "Jumlah Rumah Sakit"
  ) +
  theme_minimal() +
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1),
    plot.title = element_text(hjust = 0.5, face = "bold"),
    legend.position = "none"
  )