obese_proportions_2020.R

#### Author: M Samuel
#### Date: 24th March 2022
#### This script calculates proportions of people who had a BMI >= 30 and (univariate) differences between groups using chi squared.


## Specify libraries
library(pacman)
library(tidyverse)
library(Hmisc)
library(here)
library(arrow)
library(purrr)
library(broom)
library(data.table)
library(forcats)
library(rstatix)
library(janitor)

## Read in files  >>> Change PATH!!

# check working directory:  getwd()

BMI_complete_categories <- read_feather (here::here ("output/data", "BMI_complete_median.feather"))


BMI_complete_categories <- BMI_complete_categories %>% 
  dplyr::ungroup() %>%
  dplyr::filter (year == "2020") %>%
  dplyr::mutate (imd = as.factor(imd)) %>%
  dplyr::mutate (imd = fct_relevel(imd, "1", "2", "3", "4", "5")) %>%
  dplyr::mutate(age_group = as.factor(age_group)) %>%
  dplyr::mutate(age_group = fct_relevel(age_group, "18-39", "40-65", "65-80", "80+"))


BMI_complete_categories$obese[BMI_complete_categories$obese == 1] <- "TRUE"
BMI_complete_categories$obese[BMI_complete_categories$obese == 0] <- "FALSE"




BMI_complete_categories_2020_DT <- data.table(BMI_complete_categories)



# develop a vector of explanatory variables
explanatory_vars <- c("sex", 
                      "age_group", 
                      "region", 
                      "imd", 
                      "ethnic_no_miss", 
                      "eth_group_16",
                      "comorbid_learning_disability",     
                      "comorbid_depression",            
                      "comorbid_dementia",               
                      "comorbid_psychosis_schiz_bipolar",
                      "comorbid_diabetes_type",
                      "comorbid_diabetes_t1",           
                      "comorbid_diabetes_t2",             
                      "comorbid_asthma",                
                      "comorbid_COPD",                   
                      "comorbid_stroke_and_TIA",         
                      "comorbid_chronic_cardiac", 
                      "comorbid_hypertension",           
                      "comorbid_all_cancer")



N_obese <- BMI_complete_categories_2020_DT[, .N, ]
obese_table<- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N),] 

obese_table <- dplyr::mutate (obese_table, N = N_obese) %>%
  dplyr::mutate(proportion=n_obese/N) %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(group = "all", .before=1) %>%
  dplyr::mutate(variable = "all", .before=1)




##  Develop count table and chi squared test of difference between population long hand for each variable to reduce memory


#########  AGE GROUPS #######

#1.  count by age_group
N_age_group <- BMI_complete_categories_2020_DT[, .N, by="age_group"]
obese_age_group <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="age_group"] 
obese_age_group <-   obese_age_group[order(age_group)]

obese_age_group <- dplyr::left_join(obese_age_group, N_age_group)
obese_age_group <- obese_age_group %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_age_group <- obese_age_group %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "age_group", .before=1) %>%
  dplyr::rename(group = age_group)

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
age_group_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(age_group, obese) %>%
  select(-1) %>% 
  chisq_test() 
age_group_chisq <- dplyr::mutate (age_group_chisq, variable = "age_group") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_age_group <- obese_age_group %>%
  dplyr::left_join(age_group_chisq, by = "variable")





###### SEX

#1.  count by sex
N_sex <- BMI_complete_categories_2020_DT[, .N, by="sex"]
obese_sex <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="sex"] 
obese_sex <-   obese_sex[order(sex)]

obese_sex <- dplyr::left_join(obese_sex, N_sex)
obese_sex <- obese_sex %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_sex <- obese_sex %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "sex", .before=1) %>%
  dplyr::rename(group = sex)

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
sex_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(sex, obese) %>%
  select(-1) %>% 
  chisq_test() 
sex_chisq <- dplyr::mutate (sex_chisq, variable = "sex") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_sex <- obese_sex %>%
  dplyr::left_join(sex_chisq, by = "variable")


#########  region #######

