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
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
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
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