Analisis_Gempa_di_Pulau_jawa
220605110044_Ridho Aulia Rahman_Teknik Informatika_UINMALIKI
2022-10-15
Analisis Gempa yang terjadi pulau jawa dari kurun abad ke-19-sampai tahun 2022
library(readxl)
earthquake <- read_excel("C:/Belajar/Kalkulus/Data_3.xlsx")
earthquake## # A tibble: 16 × 43
## ID Earthq…¹ Flag …² Year Month Day Focal…³ EQ Pr…⁴ Mw Ma…⁵ Ms Ma…⁶ Mb Ma…⁷
## <chr> <chr> <dbl> <chr> <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 7181 Tsunami 2006 07 17 34 7.7 7.7 7.2 6.1
## 2 8386 Tsunami 2009 09 02 46 7 7 7 6.8
## 3 6071 Tsunami 1823 09 09 150 6.8 NA 6.8 NA
## 4 9909 <NA> 2012 06 04 50 5.8 5.8 5.2 5.9
## 5 6149 Tsunami 1863 03 16 NA NA NA NA NA
## 6 3254 <NA> 1924 11 12 NA NA NA NA NA
## 7 10038 <NA> 2013 07 08 60 5.7 5.7 NA 5.5
## 8 3064 <NA> 1916 09 11 100 7.3 NA 7.3 NA
## 9 6832 Tsunami 1994 06 03 26 6.6 6.6 6.4 6.1
## 10 10108 Tsunami 1994 06 04 11 6.5 6.5 6.3 6
## 11 5385 Tsunami 1994 06 02 18 7.8 7.8 7.2 5.7
## 12 10287 <NA> 2017 12 15 90 6.5 6.5 NA NA
## 13 4683 <NA> 1975 02 09 27 5.6 NA 5.6 5.2
## 14 4851 <NA> 1979 11 02 62 6.1 NA 6.1 NA
## 15 6833 Tsunami 1994 06 05 17 6.1 6.1 6.3 5.7
## 16 7533 <NA> 2007 08 08 289 7.5 7.5 NA NA
## # … with 33 more variables: `Ml Magnitude` <lgl>, `MFA Magnitude` <lgl>,
## # `Unknown Magnitude` <lgl>, Intensity <dbl>, Country <chr>, State <lgl>,
## # `Location name` <chr>, `Region code` <dbl>, `Earthquake : Deaths` <dbl>,
## # `Earthquake : Deaths Description` <chr>, `Earthquake : Missing` <lgl>,
## # `Earthquake : Missing Description` <lgl>, `Earthquake : Injuries` <dbl>,
## # `Earthquake : Injuries Description` <chr>,
## # `Earthquake : Damage (in M$)` <dbl>, …Untuk melihat struktur data kita digunakan fungsi str
str(earthquake)## tibble [16 × 43] (S3: tbl_df/tbl/data.frame)
## $ ID Earthquake : chr [1:16] "7181" "8386" "6071" "9909" ...
## $ Flag Tsunami : chr [1:16] "Tsunami" "Tsunami" "Tsunami" NA ...
## $ Year : num [1:16] 2006 2009 1823 2012 1863 ...
## $ Month : chr [1:16] "07" "09" "09" "06" ...
## $ Day : chr [1:16] "17" "02" "09" "04" ...
## $ Focal Depth : num [1:16] 34 46 150 50 NA NA 60 100 26 11 ...
## $ EQ Primary : num [1:16] 7.7 7 6.8 5.8 NA NA 5.7 7.3 6.6 6.5 ...
## $ Mw Magnitude : num [1:16] 7.7 7 NA 5.8 NA NA 5.7 NA 6.6 6.5 ...
## $ Ms Magnitude : num [1:16] 7.2 7 6.8 5.2 NA NA NA 7.3 6.4 6.3 ...
## $ Mb Magnitude : num [1:16] 6.1 6.8 NA 5.9 NA NA 5.5 NA 6.1 6 ...
## $ Ml Magnitude : logi [1:16] NA NA NA NA NA NA ...
## $ MFA Magnitude : logi [1:16] NA NA NA NA NA NA ...
## $ Unknown Magnitude : logi [1:16] NA NA NA NA NA NA ...
## $ Intensity : num [1:16] NA NA NA NA NA NA NA NA NA NA ...
## $ Country : chr [1:16] "INDONESIA" "INDONESIA" "INDONESIA" "INDONESIA" ...
## $ State : logi [1:16] NA NA NA NA NA NA ...
## $ Location name : chr [1:16] "INDONESIA: JAVA" "INDONESIA: JAVA" "INDONESIA: JAVA" "INDONESIA: JAVA" ...
## $ Region code : num [1:16] 60 60 60 60 60 60 60 60 60 60 ...
## $ Earthquake : Deaths : num [1:16] NA 81 NA NA NA 60 NA NA NA NA ...
## $ Earthquake : Deaths Description : chr [1:16] NA "Some (~51 to 100 deaths)" NA NA ...
## $ Earthquake : Missing : logi [1:16] NA NA NA NA NA NA ...
## $ Earthquake : Missing Description : logi [1:16] NA NA NA NA NA NA ...
## $ Earthquake : Injuries : num [1:16] NA 1297 NA 2 NA ...
## $ Earthquake : Injuries Description : chr [1:16] NA "Very Many (~1001 or more deaths)" NA "Few (~1 to 50 deaths)" ...
## $ Earthquake : Damage (in M$) : num [1:16] NA NA NA NA NA NA NA NA NA NA ...
## $ Earthquake : Damage Description : chr [1:16] NA "SEVERE (~>$5 to $24 million)" NA "LIMITED (roughly corresponding to less than $1 million)" ...
## $ Earthquakes : Houses destroyed : num [1:16] NA NA NA NA NA NA 124 340 NA NA ...
## $ Earthquakes : Houses destroyed Description : chr [1:16] NA NA NA NA ...
## $ Earthquakes : Houses damaged : chr [1:16] NA "54000" NA "23" ...
## $ Earthquakes : Houses damaged Description : chr [1:16] NA "Very Many (~1001 or more houses)" NA "Few (~1 to 50 houses)" ...
## $ Total Effects : Deaths : num [1:16] 802 81 NA NA NA NA NA NA NA NA ...
## $ Total Effects : Deaths Description : chr [1:16] "Many (~101 to 1000 deaths)" "Some (~51 to 100 deaths)" NA NA ...
## $ Total Effects : Missing : logi [1:16] NA NA NA NA NA NA ...
## $ Total Effects : Missing Description : logi [1:16] NA NA NA NA NA NA ...
## $ Total Effects : Injuries : num [1:16] 498 1297 NA 2 NA ...
## $ Total Effects : Injuries Description : chr [1:16] "Many (~101 to 1000 deaths)" "Very Many (~1001 or more deaths)" NA "Few (~1 to 50 deaths)" ...
## $ Total Effects : Damages in million Dollars : chr [1:16] "55.000" NA NA NA ...
## $ Total Effects : Damage Description : chr [1:16] "EXTREME (~$25 million or more)" "SEVERE (~>$5 to $24 million)" NA "LIMITED (roughly corresponding to less than $1 million)" ...
## $ Total Effects : Houses Destroyed : num [1:16] 1623 NA NA NA NA ...
## $ Total Effects : Houses Destroyed Description: chr [1:16] "Very Many (~1001 or more houses)" NA NA NA ...
## $ Total Effects : Houses Damaged : num [1:16] NA 54000 NA 23 NA NA NA NA NA NA ...
## $ Total Effects : Houses Damaged Description : chr [1:16] NA "Very Many (~1001 or more houses)" NA "Few (~1 to 50 houses)" ...
## $ Coordinates : chr [1:16] "-9.254, 107.411" "-7.782, 107.297" "-6.5, 108.5" "-7.692, 106.371" ...Melihat 6 data teratas yang kita punya digunakan fungsi head
head(earthquake)## # A tibble: 6 × 43
## ID Earthqu…¹ Flag …² Year Month Day Focal…³ EQ Pr…⁴ Mw Ma…⁵ Ms Ma…⁶ Mb Ma…⁷
## <chr> <chr> <dbl> <chr> <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 7181 Tsunami 2006 07 17 34 7.7 7.7 7.2 6.1
## 2 8386 Tsunami 2009 09 02 46 7 7 7 6.8
## 3 6071 Tsunami 1823 09 09 150 6.8 NA 6.8 NA
## 4 9909 <NA> 2012 06 04 50 5.8 5.8 5.2 5.9
## 5 6149 Tsunami 1863 03 16 NA NA NA NA NA
## 6 3254 <NA> 1924 11 12 NA NA NA NA NA
## # … with 33 more variables: `Ml Magnitude` <lgl>, `MFA Magnitude` <lgl>,
## # `Unknown Magnitude` <lgl>, Intensity <dbl>, Country <chr>, State <lgl>,
## # `Location name` <chr>, `Region code` <dbl>, `Earthquake : Deaths` <dbl>,
## # `Earthquake : Deaths Description` <chr>, `Earthquake : Missing` <lgl>,
## # `Earthquake : Missing Description` <lgl>, `Earthquake : Injuries` <dbl>,
## # `Earthquake : Injuries Description` <chr>,
## # `Earthquake : Damage (in M$)` <dbl>, …Melihat rangkuman statistik data gempa di pulau jawa pada kurun 1800-sekarang
summary(earthquake)## ID Earthquake Flag Tsunami Year Month
## Length:16 Length:16 Min. :1823 Length:16
## Class :character Class :character 1st Qu.:1962 Class :character
## Mode :character Mode :character Median :1994 Mode :character
## Mean :1970
## 3rd Qu.:2008
## Max. :2017
##
## Day Focal Depth EQ Primary Mw Magnitude
## Length:16 Min. : 11.00 Min. :5.600 Min. :5.700
## Class :character 1st Qu.: 26.25 1st Qu.:6.100 1st Qu.:6.200
## Mode :character Median : 48.00 Median :6.550 Median :6.550
## Mean : 70.00 Mean :6.643 Mean :6.720
## 3rd Qu.: 83.00 3rd Qu.:7.225 3rd Qu.:7.375
## Max. :289.00 Max. :7.800 Max. :7.800
## NA's :2 NA's :2 NA's :6
## Ms Magnitude Mb Magnitude Ml Magnitude MFA Magnitude
## Min. :5.200 Min. :5.200 Mode:logical Mode:logical
## 1st Qu.:6.200 1st Qu.:5.700 NA's:16 NA's:16
## Median :6.400 Median :5.900
## Mean :6.491 Mean :5.889
## 3rd Qu.:7.100 3rd Qu.:6.100
## Max. :7.300 Max. :6.800
## NA's :5 NA's :7
## Unknown Magnitude Intensity Country State
## Mode:logical Min. :7.00 Length:16 Mode:logical
## NA's:16 1st Qu.:7.25 Class :character NA's:16
## Median :7.50 Mode :character
## Mean :7.50
## 3rd Qu.:7.75
## Max. :8.00
## NA's :14
## Location name Region code Earthquake : Deaths
## Length:16 Min. :60 Min. : 1.0
## Class :character 1st Qu.:60 1st Qu.: 4.0
## Mode :character Median :60 Median :30.0
## Mean :60 Mean :35.2
## 3rd Qu.:60 3rd Qu.:60.0
## Max. :60 Max. :81.0
## NA's :11
## Earthquake : Deaths Description Earthquake : Missing
## Length:16 Mode:logical
## Class :character NA's:16
## Mode :character
##
##
##
##
## Earthquake : Missing Description Earthquake : Injuries
## Mode:logical Min. : 1.0
## NA's:16 1st Qu.: 2.0
## Median : 36.0
## Mean : 307.2
## 3rd Qu.: 200.0
## Max. :1297.0
## NA's :11
## Earthquake : Injuries Description Earthquake : Damage (in M$)
## Length:16 Min. :16
## Class :character 1st Qu.:16
## Mode :character Median :16
## Mean :16
## 3rd Qu.:16
## Max. :16
## NA's :15
## Earthquake : Damage Description Earthquakes : Houses destroyed
## Length:16 Min. :124
## Class :character 1st Qu.:178
## Mode :character Median :232
## Mean :232
## 3rd Qu.:286
## Max. :340
## NA's :14
## Earthquakes : Houses destroyed Description Earthquakes : Houses damaged
## Length:16 Length:16
## Class :character Class :character
## Mode :character Mode :character
##
##
##
##
## Earthquakes : Houses damaged Description Total Effects : Deaths
## Length:16 Min. : 1.0
## Class :character 1st Qu.: 10.5
## Mode :character Median : 55.5
## Mean :192.7
## 3rd Qu.:198.8
## Max. :802.0
## NA's :10
## Total Effects : Deaths Description Total Effects : Missing
## Length:16 Mode:logical
## Class :character NA's:16
## Mode :character
##
##
##
##
## Total Effects : Missing Description Total Effects : Injuries
## Mode:logical Min. : 1.0
## NA's:16 1st Qu.: 19.0
## Median : 200.0
## Mean : 351.0
## 3rd Qu.: 460.5
## Max. :1297.0
## NA's :9
## Total Effects : Injuries Description
## Length:16
## Class :character
## Mode :character
##
##
##
##
## Total Effects : Damages in million Dollars Total Effects : Damage Description
## Length:16 Length:16
## Class :character Class :character
## Mode :character Mode :character
##
##
##
##
## Total Effects : Houses Destroyed Total Effects : Houses Destroyed Description
## Min. : 124.0 Length:16
## 1st Qu.: 286.0 Class :character
## Median : 920.0 Mode :character
## Mean : 896.8
## 3rd Qu.:1530.8
## Max. :1623.0
## NA's :12
## Total Effects : Houses Damaged Total Effects : Houses Damaged Description
## Min. : 23 Length:16
## 1st Qu.: 2424 Class :character
## Median : 4825 Mode :character
## Mean :19616
## 3rd Qu.:29413
## Max. :54000
## NA's :13
## Coordinates
## Length:16
## Class :character
## Mode :character
##
##
##
## Dari data diatas kita dapat mengetahui dalam kurun waktu tersebut terjadi 16 gempa di pulau jawa dengan magnitudo tertinggi yaitu 7,8 SR, dan kematian terbanyak 81 jiwa.
Agar lebih mudah dipahami data akan ditampilkan dalam bentuk sliceplot
plot(earthquake$Year,xlab = "Frekuensi Gempa di pulau jawa",ylim = c(1800,2050))
Plot diatas menunjukkan data frekuensi gempa di pulau jawa dalam kurun waktu tersebut
plot(earthquake$Year,earthquake$`Mw Magnitude`)
Data yang terekam hanya >1950 disebabkan tidak ditemukannya catatan mendetail tentang gempa yang terjadi pada waktu itu