#1.  count by region
N_region <- BMI_complete_categories_2020_DT[, .N, by="region"]
obese_region <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="region"] 
obese_region <-   obese_region[order(region)]

obese_region <- dplyr::left_join(obese_region, N_region)
obese_region <- obese_region %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_region <- obese_region %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "region", .before=1) %>%
  dplyr::rename(group = region)

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
region_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(region, obese) %>%
  select(-1) %>% 
  chisq_test() 
region_chisq <- dplyr::mutate (region_chisq, variable = "region") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_region <- obese_region %>%
  dplyr::left_join(region_chisq, by = "variable")

#########  imd #######

#1.  count by imd
N_imd <- BMI_complete_categories_2020_DT[, .N, by="imd"]
obese_imd <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="imd"] 
obese_imd <-   obese_imd[order(imd)]

obese_imd <- dplyr::left_join(obese_imd, N_imd)
obese_imd <- obese_imd %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_imd <- obese_imd %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "imd", .before=1) %>%
  dplyr::rename(group = imd)

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
imd_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(imd, obese) %>%
  select(-1) %>% 
  chisq_test() 
imd_chisq <- dplyr::mutate (imd_chisq, variable = "imd") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_imd <- obese_imd %>%
  dplyr::left_join(imd_chisq, by = "variable")

#########  ethnic_no_miss #######

#1.  count by ethnic_no_miss
N_ethnic_no_miss <- BMI_complete_categories_2020_DT[, .N, by="ethnic_no_miss"]
obese_ethnic_no_miss <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="ethnic_no_miss"] 
obese_ethnic_no_miss <-   obese_ethnic_no_miss[order(ethnic_no_miss)]

obese_ethnic_no_miss <- dplyr::left_join(obese_ethnic_no_miss, N_ethnic_no_miss)
obese_ethnic_no_miss <- obese_ethnic_no_miss %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_ethnic_no_miss <- obese_ethnic_no_miss %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "ethnic_no_miss", .before=1) %>%
  dplyr::rename(group = ethnic_no_miss)

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
ethnic_no_miss_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(ethnic_no_miss, obese) %>%
  select(-1) %>% 
  chisq_test() 
ethnic_no_miss_chisq <- dplyr::mutate (ethnic_no_miss_chisq, variable = "ethnic_no_miss") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_ethnic_no_miss <- obese_ethnic_no_miss %>%
  dplyr::left_join(ethnic_no_miss_chisq, by = "variable")

#########  eth_group_16 #######

#1.  count by eth_group_16
N_eth_group_16 <- BMI_complete_categories_2020_DT[, .N, by="eth_group_16"]
obese_eth_group_16 <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="eth_group_16"] 
obese_eth_group_16 <-   obese_eth_group_16[order(eth_group_16)]

obese_eth_group_16 <- dplyr::left_join(obese_eth_group_16, N_eth_group_16)
obese_eth_group_16 <- obese_eth_group_16 %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_eth_group_16 <- obese_eth_group_16 %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "eth_group_16", .before=1) %>%
  dplyr::rename(group = eth_group_16)

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
eth_group_16_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(eth_group_16, obese) %>%
  select(-1) %>% 
  chisq_test() 
eth_group_16_chisq <- dplyr::mutate (eth_group_16_chisq, variable = "eth_group_16") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_eth_group_16 <- obese_eth_group_16 %>%
  dplyr::left_join(eth_group_16_chisq, by = "variable")

########  comorbid_learning_disability #######

#1.  count by comorbid_learning_disability
N_comorbid_learning_disability <- BMI_complete_categories_2020_DT[, .N, by="comorbid_learning_disability"]
obese_comorbid_learning_disability <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_learning_disability"] 
obese_comorbid_learning_disability <-   obese_comorbid_learning_disability[order(comorbid_learning_disability)]

