diff --git a/R/Global.R b/R/Global.R
index 06fad47..1e6ee0b 100644
--- a/R/Global.R
+++ b/R/Global.R
@@ -2,6 +2,28 @@
 # All functions for plateflow
 ################################################################################
 
+library(shiny)
+library(shinydashboard)
+library(shinyjs)
+library(shinyAce)
+library(shinydashboard)
+library(shinycssloaders)
+library(shinyBS)
+library(purrr)
+library(gslnls)
+library(tidyverse)
+library(ggplot2)
+library(reshape2)
+library(openxlsx)
+library(DT)
+library(ggpubr)
+library(gridExtra)
+library(drc)
+library(twopartm)
+library(car)
+library(dplyr)
+library(scales)
+
 
 #' Levenberg Marquard fit of 4 pl
 #'
@@ -579,7 +601,7 @@ LinPotTab <- function(circles, Lim, PureErrFlag) {
     SEsPure <- sqrt(diag(V_V) * meanPureErr)
   }
 
-  log_pot_delta <- deltaMethod(modAB, "isSample/log_dose")
+  log_pot_delta <- car::deltaMethod(modAB, "isSample/log_dose")
   if (PureErrFlag) {
     V_ <- log_pot_delta$SE^2 / SU_modAB$sigma^2
     V_p <- V_ * meanPureErr
@@ -658,7 +680,7 @@ ANOVAlintests <- function(ro_new, circles, Lim, PureErrFlag) {
   modA <- lm(readout ~ log_dose + isSample, circ_Al)
   modB <- lm(readout ~ log_dose + isSample, circ_Bl)
 
-  log_pot_delta <- deltaMethod(modAB, "isSample/log_dose")
+  log_pot_delta <- car::deltaMethod(modAB, "isSample/log_dose")
   CI_log_low <- log_pot_delta$Estimate - qt(0.975, df.residual(modAB)) * log_pot_delta$SE
   CI_log_up <- log_pot_delta$Estimate + qt(0.975, df.residual(modAB)) * log_pot_delta$SE
   ExpLinPot <- exp(c(log_pot_delta$Estimate, CI_log_low, CI_log_up))