Implementing R-code for a Lack of Fit test

Hi,

Is it possible to include a R-code for a Lack of Fit test in the existing R-code of a analysis in Quality Control>>Design of Experiments>>Analyse Desing>>Response surface design?

Tank you in advance.

Current code:

jaspQualityControl::doeAnalysis(

     data = NULL,

     version = "0.97.1",

     blocksFactorial = NULL,

     blocksResponseSurface = NULL,

     codeFactors = TRUE,

     codeFactorsManualTable = list(optionKey = "predictors", types = list("scale", "scale", "scale"), value = list(list(highValue = 1, lowValue = -1, predictors = "Temperatuur"), list(highValue = 1, lowValue = -1, predictors = "Druk"), list(highValue = 1, lowValue = -1, predictors = "Holdtime"))),

     codeFactorsMethod = "automatic",

     continuousFactorsFactorial = NULL,

     continuousFactorsResponseSurface = ~ Temperatuur + Druk + Holdtime,

     contourSurfacePlot = TRUE,

     contourSurfacePlotResponseDivision = 5,

     contourSurfacePlotType = "surfacePlot",

     contourSurfacePlotVariables = ~ Temperatuur + Holdtime,

     covariates = NULL,

     dependentFactorial = NULL,

     dependentResponseSurface = ~ `Gehalte PAKs`,

     designType = "responseSurfaceDesign",

     fixedFactorsFactorial = NULL,

     fixedFactorsResponseSurface = NULL,

     fourInOneResidualPlot = TRUE,

     highestOrder = TRUE,

     histogramBinWidthType = "sturges",

     histogramManualNumberOfBins = 30,

     modelTerms = ~ Temperatuur + Holdtime,

     normalEffectsPlot = FALSE,

     optimizationPlot = FALSE,

     optimizationPlotCustomParameterValues = list(optionKey = "variable", types = list("scale", "scale", "scale"), value = list(list(value = "", variable = "Temperatuur"), list(value = "", variable = "Druk"), list(value = "", variable = "Holdtime"))),

     optimizationPlotCustomParameters = FALSE,

     optimizationPlotPredictionType = "response",

     optimizationSolutionTable = TRUE,

     order = 2,

     plotFitted = FALSE,

     plotHeight = 320,

     plotHist = FALSE,

     plotNorm = FALSE,

     plotPareto = TRUE,

     plotRunOrder = FALSE,

     plotWidth = 480,

     responseOptimizerManualBounds = FALSE,

     responseOptimizerManualTarget = FALSE,

     responsesResponseOptimizer = list(optionKey = "variable", types = list(), value = list()),

     rsmPredefinedModel = FALSE,

     rsmPredefinedTerms = "linearAndSquared",

     squaredTerms = ~ Temperatuur + Holdtime,

     squaredTermsCoded = TRUE,

     stepwiseMethodFactorial = "enter",

     stepwiseMethodResponseSurface = "enter",

     sumOfSquaresType = "type3",

     surfacePlotHorizontalRotation = 330,

     surfacePlotVerticalRotation = 20,

     tableAlias = FALSE,

     tableEquation = TRUE)