Difference between revisions of "Worksheets/Week6"
Jump to navigation
Jump to search
Kevin Dunn (talk | contribs) (→Part 1) |
Kevin Dunn (talk | contribs) (→Part 1) |
||
| Line 2: | Line 2: | ||
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. | ||
* A: enzyme strength: -1 == 20%; +1 == 30% | |||
* B: feed concentration: -1 == 5%; +1 == 15% | |||
* C: mixer type: -1 = R mixer; +1 = W mixer | |||
<html><div data-datacamp-exercise data-lang="r" data-height="auto"> | <html><div data-datacamp-exercise data-lang="r" data-height="auto"> | ||
| Line 15: | Line 19: | ||
y <- ... | y <- ... | ||
model.stability <- lm(y ~ A*B*C) | model.stability <- lm(y ~ A*B*C) | ||
summary(model.stability.half) | summary(model.stability.half) | ||
Revision as of 20:41, 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.
- A: enzyme strength: -1 == 20%; +1 == 30%
- B: feed concentration: -1 == 5%; +1 == 15%
- C: mixer type: -1 = R mixer; +1 = W mixer
# 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 <- ...
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)