Analisis Luas Panen, Produktivitas, dan Produksi Padi Menurut Kabupaten/Kota di Provinsi Jawa Timur, tahun 2020, dengan menampilkan visualisasi boxplot, histogram, dan scatterplot

library(readxl)
Data_1 <- read_excel("C:/Belajar/Kalkulus/Data.1.xlsx", 
    sheet = "Lembar1", range = "A3:D45")

Kita dapat menggunakan head untuk memanggil 6 data pertama, tail untuk memanggil 6 data terakhir, summary untuk melihat nilai minimum, median, mean, max , quartil dll.

head(Data_1)
## # A tibble: 6 × 4
##   `Kabupaten/Kota\r\nRegency/Municipality` `Luas Lahan (ha)` Produktiv…¹ Produ…²
##   <chr>                                                <dbl>       <dbl>   <dbl>
## 1 <NA>                                                   NA         NA       NA 
## 2 Kabupaten/Regency                                      NA         NA       NA 
## 3 Pacitan                                             19360.        43.1  83378.
## 4 Ponorogo                                            67976.        58.4 396816.
## 5 Trenggalek                                          19834.        55.1 109241.
## 6 Tulungagung                                         36685.        60.9 223383.
## # … with abbreviated variable names
## #   ¹​`Produktivitas (ku/ha) Productivity (qu/ha)`, ²​`Produksi [ton]`
tail(Data_1)
## # A tibble: 6 × 4
##   `Kabupaten/Kota\r\nRegency/Municipality` `Luas Lahan (ha)` Produktiv…¹ Produ…²
##   <chr>                                                <dbl>       <dbl>   <dbl>
## 1 Pasuruan                                             1663.        60.4  1.00e4
## 2 Mojokerto                                             733.        63.4  4.65e3
## 3 Madiun                                               2279.        62.4  1.42e4
## 4 Surabaya                                             1884.        57.3  1.08e4
## 5 Batu                                                  771.        66.8  5.15e3
## 6 Jawa Timur                                        1761882.        56.9  1.00e7
## # … with abbreviated variable names
## #   ¹​`Produktivitas (ku/ha) Productivity (qu/ha)`, ²​`Produksi [ton]`
str(Data_1)
## tibble [42 × 4] (S3: tbl_df/tbl/data.frame)
##  $ Kabupaten/Kota
## Regency/Municipality    : chr [1:42] NA "Kabupaten/Regency" "Pacitan" "Ponorogo" ...
##  $ Luas Lahan (ha)                           : num [1:42] NA NA 19360 67976 19834 ...
##  $ Produktivitas (ku/ha) Productivity (qu/ha): num [1:42] NA NA 43.1 58.4 55.1 ...
##  $ Produksi [ton]                            : num [1:42] NA NA 83378 396816 109241 ...
summary(Data_1)
##  Kabupaten/Kota\r\nRegency/Municipality Luas Lahan (ha)    
##  Length:42                              Min.   :    733.3  
##  Class :character                       1st Qu.:  19597.1  
##  Mode  :character                       Median :  41212.6  
##                                         Mean   :  90352.9  
##                                         3rd Qu.:  67222.3  
##                                         Max.   :1761882.0  
##                                         NA's   :3          
##  Produktivitas (ku/ha) Productivity (qu/ha) Produksi [ton]    
##  Min.   :43.07                              Min.   :    4650  
##  1st Qu.:52.67                              1st Qu.:   92053  
##  Median :57.38                              Median :  225346  
##  Mean   :57.49                              Mean   :  513969  
##  3rd Qu.:61.78                              3rd Qu.:  403401  
##  Max.   :67.64                              Max.   :10022387  
##  NA's   :3                                  NA's   :3
Data_1
## # A tibble: 42 × 4
##    `Kabupaten/Kota\r\nRegency/Municipality` `Luas Lahan (ha)` Produkti…¹ Produ…²
##    <chr>                                                <dbl>      <dbl>   <dbl>
##  1 <NA>                                                   NA        NA       NA 
##  2 Kabupaten/Regency                                      NA        NA       NA 
##  3 Pacitan                                             19360.       43.1  83378.
##  4 Ponorogo                                            67976.       58.4 396816.
##  5 Trenggalek                                          19834.       55.1 109241.
##  6 Tulungagung                                         36685.       60.9 223383.
##  7 Blitar                                              35172.       57.4 201804.
##  8 Kediri                                              36470.       60.0 218625.
##  9 Malang                                              48475.       58.4 283097.
## 10 Lumajang                                            55725.       52.9 294865.
## # … with 32 more rows, and abbreviated variable names
## #   ¹​`Produktivitas (ku/ha) Productivity (qu/ha)`, ²​`Produksi [ton]`

Dari data yang disadur dari BPS, Survei Kerangka Sampel Area (KSA)/BPS-Statistics Indonesia, Area Sampling Frame (ASF) Survey dapat kita lihat Luas Panen, Produktivitas, dan Produksi Padi Menurut Kabupaten/Kota di Provinsi Jawa Timur, 2020

  1. Membuat Boxplot Produktivitas lahan
summary(Data_1$`Luas Lahan (ha)`)
##      Min.   1st Qu.    Median      Mean   3rd Qu.      Max.      NA's 
##     733.3   19597.1   41212.6   90352.9   67222.3 1761882.0         3
summary(Data_1$`Produksi [ton]`)
##     Min.  1st Qu.   Median     Mean  3rd Qu.     Max.     NA's 
##     4650    92053   225346   513969   403401 10022387        3
boxplot(Data_1$`Luas Lahan (ha)`)

boxplot(Data_1$`Produksi [ton]`)

Menampilkan produktivitas lahan jawa timur 2020 dengan boxplot

  1. Membuat Histogram
hist(Data_1$`Produktivitas (ku/ha) Productivity (qu/ha)`)

hist(Data_1$`Luas Lahan (ha)`)

hist(Data_1$`Produksi [ton]`)

Dapat dilihat dari histogram diatas bahwa rata rata kabupaten/kota memiliki lahan produktif sekitar 60 ku/ha, dan rata rata luas lahan < 500000 ha , dan untuk produksi 10022386.98 ton

  1. Membuat Scatterplot

    plot(Data_1$`Luas Lahan (ha)`,Data_1$`Produktivitas (ku/ha) Productivity (qu/ha)`,main = "Scatterplot antara luas lahan dan produktivitas pada kota dan kabupaten Provinsi Jawa Timur 2020")

    Disana dapat terlihat hubungan antara luas lahan dan produktivitas pada kabupaten dan kota Provinsi Jawa Timur tahun 2022

Sumber :https://jatim.bps.go.id/statictable/2021/09/06/2237/luas-panen-produktivitas-dan-produksi-padi-menurut-kabupaten-kota-di-provinsi-jawa-timur-2019-dan-2020.html