obese_comorbid_learning_disability <- dplyr::left_join(obese_comorbid_learning_disability, N_comorbid_learning_disability)
obese_comorbid_learning_disability <- obese_comorbid_learning_disability %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_learning_disability <- obese_comorbid_learning_disability %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_learning_disability", .before=1) %>%
  dplyr::rename(group = comorbid_learning_disability) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_learning_disability_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_learning_disability, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_learning_disability_chisq <- dplyr::mutate (comorbid_learning_disability_chisq, variable = "comorbid_learning_disability") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_learning_disability <- obese_comorbid_learning_disability %>%
  dplyr::left_join(comorbid_learning_disability_chisq, by = "variable")



#########  comorbid_depression #######

#1.  count by comorbid_depression
N_comorbid_depression <- BMI_complete_categories_2020_DT[, .N, by="comorbid_depression"]
obese_comorbid_depression <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_depression"] 
obese_comorbid_depression <-   obese_comorbid_depression[order(comorbid_depression)]

obese_comorbid_depression <- dplyr::left_join(obese_comorbid_depression, N_comorbid_depression)
obese_comorbid_depression <- obese_comorbid_depression %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_depression <- obese_comorbid_depression %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_depression", .before=1) %>%
  dplyr::rename(group = comorbid_depression) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_depression_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_depression, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_depression_chisq <- dplyr::mutate (comorbid_depression_chisq, variable = "comorbid_depression") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_depression <- obese_comorbid_depression %>%
  dplyr::left_join(comorbid_depression_chisq, by = "variable")


#########  comorbid_dementia #######

#1.  count by comorbid_dementia
N_comorbid_dementia <- BMI_complete_categories_2020_DT[, .N, by="comorbid_dementia"]
obese_comorbid_dementia <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_dementia"] 
obese_comorbid_dementia <-   obese_comorbid_dementia[order(comorbid_dementia)]

obese_comorbid_dementia <- dplyr::left_join(obese_comorbid_dementia, N_comorbid_dementia)
obese_comorbid_dementia <- obese_comorbid_dementia %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_dementia <- obese_comorbid_dementia %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_dementia", .before=1) %>%
  dplyr::rename(group = comorbid_dementia) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_dementia_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_dementia, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_dementia_chisq <- dplyr::mutate (comorbid_dementia_chisq, variable = "comorbid_dementia") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_dementia <- obese_comorbid_dementia %>%
  dplyr::left_join(comorbid_dementia_chisq, by = "variable")


#########  comorbid_psychosis_schiz_bipolar #######

#1.  count by comorbid_psychosis_schiz_bipolar
N_comorbid_psychosis_schiz_bipolar <- BMI_complete_categories_2020_DT[, .N, by="comorbid_psychosis_schiz_bipolar"]
obese_comorbid_psychosis_schiz_bipolar <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_psychosis_schiz_bipolar"] 
obese_comorbid_psychosis_schiz_bipolar <-   obese_comorbid_psychosis_schiz_bipolar[order(comorbid_psychosis_schiz_bipolar)]

obese_comorbid_psychosis_schiz_bipolar <- dplyr::left_join(obese_comorbid_psychosis_schiz_bipolar, N_comorbid_psychosis_schiz_bipolar)
obese_comorbid_psychosis_schiz_bipolar <- obese_comorbid_psychosis_schiz_bipolar %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_psychosis_schiz_bipolar <- obese_comorbid_psychosis_schiz_bipolar %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_psychosis_schiz_bipolar", .before=1) %>%
  dplyr::rename(group = comorbid_psychosis_schiz_bipolar) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_psychosis_schiz_bipolar_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_psychosis_schiz_bipolar, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_psychosis_schiz_bipolar_chisq <- dplyr::mutate (comorbid_psychosis_schiz_bipolar_chisq, variable = "comorbid_psychosis_schiz_bipolar") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_psychosis_schiz_bipolar <- obese_comorbid_psychosis_schiz_bipolar %>%
  dplyr::left_join(comorbid_psychosis_schiz_bipolar_chisq, by = "variable")

########################

#########  comorbid_diabetes_t1 #######

#1.  count by comorbid_diabetes_t1
N_comorbid_diabetes_t1 <- BMI_complete_categories_2020_DT[, .N, by="comorbid_diabetes_t1"]
obese_comorbid_diabetes_t1 <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_diabetes_t1"] 
obese_comorbid_diabetes_t1 <-   obese_comorbid_diabetes_t1[order(comorbid_diabetes_t1)]

