In [1]:
library(gdata)
library(ggplot2)
library(grid)
library(gridExtra)


gdata: read.xls support for 'XLS' (Excel 97-2004) files ENABLED.

gdata: read.xls support for 'XLSX' (Excel 2007+) files ENABLED.

Attaching package: ‘gdata’

The following object is masked from ‘package:stats’:

    nobs

The following object is masked from ‘package:utils’:

    object.size

The following object is masked from ‘package:base’:

    startsWith


Attaching package: ‘gridExtra’

The following object is masked from ‘package:gdata’:

    combine


In [3]:
source("C.2.5- Within_regions_generic_plot.R")

In [4]:
DT <- read.table("../Data/All_data.txt")
ratio_align <- read.table("../Data/Alignment_ratios_within_regions_across_diseases_wt_sims_patients_metrs_burdens.txt")

In [5]:
DT$Dis_lab <- DT$Disease
levels(DT$Dis_lab)


  1. 'All'
  2. 'Cardiovascular and circulatory diseases'
  3. 'Chronic respiratory diseases'
  4. 'Cirrhosis of the liver'
  5. 'Congenital anomalies'
  6. 'Diabetes, urinary diseases and male infertility'
  7. 'Diarrhea, lower respiratory infections, meningitis, and other common infectious diseases'
  8. 'Digestive diseases (except cirrhosis)'
  9. 'Gynecological diseases'
  10. 'Hemoglobinopathies and hemolytic anemias'
  11. 'Hepatitis'
  12. 'HIV/AIDS'
  13. 'Leprosy'
  14. 'Malaria'
  15. 'Maternal disorders'
  16. 'Mental and behavioral disorders'
  17. 'Musculoskeletal disorders'
  18. 'Neglected tropical diseases excluding malaria'
  19. 'Neonatal disorders'
  20. 'Neoplasms'
  21. 'Neurological disorders'
  22. 'Nutritional deficiencies'
  23. 'Oral disorders'
  24. 'Sense organ diseases'
  25. 'Sexually transmitted diseases excluding HIV'
  26. 'Skin and subcutaneous diseases'
  27. 'Sudden infant death syndrome'
  28. 'Tuberculosis'

In [6]:
#Disease labels for plot
levels(DT$Dis_lab) <- c("",
               "Cardiovasc.\nand circulatory",
               "Chronic\nrespiratory",
               "Cirrhosis",
               "Congenital\nanomalies",
               "Diabetes, urinary\nmale infertility",
               "Common\ninfect. dis.",
               "Digestive",
               "Gynecol.",
               "Hemoglob. and\nhemolytic anemia",
               "Hepatitis",
               "HIV",
               "Leprosy",
               "Malaria",
               "Maternal\ndisorders",
               "Mental and\nbehavioral",
               "Musculosk.",
               "Neglected trop.",
               "Neonatal\ndisorders",
               "Neoplasms",
               "Neurological",
               "Nutritional",
               "Oral",
               "Sense organ",
               "STD",
               "Skin and\nsubcutaneous",
               "Sudden infant death",
               "Tuberculosis")

In [7]:
#Region labels
DT$regs_lab <- DT$Region
levels(DT$regs_lab)


  1. 'All'
  2. 'Central Europe, Eastern Europe, and Central Asia'
  3. 'High-income'
  4. 'Latin America and Caribbean'
  5. 'Non-HI'
  6. 'North Africa and Middle East'
  7. 'South Asia'
  8. 'Southeast Asia, East Asia and Oceania'
  9. 'Sub-Saharian Africa'

In [8]:
levels(DT$regs_lab) <-  c("World",
                          "Eastern Europe and Central Asia",
                          "High-income countries",
                          "Latin America and Caribbean",
                          "Non-high-income countries",
                          "North Africa and Middle East",
                          "South Asia",
                          "Southeast Asia, East Asia and Oceania",
                          "Sub-Saharian Africa")

In [9]:
L <- Within_2Regions_plot(metr_burden = "yld",metr_res = "RCTs",region1 = "High-income",region2="Sub-Saharian Africa")

In [10]:
options(repr.plot.width=10, repr.plot.height=5)
grid.arrange(L[[1]],L[[2]], ncol=2)



In [12]:
options(repr.plot.width=5, repr.plot.height=5)
Within_Region_plot(metr_burden = "yld",metr_res = "Patients",region = "South Asia")



In [ ]: