diff --git a/scripts/data_analysis/0_before_after_datmod_prep.R b/scripts/data_analysis/0_before_after_datmod_prep.R
index f557e12..ef8886d 100644
--- a/scripts/data_analysis/0_before_after_datmod_prep.R
+++ b/scripts/data_analysis/0_before_after_datmod_prep.R
@@ -2,6 +2,12 @@
 temp.filename <- paste0(site, '_mod_dat.csv')
 dat <- read.csv(file.path('data_cached', temp.filename))
 
+# change infinite values to 365, representing a full year since application
+dat<-do.call(data.frame,lapply(dat, function(x) replace(x, is.infinite(x),365)))
+                               
+# change 'NaN' values to NA to omit during data processing
+dat<-do.call(data.frame,lapply(dat, function(x) replace(x, "NaN",NA)))
+                               
 # define predictors
 load('data_cached/modvars.Rdata')
 dat.mod <- dat[,predictors]
@@ -20,6 +26,9 @@ drop.predictors <- caret::findCorrelation(predictors.cor, cutoff = 0.95, verbose
 
 predictors.keep <- c(names.cor[-drop.predictors], 'frozen')
 
+# Change all 0 values to 0.0000001, so we can take the log10 of those vlaues
+dat.mod[dat.mod==0] <- 0.0000001
+                               
 # log transform response vars
 dat.mod[,responses] <- log10(dat.mod[,responses])
 sums <- colSums(dat.mod[,responses])