obese_comorbid_diabetes_t1 <- dplyr::left_join(obese_comorbid_diabetes_t1, N_comorbid_diabetes_t1)
obese_comorbid_diabetes_t1 <- obese_comorbid_diabetes_t1 %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_diabetes_t1 <- obese_comorbid_diabetes_t1 %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_diabetes_t1", .before=1) %>%
  dplyr::rename(group = comorbid_diabetes_t1) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_diabetes_t1_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_diabetes_t1, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_diabetes_t1_chisq <- dplyr::mutate (comorbid_diabetes_t1_chisq, variable = "comorbid_diabetes_t1") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_diabetes_t1 <- obese_comorbid_diabetes_t1 %>%
  dplyr::left_join(comorbid_diabetes_t1_chisq, by = "variable")


#########  comorbid_diabetes_t2 #######

#1.  count by comorbid_diabetes_t2
N_comorbid_diabetes_t2 <- BMI_complete_categories_2020_DT[, .N, by="comorbid_diabetes_t2"]
obese_comorbid_diabetes_t2 <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_diabetes_t2"] 
obese_comorbid_diabetes_t2 <-   obese_comorbid_diabetes_t2[order(comorbid_diabetes_t2)]

obese_comorbid_diabetes_t2 <- dplyr::left_join(obese_comorbid_diabetes_t2, N_comorbid_diabetes_t2)
obese_comorbid_diabetes_t2 <- obese_comorbid_diabetes_t2 %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_diabetes_t2 <- obese_comorbid_diabetes_t2 %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_diabetes_t2", .before=1) %>%
  dplyr::rename(group = comorbid_diabetes_t2) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_diabetes_t2_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_diabetes_t2, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_diabetes_t2_chisq <- dplyr::mutate (comorbid_diabetes_t2_chisq, variable = "comorbid_diabetes_t2") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_diabetes_t2 <- obese_comorbid_diabetes_t2 %>%
  dplyr::left_join(comorbid_diabetes_t2_chisq, by = "variable")



#########  comorbid_asthma #######

#1.  count by comorbid_asthma
N_comorbid_asthma <- BMI_complete_categories_2020_DT[, .N, by="comorbid_asthma"]
obese_comorbid_asthma <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_asthma"] 
obese_comorbid_asthma <-   obese_comorbid_asthma[order(comorbid_asthma)]

obese_comorbid_asthma <- dplyr::left_join(obese_comorbid_asthma, N_comorbid_asthma)
obese_comorbid_asthma <- obese_comorbid_asthma %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_asthma <- obese_comorbid_asthma %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_asthma", .before=1) %>%
  dplyr::rename(group = comorbid_asthma) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_asthma_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_asthma, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_asthma_chisq <- dplyr::mutate (comorbid_asthma_chisq, variable = "comorbid_asthma") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_asthma <- obese_comorbid_asthma %>%
  dplyr::left_join(comorbid_asthma_chisq, by = "variable")


#########  comorbid_COPD #######

#1.  count by comorbid_COPD
N_comorbid_COPD <- BMI_complete_categories_2020_DT[, .N, by="comorbid_COPD"]
obese_comorbid_COPD <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_COPD"] 
obese_comorbid_COPD <-   obese_comorbid_COPD[order(comorbid_COPD)]

obese_comorbid_COPD <- dplyr::left_join(obese_comorbid_COPD, N_comorbid_COPD)
obese_comorbid_COPD <- obese_comorbid_COPD %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_COPD <- obese_comorbid_COPD %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_COPD", .before=1) %>%
  dplyr::rename(group = comorbid_COPD) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_COPD_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_COPD, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_COPD_chisq <- dplyr::mutate (comorbid_COPD_chisq, variable = "comorbid_COPD") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_COPD <- obese_comorbid_COPD %>%
  dplyr::left_join(comorbid_COPD_chisq, by = "variable")

#########  comorbid_stroke_and_TIA #######

