rsmDesign: Generate a response surface design (i.e. central composite...

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

View source: R/rsm_n.r

Description

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

Usage

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

Arguments

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’.

Details

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).

Value

rsmDesign returns an object of class facDesign.

Note

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.

Author(s)

Thomas Roth [email protected]

See Also

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
#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

Example output

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

qualityTools documentation built on May 30, 2017, 1:43 a.m.