mf0412
. One of the factors has 3 levels and the
other three have two levels each. The design consists of 12 runs and
is a 1/2 fraction of the full factorial. The response is the percent of
DNA liberated compared to a gold standard.
Experimental data provided by W. Keating of the Becton Dickinson Research Center, Research Triangle Park, NC. Used by permission of Becton Dickinson.
# This design is already available in S-PLUS # under the name lysis.df. The following commands # were used to create the data frame: lysis.names<- list(temp=c(25,60,95),sonic=c(0,5), edta=c(0,1),bme=c(0,5)) lysis.design <-design.digest('mf0412', lysis.names) lysis.gpct <- c(37.2,32.3,51.4,40.8,74.9,74.4,61.3, 86.3,87.4,100.2,139.3,153.2) lysis.df <- cbind(lysis.design,gpct=lysis.gpct) # sample analysis summary(lysis.df) lysis.fac <- fac.aov(lysis.df) summary(lysis.fac) pareto(lysis.fac) qqnorm(lysis.fac) tfiplot(lysis.fac,~temp:sonic) acplot(lysis.fac) acplot(lysis.fac,.4,5) ebplot(lysis.fac) lysis.rmod <- update(lysis.fac,. ~temp*sonic) pareto(lysis.rmod) plot(lysis.rmod)