cordexch

协调交换

语法

dCE = cordexch (nfactors nruns) (dCE, X) = cordexch (nfactors nruns) (dCE, X) = cordexch (nfactors, nruns。”模型”) (dCE, X) = cordexch (……”的名字”,价值)

描述

dCE = cordexch (nfactors nruns)使用一个coordinate-exchange算法来生成一个D最优设计dCE与nruns(行运行dCE)线性相加模型nfactors因素(的列dCE)。模型包括一个常数项。

(dCE, X) = cordexch (nfactors nruns)还返回相关的设计矩阵X列的模型项评估在每个治疗(行)dCE。

(dCE, X) = cordexch (nfactors, nruns。”模型”)使用线性回归模型中指定模型。模型是下列之一:

“线性”——常数和线性。这是默认的。
“互动”常数,线性和交互
“二次”常数,线性、交互和平方项
“purequadratic”常数,线性,平方项

列的顺序X对于一个完整的二次模型n条件是:

常数项
线性条件为1、2、…,n
交互条件(1、2),(1,3),…(1,n),(2,3),…,(n- 1,n)
平方项的1、2、…n

其他模型使用这些术语的一个子集,在相同的顺序。

另外,模型可以是一个多项式矩阵指定任意的顺序。在这种情况下,模型应该为每个因子和一行一列模型中的每一项。在任何的行条目模型是权力因素的列。例如,如果一个模型的因素X1,X2,X3,然后一行[0 1 2]在模型指定这个词(X1。^ 0) * (X2 ^ 1)。* (X3。^ 2)。连续的0模型指定一个常数项,可以省略。

(dCE, X) = cordexch (……”的名字”,价值)指定一个或多个可选的名称/值对的设计。有效的参数及其值在下表中列出。指定的名字在单引号。

的名字	价值
`界限`	上下边界为每一个因素,作为一个指定`2`——- - - - - -`nfactors`矩阵。另外,这个值可以是一个包含细胞数组`nfactors`元素,每个元素指定容许值的向量对应的因素。
`分类`	分类预测的指标。
`显示`	要么`“上”`或`“关闭”`来控制迭代计数器的显示。默认值是`“上”`。
`excludefun`	处理函数,排除了不受欢迎的。如果函数是f,它必须支持的语法金宝appb=f(年代),年代是一个矩阵的治疗`nfactors`列和b是一个向量的布尔值的行数一样吗年代。b(我如果方法应该排除)是正确的我th行年代。
`初始化`	初步设计的`nruns`——- - - - - -`nfactors`矩阵。默认是一个随机选择的点集。
`水平`	向量的数量为每个因素水平。不使用时`界限`被指定为一个细胞数组。
`麦克斯特`	最大迭代次数。默认值是`10`。
`试`	的次数来生成一个设计从一个新的起点。算法使用随机点为每个试,除了可能是第一个。默认值是`1`。
`选项`	结构,指定是否要并行运行,并指定随机流或流。并行计算需要并行计算工具箱™。创建`选项`结构`statset`。结构方面: `UseParallel`——设置为`真正的`并行计算。默认是`假`。 `UseSubstreams`——设置为`真正的`以可再生的方式并行计算。默认是`假`。计算重复性良好`流`一种允许substreams:`“mlfg6331_64”`或`“mrg32k3a”`。 `流`——一个`RandStream`这样的对象的对象或单元阵列。如果你不指定`流`,`cordexch`使用默认的河流或溪流。如果你选择指定`流`,除了在使用一个对象 `UseParallel`是`真正的` `UseSubstreams`是`假` 在这种情况下,使用一个单元阵列大小相同的平行池。

例子

假设你想要一个设计参数估计在接下来的三因子,seven-term交互模型:

$y = β_{0} + β_{1} x {}_{1}+ β_{2} x {}_{2}+ β_{3} x {}_{3}+ β_{12} x {}_{1}x {}_{2}+ β_{13} x {}_{1}x {}_{3}+ β_{23} x {}_{2}x {}_{3}+ ε$

使用cordexch生成一个D最优设计,七分:

nfactors = 3;nruns = 7;(dCE, X) = cordexch (nfactors nruns,“互动”,“尝试”,10)dCE = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 X = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

设计矩阵的列X条件评估模型在每一行的设计dCE。条款中出现顺序从左到右:常数项,线性条件(1、2、3),交互条款(12日,13日,23)。使用X以适应模型中描述线性回归在设计点,响应数据测量dCE。

算法

这两个cordexch和rowexch使用迭代搜索算法。他们通过逐步改变一个初始设计矩阵X增加D= |X^TX|每一步。内置在两种算法,随机选择初始设计和选择的增量变化。结果,两种算法可能返回本地,但不是在全球范围内,D最优设计。多次运行每个算法,选择最好的为您的最终设计结果。两个函数有一个“尝试”参数,自动重复和比较。

与row-exchange算法所使用的rowexch,cordexch不使用一组候选人。(或者更确切地说,候选人设置是整个设计空间。)在每一步,coordinate-exchange算法交流单个元素X用一个新的元素在设计空间的邻近点评估。没有一组候选人降低要求记忆,但规模较小的搜索意味着coordinate-exchange算法更有可能被困在一个局部最小值。