rsmDesign: Generate a response surface design (i.e. central composite...
In qualityTools: Statistical Methods for Quality Science

Description Usage Arguments Details Value Note Author(s) See Also Examples

Generates a response surface design containing a cube, centerCube, star and centerStar portion.

1 2	rsmDesign(k = 3, p = 0, alpha = "rotatable", blocks = 1, cc = 1, cs = 1, fp = 1, sp = 1, faceCentered = FALSE)

`k`	integer value giving the number of factors. By default `k` has the value ‘3’.
`p`	integer value giving the number of additional factors in the response surface design by aliasing effects. By default `p` has the value ‘0’.
`alpha`	character string - alpha shoud be “rotatable” (default), “orthogonal” or “both”. If “both” values for cc and cs will be DISCARDED.
`blocks`	integer value giving the number of blocks in the response surface design. By default `blocks` has the value ‘1’.
`cc`	integer value giving the number of centerpoints (per block) in the cube portion (i.e. the factorial 2^k design) of the response surface design. Replications of centerpoints in the cube portion can be set with this.By default `cc` has the value ‘1’.
`cs`	integer value giving the number of centerpoints in the star portion (alpha) of the response surface design. Replications of centerpoints in the star portion can be set with this. By default `cs` has the value ‘1’.
`fp`	integer value giving the number of replications per factorial point (i.e. corner points e.g. (-1,1)). By default `fp` has the value ‘1’.
`sp`	integer value giving the number of replications per star point (i.e. alpha). By default `sp` has the value ‘1’.
`faceCentered`	logical value wheter to use a faceCentered response surface design or not (i.e. alpha = 1). By default `faceCentered` has the value ‘FALSE’.

Generated designs consist of a cube, centerCube, star and a centerStar portion. The replication structure can be set with the parameters cc (centerCube), cs (centerStar), fp (factorialPoints) and sp (starPoints).

rsmDesign returns an object of class facDesign.

For an example in context which shows the usage of the function rsmDesign() please read the vignette for the package qualityTools at http://www.r-qualitytools.org/html/Improve.html.

Thomas Roth thomas.roth@tu-berlin.de

facDesign
fracDesign
fracChoose
pbDesign
rsmChoose
http://www.r-qualitytools.org/html/Improve.html

#central composite design for 2 factors with 2 blocks, alpha = 1.41, 
#5 centerpoints in the cube portion and 3 centerpoints in the star portion:
rsmDesign(k = 2, blocks = 2, alpha = sqrt(2),cc = 5, cs = 3)

#central composite design with both, orthogonality and near rotatability 
rsmDesign(k = 2, blocks = 2, alpha = "both")

#central composite design with
#2 centerpoints in the factorial portion of the design i.e 2
#1 centerpoint int the star portion of the design i.e. 1
#2 replications per factorial point i.e. 2^3*2 = 16
#3 replications per star points 3*2*3 = 18
#makes a total of 37 factor combinations
rsdo = rsmDesign(k = 3, blocks = 1, alpha = 2, cc = 2, cs = 1, fp = 2, sp = 3)
nrow(rsdo) #37

Loading required package: Rsolnp
Loading required package: MASS

Attaching package: 'qualityTools'

The following object is masked from 'package:stats':

    sigma

   StandOrd RunOrder Block      A      B  y
4         4        1     1  1.000  1.000 NA
3         3        2     1 -1.000  1.000 NA
2         2        3     1  1.000 -1.000 NA
5         5        4     1  0.000  0.000 NA
9         9        5     1  0.000  0.000 NA
6         6        6     1  0.000  0.000 NA
7         7        7     1  0.000  0.000 NA
8         8        8     1  0.000  0.000 NA
1         1        9     1 -1.000 -1.000 NA
13       13       10     2  0.000  1.414 NA
11       11       11     2  1.414  0.000 NA
16       16       12     2  0.000  0.000 NA
14       14       13     2  0.000  0.000 NA
10       10       14     2 -1.414  0.000 NA
15       15       15     2  0.000  0.000 NA
12       12       16     2  0.000 -1.414 NA
Warning messages:
1: In `[<-`(`*tmp*`, i, value = <S4 object of class "doeFactor">) :
  implicit list embedding of S4 objects is deprecated
2: In `[<-`(`*tmp*`, i, value = <S4 object of class "doeFactor">) :
  implicit list embedding of S4 objects is deprecated
   StandOrd RunOrder Block      A      B  y
3         3        1     1 -1.000  1.000 NA
7         7        2     1  0.000  0.000 NA
2         2        3     1  1.000 -1.000 NA
6         6        4     1  0.000  0.000 NA
4         4        5     1  1.000  1.000 NA
5         5        6     1  0.000  0.000 NA
1         1        7     1 -1.000 -1.000 NA
13       13        8     2  0.000  0.000 NA
9         9        9     2  1.414  0.000 NA
12       12       10     2  0.000  0.000 NA
8         8       11     2 -1.414  0.000 NA
14       14       12     2  0.000  0.000 NA
11       11       13     2  0.000  1.414 NA
10       10       14     2  0.000 -1.414 NA
Warning messages:
1: In `[<-`(`*tmp*`, i, value = new("doeFactor")) :
  implicit list embedding of S4 objects is deprecated
2: In `[<-`(`*tmp*`, i, value = new("doeFactor")) :
  implicit list embedding of S4 objects is deprecated
Warning messages:
1: In `[<-`(`*tmp*`, i, value = new("doeFactor")) :
  implicit list embedding of S4 objects is deprecated
2: In `[<-`(`*tmp*`, i, value = new("doeFactor")) :
  implicit list embedding of S4 objects is deprecated
3: In `[<-`(`*tmp*`, i, value = new("doeFactor")) :
  implicit list embedding of S4 objects is deprecated
[1] 37