中心复合设计CCD

中心复合设计CCD

一、描述

这个功能可以灵活地选择复制，走样的预测和小数块，并选择的轴或中心点产生了广泛的中心复合设计。

二、用法

ccd(basis, generators, blocks = "Block", n0 = 4, alpha = "orthogonal", wbreps = 1, bbreps = 1, randomize = TRUE, inscribed = FALSE, coding)

三、参数

四、详细信息

中央复合设计的CCD响应面勘探中使用的流行的设计。他们被堵塞的设计组成的至少一种的多维数据集的块（两水平因子或部分因子，附加中心点），和至少一个星状块（点沿每个轴的位置-α+α-），加上中心点。一切都被假定为规模的编码的多维数据集的设计部分，其中为每个变量的值-1和1，中心点0。使用的编码参数，如果需要的话，编码可以提供。

基础参数确定一个基本的设计，用来创建立方体块。例如，basis = ~A+B+C会生成一个基本的8次试验的设计。如果你需要额外的变量，可使用generators，例如，generators = c(D~-A*B, E~B*C)加在上面定义关系会产生一个5因素设计I = ABD = BCE = ACDE 。

如果你想立方体点分为零碎块，给块参数公式。例如，basis =A+B+C+D+E

generators =F~-A*C*D, 和blocks =Day ~c(A*B*C,C*D*E)。阻挡变量将被命名

为“Day”，将其分为4块8和32运行的基本设计（半部分6因素）运行每一个的基础上，结合标志的A * B * C 和C * D * E 。请注意，这一天将混淆与发电机组，及其相互作用，以及所有这些别名：Day = ABC = CDE = ABDE = -BDF = -ADCF = -BCEF = -AEF。对于每个模块，我们将添加N0（N0 [1 ]）的中心点。

A星块包含N0（N0 [2]）的中心点，在+/- alpha在每个坐标轴加点。你可以指定你想要的字母，或选择以实现正交块（默认）或可旋转。可能的设计点值有限制的情况下，你可能要指定inscribed = TRUE。这将缩减本设计所以没有编码值超过1 。

wbreps值大于1的原因的立方体点或星点在同一个块被复制。bbreps值大于1，导致了额外的块被添加到设计中。默认情况下，设计是随机的，所以每块的运行顺序是随机的。块的出现顺序是随机的。

如果基础是一个整数，而不是一个公式，默认的环境变量为x1，x2 ， ...用于，例如，basis = 3是相当于basis = ~x1+x2+x3的。您可以指定使用相同的符号的或分数块。您可能对向的左手侧上的变量的基础是一个公式时，这些变量被添加到返回的设计。例如，Yield+ProcessTime ~ A+B+C会生成一个设计A，B，C ，以及在NA 中额外的变量因素，产量和初始处理时间。

五、值

数据框所产生的设计，或如果指定编码的coded.data对象。阻塞变量将是一个因素，所有其他变量将是数字。

六、例子

library(rsm)

### Generate a standard 3-variable CCD with 12 runs in each block

des = ccd(3, n0=c(4,6), coding = list(x1 ~ (Temp - 150)/10,

x2 ~ (Pres - 50)/5, x3 ~ Feedrate - 4))

decode.data(des)

### Same as above, except make the design rotatable,

### and inscribed so that no coded value exceeds 1

des2 = ccd(3, n0=c(4,6), alpha = "rotatable", inscribed = TRUE,

coding = list(x1 ~ (Temp - 150)/10,

x2 ~ (Pres - 50)/5, x3 ~ Feedrate - 4))

decode.data(des2)

### Generate a 5-variable design in 2 blocks. The cube block has 16 runs ### This design will have alpha = 2; it is both orthogonal and rotatable ccd(~ x1 + x2 + x3 + x4, x5 ~ x1 * x2 * x3 * x4, n0 = c(6,1))

### Generate a 5-variable design with 5 blocks:

### 4 blocks of 8 cube points each, and 1 block with star points

### You'll get alpha = 2.366; if you add alpha = "rot", you'll get 2.378 des = ccd(y1 + y2 ~ A + B + C + D + E, , Shift ~ c(-A*B*C, C*D*E), n0=c(2,4))