Creating Bar plot with Standard Error Bars in RStudio

I currently have a .csv containing water quality data taken from 22 rainfall events. For each event, I collected pH, TDS, EC, TSS, DOC, NH4, NO3, NO2, PO4, Zn, Fe, and Cl measurements at 12 sites (identified as "C1Up", "C1in1", "C1In2", "C1Out", "C1Down", "C2In", "C2Out", "C3In", and "C3Out". The format of my date file is:

Values showing NA indicate that that parameter was not measured at that site on that day.

EventNo,Date,Year,Site,pH ,TDS,Ec,Turbidity (NTU),TSS (g/L),DOC (mg/L),Ammonia,Nitrate,Nitrite,Phosphate,Zinc ,Iron,Chloride
1,05-Aug,2022,C1Up,7.84,224,447,1.64,NA,NA,0.78,0.172,0.004,0.14,0.67,0.17,16.0
1,05-Aug,2022,C1In1,6.24,39,78,7.25,NA,NA,4.32,0.161,0.030,0.23,0.82,0.14,14.0
1,05-Aug,2022,C1In2,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,
1,05-Aug,2022,C1Out,7.51,165,330,5.63,NA,NA,NA,0.54,0.155,0.013,0.11,0.33,0.09,3.4
1,05-Aug,2022,C1Down,8.25,230,460,2.06,NA,NA,0.78,0.434,0.029,0.13,0.08,0.12,7.4
1,05-Aug,2022,C2In,8.52,45,90,11.1,NA,NA,0.24,0.131,0.007,0.15,0.09,0.04,2.7
1,05-Aug,2022, C2Out,7.85,42,85,3.95,NA,NA,0.36,0.178,0.009,0.16,0.10,0.17,22.0

I want to use this data to create a bar plot with means and standard error bars, showing the water quality data among sites. I have provided a sample plot below:

Does anyone have any advice as to how to code this in R?

TIA

Hi, @brant , function geom_errorbar() and geom_text() in ggplot2 maybe helpful. You coluld learn it.