#THIS IS AN R CODE SCRIPT


# Load required packages
library(rstudioapi)
library(readxl)
library(meta)


# Interactively select the Excel file
filepath <- "E:/SAVE R FILES HERE/Meta-analysis/Pip-tazo vs Cefepime/Extracted data.xlsx"

# Get sheet names from the Excel file
sheetsnames <- excel_sheets(filepath)
print(sheetsnames)

# Read all sheets into a named list of data frames
ma <- lapply(sheetsnames, function(x) {
  as.data.frame(read_excel(filepath, sheet = x))
})

# Assign sheet names to the list
names(ma) <- sheetsnames

# Select the outcome sheet
data <- ma$Mortality
print(data)

# Run meta-analysis for binary outcome using Risk Ratio (RR)
meta <- metabin(events_interv, n_interv, 
                events_ctrl, n_ctrl,
                studlab = Authors,
                data = data,
                sm = "RR",                
                method = "MH",            
                method.tau = "REML",      
                method.random.ci = "classic",
                random = TRUE,           
                common = FALSE,           
                prediction = FALSE)       
# Display meta-analysis results
summary(meta)


png(
  filename = "E:/SAVE R FILES HERE/Meta-analysis/Pip-tazo vs Cefepime/SUCCESS FOREST PLOT.png",
  width = 5000,   # wide enough for long study names
  height = 1500,
  res = 300
)


forest(meta,
       layout = "Revman",
       random = TRUE,
       common = FALSE,
       test.overall.random = TRUE,
       pooled.events = TRUE,
       label.e = "Pip-Tazo",
       label.c = "Cefepime",
       label.left = "Favors Pip-Tazo",
       label.right = "Favors Cefepime",
       leftcols = c("studlab",
                    "event.e","n.e",
                    "event.c","n.c",
                    "w.random",
                    "effect","ci"),
       leftlabs = c("Study", 
                    NA, NA, 
                    NA, NA,
                    "Weight (%)", 
                    "Effect size", "95% CI"),
       
       ff.xlab = "bold",
       xlim = c(0.1, 80),
       col.square = "darkblue",
       col.square.lines = "black",
       col.diamond.random = "white",
       col.diamond.lines.random = "black",
       colgap = "5mm",
       colgap.studlab = "10mm",
       colgap.forest.left = "5mm",
       just = "center")


dev.off()