Дробные генераторы факториального проекта
generators = fracfactgen(terms)
generators = fracfactgen(terms,k)
generators = fracfactgen(terms,k,R)
generators = fracfactgen(terms,k,R,basic)
generators = fracfactgen(terms)
использует алгоритм Франклина-Бейли, чтобы найти генераторы для наименьшего двухуровневого дробно-факториального проекта для оценки терминов линейной модели, заданных terms
. terms
- вектор символов или строковый скаляр, состоящий из слов, образованных из 52
буквы с учетом регистра a
- Z
, разделенный пространствами. Использование 'a'-'z'
для первых 26 факторов и, при необходимости, 'A'-'Z'
для остальных факторов. Для примера, terms = 'a b c ab ac'
. Однобуквенные слова указывают на основные эффекты, которые необходимо оценить; слова с несколькими буквами указывают на взаимодействие. Другой способ terms
является m -by - n матрицей 0
s и 1
s, где m - количество модельных членов, которые будут оценены, и n - количество факторов. Для примера, если terms
содержит строки [0 1 0 0]
и [1 0 0 1]
, затем коэффициент b
и взаимодействие между факторами a
и d
включены в модель. generators
- массив ячеек из векторов символов с одним генератором на камеру. Передайте generators
кому fracfact
получить дробно-факториальный проект и соответствующий смешивающий шаблон.
generators = fracfactgen(terms,k)
возвращает генераторы для двухуровневого дробно-факториального проекта с 2k
-запуски, если возможно. Если k
является []
, fracfactgen
находит наименьший проект.
generators = fracfactgen(terms,k,R)
находит проект с разрешением R
, по возможности. Разрешение по умолчанию является 3
.
Проект resolution R является такой, в которой никакое n-факторное взаимодействие не путается с любым другим эффектом, содержащим менее R- n факторов. Таким образом, проект разрешения III не путает основные эффекты друг с другом, но может смешивать их с двухсторонними взаимодействиями, в то время как проект разрешения IV не путает ни основные эффекты, ни двухсторонние взаимодействия, но может мешать двухсторонним взаимодействиям друг с другом.
Если fracfactgen
не может найти проект в запрошенном разрешении, он пытается найти проект с более низким разрешением, достаточный для калибровки модели. Если он успешен, он возвращает генераторы для проекта с более низким разрешением вместе с предупреждением. Если это не удается, возвращается ошибка.
generators = fracfactgen(terms,k,R,basic)
также принимает вектор basic
определение индексов факторов, которые должны рассматриваться как основные. Эти факторы получают полное факториальное лечение в проекте. По умолчанию включают факторы, которые являются частью взаимодействия высшего порядка в terms
.
Предположим, что вы хотите определить эффекты четырех двухуровневых факторов, для которых могут быть двухсторонние взаимодействия. Полный факториальный проект потребует 24 = 16 запуски. fracfactgen
функция находит генераторы для разрешения IV (разделения основных эффектов) дробно-факториального проекта, которая требует только 23 = 8 запуски:
generators = fracfactgen('a b c d',3,4) generators = 'a' 'b' 'c' 'abc'
Более экономичный проект и соответствующий смешивающий шаблон возвращаются fracfact
:
[dfF,confounding] = fracfact(generators) dfF = -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 confounding = 'Term' 'Generator' 'Confounding' 'X1' 'a' 'X1' 'X2' 'b' 'X2' 'X3' 'c' 'X3' 'X4' 'abc' 'X4' 'X1*X2' 'ab' 'X1*X2 + X3*X4' 'X1*X3' 'ac' 'X1*X3 + X2*X4' 'X1*X4' 'bc' 'X1*X4 + X2*X3' 'X2*X3' 'bc' 'X1*X4 + X2*X3' 'X2*X4' 'ac' 'X1*X3 + X2*X4' 'X3*X4' 'ab' 'X1*X2 + X3*X4'
Смешанный шаблон показывает, например, что двухстороннее взаимодействие между X1
и X2
смущается двухсторонним взаимодействием между X3
и X4
.
[1] Бокс, Г. Э. П., У. Г. Хантер и Дж. С. Хантер. Статистика для экспериментаторов. Hoboken, NJ: Wiley-Interscience, 1978.