#1.  count by comorbid_stroke_and_TIA
N_comorbid_stroke_and_TIA <- BMI_complete_categories_2020_DT[, .N, by="comorbid_stroke_and_TIA"]
obese_comorbid_stroke_and_TIA <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_stroke_and_TIA"] 
obese_comorbid_stroke_and_TIA <-   obese_comorbid_stroke_and_TIA[order(comorbid_stroke_and_TIA)]

obese_comorbid_stroke_and_TIA <- dplyr::left_join(obese_comorbid_stroke_and_TIA, N_comorbid_stroke_and_TIA)
obese_comorbid_stroke_and_TIA <- obese_comorbid_stroke_and_TIA %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_stroke_and_TIA <- obese_comorbid_stroke_and_TIA %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_stroke_and_TIA", .before=1) %>%
  dplyr::rename(group = comorbid_stroke_and_TIA) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_stroke_and_TIA_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_stroke_and_TIA, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_stroke_and_TIA_chisq <- dplyr::mutate (comorbid_stroke_and_TIA_chisq, variable = "comorbid_stroke_and_TIA") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_stroke_and_TIA <- obese_comorbid_stroke_and_TIA %>%
  dplyr::left_join(comorbid_stroke_and_TIA_chisq, by = "variable")

#########  comorbid_chronic_cardiac #######

#1.  count by comorbid_chronic_cardiac
N_comorbid_chronic_cardiac <- BMI_complete_categories_2020_DT[, .N, by="comorbid_chronic_cardiac"]
obese_comorbid_chronic_cardiac <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_chronic_cardiac"] 
obese_comorbid_chronic_cardiac <-   obese_comorbid_chronic_cardiac[order(comorbid_chronic_cardiac)]

obese_comorbid_chronic_cardiac <- dplyr::left_join(obese_comorbid_chronic_cardiac, N_comorbid_chronic_cardiac)
obese_comorbid_chronic_cardiac <- obese_comorbid_chronic_cardiac %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_chronic_cardiac <- obese_comorbid_chronic_cardiac %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_chronic_cardiac", .before=1) %>%
  dplyr::rename(group = comorbid_chronic_cardiac) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_chronic_cardiac_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_chronic_cardiac, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_chronic_cardiac_chisq <- dplyr::mutate (comorbid_chronic_cardiac_chisq, variable = "comorbid_chronic_cardiac") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_chronic_cardiac <- obese_comorbid_chronic_cardiac %>%
  dplyr::left_join(comorbid_chronic_cardiac_chisq, by = "variable")


#########  comorbid_hypertension #######

#1.  count by comorbid_hypertension
N_comorbid_hypertension <- BMI_complete_categories_2020_DT[, .N, by="comorbid_hypertension"]
obese_comorbid_hypertension <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_hypertension"] 
obese_comorbid_hypertension <-   obese_comorbid_hypertension[order(comorbid_hypertension)]

obese_comorbid_hypertension <- dplyr::left_join(obese_comorbid_hypertension, N_comorbid_hypertension)
obese_comorbid_hypertension <- obese_comorbid_hypertension %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_hypertension <- obese_comorbid_hypertension %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_hypertension", .before=1) %>%
  dplyr::rename(group = comorbid_hypertension) %>%
  dplyr::mutate(group = as.character(group)) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_hypertension_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_hypertension, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_hypertension_chisq <- dplyr::mutate (comorbid_hypertension_chisq, variable = "comorbid_hypertension") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_hypertension <- obese_comorbid_hypertension %>%
  dplyr::left_join(comorbid_hypertension_chisq, by = "variable")


#########  comorbid_all_cancer #######

#1.  count by comorbid_all_cancer
N_comorbid_all_cancer <- BMI_complete_categories_2020_DT[, .N, by="comorbid_all_cancer"]
obese_comorbid_all_cancer <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="comorbid_all_cancer"] 
obese_comorbid_all_cancer <-   obese_comorbid_all_cancer[order(comorbid_all_cancer)]

