Sample rainfall data (you can replace later with your real data)

data <- data.frame( Date = as.Date(c(“2025-01-01”,“2025-01-02”,“2025-03-10”,“2025-03-11”,“2025-06-20”,“2025-06-21”,“2025-09-15”,“2025-09-16”)), Rainfall = c(3.2, 0.0, 1.5, 0.2, 18.4, 12.7, 8.9, 10.3), Season = c(“Winter”,“Winter”,“Summer”,“Summer”,“Monsoon”,“Monsoon”,“Autumn”,“Autumn”) )

print(data)

plot(data\(Date, data\)Rainfall, type=“l”, xlab=“Date”, ylab=“Rainfall (mm)”, main=“Daily Rainfall Trend”)

Calculate average rainfall per season

avg_season <- aggregate(Rainfall ~ Season, data, mean) print(avg_season)

barplot(avg_season\(Rainfall, names.arg = avg_season\)Season, xlab = “Season”, ylab = “Average Rainfall (mm)”, main = “Average Rainfall by Season”)

boxplot(Rainfall ~ Season, data = data, xlab = “Season”, ylab = “Rainfall (mm)”, main = “Rainfall Distribution by Season”)

hist(data$Rainfall, xlab=“Rainfall (mm)”, main=“Rainfall Frequency Distribution”)