中心复合设计CCD

1、中心复合设计CCD一、描述这个功能可以灵活地选择复制，走样的预测和小数块，并选择的轴或中心点产生了广泛 的中心复合设计。二、用法ccd（basis, generators, blocks = Block, n0 = 4, alpha = orthogonal,wbreps = 1, bbreps = 1, randomize= TRUE, inscribed = FALSE, coding）三、参数basis表示公式，或一个给变量的数目的整数。如果这个公式有一个左手侧，即名 为变量所附的设计，并初始化为NA。generators表示可选的公式或公式列表生成别名变量blocks表示一个字符串或公

2、式。如果一个字符串，它是名称的blcoking因素，如果 公式中，左手侧是用来作为名称的blcoking因素，和式（次）的右手侧上的 用于生成单独的零碎块。n0表示整数或载体给中心点的数量为两个整数。如果只有一个值，这是在每个 块的数目中心pouints的。如果给定两个值，N01指定的数量中心点立方体 块，N02指定中心点星状块数。alpha表示如果是数字在星状点的位置，也可能是 个字符串，它是否符 合“正交”，星点定位设计正交阻止它是否符合“可旋转”，星点 的选择，使设计有可旋转性。默认值是产生正交阻塞的设计。需要 注意的是面心的设计是通过以下方式获得指定a =1。wbreps表示括号（）块

3、内复制。如果这是一个长度为2的向量，然后使用单独的数 字“多维数据集”和“轴块。bbreps表示括号（）之间的块的复制（即，每个块的重复数）f这是一个长度为2 的向量，然后使用单独的数字“多维数据集”和“轴块。randomize表示布尔值确定是否随机化设计。每个块是独立随机化。inscribed表示布尔值。如果设置为FALSE，多维数据集的点是在+ / - 1轴点是至少为 1。如果TRUE，整个DESGN的被缩放下来，该轴点是在+/ - 1和立方体点 是在空间内的位置.coding表示可选的公式列表。如果设置此参数，表示一个coded.data对象返回与 指定的编码。四、详细信息中央复合设计的

4、CCD响应面勘探中使用的流行的设计。他们被堵塞的设计组成的 至少一种的多维数据集的块（两水平因子或部分因子，附加中心点），和至少一个星状块 （点沿每个轴的位置-a +a -），加上中心点。一切都被假定为规模的编码的多维数据集 的设计部分，其中为每个变量的值-1和1，中心点0。使用的编码参数，如果需要的话， 编码可以提供。基础参数确定一个基本的设计，用来创建立方体块。例如，basis =A+B+C会生 成一个基本的8次试验的设计。如果你需要额外的变量，可使用generators，例如， generators = c（D-A*B, EB*C）加在上面定义关系会产生一个5因素设计I = ABD =B

5、CE = 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（ N

6、0 2）的中心点，在+/- alpha在每个坐标轴加点。你可以指 定你想要的字母，或选择以实现正交块（默认）或可旋转。可能的设计点值有限制的情 况下，你可能要指定inscribed = TRUE。这将缩减本设计所以没有编码值超过1。 wbreps值大于1的原因的立方体点或星点在同一个块被复制。bbreps值大于1，导致 了额外的块被添加到设计中。默认情况下，设计是随机的，所以每块的运行顺序是随机 的。块的出现顺序是随机的。如果基础是一个整数，而不是一个公式，默认的环境变量为x1, x2，.用于，例 如，basis = 3是相当于basis = x1+x2+x3的。您可以指定使用相同的符号的或

7、分数块。您可能对向的左手侧上的变量的基础是一个公式时，这些变量被添加到返回的 设计。例如，Yield+ProcessTime A+B+C会生成一个设计A，B，C，以及在NA 中额外的变量因素，产量和初始处理时间。五、值数据框所产生的设计，或如果指定编码的coded.data对象。阻塞变量将是一个因素，所 有其他变量将是数字。六、例子library(rsm)# Generate a standard 3-variable CCD with 12 runs in each blockdes = ccd(3, n0=c(4,6), coding = list(x1 (Temp - 150)/10,

8、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 1des2 = 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 rotatableccd( 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 poin