plotResiduals

Класс: LinearMixedModel

Постройте остаточные значения линейной модели смешанных эффектов

расширьте все на странице

Синтаксис

plotResiduals(lme,plottype)

plotResiduals(lme,plottype,Name,Value)

h = plotResiduals(___)

Описание

пример

plotResiduals(lme,plottype) строит необработанные условные остаточные значения линейной модели lme смешанных эффектов в графике типа, заданного plottype.

пример

plotResiduals(lme,plottype,Name,Value) также строит остаточные значения линейной модели lme смешанных эффектов с дополнительными опциями, заданными одним или несколькими аргументами пары "имя-значение". Например, можно задать остаточный тип, чтобы построить.

plotResiduals также принимает некоторые другие аргументы пары "имя-значение", которые задают свойства первичной линии в графике. Для тех пар "имя-значение" смотрите plot.

h = plotResiduals(___) возвращает указатель, h, к линиям или закрашенным фигурам в графике остаточных значений.

Входные параметры

развернуть все

`lme` — Линейная модель смешанных эффектов
`LinearMixedModel` объект

Линейная модель смешанных эффектов в виде LinearMixedModel объект создал использование fitlme или fitlmematrix.

`plottype` — Тип остаточного графика
`'histogram'` (значение по умолчанию) | `'caseorder'` | `'fitted'` | `'lagged'` | `'probability'` | `'symmetry'`

Тип остаточного графика в виде одного из следующих.

`'histogram'`	Значение по умолчанию. Гистограмма остаточных значений
`'caseorder'`	Остаточные значения по сравнению со случаем (строка) порядок
`'fitted'`	Остаточные значения по сравнению с подходящими значениями
`'lagged'`	Остаточные значения по сравнению с изолированной невязкой (r (t) по сравнению с r (t – 1))
`'probability'`	График нормального распределения
`'symmetry'`	График симметрии

Пример: plotResiduals(lme,'lagged')

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

`'ResidualType'` — Остаточный тип
`'Raw'` (значение по умолчанию) | `'Pearson'` | `'Standardized'`

Остаточный тип, заданный разделенной запятой парой, состоящей из ResidualType и одно из следующих.

Остаточный тип Условное выражение Крайний

Остаточный тип	Условное выражение	Крайний
`'Raw'`	$r_{i}^{C} = {[y - X \hat{β} - Z \hat{b}]}_{i}$	$r_{i}^{M} = {[y - X \hat{β}]}_{i}$
`'Pearson'`	$p r_{i}^{C} = \frac{r_{i}^{C}}{{\sqrt{[{\hat{V a r}}_{y, b} (y - X β - Z b)]}}_{i i}}$	$p r_{i}^{M} = \frac{r_{i}^{M}}{\sqrt{{[{\hat{V a r}}_{y} (y - X β)]}_{i i}}}$
`'Standardized'`	$s t_{i}^{C} = \frac{r_{i}^{C}}{\sqrt{{[{\hat{V a r}}_{y} (r^{C})]}_{i i}}}$	$s t_{i}^{M} = \frac{r_{i}^{M}}{\sqrt{{[{\hat{V a r}}_{y} (r^{M})]}_{i i}}}$

'Raw'

$r_{i}^{C} = {[y - X \hat{β} - Z \hat{b}]}_{i}$

$r_{i}^{M} = {[y - X \hat{β}]}_{i}$

'Pearson'

$p r_{i}^{C} = \frac{r_{i}^{C}}{{\sqrt{[{\hat{V a r}}_{y, b} (y - X β - Z b)]}}_{i i}}$

$p r_{i}^{M} = \frac{r_{i}^{M}}{\sqrt{{[{\hat{V a r}}_{y} (y - X β)]}_{i i}}}$

'Standardized'

$s t_{i}^{C} = \frac{r_{i}^{C}}{\sqrt{{[{\hat{V a r}}_{y} (r^{C})]}_{i i}}}$

$s t_{i}^{M} = \frac{r_{i}^{M}}{\sqrt{{[{\hat{V a r}}_{y} (r^{M})]}_{i i}}}$

Для получения дополнительной информации об условных и крайних остаточных значениях и остаточных отклонениях, смотрите Definitions в конце этой страницы.

Пример: 'ResidualType','Standardized'

Выходные аргументы

развернуть все

`h` — Обработайте к остаточному графику
указатель

Обработайте к остаточному графику, возвращенному как указатель.

Примеры

развернуть все

Исследуйте остаточные значения