obese_comorbid_all_cancer <- dplyr::left_join(obese_comorbid_all_cancer, N_comorbid_all_cancer)
obese_comorbid_all_cancer <- obese_comorbid_all_cancer %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_comorbid_all_cancer <- obese_comorbid_all_cancer %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "comorbid_all_cancer", .before=1) %>%
  dplyr::rename(group = comorbid_all_cancer) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
comorbid_all_cancer_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(comorbid_all_cancer, obese) %>%
  select(-1) %>% 
  chisq_test() 
comorbid_all_cancer_chisq <- dplyr::mutate (comorbid_all_cancer_chisq, variable = "comorbid_all_cancer") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_comorbid_all_cancer <- obese_comorbid_all_cancer %>%
  dplyr::left_join(comorbid_all_cancer_chisq, by = "variable")

## PRECOVID OBESE

#1.  count by precovid_obese_flag
N_precovid_obese_flag <- BMI_complete_categories_2020_DT[, .N, by="precovid_obese_flag"]
obese_precovid_obese_flag <- BMI_complete_categories_2020_DT[obese=="TRUE", .(n_obese= .N), by="precovid_obese_flag"] 
obese_precovid_obese_flag <-   obese_precovid_obese_flag[order(precovid_obese_flag)]

obese_precovid_obese_flag <- dplyr::left_join(obese_precovid_obese_flag, N_precovid_obese_flag)
obese_precovid_obese_flag <- obese_precovid_obese_flag %>%
  dplyr::mutate(proportion=n_obese/N) 

# 2. calculate confidence interval of propotions
obese_precovid_obese_flag <- obese_precovid_obese_flag %>%
  dplyr::mutate(lower_limit = (proportion - ((proportion*(1-proportion))/N*1.96))) %>%   # confidence interval of proporion
  dplyr::mutate(upper_limit = (proportion + ((proportion*(1-proportion))/N*1.96))) %>%
  dplyr::mutate(across(where(is.numeric), round, 4)) %>%
  dplyr::mutate(variable = "precovid_obese_flag", .before=1) %>%
  dplyr::rename(group = precovid_obese_flag) %>%
  dplyr::mutate(group = as.character(group))

#.... confidence interval proportion: (((proportion(1-proportion)/N))^0.5) * 1.96

# 3. chisq test
precovid_obese_flag_chisq <- as_tibble(BMI_complete_categories) %>%
  tabyl(precovid_obese_flag, obese) %>%
  select(-1) %>% 
  chisq_test() 
precovid_obese_flag_chisq <- dplyr::mutate (precovid_obese_flag_chisq, variable = "precovid_obese_flag") %>%
  dplyr::select("variable", "p", "method")

# 4.  Final table 
obese_precovid_obese_flag <- obese_precovid_obese_flag %>%
  dplyr::left_join(precovid_obese_flag_chisq, by = "variable")



obese_table <- obese_table %>%
  bind_rows(obese_age_group) %>%
  bind_rows (obese_sex) %>%
  bind_rows (obese_precovid_obese_flag) %>%
  bind_rows (obese_region) %>%
  bind_rows (obese_imd) %>%
  bind_rows (obese_ethnic_no_miss) %>%
  bind_rows (obese_eth_group_16) %>%   
  bind_rows (obese_comorbid_hypertension) %>%
  bind_rows (obese_comorbid_diabetes_t1) %>%
  bind_rows (obese_comorbid_diabetes_t2) %>%
  bind_rows (obese_comorbid_asthma) %>%
  bind_rows (obese_comorbid_COPD) %>%
  bind_rows (obese_comorbid_learning_disability) %>%
  bind_rows (obese_comorbid_depression) %>%
  bind_rows (obese_comorbid_psychosis_schiz_bipolar) %>%    
  bind_rows (obese_comorbid_dementia) %>%
  bind_rows (obese_comorbid_stroke_and_TIA) %>%
  bind_rows (obese_comorbid_chronic_cardiac) %>%
  bind_rows (obese_comorbid_all_cancer)

obese_table

write.csv (obese_table, here::here ("output/data","proportion_obese_2020.csv"))