Генерация набора кандидатов
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 функция автоматически генерирует набор кандидатов, используя candgen, а затем создает D -оптимальный проект из этого набора кандидатов, используяcandexch. Звонить candexch отдельно, чтобы задать свой собственный набор кандидатов для алгоритма обмена строками.
Следующий пример использует 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