D -оптимальное увеличение
dCE2 = daugment(dCE,mruns)
[dCE2,X] = daugment(dCE,mruns)
[dCE2,X] = daugment(dCE,mruns,model)
[dCE2,X] = daugment(...,param1,val1,param2,val2,...)
dCE2 = daugment(dCE,mruns) использует алгоритм координатного обмена, чтобы D -optimally add mruns выполняется в существующем экспериментальном проекте dCE для линейной аддитивной модели.
[dCE2,X] = daugment(dCE,mruns) также возвращает матрицу проекта X связана с дополненным проектом.
[dCE2,X] = daugment(dCE,mruns, использует линейную регрессионую модель, заданную в model)model. model является одним из следующих:
'linear' - Постоянные и линейные условия. Это значение по умолчанию.
'interaction' - Константа, линейная и условия взаимодействия
'quadratic' - Константа, линейный, взаимодействие и квадратные условия
'purequadratic' - Постоянные, линейные и квадратные условия
Порядок столбцов X для полной квадратичной модели с n терминами:
Постоянный член
Линейные условия в порядок 1 , 2,..., n
Условия взаимодействия по порядку (1, 2), (1, 3),..., (1, n), (2, 3),..., (n - 1, n)
Квадратные условия в порядок 1 , 2,..., n
Другие модели используют подмножество этих терминов в том же порядке.
Другой способ model может быть матрицей, задающей полиномиальные условия произвольного порядка. В этом случае model должен иметь по одному столбцу для каждого коэффициента и по одной строке для каждого члена в модели. Значения в любой строке model являются степенями для факторов в столбцах. Для примера, если в модели есть факторы X1, X2, и X3, затем строка [0 1 2] в model задает термин (X1.^0).*(X2.^1).*(X3.^2). Строка всех нулей в model задает постоянный термин, который может быть опущен.
[dCE2,X] = daugment(..., задает дополнительные пары параметр/значение для проекта. Действительные параметры и их значения перечислены в следующей таблице.param1,val1,param2,val2,...)
| Параметр | Значение |
|---|---|
'bounds' | Нижняя и верхняя границы для каждого фактора, заданные как |
'categorical' | Индексы категориальных предикторов. |
'display' | Либо |
'excludefun' | Указатель на функцию, которая исключает нежелательные запусков. Если функция f, она должна поддерживать синтаксис b = f (S), где S является матрицей обработок с |
'init' | Первоначальный проект как |
'levels' | Вектор количества уровней для каждого фактора. |
'maxiter' | Максимальное количество итераций. Значение по умолчанию является |
'options' | Значение является структурой, которая содержит опции, определяющие, вычислять ли несколько попыток параллельно, и определяющие, как использовать случайные числа при генерации начальных точек для попытки. Создайте структуру опций с
|
'tries' | Количество попыток сгенерировать проект из новой начальной точки. Алгоритм использует случайные точки для каждой попытки, кроме, возможно, первой. Значение по умолчанию является |
Примечание
daugment функция увеличивает существующий проект с помощью алгоритма координат-обменов; а 'start' параметром candexch функция обеспечивает ту же функциональность с помощью алгоритма обмена строками.
Следующий восьмиступенчатый проект подходит для оценки основных эффектов в четырехфакторной модели:
dCEmain = cordexch(4,8)
dCEmain =
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Чтобы оценить шесть членов взаимодействия в модели, дополните проект восемью дополнительными запусками:
dCEinteraction = daugment(dCEmain,8,'interaction')
dCEinteraction =
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 1 1 -1
1 1 -1 -1
1 -1 1 1
1 1 1 -1Дополненный проект является полным факториалом, с исходными восемью запускается в первых восьми строках.