Climate_Assessment.R

#setting the working directory
setwd('C:/Users/PC/Desktop/Pula')
#getting the working directory
getwd()

## [1] "C:/Users/PC/Desktop/Pula"

dir()

##  [1] "~$daily_rainfall.xlsx"       "~$Data.xlsx"                
##  [3] "Bar graph.png"               "Climate_Assessment.html"    
##  [5] "Climate_Assessment.R"        "Climate_Assessment.spin.R"  
##  [7] "Climate_Assessment.spin.Rmd" "daily_rainfall.csv"         
##  [9] "daily_rainfall.xlsx"         "Data.xlsx"                  
## [11] "ee-chart.csv"                "ee-chart.png"               
## [13] "Farmer's_interview.csv"      "Farmer's_interview.xlsx"    
## [15] "Hospital_Ward.cpg"           "Hospital_Ward.dbf"          
## [17] "Hospital_Ward.prj"           "Hospital_Ward.qpj"          
## [19] "Hospital_Ward.shp"           "Hospital_Ward.shx"          
## [21] "Kimimini_subcounty.cpg"      "Kimimini_subcounty.dbf"     
## [23] "Kimimini_subcounty.prj"      "Kimimini_subcounty.qpj"     
## [25] "Kimimini_subcounty.shp"      "Kimimini_subcounty.shx"     
## [27] "kimimini_wards.cpg"          "kimimini_wards.dbf"         
## [29] "kimimini_wards.prj"          "kimimini_wards.qpj"         
## [31] "kimimini_wards.shp"          "kimimini_wards.shx"         
## [33] "Map.png"                     "Monthly_rainfall.csv"       
## [35] "Monthly_rainfall.xlsx"       "Project.qgz"                
## [37] "rainfall.png"                "Rainfall_data.csv"          
## [39] "Rplot.png"                   "Sub_county.cpg"             
## [41] "Sub_county.dbf"              "Sub_county.prj"             
## [43] "Sub_county.qpj"              "Sub_county.shp"             
## [45] "Sub_county.shx"              "Towns.cpg"                  
## [47] "Towns.dbf"                   "Towns.prj"                  
## [49] "Towns.qpj"                   "Towns.shp"                  
## [51] "Towns.shx"                   "Trans_Nzoia County.cpg"     
## [53] "Trans_Nzoia County.dbf"      "Trans_Nzoia County.prj"     
## [55] "Trans_Nzoia County.qpj"      "Trans_Nzoia County.shp"     
## [57] "Trans_Nzoia County.shx"

#calling the libraries
library(sf)

## Linking to GEOS 3.9.1, GDAL 3.2.1, PROJ 7.2.1

library(sp)
library(ggplot2)
#importing the rainfall data
Rainfall<-st_read("Monthly_rainfall.csv")

## Reading layer `Monthly_rainfall' from data source 
##   `C:\Users\PC\Desktop\Pula\Monthly_rainfall.csv' using driver `CSV'

## Warning: no simple feature geometries present: returning a data.frame or tbl_df

#changing the average column from character to numeric
Rainfall<-transform(Rainfall, Average=as.numeric(Average), 
                    Precipitation=as.numeric(Precipitation))

#Performing descriptive analysis(calculating the mean,median,min,
#max, standard deviation and quartile)

#maximum value
max(Rainfall$Precipitation)

## [1] 293.902

#minimum value
min(Rainfall$Precipitation)

## [1] 91.349

#mean value
mean(Rainfall$Precipitation)

## [1] 160.8124

#standard deviation
sd(Rainfall$Precipitation)

## [1] 78.57637

#using the summary function
summary(Rainfall)

##     Month           Precipitation       Average     
##  Length:5           Min.   : 91.35   Min.   :3.044  
##  Class :character   1st Qu.:119.63   1st Qu.:3.988  
##  Mode  :character   Median :139.32   Median :4.644  
##                     Mean   :160.81   Mean   :5.297  
##                     3rd Qu.:159.86   3rd Qu.:5.329  
##                     Max.   :293.90   Max.   :9.480