TUGAS-MASUK-2.R

#Nama : Fairuz Putri Fahriyyah
#NIM : 2507016022
#Kelas : 1 A
#TUGAS MASUK 2 PRAKTIKUM KOMPUTASI STATISTIKA

Data <- read.delim(file.choose(), header = TRUE, sep = "\t", na.strings = "NA")
Data

##                Provinsi   TGM
## 1                  Aceh 59.23
## 2        Sumatera Utara 59.36
## 3        Sumatera Barat 59.42
## 4                  Riau 58.59
## 5                 Jambi 57.77
## 6      Sumatera Selatan 60.86
## 7              Bengkulu 57.74
## 8               Lampung 58.79
## 9  Kep. Bangka Belitung 58.00
## 10            Kep. Riau 59.33
## 11          DKI Jakarta 57.23
## 12           Jawa Barat 59.19
## 13          Jawa Tengah 57.11
## 14        DI Yogyakarta 55.44
## 15           Jawa Timur 58.86
## 16               Banten 56.53
## 17                 Bali    NA
## 18  Nusa Tenggara Barat 61.19
## 19  Nusa Tenggara Timur 62.05
## 20     Kalimantan Barat 59.85
## 21    Kalimantan Tengah 56.24
## 22   Kalimantan Selatan 57.10
## 23     Kalimantan Timur 57.86
## 24     Kalimantan Utara 58.89
## 25       Sulawesi Utara 58.75
## 26      Sulawesi Tengah 59.51
## 27     Sulawesi Selatan 59.84
## 28    Sulawesi Tenggara 58.64
## 29            Gorontalo 59.20
## 30       Sulawesi Barat 48.63
## 31               Maluku 57.38
## 32         Maluku Utara 60.66
## 33          Papua Barat 50.77
## 34     Papua Barat Daya    NA
## 35                Papua 51.12
## 36        Papua Selatan    NA
## 37         Papua Tengah    NA
## 38     Papua Pegunungan    NA

length(Data$TGM)

## [1] 38

min(Data$TGM, na.rm = TRUE)

## [1] 48.63

max(Data$TGM, na.rm = TRUE)

## [1] 62.05

range(Data$TGM, na.rm = TRUE)

## [1] 48.63 62.05

max(Data$TGM, na.rm = TRUE) - min(Data$TGM, na.rm = TRUE)

## [1] 13.42

mean(Data$TGM, na.rm = TRUE)

## [1] 57.91303

median(Data$TGM, na.rm = TRUE)

## [1] 58.75

var(Data$TGM, na.rm = TRUE)

## [1] 8.315253

sd(Data$TGM, na.rm = TRUE)

## [1] 2.883618

quantile(Data$TGM, na.rm = TRUE)

##    0%   25%   50%   75%  100% 
## 48.63 57.23 58.75 59.36 62.05

#Histogram
library(ggplot2)

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

ggplot(Data, aes(x = TGM)) +
  geom_histogram(
    fill = "skyblue",
    color = "pink",
    na.rm = TRUE
  ) +
  labs(
    title = "Histogram Tingkat Kegemaran Membaca",
    x = "Tingkat Kegemaran Membaca",
    y = "Frekuensi"
  )

## `stat_bin()` using `bins = 30`. Pick better value `binwidth`.

# Boxplot
ggplot(Data, aes(y = TGM)) +
  geom_boxplot(
    fill = "violet",
    color = "pink",
    na.rm = TRUE
  ) +
  labs(
    title = "Boxplot Tingkat Kegemaran Membaca",
    y = "Tingkat Kegemaran Membaca"
  )

Data_Kalimantan <- Data[grep("^Kalimantan", Data$Provinsi), ]
Data_Kalimantan

##              Provinsi   TGM
## 20   Kalimantan Barat 59.85
## 21  Kalimantan Tengah 56.24
## 22 Kalimantan Selatan 57.10
## 23   Kalimantan Timur 57.86
## 24   Kalimantan Utara 58.89

# TINGKAT KEGEMARAN MEMBACA DI KALIMANTAN
ggplot(Data_Kalimantan,
       aes(x = Provinsi, y = TGM, fill = Provinsi)) +
  geom_col(
    color = "pink",
    na.rm = TRUE
  ) +
  labs(
    title = "Tingkat Kegemaran Membaca Provinsi di Kalimantan",
    x = "Provinsi",
    y = "Tingkat Kegemaran Membaca"
  )