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-оптимального добавления 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Увеличенный дизайн является полным факториалом, с первоначальными восемью прогонами в первых восьми рядах.