diff --git a/R/doeAnalysis.R b/R/doeAnalysis.R
index 3425ba8f..cd7b28ca 100644
--- a/R/doeAnalysis.R
+++ b/R/doeAnalysis.R
@@ -108,89 +108,6 @@ doeAnalysis <- function(jaspResults, dataset, options, ...) {
   return(deps)
 }
 
-# dataset <- read.csv("C:/Users/Jonee/Google Drive/JASP/SKF Six Sigma/JASP Data Library/4_3_AnalyzeDesign/ResponseSurfaceDesignAnalysis.csv")
-# dataset <- dataset[4:7]
-# options <- list()
-# options[["continuousFactorsResponseSurface"]] <- c("Inlet_feeding")
-#
-# options[["dependent"]] <- "Vdk"
-# options[["fixedFactorsResponseSurface"]] <- c("Time", "Oil_temperature")
-#
-# options$modelTerms <- list(list(components = "Inlet_feeding"), list(components = "Time"),
-#                           list(components = "Oil_temperature"),
-#                           list(components = c("Inlet_feeding", "Time")),
-#                           list(components = c("Time", "Oil_temperature")),
-#                           list(components = c("Inlet_feeding", "Oil_temperature")))
-# dataset$Time <- as.factor(dataset$Time)
-# dataset$Oil_temperature <- as.factor(dataset$Oil_temperature)
-# colnames(dataset) <- jaspBase::encodeColNames(colnames(dataset))
-#
-# options$modelTerms <- list(list(components = "Inlet_feeding"), list(components = "Time"),
-#                            list(components = "Oil_temperature"), list(components = c("Inlet_feeding",
-#                                                                                     "Time")), list(components = c("Time", "Oil_temperature"
-#                                                                                     )), list(components = c("Inlet_feeding", "Oil_temperature")),
-#                            list(components = c("Inlet_feeding", "Time", "Oil_temperature"
-#                            )))
-# options$squaredTerms <- list(components = c("Oil_temperature"))
-#
-# testModel <- lm(Vdk ~ Inlet_feeding + Time + Oil_temperature + I(Oil_temperature^2), data = dataset)
-# testModel2 <- lm(Vdk ~ Inlet_feeding + Time + Oil_temperature, data = dataset)
-# car::Anova(testModel)
-# car::Anova(testModel2)
-
-#
-#
-# options[["continuousFactorsResponseSurface"]] <- c("Inlet_feeding")
-# options$modelTerms <- list(list(components = "Inlet_feeding"))
-
-# dataset <- read.csv("C:/Users/Jonee/Google Drive/JASP/SKF Six Sigma/JASP Data Library/4_3_AnalyzeDesign/FactorialDesignAnalysis.csv")
-# options <- list()
-# dataset <- dataset[5:8]
-# dataset[1] <- as.factor(dataset[[1]])
-# dataset[2] <- as.factor(dataset[[2]])
-# dataset[3] <- as.factor(dataset[[3]])
-# options[["continuousFactorsResponseSurface"]] <- NULL
-# options[["dependent"]] <- "Yield"
-# options[["fixedFactors"]] <- c("Exposure_time", "Develop_time", "Mask_dimension")
-# options$modelTerms <- list(list(components = "Exposure_time"),
-#                            list(components = "Develop_time"),
-#                            list(components = "Mask_dimension"),
-#                            list(components = c("Mask_dimension", "Exposure_time")))
-#
-# dataset <- read.csv("C:/Users/Jonee/Google Drive/JASP/SKF Six Sigma/JASP Data Library/4_3_AnalyzeDesign/FactorialDesignAnalysis.csv")
-# options <- list()
-# dataset <- dataset[5:8]
-# dataset[1] <- as.factor(dataset[[1]])
-# dataset[3] <- as.factor(dataset[[3]])
-# options[["continuousFactorsResponseSurface"]] <- NULL
-# options[["dependent"]] <- "Yield"
-# options[["fixedFactorsFactorial"]] <- c("Exposure_time", "Mask_dimension")
-# options[["covariates"]] <- c("Develop_time")
-# options$modelTerms <- list(list(components = "Exposure_time"),
-#                            list(components = "Mask_dimension"),
-#                            list(components = c("Mask_dimension", "Exposure_time")))
-#
-#
-# dataset <- read.csv("C:/Users/Jonee/Google Drive/JASP/SKF Six Sigma/Datasets/DOE_FAC_withBlocks.csv", sep = ",")
-# options <- list()
-# dataset <- dataset[4:8]
-# dataset[1] <- as.factor(dataset[[1]])
-# dataset[2] <- as.factor(dataset[[2]])
-# dataset[3] <- as.factor(dataset[[3]])
-# dataset[4] <- as.factor(dataset[[4]])
-# options[["continuousFactorsResponseSurface"]] <- NULL
-# options[["dependent"]] <- "Response"
-# options[["fixedFactorsFactorial"]] <- c("A", "B", "C")
-# options[["blocks"]] <- c("Blocks")
-# options$modelTerms <- list(list(components = "A"),
-#                            list(components = "B"),
-#                            list(components = "C"),
-#                            list(components = c("A", "B")),
-#                            list(components = c("A", "C")),
-#                            list(components = c("B", "C")),
-#                            list(components = c("A", "B", "C")))
-
-
 .scaleDOEvariable <- function(x){2*(x-min(x))/(max(x)-min(x))-1}
 
 .doeAnalysisMakeState <- function(jaspResults, dataset, options, continuousPredictors, discretePredictors, blocks, covariates, dependent, ready) {