Формирование набора кандидатов
dC = candgen(nfactors,'model')
[dC,C] = candgen(nfactors,'model')
[...] = candgen(nfactors,'model','Name',value)
dC = candgen(nfactors,' генерирует набор кандидатов model')dC обработок, подходящих для оценки параметров в model с nfactors факторы. dC имеет nfactors столбцы и одна строка для каждого возможного лечения. model является одним из следующих:
'linear' - Постоянный и линейный члены. Это значение по умолчанию.
'interaction' - Постоянные, линейные и интерактивные термины
'quadratic' - постоянные, линейные, интерактивные и квадратные члены;
'purequadratic' - Постоянные, линейные и квадратные члены
В качестве альтернативы, model может быть матрицей, задающей многочлены произвольного порядка. В этом случае model должен иметь один столбец для каждого фактора и одну строку для каждого термина в модели. Записи в любой строке model являются полномочиями для факторов в столбцах. Например, если модель имеет коэффициенты X1, X2, и X3, затем ряд [0 1 2] в model задает термин (X1.^0).*(X2.^1).*(X3.^2). Строка всех нулей в model указывает постоянный член, который может быть опущен.
[dC,C] = candgen(nfactors,' также возвращает матрицу проектирования model')C оценивается на обработках в dC. Порядок столбцов C для полной квадратичной модели с n членами:
Постоянный член
Линейные члены порядка 1, 2,..., n
Члены взаимодействия в порядке (1, 2), (1, 3),..., (1, n), (2, 3),..., (n - 1, n)
Квадратные члены в порядке 1, 2,..., n
Другие модели используют подмножество этих терминов в том же порядке.
Проход C кому candexch для формирования D-оптимальной конструкции с использованием алгоритма обмена координатами.
[...] = candgen(nfactors,' указывает одну или несколько дополнительных пар имя/значение для проекта. Допустимые параметры и их значения перечислены в следующей таблице. Определить model','Name',value)Name внутри одиночных кавычек.
| Имя | Стоимость |
|---|---|
bounds | Нижняя и верхняя границы для каждого коэффициента, указанного как |
categorical | Индексы категориальных предикторов. |
levels | Вектор количества уровней для каждого фактора. |
В следующем примере используется rowexch для создания пятиэтапной конструкции для двухфакторной чистой квадратичной модели с использованием набора кандидатов, который производится внутри:
dRE1 = rowexch(2,5,'purequadratic','tries',10)
dRE1 =
-1 1
0 0
1 -1
1 0
1 1То же самое можно сделать с помощью candgen и candexch последовательно:
[dC,C] = candgen(2,'purequadratic') % Candidate set, C
dC =
-1 -1
0 -1
1 -1
-1 0
0 0
1 0
-1 1
0 1
1 1
C =
1 -1 -1 1 1
1 0 -1 0 1
1 1 -1 1 1
1 -1 0 1 0
1 0 0 0 0
1 1 0 1 0
1 -1 1 1 1
1 0 1 0 1
1 1 1 1 1
treatments = candexch(C,5,'tries',10) % Find D-opt subset
treatments =
2
1
7
3
4
dRE2 = dC(treatments,:) % Display design
dRE2 =
0 -1
-1 -1
-1 1
1 -1
-1 0