qualityTools (version 1.55)

starDesign: Axial Design

Description

starDesign is a function to create the star portion of a response surface design. The starDesign function can be used to create a star portion of a response surface design for a sequential assembly strategy. One can either specify k and p and alpha and cs and cc OR simply simply pass an object of class facDesign to the data. In the latter an object of class facDesign otherwise a list containing the axial runs and centerpoints is returned.

Usage

starDesign(k, p = 0, alpha = c("both", "rotatable", "orthogonal"), cs, cc, data)

Arguments

k

integer value giving number of factors.

p

integer vaue giving the number of factors via aliasing. By default set to ‘0’.

alpha

if no numeric value is given defaults to “both” i.e. “orthogonality” and “rotatibility” which can be set as character strings too.

cs

integer value giving the number of centerpoints in the star portion of the design.

cc

integer value giving the number of centerpoints in the cube portion of the design.

data

optional - an object of class facDesign.

Value

starDesign returns a facDesign object if an object of class facDesign is given or a list containing entries for axial runs and center points in the cube and the star portion of a design.

References

MYERS, R.H.; MONTGOMERY, D.C.; ANDERSON-COOK, C.M.: Response Surface Methodology. New Jersey: Wiley,2009.

See Also

facDesign for 2^k factorial designs fracDesign for 2^k-p fractional factorial designs rsmDesign for response surface designs mixDesign for mixture designs http://www.r-qualitytools.org

Examples

#Example 1 - sequential assembly
#factorial design with one center point in the cube portion
fdo = facDesign(k = 3, centerCube = 1)  
fdo

#set the response via generic response method
response(fdo) = 1:9 
 
#sequential assembly of a response surface design (rsd)
rsd = starDesign(data = fdo)  
rsd

#Example 2 - returning a list
starDesign(k = 3, cc = 2, cs = 2, alpha = "orthogonal")
starDesign(k = 3, cc = 2, cs = 2, alpha = "rotatable")
starDesign(k = 3, cc = 2, cs = 2, alpha = "both")