Proyek_ADV

#memasukkan data
library(readxl)
library(tidyverse)

## Warning: package 'tidyverse' was built under R version 4.3.3

## Warning: package 'tidyr' was built under R version 4.3.3

## Warning: package 'purrr' was built under R version 4.3.3

## Warning: package 'dplyr' was built under R version 4.3.3

## Warning: package 'forcats' was built under R version 4.3.3

## Warning: package 'lubridate' was built under R version 4.3.3

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ ggplot2   3.4.4     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── 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

ternak = read_xlsx("C:/Users/ENGGAR/Downloads/Jumlah Populasi Ternak Per Kecamatan Di Kabupaten Garut -  Tahun 2023/Jumlah Populasi Ternak Per Kecamatan Di Kabupaten Garut -  Tahun 2023.xlsx")
head(ternak)

## # A tibble: 6 × 16
##      id kode_provinsi nama_provinsi kode_kabupaten nama_kabupaten kode_kecamatan
##   <dbl>         <dbl> <chr>                  <dbl> <chr>                   <dbl>
## 1     1            32 JAWA BARAT              3205 GARUT                 3205130
## 2     2            32 JAWA BARAT              3205 GARUT                 3205230
## 3     3            32 JAWA BARAT              3205 GARUT                 3205150
## 4     4            32 JAWA BARAT              3205 GARUT                 3205290
## 5     5            32 JAWA BARAT              3205 GARUT                 3205030
## 6     6            32 JAWA BARAT              3205 GARUT                 3205011
## # ℹ 10 more variables: kecamatan <chr>, sapi_potong <dbl>, sapi_perah <dbl>,
## #   kerbau <dbl>, domba <dbl>, kambing <dbl>, ayam_ras_pedaging <dbl>,
## #   ayam_ras_petelur <dbl>, ayam_buras <dbl>, itik <dbl>

#memilih data yang akan digunakan
ternak = ternak %>% rename(ayam = ayam_buras, sapi = sapi_potong) %>% select(sapi, kerbau, domba, kambing, ayam, itik) %>% gather(key = "ternak", value = "jumlah")
head(ternak)

## # A tibble: 6 × 2
##   ternak jumlah
##   <chr>   <dbl>
## 1 sapi     79.4
## 2 sapi    120. 
## 3 sapi     74.5
## 4 sapi    582. 
## 5 sapi   1099. 
## 6 sapi    685.

#anova
ternak_anova = aov(jumlah ~  ternak, data = ternak)
summary(ternak_anova)

##              Df    Sum Sq   Mean Sq F value Pr(>F)    
## ternak        5 5.283e+10 1.057e+10   438.3 <2e-16 ***
## Residuals   246 5.929e+09 2.410e+07                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

#Pie Chart
ggplot(ternak, aes(x = "", y = jumlah, fill = ternak)) +
  geom_bar(stat = "identity", width = 1) +
  coord_polar("y", start = 0) +
  labs(title = "Persentase Jenis Ternak") +
  theme_void()

#Bar Chart
ggplot(ternak, aes(x = ternak, y = jumlah)) +
  geom_bar(stat = "identity", fill = "skyblue") +
  labs(title = "Boxplot of Jumlah Ternak", x = "Ternak", y = "Jumlah") +
  theme_minimal()