TM2.Komstat_StellaBunganAnjani

data <- read.delim(file.choose(),
                   header=TRUE,
                   sep="\t",
                   dec=",")
names(data) <- c("Provinsi","Merokok")
data <- data[data$Provinsi != "INDONESIA",]
data

##                Provinsi Merokok
## 1                  ACEH   28.19
## 2        SUMATERA UTARA   26.41
## 3        SUMATERA BARAT   30.48
## 4                  RIAU   28.88
## 5                 JAMBI   30.29
## 6      SUMATERA SELATAN   31.02
## 7              BENGKULU   32.03
## 8               LAMPUNG   34.25
## 9  KEP. BANGKA BELITUNG   27.71
## 10            KEP. RIAU   24.67
## 11          DKI JAKARTA   25.17
## 12           JAWA BARAT   32.07
## 13          JAWA TENGAH   28.48
## 14        DI YOGYAKARTA   23.29
## 15           JAWA TIMUR   28.11
## 16               BANTEN   31.19
## 17                 BALI   19.02
## 18  NUSA TENGGARA BARAT   32.32
## 19  NUSA TENGGARA TIMUR   28.55
## 20     KALIMANTAN BARAT   27.80
## 21    KALIMANTAN TENGAH   27.13
## 22   KALIMANTAN SELATAN   23.21
## 23     KALIMANTAN TIMUR   22.42
## 24     KALIMANTAN UTARA   25.83
## 25       SULAWESI UTARA   25.81
## 26      SULAWESI TENGAH   28.48
## 27     SULAWESI SELATAN   24.05
## 28    SULAWESI TENGGARA   25.42
## 29            GORONTALO   30.72
## 30       SULAWESI BARAT   26.31
## 31               MALUKU   27.59
## 32         MALUKU UTARA   29.08
## 33          PAPUA BARAT   27.76
## 34     PAPUA BARAT DAYA   24.94
## 35                PAPUA   22.02
## 36        PAPUA SELATAN   31.62
## 37         PAPUA TENGAH   23.98
## 38     PAPUA PEGUNUNGAN   26.52
## 39                           NA
## 40                           NA

data2 <- data[order(data$Merokok),]
par(mar=c(4, 6, 2, 8))
barplot(data2$Merokok,
        names.arg=data2$Provinsi,
        horiz=TRUE,
        las=1,
        cex.names=0.6,
        main="Persentase Penduduk yang Merokok Tahun 2025",
        xlab="Persentase (%)",
        ylab="Provinsi")

boxplot(data$Merokok,
        main="Boxplot Persentase Penduduk yang Merokok Tahun 2025",
        ylab="Persentase (%)")

Kalimantan <- subset(data,
                     Provinsi %in% c("KALIMANTAN BARAT",
                                     "KALIMANTAN TENGAH",
                                     "KALIMANTAN SELATAN",
                                     "KALIMANTAN TIMUR",
                                     "KALIMANTAN UTARA"))
Kalimantan

##              Provinsi Merokok
## 20   KALIMANTAN BARAT   27.80
## 21  KALIMANTAN TENGAH   27.13
## 22 KALIMANTAN SELATAN   23.21
## 23   KALIMANTAN TIMUR   22.42
## 24   KALIMANTAN UTARA   25.83

barplot(Kalimantan$Merokok,
        names.arg=Kalimantan$Provinsi,
        las=2,
        main="Persentase Penduduk yang Merokok di Pulau Kalimantan",
        xlab="Provinsi",
        ylab="Persentase (%)")