library(readr)
library(Hmisc)
## 
## Attaching package: 'Hmisc'
## The following objects are masked from 'package:base':
## 
##     format.pval, units
#INPUT DATA KEDALAM R
setwd("D:/2024/Datathon/Pengolahan Data/IKL 2023/")
df <- read_csv("Data_Merge_2023.csv")
## Rows: 5888 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (12): GID_4, GID_0, COUNTRY, GID_1, NAME_1, GID_2, NAME_2, GID_3, NAME_3...
## dbl (12): CC_4, lat, long, AOD_2023mean, NO2_2023mean, CO_2023mean, SO2_2023...
## lgl  (1): VARNAME_4
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
dataku <- df[,17:25]


#MELIHAT SUMMARY DARI DATA
summary(dataku)
##   AOD_2023mean     NO2_2023mean      CO_2023mean       SO2_2023mean   
##  Min.   :0.1670   Min.   :0.03869   Min.   :0.01658   Min.   :0.1642  
##  1st Qu.:0.3140   1st Qu.:0.12358   1st Qu.:0.40389   1st Qu.:0.4498  
##  Median :0.3576   Median :0.19454   Median :0.53194   Median :0.5052  
##  Mean   :0.3792   Mean   :0.22630   Mean   :0.55086   Mean   :0.5061  
##  3rd Qu.:0.4356   3rd Qu.:0.26171   3rd Qu.:0.71332   3rd Qu.:0.5592  
##  Max.   :0.6929   Max.   :0.97025   Max.   :0.99310   Max.   :0.9525  
##  NDVI_2023mean     NTL_2023mean        LST_2023mean     TSS_2023mean   
##  Min.   :0.1281   Min.   :0.0004518   Min.   :0.2288   Min.   :0.2546  
##  1st Qu.:0.2258   1st Qu.:0.0037260   1st Qu.:0.5593   1st Qu.:0.7001  
##  Median :0.3084   Median :0.0078841   Median :0.6471   Median :0.7162  
##  Mean   :0.2981   Mean   :0.0176532   Mean   :0.6486   Mean   :0.7186  
##  3rd Qu.:0.3667   3rd Qu.:0.0188033   3rd Qu.:0.7323   3rd Qu.:0.7381  
##  Max.   :0.5855   Max.   :0.1894942   Max.   :0.9874   Max.   :0.8501  
##   IKL_2023mean     
##  Min.   :-0.45750  
##  1st Qu.:-0.09348  
##  Median : 0.04664  
##  Mean   : 0.06077  
##  3rd Qu.: 0.19809  
##  Max.   : 0.66623
hist.data.frame(dataku)