Класс: GeneralizedLinearMixedModel
Постройте остаточные значения обобщенной линейной модели смешанных эффектов
plotResiduals(
строит условные остаточные значения glme
,plottype
,Name,Value
)glme
использование дополнительных опций задано одним или несколькими Name,Value
парные аргументы. Например, можно задать, чтобы построить остаточные значения Пирсона.
возвращает указатель, h
= plotResiduals(___)h
, к линиям или закрашенным фигурам в графике остаточных значений.
glme
— Обобщенная линейная модель смешанных эффектовGeneralizedLinearMixedModel
объектОбобщенная линейная модель смешанных эффектов в виде GeneralizedLinearMixedModel
объект. Для свойств и методов этого объекта, смотрите GeneralizedLinearMixedModel
.
plottype
— Тип остаточного графика'histogram'
(значение по умолчанию) | 'caseorder'
| 'fitted'
| 'lagged'
| 'probability'
| 'symmetry'
Тип остаточного графика в виде одного из следующих.
Значение | Описание |
---|---|
'histogram' | Гистограмма остаточных значений |
'caseorder' | Остаточные значения по сравнению с порядком случая. Порядок случая совпадает с порядком строк, используемым во входных данных tbl при подборе кривой использованию модели fitglme . |
'fitted' | Остаточные значения по сравнению с подходящими значениями |
'lagged' | Остаточные значения по сравнению с изолированной невязкой (r (t) по сравнению с r (t – 1)) |
'probability' | График нормального распределения |
'symmetry' | График симметрии |
Пример: plotResiduals(glme,'lagged')
Задайте дополнительные разделенные запятой пары Name,Value
аргументы. Name
имя аргумента и Value
соответствующее значение. Name
должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN
.
'ResidualType'
— Остаточный тип'raw'
(значение по умолчанию) | 'Pearson'
Остаточный тип, заданный разделенной запятой парой, состоящей из ResidualType
и одно из следующих.
Остаточный тип | Формула |
---|---|
'raw' |
|
'Pearson' |
|
В каждом из этих уравнений:
yi является i th элемент n-by-1 вектор отклика, y, где i = 1..., n.
g-1 является обратной функцией ссылки для модели.
xi T является i th, строка фиксированных эффектов проектируют матричный X.
zi T является i th, строка случайных эффектов проектируют матричный Z.
δi является i th значение смещения.
σ 2 является дисперсионным параметром.
wi является i th вес наблюдения.
vi является термином отклонения для i th наблюдение.
μi является средним значением ответа для i th наблюдение.
и ориентировочные стоимости β и b.
Необработанные остаточные значения обобщенной линейной модели смешанных эффектов имеют непостоянное отклонение. Остаточные значения Пирсона, как ожидают, будут иметь приблизительно постоянное отклонение и обычно используются для анализа.
Пример: 'ResidualType','Pearson'
h
— Обработайте к остаточному графикуОбработайте к остаточному графику, возвращенному как графический объект. Можно использовать запись через точку, чтобы изменить определенные значения свойств объекта, включая цвет поверхности для гистограммы, и стиль маркера и цвет для scatterplot. Для получения дополнительной информации см. Доступ к значениям свойств.
Загрузите выборочные данные.
load mfr
Эти симулированные данные от компании-производителя, которая управляет 50 фабриками во всем мире с каждой фабрикой, запускающей процесс пакетной обработки, чтобы создать готовое изделие. Компания хочет сократить число дефектов в каждом пакете, таким образом, это разработало новый производственный процесс. Чтобы протестировать эффективность нового процесса, компания выбрала 20 своих фабрик наугад, чтобы участвовать в эксперименте: Десять фабрик реализовали новый процесс, в то время как другие десять продолжали запускать старый процесс. На каждой из этих 20 фабрик компания запустила пять пакетов (для в общей сложности 100 пакетов) и записала следующие данные:
Отметьте, чтобы указать, использовал ли пакет новый процесс (newprocess
)
Время вычислений для каждого пакета, в часах (time
)
Температура пакета, в градусах Цельсия (temp
)
Категориальная переменная, указывающая на поставщика (A
B
, или C
) из химиката, используемого в пакете (supplier
)
Количество дефектов в пакете (defects
)
Данные также включают time_dev
и temp_dev
, которые представляют абсолютное отклонение времени и температуры, соответственно, из стандарта процесса 3 часов на уровне 20 градусов Цельсия.
Подбирайте обобщенную линейную модель смешанных эффектов использование newprocess
, time_dev
, temp_dev
, и supplier
как предикторы фиксированных эффектов. Включайте термин случайных эффектов для точки пересечения, сгруппированной factory
, с учетом качественных различий, которые могут существовать из-за специфичных для фабрики изменений. Переменная отклика defects
имеет распределение Пуассона, и соответствующая функция ссылки для этой модели является журналом. Используйте подходящий метод Лапласа, чтобы оценить коэффициенты. Задайте фиктивную переменную, кодирующую как 'effects'
, таким образом, фиктивные переменные коэффициенты суммируют к 0.
Количество дефектов может быть смоделировано с помощью распределения Пуассона:
Это соответствует обобщенной линейной модели смешанных эффектов
где
количество дефектов, наблюдаемых в пакете, произведенном фабрикой во время пакета .
среднее количество дефектов, соответствующих фабрике (где ) во время пакета (где ).
, , и измерения для каждой переменной, которые соответствуют фабрике во время пакета . Например, указывает ли пакет, произведенный фабрикой во время пакета используемый новый процесс.
и фиктивные переменные, которые используют эффекты (сумма к нулю) кодирование, чтобы указать ли компания C
или B
, соответственно, предоставленный химикаты процесса для пакета производятся фабрикой во время пакета .
точка пересечения случайных эффектов для каждой фабрики это составляет специфичное для фабрики изменение по качеству.
glme = fitglme(mfr,'defects ~ 1 + newprocess + time_dev + temp_dev + supplier + (1|factory)','Distribution','Poisson','Link','log','FitMethod','Laplace','DummyVarCoding','effects');
Создайте диагностические графики с помощью остаточных значений Пирсона, чтобы протестировать предположения модели.
Постройте гистограмму, чтобы визуально подтвердить, что среднее значение остаточных значений Пирсона равно 0. Если модель будет правильна, мы ожидаем, что остаточные значения Пирсона будут сосредоточены в 0.
plotResiduals(glme,'histogram','ResidualType','Pearson')
Гистограмма показывает, что остаточные значения Пирсона сосредоточены в 0.
Постройте остаточные значения Пирсона по сравнению с подходящими значениями, чтобы проверять на знаки непостоянного отклонения среди остаточных значений (heteroscedasticity). Мы ожидаем, что условное выражение остаточные значения Пирсона будет иметь постоянное отклонение. Поэтому график условного выражения, остаточные значения Пирсона по сравнению с условным выражением соответствовали значениям, не должен показывать, что любая систематическая зависимость от условного выражения соответствовала значениям.
plotResiduals(glme,'fitted','ResidualType','Pearson')
График не показывает систематическую зависимость от подходящих значений, таким образом, нет никаких знаков непостоянного отклонения среди остаточных значений.
Постройте остаточные значения Пирсона по сравнению с изолированными остаточными значениями, чтобы проверять на корреляцию среди остаточных значений. Условное предположение независимости в GLME подразумевает, что условное выражение остаточные значения Пирсона является приблизительно некоррелированым.
plotResiduals(glme,'lagged','ResidualType','Pearson')
Нет никакого шаблона к графику, таким образом, нет никаких знаков корреляции среди остаточных значений.
fitglme
| fitted
| GeneralizedLinearMixedModel
| plot
| residuals
У вас есть модифицированная версия этого примера. Вы хотите открыть этот пример со своими редактированиями?
1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.
2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.
3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.
4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.
5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.