daewr
daewr
daewr
September 9, 2023
Type Package
Title Design and Analysis of Experiments with R
Version 1.2-11
Date 2023-09-04
Maintainer John Lawson <lawsonjsl7net@gmail.com>
Description Contains Data frames and functions used in the book ``Design and Analysis of Experi-
ments with R'', Lawson(2015) ISBN-13:978-1-4398-6813-3.
License GPL-2
Depends R (>= 3.5.0)
Encoding UTF-8
LazyLoad true
LazyData true
Imports stringi,stats,graphics,grDevices,lattice
RoxygenNote 7.2.3
NeedsCompilation no
Author John Lawson [aut, cre],
Gerhard Krennrich [aut],
Ruben Amoros [ctr]
Repository CRAN
Date/Publication 2023-09-09 08:20:07 UTC
R topics documented:
daewr-package . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
Altscreen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
antifungal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Apo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
apple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
arso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
augm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1
2 R topics documented:
Bdish . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Bff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
bha . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
BIBsize . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
bioequiv . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
bioeqv . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
blood . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
BoxM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
BPmonitor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
bread . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
cakeb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
cement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
chem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
chipman . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
COdata . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
colormap . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
connector . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
cont . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
cpipe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
culture . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
dairy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
DefScreen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
drug . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
EEw1s1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
EEw1s2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
EEw1s3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
EEw2s1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
EEw2s2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
EEw2s3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
EEw3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
eptaxr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
eptaxs2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
eptaxyb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Fcrit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
fhstep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
fhstepDS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
firstm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
FitDefSc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
fnextrm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
fntrmDS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Fpower . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Fpower1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Fpower2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
fullnormal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
gagerr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
gapstat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Gaptest . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
gear . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
R topics documented: 3
halfnorm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
hardwood . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
HierAFS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
ihstep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
inject . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
interleave . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
LenthPlot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
LGB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
LGBc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
mod . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
ModelRobust . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
MPV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
Naph . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
OptPB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
pastry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
PBDes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
pest . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
pesticide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
plasma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
polvdat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
polymer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
prodstd . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
qsar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
Rations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
rcb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
residue . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
rubber . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
sausage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
Smotor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
soup . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
soupmx . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
splitPdes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
SPMPV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
stdord . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
strung . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
strungtile . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
sugarbeet . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
taste . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
teach . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
Tet . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
tile . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
Treb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
Tukey1df . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
vci . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
vinyl . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
virus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
volt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
web . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
4 Altscreen
WeldS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
Index 77
daewr-package Data frames and functions for Design and Analysis of Experiments
with R
Description
This package contains the data sets and functions from the book Design and Analysis of Experi-
ments with R published by CRC in 2013.
Author(s)
John Lawson
Maintainer: John Lawson <lawsonjsl7net@gmail.com>
References
J. Lawson, Design and Analysis of Experiments with R, CRC 2013.
Description
Recalls Jones and Montgomery’s 16 run screening designs from data frames
Usage
Altscreen(nfac, randomize=FALSE)
Arguments
nfac input- an integer
randomize input - logical
Value
a data frame containing the alternate screening design
Author(s)
John Lawson
antifungal 5
References
Description
Data from the Two-period crossover study of an antifungal agent in chapter 9 of Design and Analysis
of Experiments with R
Usage
data(antifungal)
Format
Source
Examples
data(antifungal)
6 apple
Description
Data from the apolipoprotein survey variance component study of Chapter 5 in Design and Analysis
of Experiments with R
Usage
data(Apo)
Format
A data frame with 30 observations on the following 2 variables.
lab a factor with levels A B C D
conc a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(Apo)
Description
Data from the confounded apple slice browning experiment in chapter 7 of Design and Analysis of
Experiments with R
Usage
data(apple)
Format
A data frame with 24 observations on the following 4 variables.
Block a factor with levels 1 2 3 4
A a factor with levels 0 1 2 3
B a factor with levels 0 1 2
rating a numeric vector containing the response
arso 7
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(apple)
Description
Data from the 2(7−3) arsenic removal experiment in chapter 6 of Design and Analysis of Experi-
ments with R
Usage
data(arso)
Format
A data frame with 8 observations on the following 8 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(arso)
8 augm
Description
Data from the 2(7−3) arsenic removal experiment augmented with mirror image in chapter 6 of
Design and Analysis of Experiments with R
Usage
data(augm)
Format
Source
Examples
data(augm)
Bdish 9
Description
Data from the Confounded Block Dishwashing Experiment in chapter 7 of Design and Analysis of
Experiments with R
Usage
data(Bdish)
Format
A data frame with 16 observations on the following 5 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(Bdish)
Description
Data from the Confounded block fractional factorial mouse growth experiment in chapter 7 of
Design and Analysis of Experiments with R
Usage
data(Bff)
10 bha
Format
A data frame with 16 observations on the following 5 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(Bff)
Description
Data from the mouse liver enzyme experiment in chapter 4 of Design and Analysis of Experiments
with R
Usage
data(bha)
Format
A data frame with 16 observations on the following 4 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(bha)
Description
This function computes the number of blocks, treatment frequency and lambda for a potential BIB
design
Usage
BIBsize(t,k)
Arguments
t input - number of levels of the treatment factor
k input - blocksize or number of experimental units per block
Value
a list containing the b=number of blocks, r=number of treatment replicates and lambda for a poten-
tial BIB design with t levels of treatment factor and blocksize k.
Author(s)
John Lawson
Description
Data from the extra-period crossover bioequivalence study in chapter 9 of Design and Analysis of
Experiments with R
Usage
data(bioequiv)
12 bioeqv
Format
A data frame with 108 observations on the following 5 variables.
Group a factor with levels 1 2
Subject a factor with levels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 21 23 24 25 26 27 28
30 31 32 33 34 35 36 120 122 129
Period a factor with levels 1 2 3
Treat a factor with levels A B
Carry a factor with levels none A B
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(bioequiv)
Description
Data from the Latin Square bioequivalence experiment in chapter 4 of Design and Analysis of
Experiments with R
Usage
data(bioeqv)
Format
A data frame with 9 observations on the following 4 variables.
Period a factor with levels 1 2 3
Subject a factor with levels 1 2 3
Treat a factor with levels A B C
AUC a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(bioeqv)
blood 13
Description
Data from the Variance component study of calcium in blood serum in chapter 5 of Design and
Analysis of Experiments with R
Usage
data(blood)
Format
Source
Examples
data(blood)
Description
Data from Box and Meyer’s unreplicated 24 in chapter 3 of Design and Analysis of Experiments
with R
Usage
data(BoxM)
14 BPmonitor
Format
A data frame with 16 observations on the following 4 variables.
A a numeric vector containing the coded (-1,1) levels of factor A
B a numeric vector containing the coded (-1,1) levels of factor B
C a numeric vector containing the coded (-1,1) levels of factor C
D a numeric vector containing the coded (-1,1) levels of factor D
y a numeric vector containing the response
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
References
Box, G. E. P. "George’s Column", Quality Engineering, Vol. 3, pp. 405-410.
Examples
data(BoxM)
Description
Data from the blood pressure monitor experiment experiment in Chapter 7 of Design and Analysis
of Experiments with R
Usage
data(BPmonitor)
Format
A data frame with 12 observations on the following 3 variables.
Block a factor with levels 1 2 3 4 5 6
Treatment a factor with levels "P" "A" "B" "C"
pressure a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(BPmonitor)
bread 15
Description
Data from the bread rise experiment in chapter 2 of Design and Analysis of Experiments with R
Usage
data(bread)
Format
Source
Examples
data(bread)
Description
Data from the Split-Plot response surface for cake baking experiment in chapter 10 of Design and
Analysis of Experiments with R
Usage
data(cakeb)
16 cement
Format
A data frame with 11 observations on the following 6 variables.
Ovenrun a factor with levels 1 2 3 4
x1 a numeric vector
x2 a numeric vector
y a numeric vector
x1sq a numeric vector
x2sq a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(cakeb)
Description
Data from the CCD design for cement workability experiment in chapter 10 of Design and Analysis
of Experiments with R
Usage
data(cement)
Format
A data frame with 20 observations on the following 4 variables.
Block a factor with levels 1 2
x1 a numeric vector
x2 a numeric vector
x3 a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(cement)
chem 17
Description
Data from the Chemical process experiment in chapter 3 of Design and Analysis of Experiments
with R
Usage
data(chem)
Format
A data frame with 16 observations on the following 4 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(chem)
Description
Data from the Williams’ crossover design for sprinting experiment in chapter 9 of Design and
Analysis of Experiments with R
Usage
data(chipman)
18 COdata
Format
A data frame with 36 observations on the following 5 variables.
Square a factor with levels 1 2
Group a factor with levels 1 2 3
Subject a factor with levels 1 2 3 4 5 6 7 8 9 10 11 12
Period a factor with levels 1 2 3
Treat a factor with levels 1 2 3
Carry a factor with levels 0 1 2 3
Time a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(chipman)
Description
Data from the CO emissions experiment in chapter 3 of Design and Analysis of Experiments with
R
Usage
data(COdata)
Format
A data frame with 18 observations on the following 3 variables.
Eth a factor with levels 0.1 0.2 0.3
Ratio a factor with levels 14 15 16
CO a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(COdata)
colormap 19
Description
This function makes a colormap of the correlations of a design matrix stored in the data frame
design
Usage
colormap(design, mod)
Arguments
design input - a data frame containing columns of the numeric factor levels
mod input - a number indicationg the model for the colormap 1 = linear model con-
taining only the terms in the dataframe 2 = linear model plus two factor interac-
tions 3 = linear model plus 2 and 3 factor interactions 4 = linear model plus 2,
3, and 4 factor interactions
Author(s)
John Lawson
Examples
## The function is currently defined as
function(design,mod) {
##################### Inputs ###########################################
# design - a data frame containing columns of the numeric factor levels
# mod - the model for the color plot of correlations
# 1 = Linear model containing only the terms in the data frame
# 2 = Linear model plus two factor interactions
# 3 = Linear model plus 2 and 3 factor interactions
# 4 = Linear model plus 2, 3 and 4 factor interactions
########################################################################
y<-runif(nrow(design),0,1)
if(mod==1) {test <- model.matrix(lm(y~(.),data=design))}
if(mod==2) {test <- model.matrix(lm(y~(.)^2,data=design))}
if(mod==3) {test <- model.matrix(lm(y~(.)^3,data=design))}
if(mod==4) {test <- model.matrix(lm(y~(.)^4,data=design))}
names<-colnames(test)
names<-gsub(':','',names)
names<-gsub('1','',names)
colnames(test)<-names
cmas<-abs(cor(test[,ncol(test):2]))
cmas<-cmas[c((ncol(cmas)):1), ]
rgb.palette <- colorRampPalette(c("white", "black"), space = "rgb")
levelplot(cmas, main="Map of absolute correlations", xlab="", ylab="", col.regions=rgb.palette(120),
cuts=100, at=seq(0,1,0.01),scales=list(x=list(rot=90))) }
20 cont
Description
Data from the Single Array Experiment with an Elastometric Connector in Chapter 12 of Design
and Analysis of Experiments with R. The control and noise factors are in coded levels.
Usage
data(connector)
Format
A data frame with 32 observations on the following 8 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(connector)
cont Control factor array and summary statistics for controller circuit de-
sign experiment
Description
Data from the control factor array and summary statistics for controller circuit design experiment
in chapter 12 of Design and Analysis of Experiments with R
Usage
data(cont)
cpipe 21
Format
A data frame with 18 observations on the following 6 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
F a numeric vector
lns2 a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(cont)
Description
Data from the Split-plot response surface for ceramic pipe experiment in chapter 10 of Design and
Analysis of Experiments with R
Usage
data(cpipe)
Format
A data frame with 48 observations on the following 6 variables.
WP a factor with levels 1 2 3 4 5 6 7 8 9 10 11 12
A a numeric vector
B a numeric vector
P a numeric vector
Q a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(cpipe)
22 culture
Description
Data from the paecilomyces variotii culture experiment experiment in chapter 6 of Design and
Analysis of Experiments with R
Usage
data(culture)
Format
Source
Examples
data(culture)
dairy 23
Description
Data from the Repeated measures study with dairy cow diets in chapter 9 of Design and Analysis
of Experiments with R (compact format)
Usage
data(dairy)
Format
A data frame with 120 observations on the following 5 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(dairy)
Description
Recalls Jones and Nachtsheim’s Definitive screening designs for 3-level factors and 3-level factors
with added 2-level categorical factors.
Usage
DefScreen(m, c=0, center=0, randomize=FALSE)
24 drug
Arguments
m input- an integer, the m=number of 3-level factors
c input- an integer, the m=number of 2-level categorical factors, default is zero if
not supplied
center input- an integer, the number of extra center points. This must be zero when c>0
randomize input - logical
Value
a data frame containing the definitive screening design with 3-level factors first followed by 2-level
factors.
Author(s)
John Lawson
References
Jones, B. and Nachtsheim, C. J. (2011) "A Class of Three Level Designs for Definitive Screening
in the Presence of Second-Order Effects", Journal of Quality Technology, Vol 43, No. 1, 2011, pp
1-15. Jones, B. and Nachtsheim, C. J. (2013) "Definitive Screening Designs with Added Two-Level
Categorical Factors", Journal of Quality Technology, Vol 44, No. 2, 2013, pp. 121-129.
Description
Data from rat behavior experiment in Chapter 4 of Design and Analysis of Experiments with R
Usage
data(drug)
Format
A data frame with 50 observations on the following 3 variables.
rat a factor with levels 1 2 3 4 5 6 7 8 9 10
dose a factor with levels 0.0 0.5 1.0 1.5 2.0
rate a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(drug)
EEw1s1 25
Description
Recalls Jones and Goos JQT Estimation Equivalent Response Surface Designs for 1 whole plot
factor and 1 sub-plot factor from a catalog
Usage
EEw1s1(des, randomize=FALSE)
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Jones, B. and Goos, P.(2012) "An Algorithm for Finding D-Efficient Equivalent-Estimation Second-
Order Split Plot Designs", Journal of Quality Technology, Vol 44, No. 4, pp281-303, 2012.
Description
Recalls Jones and Goos JQT Estimation Equivalent Response Surface Designs for 1 whole plot
factor and 2 sub-plot factors from a catalog
Usage
EEw1s2(des, randomize=FALSE)
26 EEw1s3
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Jones, B. and Goos, P.(2012) "An Algorithm for Finding D-Efficient Equivalent-Estimation Second-
Order Split Plot Designs", Journal of Quality Technology, Vol 44, No. 4, pp281-303, 2012.
Description
Recalls Jones and Goos JQT Estimation Equivalent Response Surface Designs for 1 whole plot
factor and 3 sub-plot factors from a catalog
Usage
EEw1s3(des, randomize=FALSE)
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Jones, B. and Goos, P.(2012) "An Algorithm for Finding D-Efficient Equivalent-Estimation Second-
Order Split Plot Designs", Journal of Quality Technology, Vol 44, No. 4, pp281-303, 2012.
EEw2s1 27
Description
Recalls Jones and Goos JQT Estimation Equivalent Response Surface Designs for 2 whole plot
factors and 1 sub-plot factor from a catalog
Usage
EEw2s1(des, randomize=FALSE)
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Jones, B. and Goos, P.(2012) "An Algorithm for Finding D-Efficient Equivalent-Estimation Second-
Order Split Plot Designs", Journal of Quality Technology, Vol 44, No. 4, pp281-303, 2012.
Description
Recalls Jones and Goos JQT Estimation Equivalent Response Surface Designs for 2 whole plot
factors and 1 sub-plot factor from a catalog
Usage
EEw2s2(des, randomize=FALSE)
28 EEw2s3
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Jones, B. and Goos, P.(2012) "An Algorithm for Finding D-Efficient Equivalent-Estimation Second-
Order Split Plot Designs", Journal of Quality Technology, Vol 44, No. 4, pp281-303, 2012.
Description
Recalls Jones and Goos JQT Estimation Equivalent Response Surface Designs for 2 whole plot
factors and 1 sub-plot factor from a catalog
Usage
EEw2s3(des, randomize=FALSE)
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Jones, B. and Goos, P.(2012) "An Algorithm for Finding D-Efficient Equivalent-Estimation Second-
Order Split Plot Designs", Journal of Quality Technology, Vol 44, No. 4, pp281-303, 2012.
EEw3 29
Description
Recalls Jones and Goos JQT Estimation Equivalent Response Surface Designs for 3 whole plot
factors and 1-2 sub-plot factors from a catalog
Usage
EEw3(des, randomize=FALSE)
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Jones, B. and Goos, P.(2012) "An Algorithm for Finding D-Efficient Equivalent-Estimation Second-
Order Split Plot Designs", Journal of Quality Technology, Vol 44, No. 4, pp281-303, 2012.
eptaxr Single array and raw response for silicon layer growth experiment
Description
Data from the single array and raw response for silicon layer growth experiment in chapter 12 of
Design and Analysis of Experiments with R
Usage
data(eptaxr)
30 eptaxs2
Format
A data frame with 64 observations on the following 9 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
H a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(eptaxr)
eptaxs2 Control array and variance of response for silicon layer growth exper-
iment
Description
Data from the control array and variance of response for silicon layer growth experiment in chapter
12 of Design and Analysis of Experiments with R
Usage
data(eptaxs2)
Format
A data frame with 16 observations on the following 9 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
H a numeric vector
s2 a numeric vector
eptaxyb 31
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(eptaxs2)
eptaxyb Control array and mean response for silicon layer growth experiment
Description
Data from the control array and mean response for silicon layer growth experiment in chapter 12 of
Design and Analysis of Experiments with R
Usage
data(eptaxyb)
Format
A data frame with 16 observations on the following 9 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
H a numeric vector
ybar a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(eptaxyb)
32 fhstep
Description
Usage
Arguments
Value
Author(s)
John Lawson
Description
This function performs a single step of a hierarchical forward stepwise regression by entering ad-
ditional term(s) to a model already created by ihstep or fhstep. If an interaction or quadratic term
is entered first, the parent main effects are also entered into the model. This function is called by
HierAFS.R
Usage
fhstep(y,des,m,c,prvm)
fhstepDS 33
Arguments
y input - this is a data frame containing a single numeric column of response data.
des input - this is a data frame containing the numeric columns of the candidate
independent variables. The column names of des are of length 1 i.e., letters of
the alphabet. The m three-level factors always preceed the c two-level factors in
the design.
m input - this is an integer equal to the number of three-level factors in the design
c input - this is an integer equal to the number of two-level factors in the design.
Note m+c must be equal to the number of columns of des.
prvm input - this is a vector of text names of the terms in the model. This is created as
the value resulting from running ihstep or fhstep.
Value
returned vector of terms entered in the model at this step.
Author(s)
John Lawson
Description
This function performs a single step of a forward stepwise regression by entering an additional
2nd order term to a model already created by FitDefSc.R or fhstepDS.R This function is called by
FitDefSc.R
Usage
fhstepDS(y,des,m,c,prvm)
Arguments
y input - this is a data frame containing a single numeric column of response data.
des input - this is a data frame containing the numeric columns of the candidate
independent variables. The column names of des are of length 1 i.e., letters of
the alphabet. The m three-level factors always preceed the c two-level factors in
the design.
m input - this is an integer equal to the number of three-level factors in the design
c input - this is an integer equal to the number of two-level factors in the design.
Note m+c must be equal to the number of columns of des.
prvm input - this is a vector of text names of the terms in the model. This is created as
the value resulting from running ihstep or fhstep.
34 firstm
Value
Author(s)
John Lawson
firstm Find first term to enter forward stepwise regression that preserves
model hierarchy
Description
This function finds the first term to enter a hierarchical forward stepwise regression. If the term is
an interaction or quadratic term, the parent main effects are also included. This function is called
by ihstep.R
Usage
firstm(y,des)
Arguments
y input - this is a data frame containing a single numeric column of response data.
des input - this is a data frame containing the numeric columns of the candidate
independent variables. The column names of des are of length 1. The m three-
level factors always preceed the c two-level factors in the design.
Value
Author(s)
John Lawson
FitDefSc 35
Description
This function performs fits a model to a Definitive Screeing Design by first restricting main effects
to the smallest main effects and those significant at at least the .20 level in a main effects model.
Next forward stepwise selection is used to enter 2 factor interactions and quadratic effects.
Usage
FitDefSc(y,design,alpha=.05)
Arguments
y input - this is a vector containing a single numeric column of response data.
design input - this is a data frame containing the numeric columns of the candidate
independent variables created by the DefScreen function with only numerical
factors i.e. c=0. The factor names or colnames(design) should always be of
length 1 (for example letters of the alphabet "A", "B", etc.)
alpha input - alpha to enter in the forward stepwise regression with second order can-
didates should be between 0.05 and 0.20
Author(s)
John Lawson
fnextrm Find first term to enter forward stepwise regression that preserves
model hierarchy
Description
This function finds the first term to enter a hierarchical forward stepwise regression. If the term is
an interaction or quadratic term, the parent main effects are also included. This function is called
by ihstep.R
Usage
fnextrm(y,des,prvm)
36 fntrmDS
Arguments
y input - this is a data frame containing a single numeric column of response data.
des input - this is a data frame containing the numeric columns of the candidate
independent variables. The column names of des are of length 1. The m three-
level factors always preceed the c two-level factors in the design.
prvm input - this is a vector of text names of the terms in the model. This is created as
the value resulting from running ihstep or fhstep.
Value
returned vector of terms to be entered in the model at the next step.
Author(s)
John Lawson
fntrmDS Find first term to enter forward stepwise regression that preserves
model hierarchy
Description
This function finds the first term to enter a hierarchical forward stepwise regression. If the term is
an interaction or quadratic term, the parent main effects are also included. This function is called
by ihstep.R
Usage
fntrmDS(y,des,prvm)
Arguments
y input - this is a data frame containing a single numeric column of response data.
des input - this is a data frame containing the numeric columns of the candidate
independent variables. The column names of des are of length 1. The m three-
level factors always preceed the c two-level factors in the design.
prvm input - this is a vector of text names of the terms in the model. This is created as
the value resulting from running ihstep or fhstep.
Value
returned vector of terms to be entered in the model at the next step.
Author(s)
John Lawson
Fpower 37
Description
Calculates the power for the non-central F-distribution
Usage
Fpower(alpha, nu1, nu2, nc)
Arguments
alpha input - critical value alpha
nu1 input - degrees of freedom for numerator
nu2 input - degrees of freedom for denominator
nc input - noncentrality parameter
Value
probability of exceeding fcrit(alpha, nu1,nu2) with the non-central F-distribution with nu1 and nu2
degrees of freedom and noncentrality parameter nc
Author(s)
John Lawson
Description
Calculates the power for one-way ANOVA
Usage
Fpower1(alpha,nlev,nreps,Delta,sigma)
Arguments
alpha input - significance level of the F-test.
nlev input - the number of levels of the factor
nreps input - the number of replicates in each level of the factor.
Delta input - the size of a practical difference in two cell means.
sigma input - the standard deviation of the experimental error.
38 Fpower2
Value
probability of exceeding fcrit(alpha, nu1,nu2) with the non-central F-distribution with nu1 and nu2
degrees of freedom and noncentrality parameter nc
Author(s)
John Lawson
Description
Usage
Fpower2(alpha,nlev,nreps,Delta,sigma)
Arguments
Value
probability of exceeding fcrit(alpha, nu1,nu2) with the non-central F-distribution with nu1 and nu2
degrees of freedom and noncentrality parameter nc
Author(s)
John Lawson
fullnormal 39
fullnormal This function makes a full normal plot of the elements of the vector
called effects
Description
This function makes a full normal plot of the elements of the vector called effects
Usage
fullnormal(effects, labs, alpha = 0.05, refline = "TRUE")
Arguments
effects input - vector of effects to be plotted
labs input - vector of labels of the effects to be plotted
alpha input - alpha level for labeling of significant effects using Lenth statistic
refline input - logical variable that indicates whether a reference line is added to the
plot (default is "TRUE")
Author(s)
John Lawson
Description
Data from the Gauge R&R Study in chapter 5 of Design and Analysis of Experiments with R
Usage
data(gagerr)
Format
A data frame with 60 observations on the following 3 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(gagerr)
gapstat This function computes the gap statistic which is used to test for an
outlier using Daniels method
Description
This function computes the gap statistic which is used to test for an outlier using Daniels method
Usage
gapstat(beta, pse)
Arguments
beta input - vector of coefficients from saturated model fit to the data
pse input - Lenth’s PSE statistic calculated from the elements of beta
Value
returned gap statistic
Author(s)
John Lawson
Description
This function uses Daniel’s Method to find an outlier in an unreplicated 2(k−p) design.
Usage
Gaptest(DesY)
gear 41
Arguments
DesY input this is a data frame containing an unreplicated 2(k−p) design. The last
variable in the data frame should be the numeric response.
Author(s)
John Lawson
References
Box, G.E.P. (1991) "George’s column: Finding bad values in factorial designs", Quality Engineer-
ing, 3, 249-254.
Description
Data from the unreplicated split-plot fractional-factorial experiment on geometric distortion of drive
gears in chapter 8 of Design and Analysis of Experiments with R
Usage
data(gear)
Format
A data frame with 16 observations on the following 6 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(gear)
42 hardwood
halfnorm This function makes a half normal plot of the elements of the vector
called effects
Description
This function makes a half normal plot of the elements of the vector called effects
Usage
halfnorm(effects, labs, alpha = 0.05, refline = "TRUE")
Arguments
effects input - vector of effects to be plotted
labs input - vector of labels of the effects to be plotted
alpha input - alpha level for labeling of significant effects using Lenth statistic
refline input - logical variable that indicates whether a reference line is added to the
plot (default is "TRUE")
Author(s)
John Lawson
Description
Data from the low grade hardwood conjoint study in chapter 6 of Design and Analysis of Experi-
ments with R
Usage
data(hardwood)
Format
A data frame with 12 observations on the following 5 variables.
Design a factor with levels "RC" "AC" "OCI" "OCII"
Price a numeric variable
Density a factor with levels "Clear" "Heavy" "Medium"
Guarantee a factor with levels "1y" "Un"
Rating a numeric vector
HierAFS 43
Source
Examples
data(hardwood)
Description
Usage
HierAFS(y,x,m,c,step)
Arguments
Value
returned data frame the first column is a factor variable containing the formula for the model fit at
each step, the second numeric column is the R-square statistic for the model fit with each formula.
Author(s)
ihstep First step in a forward stepwise regression that preserves model hier-
archy
Description
This function performs the first step of a hierarchical forward stepwise regression. If an interaction
or quadratic term is entered first, the parent main effects are also entered into the model. This
function is called by HierAFS.R
Usage
ihstep(y,des,m,c)
Arguments
y input - this is a data frame containing a single numeric column of response data.
des input - this is a data frame containing the numeric columns of the candidate
independent variables. The column names of des are of length 1 i.e., letters of
the alphabet. The m three-level factors always preceed the c two-level factors in
the design.
m input - this is an integer equal to the number of three level factors in the design
c input - this is an integer equal to the number of two level factors in the design.
Note m+c must be equal to the number of columns of des.
Value
returned vector of terms entered in the model at this step.
Author(s)
John Lawson
Description
Data from the single array for injection molding experiment in chapter 12 of Design and Analysis
of Experiments with R
Usage
data(inject)
interleave 45
Format
A data frame with 20 observations on the following 8 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
shrinkage a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(inject)
Description
interleaves two vectors
Usage
interleave(v1,v2)
Arguments
v1 input - first vector
v2 input - second vector
Value
vector
46 LenthPlot
Description
Plot of the factor effects with significance levels based on robust estimation of contrast standard
errors.
Usage
LenthPlot(obj, alpha = 0.05, plt = TRUE, limits = TRUE,
xlab = "factors", ylab = "effects", faclab = NULL, cex.fac = graphics::par("cex.lab"),
cex.axis=graphics::par("cex.axis"), adj = 1, ...)
Arguments
obj object of class lm or vector with the factor effects.
alpha numeric. Significance level used for the margin of error (ME) and simultaneous
margin of error (SME). See Lenth(1989).
plt logical. If TRUE, a spikes plot with the factor effects is displayed. Otherwise, no
plot is produced.
limits logical. If TRUE ME and SME limits are displayed and labeled.
xlab character string. Used to label the x-axis. "factors" as default.
ylab character string. Used to label the y-axis. "effects" as default.
faclab list with components idx (numeric vector) and lab (character vector). The idx
entries of effects vector (taken from obj) are labelled as lab. The rest of the
effect names are blanked. If NULL all factors are labelled using the coefficients’
name.
cex.fac numeric. Character size used for the factor labels.
cex.axis numeric. Character size used for the axis.
adj numeric between 0 and 1. Determines where to place the "ME" (margin of
error) and the "SME" (simultaneous margin of error) labels (character size of
0.9*cex.axis). 0 for extreme left hand side, 1 for extreme right hand side.
... extra parameters passed to plot.
Details
If obj is of class lm, 2*coef(obj) is used as factor effect with the intercept term removed. Other-
wise, obj should be a vector with the factor effects. Robust estimate of the contrasts standard error
is used to calculate marginal (ME) and simultaneous margin of error (SME) for the provided signif-
icance (1 - alpha) level. See Lenth(1989). Spikes are used to display the factor effects. If faclab
is NULL, factors are labelled with the effects or coefficient names. Otherwise, those faclab\$idx
factors are labelled as faclab\$lab. The rest of the factors are blanked.
LGB 47
Value
The function is called mainly for its side effect. It returns a vector with the value of alpha used, the
estimated PSE, ME and SME.
Author(s)
Ernesto Barrios. Extension provided by Kjetil Kjernsmo (2013).
References
Lenth, R. V. (1989). "Quick and Easy Analysis of Unreplicated Factorials". Technometrics Vol. 31,
No. 4. pp. 469–473.
LGB This function uses the LGB Method to detect significant effects in un-
replicated fractional factorials.
Description
This function uses the LGB Method to detect significant effects in unreplicated fractional factorials.
Usage
LGB(Beta, alpha = 0.05, rpt = TRUE, plt = TRUE, pltl = TRUE)
Arguments
Beta input - this is the numeric vector of effects or coefficients to be tested
alpha input - This is the significance level of the test
rpt input - this is a logical variable that controls whether the report is written (default
is TRUE)
plt input - this is a logical variable that controls whether a half-normal plot is made
(default is TRUE)
pltl input - this is a logical variable that controls whether the significance limit line
is drawn on the half-normal plot (default is TRUE)
Author(s)
John Lawson
References
Lawson, J., Grimshaw, S., Burt, J. (1998) "A quantitative method for identifying active contrasts
in unreplicated factorial experiments based on the half-normal plot", Computational Statistics and
Data Analysis, 26, 425-436.
48 mod
LGBc This function does the calculations for the LGB Method to detect sig-
nificant effects in unreplicated fractional factorials.
Description
This function uses the LGB Method to detect significant effects in unreplicated fractional factorials.
Usage
LGBc(Beta, alpha = 0.05, rpt = TRUE, plt = TRUE, pltl = TRUE)
Arguments
Beta input - this is the numeric vector of effects or coefficients to be tested
alpha input - This is the significance level of the test
rpt input - this is a logical variable that controls whether the report is written (default
is TRUE)
plt input - this is a logical variable that controls whether a half-normal plot is made
(default is TRUE)
pltl input - this is a logical variable that controls whether the significance limit line
is drawn on the half-normal plot (default is TRUE)
Author(s)
John Lawson
References
Lawson, J., Grimshaw, S., Burt, J. (1998) "A quantitative method for identifying active contrasts
in unreplicated factorial experiments based on the half-normal plot", Computational Statistics and
Data Analysis, 26, 425-436.
Description
Gets mod of a to base b
Usage
mod(a,b)
ModelRobust 49
Arguments
a input- an integer
b input - an integer
Value
remainder of a/b or mod(a,b)
Author(s)
John Lawson
Description
Recalls Li and Nachtsheim’s model robust factorial designs from a catalog of data frames
Usage
ModelRobust(des, randomize=FALSE)
Arguments
des input- a character variable containing the name of a design in the catalog. If left
blank, the function prints a table showing all the design names in the catalog
randomize input- a logical
Value
design
Author(s)
John Lawson
References
Li, W. and Nachtsheim, C. J. (2000) "Model Robust factorial Designs", Technometrics, Vol 42, No.
4, pp345-352, 2000.
50 Naph
Description
Data from the mixture process variable experiment with mayonnaise in chapter 11 of Design and
Analysis of Experiments with R
Usage
data(MPV)
Format
A data frame with 35 observations on the following 4 variables.
x1 a numeric vector
x2 a numeric vector
x3 a numeric vector
z1 a numeric vector
z2 a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(MPV)
Description
Data from the Yields of naphthalene black of Chapter 5 in Design and Analysis of Experiments
with R
Usage
data(Naph)
OptPB 51
Format
A data frame with 30 observations on the following 2 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(Naph)
Description
Selects the columns from a Plackett-Burman Design Produced by FrF2 that will minimize model
dependence for main effects and two factor interactions and returns the design in a data frame
Usage
OptPB(nruns, nfactors, randomize=FALSE)
Arguments
nruns input- an integer representing the number of runs in the design
nfactors input - in integer representing the number of factors in the design
randomize input - logical
Value
design
Author(s)
John Lawson
References
Fairchild, K. (2011) "Screening Designs that Minimize Model Dependence", MS Project Depart-
ment of Statistics Brigham Young University, Dec. 2011.
52 PBDes
Description
Data from the Blocked response surface design for pastry dough experiment in chapter 10 of Design
and Analysis of Experiments with R
Usage
data(pastry)
Format
A data frame with 28 observations on the following 5 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(pastry)
Description
Creates a 12, 20, or 24 run Plackett-Burman design in a data frame with numeric factor levels by
cyclically rotating the factor leves in the first row
Usage
PBDes(nruns, nfactors, randomize=FALSE)
pest 53
Arguments
nruns input- an integer representing the number of runs in the design
nfactors input - in integer representing the number of factors in the design
randomize input - logical
Value
design
Author(s)
John Lawson
References
Lawson, J. (2015) "Design and Analysis of Experiments with R page 229",CRC Press, Boca Raton,
2015.
Description
Data from the Pesticide formulation experiment in chapter 11 of Design and Analysis of Experi-
ments with R
Usage
data(pest)
Format
A data frame with 13 observations on the following 4 variables.
x1 a numeric vector
x2 a numeric vector
x3 a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(pest)
54 plasma
Description
Data from the pesticide application experiment in chapter 5 of Design and Analysis of Experiments
with R
Usage
data(pesticide)
Format
Source
Examples
data(pesticide)
Description
Data from the unreplicated split-plot 25 experiment on plasma treatment of paper in chapter 8 of
Design and Analysis of Experiments with R
Usage
data(plasma)
polvdat 55
Format
A data frame with 32 observations on the following 6 variables.
A a factor with levels -1 1
B a factor with levels -1 1
C a factor with levels -1 1
D a factor with levels -1 1
E a factor with levels -1 1
y a numeric vector containing the response
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(plasma)
Description
Data from the Polvoron mixture experiment in chapter 11 of Design and Analysis of Experiments
with R
Usage
data(polvdat)
Format
A data frame with 12 observations on the following 4 variables.
x1 a numeric vector
x2 a numeric vector
x3 a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(polvdat)
56 prodstd
Description
Data from the polymerization strength variability study in chapter 5 of Design and Analysis of
Experiments with R
Usage
data(polymer)
Format
A data frame with 120 observations on the following 5 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(polymer)
prodstd Complete control factor array and noise factor array for connector
experiment
Description
Data from the complete control factor array and noise factor array for connector experiment in
chapter 12 of Design and Analysis of Experiments with R
Usage
data(prodstd)
qsar 57
Format
A data frame with 16 observations on the following 16 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
Pof a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(prodstd)
Description
Data from the Library of substituted hydroxyphenylurea compounds in chapter 10 of Design and
Analysis of Experiments with R (compact format)
Usage
data(qsar)
Format
A data frame with 36 observations on the following 4 variables.
Compound a numeric vector
HE a numeric vector
DMz a numeric vector
S0K a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(qsar)
58 rcb
Description
Data from the cattle rations design experiment in chapter 10 of Design and Analysis of Experiments
with R
Usage
data(Rations)
Format
Source
Examples
data(Rations)
Description
Data from the generalized RCB golf driving experiment in chapter 4 of Design and Analysis of
Experiments with R
Usage
data(rcb)
residue 59
Format
A data frame with 135 observations on the following 3 variables.
id a factor with levels 1 2 3 4 5 6 7 8 9
teehgt a factor with levels 1 2 3
cdistance a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(rcb)
Description
Data from the Herbicide degradation experiment in chapter 9 of Design and Analysis of Experi-
ments with R
Usage
data(residue)
Format
A data frame with 16 observations on the following 3 variables.
soil a factor with levels "C" "P"
moisture a factor with levels "L" "H"
temp a factor with levels 10 30
X1 a numeric vector
X2 a numeric vector
X3 a numeric vector
X4 a numeric vector
X5 a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(residue)
60 sausage
Description
Data from the Rubber Elasticity Study in chapter 5 of Design and Analysis of Experiments with R
Usage
data(rubber)
Format
Source
Examples
data(rubber)
Description
Data from the Split-plot experiment on sausage casing with RCB in whole plot in chapter 7 of
Design and Analysis of Experiments with R
Usage
data(sausage)
Smotor 61
Format
A data frame with 32 observations on the following 5 variables.
Block a factor with levels 1 2
Gbatch a factor with levels 1 2 3 4
A a factor with levels -1 1
B a factor with levels -1 1
C a factor with levels -1 1
D a factor with levels -1 1
ys a numeric vector containing the response
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(sausage)
Description
Data from the single array for starting motor experiment in chapter 12 of Design and Analysis of
Experiments with R
Usage
data(Smotor)
Format
A data frame with 18 observations on the following 6 variables.
A a factor with levels 1 2
B a factor with levels 1 2 3
C a factor with levels 1 2 3
D a factor with levels 1 2 3
E a factor with levels 1 2
torque a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
62 soup
Examples
data(Smotor)
Description
Data from the dry mix soup experiment in chapter 6 of Design and Analysis of Experiments with R
Usage
data(soup)
Format
Source
Examples
data(soup)
soupmx 63
Description
Data from the dry soup mix variance component study of Chapter 5 in Design and Analysis of
Experiments with R
Usage
data(soupmx)
Format
A data frame with 12 observations on the following 2 variables.
batch a factor with levels 1 2 3 4
weight a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(soupmx)
Description
Data from the Split-plot cookie baking experiment in chapter 8 of Design and Analysis of Experi-
ments with R
Usage
data(splitPdes)
Format
A data frame with 24 observations on the following 5 variables.
short a factor with levels 100 80
trayT a factor with levels RoomT Hot
bakeT a factor with levels low mid high
batch a factor with levels 1 2
y a numeric vector
64 SPMPV
Source
Examples
data(splitPdes)
Description
Data from the Split-plot mixture process variable experiment with vinyl in chapter 10 of Design and
Analysis of Experiments with R
Usage
data(SPMPV)
Format
Source
Examples
data(SPMPV)
stdord 65
Description
Makes standard order
Usage
stdord(m)
Arguments
m input - vector length
Value
vector in standard order
Description
Data from the Repeated measures study with dairy cow diets in chapter 9 of Design and Analysis
of Experiments with R (strung out format)
Usage
data(strung)
Format
A data frame with 120 observations on the following 5 variables.
Diet a factor with levels "Barley" "Mixed" "Lupins"
Cow a factor with levels 1 2 3 4 5 6 7 8 9 10
week a factor with levels 1 2 3 4
protein a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(strung)
66 strungtile
strungtile Strung out control factor array and raw response data for Ina tile ex-
periment
Description
Data from the strung out control factor array and raw response data for Ina tile experiment in chapter
12 of Design and Analysis of Experiments with R
Usage
data(strungtile)
Format
A data frame with 16 observations on the following 16 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
H a numeric vector
AH a numeric vector
BH a numeric vector
CH a numeric vector
DH a numeric vector
EH a numeric vector
FH a numeric vector
GH a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(strungtile)
sugarbeet 67
Description
Sugarbeet data from chapter 2 of Design and Analysis of Experiments with R
Usage
data(sugarbeet)
Format
A data frame with 18 observations on the following 2 variables.
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(sugarbeet)
Description
Data from the taste test panel experiment in Chapter 7 of Design and Analysis of Experiments with
R
Usage
data(taste)
Format
A data frame with 24 observations on the following 3 variables.
Source
Examples
data(taste)
Description
Data from the teaching experiment in chapter 2 of Design and Analysis of Experiments with R
Usage
data(teach)
Format
Source
Examples
data(teach)
Tet 69
Description
Data from the Tetracycline concentration in plasma study in chapter 10 of Design and Analysis of
Experiments with R (compact format)
Usage
data(Tet)
Format
Source
Examples
data(Tet)
tile Control factor array and summary statistics for Ina tile experiment
Description
Data from the control factor array and summary statistics for Ina tile experiment in chapter 12 of
Design and Analysis of Experiments with R
Usage
data(tile)
70 Treb
Format
A data frame with 8 observations on the following 11 variables.
A a numeric vector
B a numeric vector
C a numeric vector
D a numeric vector
E a numeric vector
F a numeric vector
G a numeric vector
y1 a numeric vector
y2 a numeric vector
ybar a numeric vector
lns2 a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(tile)
Description
Data from the Box-Behnken design for trebuchet experiment in chapter 10 of Design and Analysis
of Experiments with R
Usage
data(Treb)
Format
A data frame with 15 observations on the following 4 variables.
x1 a numeric vector
x2 a numeric vector
x3 a numeric vector
y a numeric vector
Tukey1df 71
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(Treb)
Tukey1df This function performs Tukey’s single degree of freedom test for inter-
action in an unreplicated two-factor design
Description
This function performs Tukey’s single degree of freedom test for interaction in an unreplicated
two-factor design
Usage
Tukey1df(data)
Arguments
data input - this is a data frame with three variables, the first variable is a numeric
response and next two variables are factors. There should be ab lines in the data
frame where a is the number of levels of the first factor, and b is the number of
levels of the second factor.
Author(s)
John Lawson
Description
function for getting confidence intervals on variance components estimated by the method of mo-
ments
Usage
vci(confl,c1,ms1,nu1,c2,ms2,nu2)
72 vinyl
Arguments
confl input- confidence level
c1 input - linear combination coefficient of ms1 in the estimated variance compo-
nent
ms1 input - Anova mean square 1
nu1 input - Anova degrees of freedom for mean square 1
c2 input - linear combination coefficient of ms2 in the estimated variance compo-
nent
ms2 input - Anova mean square 2
nu2 input - Anova degrees of freedom for mean square 2
Value
returned delta, Lower and Upper limits
Author(s)
John Lawson
Description
Data from vinyl plasticiser formulation experiment in chapter 11 of Design and Analysis of Exper-
iments with R
Usage
data(vinyl)
Format
A data frame with 40 observations on the following 7 variables.
WP a numeric vector
x1 a numeric vector
x2 a numeric vector
x3 a numeric vector
z1 a numeric vector
z2 a numeric vector
y a numeric vector
virus 73
Source
Examples
data(vinyl)
Description
Data from the Assay of Viral Contamination experiment in chapter 3 of Design and Analysis of
Experiments with R
Usage
data(virus)
Format
y a numeric vector
Sample a factor with levels 1 2 3 4 5 6
Dilution a factor with levels 3 4 5
Source
Examples
data(virus)
74 web
Description
Data from the Volt meter experiment in chapter 3 of Design and Analysis of Experiments with R
Usage
data(volt)
Format
y a numeric vector
A a factor containing the levels (22, 32) of factor A
B a factor containing the levels (0.5, 5.0) of factor B
C a factor containing the levels (0.5, 5.0) of factor C
Source
Examples
data(volt)
Description
Data from the web page design experiment in chapter 3 of Design and Analysis of Experiments
with R
Usage
data(web)
WeldS 75
Format
A data frame with 36 observations on the following 6 variables.
A a factor with levels 1 2
B a factor with levels 1 2
C a factor with levels 1 2
D a factor with levels 1 2
visitors a numeric vector
signup a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(web)
Description
Data from the Single Array Experiment in Exercise 5 of Chapter 12 in Design and Analysis of
Experiments with R. The factors are in coded levels.
Usage
data(WeldS)
Format
A data frame with 16 observations on the following 16 variables.
D a numeric vector
H a numeric vector
G a numeric vector
A a numeric vector
F a numeric vector
GH a numeric vector
C a numeric vector
B a numeric vector
J a numeric vector
E a numeric vector
76 WeldS
AC a numeric vector
AH a numeric vector
AG a numeric vector
e1 a numeric vector
e2 a numeric vector
y a numeric vector
Source
Design and Analysis of Experiments with R, by John Lawson, CRC/Chapman Hall
Examples
data(WeldS)
Index
∗ datagen cont, 20
Altscreen, 4 cpipe, 21
BIBsize, 11 culture, 22
DefScreen, 23 dairy, 23
EEw1s1, 25 drug, 24
EEw1s2, 25 eptaxr, 29
EEw1s3, 26 eptaxs2, 30
EEw2s1, 27 eptaxyb, 31
EEw2s2, 27 gagerr, 39
EEw2s3, 28 gear, 41
EEw3, 29 hardwood, 42
Fcrit, 32 inject, 44
Fpower1, 37 MPV, 50
Fpower2, 38 Naph, 50
mod, 48 pastry, 52
ModelRobust, 49 pest, 53
OptPB, 51 pesticide, 54
PBDes, 52 plasma, 54
∗ datasets polvdat, 55
antifungal, 5 polymer, 56
Apo, 6 prodstd, 56
apple, 6 qsar, 57
arso, 7 Rations, 58
augm, 8 rcb, 58
Bdish, 9 residue, 59
Bff, 9 rubber, 60
bha, 10 sausage, 60
bioequiv, 11 Smotor, 61
bioeqv, 12 soup, 62
blood, 13 soupmx, 63
BoxM, 13 splitPdes, 63
BPmonitor, 14 SPMPV, 64
bread, 15 strung, 65
cakeb, 15 strungtile, 66
cement, 16 sugarbeet, 67
chem, 17 taste, 67
chipman, 17 teach, 68
COdata, 18 Tet, 69
connector, 20 tile, 69
77
78 INDEX
Treb, 70 cement, 16
vinyl, 72 chem, 17
virus, 73 chipman, 17
volt, 74 COdata, 18
web, 74 colormap, 19
WeldS, 75 connector, 20
∗ design cont, 20
LenthPlot, 46 cpipe, 21
∗ hplot culture, 22
colormap, 19
fullnormal, 39 daewr (daewr-package), 4
halfnorm, 42 daewr-package, 4
∗ htest dairy, 23
fhstep, 32 DefScreen, 23
fhstepDS, 33 drug, 24
firstm, 34
EEw1s1, 25
FitDefSc, 35
EEw1s2, 25
fnextrm, 35
EEw1s3, 26
fntrmDS, 36
EEw2s1, 27
gapstat, 40
EEw2s2, 27
Gaptest, 40
EEw2s3, 28
HierAFS, 43
EEw3, 29
ihstep, 44
eptaxr, 29
LGB, 47
eptaxs2, 30
LGBc, 48
eptaxyb, 31
Tukey1df, 71
vci, 71 Fcrit, 32
∗ package fhstep, 32
daewr-package, 4 fhstepDS, 33
firstm, 34
Altscreen, 4
FitDefSc, 35
antifungal, 5
fnextrm, 35
Apo, 6 fntrmDS, 36
apple, 6 Fpower, 37
arso, 7 Fpower1, 37
augm, 8 Fpower2, 38
fullnormal, 39
Bdish, 9
Bff, 9 gagerr, 39
bha, 10 gapstat, 40
BIBsize, 11 Gaptest, 40
bioequiv, 11 gear, 41
bioeqv, 12
blood, 13 halfnorm, 42
BoxM, 13 hardwood, 42
BPmonitor, 14 HierAFS, 43
bread, 15
ihstep, 44
cakeb, 15 inject, 44
INDEX 79
interleave, 45 virus, 73
volt, 74
LenthPlot, 46
LGB, 47 web, 74
LGBc, 48 WeldS, 75
mod, 48
ModelRobust, 49
MPV, 50
Naph, 50
OptPB, 51
pastry, 52
PBDes, 52
pest, 53
pesticide, 54
plasma, 54
polvdat, 55
polymer, 56
prodstd, 56
qsar, 57
Rations, 58
rcb, 58
residue, 59
rubber, 60
sausage, 60
Smotor, 61
soup, 62
soupmx, 63
splitPdes, 63
SPMPV, 64
stdord, 65
strung, 65
strungtile, 66
sugarbeet, 67
taste, 67
teach, 68
Tet, 69
tile, 69
Treb, 70
Tukey1df, 71
vci, 71
vinyl, 72