Difference between revisions of "Worksheets/Week3"
Kevin Dunn (talk | contribs) (→Part 3) |
Kevin Dunn (talk | contribs) (→Part 4) |
||
| Line 123: | Line 123: | ||
</code> | </code> | ||
</div></html> | </div></html> | ||
=== Part 5 === | |||
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. | |||
<html><div data-datacamp-exercise data-lang="r" data-height="500px"> | |||
<code data-type="sample-code"> | |||
# Define the 3 factors. This code is a template that you | |||
# can easily extend and reuse for full factorial designs: | |||
base <- c(-1, +1) | |||
design <- expand.grid(A=base, B=base, C=base) | |||
A <- design$A | |||
B <- design$B | |||
C <- design$C | |||
# Type "A", and "B" and "C" at the command prompt | |||
# to verify what these letters contain. Are they in | |||
# standard order? | |||
# The response: stability, in number of days. | |||
y <- c(40, 27, 35, 21, 41, 27, 31, 20) | |||
# 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 | |||
mod.stability <- lm(y ~ A*B*C) | |||
summary(mod.stability) | |||
# Uncomment this line if you run the code in RStudio | |||
#library(pid) | |||
# Comment this line if you run this code in RStudio | |||
source('https://yint.org/paretoPlot.R') | |||
source('https://yint.org/contourPlot.R') | |||
paretoPlot(mod.stability) | |||
contourPlot(mod.stability, 'A', 'B') | |||
contourPlot(mod.stability, 'A', 'C') | |||
</code> | |||
</div></html> | |||
base <- c(-1, +1) | |||
design <- expand.grid(A=base, B=base, C=base) | |||
A <- design$A | |||
B <- design$B | |||
C <- design$C | |||
y <- c(40, 27, 35, 21, 41, 27, 31, 20) | |||
mod.stability <- lm(y ~ A*B*C) | |||
summary(mod.stability) | |||
paretoPlot(mod.stability) | |||
contourPlot(mod.stability, 'A', 'B') | |||
contourPlot(mod.stability, 'A', 'C') | |||
Revision as of 16:38, 10 March 2019
Part 1
A factorial experiment was run to investigate the settings that minimize the production of an unwanted side product. The two factors being investigated are called A and B for simplicity.
Part 2
Continuing from above, with 2 extra experimental points:
Part 3
Your family runs a small business selling products online. The first factor of interest is whether to provide free shipping over €30 or over €50. The second factor is whether or not the purchaser must first create a profile before completing the transaction. The purchaser can still complete their transaction without creating a profile. Below are the data collected, showing the 8 experiments.
Part 4
Continuing with the case above:
Part 5
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.
base <- c(-1, +1) design <- expand.grid(A=base, B=base, C=base) A <- design$A B <- design$B C <- design$C y <- c(40, 27, 35, 21, 41, 27, 31, 20) mod.stability <- lm(y ~ A*B*C) summary(mod.stability) paretoPlot(mod.stability) contourPlot(mod.stability, 'A', 'B') contourPlot(mod.stability, 'A', 'C')