# R code to display the barplot of number of OCV agaist the count as per county 
library(readxl)
ACTIVE <- read_excel("D:/PERIS/ACTIVE.xlsx")
ACTIVE

## # A tibble: 64,162 × 57
##    cbo_id cbo          ward_id ward  consituency_id constituency countyid county
##     <dbl> <chr>          <dbl> <chr>          <dbl> <chr>           <dbl> <chr> 
##  1   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
##  2   7306 CHEER UP PR…     911 Kiuu             162 Ruiru              22 Kiambu
##  3   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
##  4   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
##  5   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
##  6   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
##  7   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
##  8   7306 CHEER UP PR…     911 Kiuu             162 Ruiru              22 Kiambu
##  9   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
## 10   7306 CHEER UP PR…     912 Mwiki            162 Ruiru              22 Kiambu
## # ℹ 64,152 more rows
## # ℹ 49 more variables: cpims_ovc_id <dbl>, ovc_names <chr>, gender <chr>,
## #   dob <chr>, date_of_birth <chr>, age <dbl>, age_at_reg <dbl>,
## #   agerange <chr>, birthcert <chr>, bcertnumber <chr>, ovcdisability <chr>,
## #   ncpwdnumber <chr>, ovchivstatus <chr>, artstatus <chr>, facility_id <dbl>,
## #   facility <chr>, facility_mfl_code <dbl>, date_of_linkage <chr>,
## #   ccc_number <chr>, duration_on_art <chr>, viral_load <dbl>, …

number_of_children <- table(ACTIVE$county)
barplot(number_of_children, xlab = "county", ylab = "count", col = "blue", main = "BARPLOT OF NUMBER OF OVC AGAINST THE COUNT")

barplot

# with or without birth certificate
OVCwith_without_berthcertficate <- table(ACTIVE$birthcert)
barplot(OVCwith_without_berthcertficate, xlab = "Birth certificate status", ylab = "count", main = "BARPLOT OF OVC WITH AND WITHOUT BIRTH CERTIFICATE", col = "red")

barplot

# with or without HIV
with_without_HIV_ovc <- table(ACTIVE$ovchivstatus)
barplot(with_without_HIV_ovc, xlab = "HIV status", ylab = "count", col = "green", main = "BARPLOT TO SHOW OVC WITH/WITHOUT HIV")

barplot

# barplot to show the age range of OVC
age_range <- table(ACTIVE$agerange)
barplot(age_range, xlab = "Age range", ylab = "count", col = "skyblue", main = "BARPLOT OF OVC AGE RANGE AGAINST COUNT")

barplot