Скрипт Open Live Script

Загрузите выборочные данные.

load('weight.mat')

weight содержит данные из продольного исследования, где 20 предметов случайным образом присвоены 4 программам подготовки, и их потеря веса зарегистрирована более чем шесть 2-недельных периодов времени. Это - симулированные данные.

Храните данные в таблице. Задайте Subject и Program как категориальные переменные.

tbl = table(InitialWeight,Program,Subject,Week,y);
tbl.Subject = categorical(tbl.Subject);
tbl.Program = categorical(tbl.Program);

Подбирайте линейную модель смешанных эффектов, где начальный вес, тип программы, неделя и взаимодействие между неделей и типом программы являются фиксированными эффектами. Точка пересечения и неделя варьируется предметом.

lme = fitlme(tbl,'y ~ InitialWeight + Program*Week + (Week|Subject)');

Постройте гистограмму необработанных остаточных значений.

plotResiduals(lme)

Figure contains an axes. The axes with title Histogram of residuals contains an object of type patch.

Постройте остаточные значения по сравнению с подходящими значениями.

plotResiduals(lme,'fitted')

Figure contains an axes. The axes with title Plot of residuals vs. fitted values contains 2 objects of type line.

Нет никакого очевидного шаблона, таким образом, нет никаких мгновенных знаков heteroscedasticity.

Создайте график нормального распределения остаточных значений.

plotResiduals(lme,'probability')

Figure contains an axes. The axes with title Normal probability plot of residuals contains 2 objects of type line.

Данные, кажется, нормальны.

Найдите номер наблюдения для данных, которые, кажется, выброс справа от графика.

find(residuals(lme)>0.25)

ans = 101

Создайте диаграмму сырых данных, Пирсона и стандартизированных остаточных значений.

r = residuals(lme);
pr = residuals(lme,'ResidualType','Pearson');
st = residuals(lme,'ResidualType','Standardized');
X = [r pr st];
boxplot(X,'labels',{'Raw','Pearson','Standardized'})

Figure contains an axes. The axes contains 21 objects of type line.

Все три диаграммы указывают на выброс на правом хвосте распределения. Диаграммы сырых данных и остаточных значений Пирсона также указывают на второй возможный выброс на левом хвосте. Найдите соответствующий номер наблюдения.

find(pr<-2)

ans = 10

Постройте необработанные остаточные значения по сравнению с изолированными остаточными значениями.

plotResiduals(lme,'lagged')

Figure contains an axes. The axes with title Plot of residuals vs. lagged residuals contains 3 objects of type line.

В графике нет никакого очевидного шаблона. Остаточные значения, кажется, не коррелируются.

Документация

plotResiduals

Синтаксис

Описание

Входные параметры

`lme` — Линейная модель смешанных эффектов
`LinearMixedModel` объект

`plottype` — Тип остаточного графика
`'histogram'` (значение по умолчанию) | `'caseorder'` | `'fitted'` | `'lagged'` | `'probability'` | `'symmetry'`

Аргументы в виде пар имя-значение

`'ResidualType'` — Остаточный тип
`'Raw'` (значение по умолчанию) | `'Pearson'` | `'Standardized'`

Выходные аргументы

`h` — Обработайте к остаточному графику
указатель

Примеры

Исследуйте остаточные значения

Смотрите также

Документация Statistics and Machine Learning Toolbox

Поддержка

Документация

plotResiduals

Синтаксис

Описание

Входные параметры

lme — Линейная модель смешанных эффектов LinearMixedModel объект

plottype — Тип остаточного графика 'histogram' (значение по умолчанию) | 'caseorder' | 'fitted' | 'lagged' | 'probability' | 'symmetry'

Аргументы в виде пар имя-значение

'ResidualType' — Остаточный тип 'Raw' (значение по умолчанию) | 'Pearson' | 'Standardized'

Выходные аргументы

h — Обработайте к остаточному графику указатель

Примеры

Исследуйте остаточные значения

Смотрите также

Документация Statistics and Machine Learning Toolbox

Поддержка

`lme` — Линейная модель смешанных эффектов
`LinearMixedModel` объект

`plottype` — Тип остаточного графика
`'histogram'` (значение по умолчанию) | `'caseorder'` | `'fitted'` | `'lagged'` | `'probability'` | `'symmetry'`

`'ResidualType'` — Остаточный тип
`'Raw'` (значение по умолчанию) | `'Pearson'` | `'Standardized'`

`h` — Обработайте к остаточному графику
указатель