Difference between revisions of "Worksheets/Week6"
Jump to navigation
Jump to search
Kevin Dunn (talk | contribs) (Created page with "=== Part 3 === 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...") |
Kevin Dunn (talk | contribs) (→Part 3) |
||
Line 1: | Line 1: | ||
=== Part | === 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. | 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. |
Revision as of 20:24, 31 March 2019
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.
# Define the 3 factors. This code is a template that you
# can easily extend and reuse for full factorial designs:
A <- c(-1, +1, -1, +1)
B <- c(-1, -1, +1, +1)
C <- A * B
# The response: stability [units=days]
y <- c...
# Linear model to predict stability from
# A: enzyme strength: -1 == 20%; +1 == 30%
# B: feed concentration: -1 == 5%; +1 == 15%
# C: mixer type: -1 = R mixer; +1 = W mixer
model.stability <- lm(y ~ A*B*C)
summary(model.stability.half)
# 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)