Worksheets/Week6

Part 1

Your group is developing a new product, but have been struggling to get the product’s stability, measured in days, to the level required. You are aiming for a stability value of 50 days or more.

• A: enzyme strength: -1 == 20%; +1 == 30%
• B: feed concentration: -1 == 5%; +1 == 15%
• C: mixer type: -1 = R mixer; +1 = W mixer

# This is the half-fraction, when C = A*B A <- c(-1, +1, -1, +1) B <- c(-1, -1, +1, +1) C <- A * B # The response: stability [units=days] y <- ... model.stability.poshalf <- lm(y ~ A*B*C) summary(model.stability.poshalf) # Uncomment this line if you run the code in RStudio #library(pid) # Comment these 2 lines if you run this code in RStudio source('https://yint.org/paretoPlot.R') source('https://yint.org/contourPlot.R') paretoPlot(model.stability.poshalf) # This is the other half-fraction, when C = -A*B A <- c(-1, +1, -1, +1) B <- c(-1, -1, +1, +1) C <- -1 * A * B # The response: stability [units=days] y <- ... model.stability.neghalf <- lm(y ~ A*B*C) summary(model.stability.neghalf)

Part 2

base <- c(-1, +1) design <- expand.grid(A=base, B=base, C=base) A <- design$A B <- design$B C <- design$C D <- A * B * C y <- c(60, 63, 70, 61, 44, 61, 94, 77) model.bio <- lm(y ~ A*B*C*D) summary(model.bio) # Uncomment this line if you run the code in RStudio #library(pid) # Comment these 2 lines if you run this code in RStudio source('https://yint.org/paretoPlot.R') source('https://yint.org/contourPlot.R') paretoPlot(model.bio)