невязки

Класс: LinearMixedModel

Невязки подходящей линейной модели смешанных эффектов

Синтаксис

R = residuals(lme)
R = residuals(lme,Name,Value)

Описание

пример

R = residuals(lme) возвращает необработанные условные невязки в подходящую линейную модель lme смешанных эффектов.

пример

R = residuals(lme,Name,Value) возвращает невязки в линейную модель lme смешанных эффектов с дополнительными опциями, заданными одним или несколькими аргументами пары Name,Value.

Например, можно задать Пирсона или стандартизированные невязки или невязки с вкладами только от фиксированных эффектов.

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

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

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

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

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

Индикатор для условных невязок, заданных как пара, разделенная запятой, состоящая из 'Conditional' и одно из следующих.

TrueВклад и от зафиксированных эффектов и от случайных эффектов (условное выражение)
FalseВклад только от фиксированных (крайних) эффектов

Пример: 'Conditional,'False'

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

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

riC=[yXβ^Zb^]i

riM=[yXβ^]i

'Pearson'

priC=riC[Var^y,b(yXβZb)]ii

priM=riM[Var^y(yXβ)]ii

'Standardized'

stiC=riC[Var^y(rC)]ii

stiM=riM[Var^y(rM)]ii

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

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

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

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

Невязки подходящей линейной модели lme смешанных эффектов возвратились как n-by-1 вектор, где n является количеством наблюдений.

Примеры

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

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

load(fullfile(matlabroot,'examples','stats','weight.mat'));

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

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

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

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

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

Вычислите подходящие значения и необработанные невязки.

F = fitted(lme);
R = residuals(lme);

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

plot(F,R,'bx')
xlabel('Fitted Values')
ylabel('Residuals')

Теперь, постройте невязки по сравнению с подходящими значениями, сгруппированными программой.

figure();
gscatter(F,R,Program)

Невязки, кажется, ведут себя так же через уровни программы как ожидалось.

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

load carbig

Сохраните переменные для миль на галлон (MPG), ускорения, лошадиной силы, цилиндров, и модельный год в таблице.

tbl = table(MPG,Acceleration,Horsepower,Cylinders,Model_Year);

Соответствуйте линейной модели смешанных эффектов для миль на галлон (MPG), с фиксированными эффектами для ускорения, лошадиной силы и цилиндров, и потенциально коррелировал случайные эффекты для прерывания и ускорения, сгруппированного модельным годом.

lme = fitlme(tbl,'MPG ~ Acceleration + Horsepower + Cylinders + (Acceleration|Model_Year)');

Вычислите условное выражение невязки Пирсона и отобразите первые пять невязок.

PR = residuals(lme,'ResidualType','Pearson');
PR(1:5)
ans = 5×1

   -0.0533
    0.0652
    0.3655
   -0.0106
   -0.3340

Вычислите крайние невязки Пирсона и отобразите первые пять невязок.

PRM = residuals(lme,'ResidualType','Pearson','Conditional',false);
PRM(1:5)
ans = 5×1

   -0.1250
    0.0130
    0.3242
   -0.0861
   -0.3006

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

load carbig

Сохраните переменные для миль на галлон (MPG), ускорения, лошадиной силы, цилиндров, и модельный год в таблице.

tbl = table(MPG,Acceleration,Horsepower,Cylinders,Model_Year);

Соответствуйте линейной модели смешанных эффектов для миль на галлон (MPG), с фиксированными эффектами для ускорения, лошадиной силы и цилиндров, и потенциально коррелировал случайные эффекты для прерывания и ускорения, сгруппированного модельным годом.

lme = fitlme(tbl,'MPG ~ Acceleration + Horsepower + Cylinders + (Acceleration|Model_Year)');

Чертите гистограмму необработанных невязок с нормальной подгонкой.

r = residuals(lme);
histfit(r)

Нормальное распределение, кажется, подходящий вариант для невязок.

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

pr = residuals(lme,'ResidualType','Pearson');
st = residuals(lme,'ResidualType','Standardized');
X = [r pr st];
figure();
boxplot(X)

Красные знаки "плюс" показывают наблюдения с невязками выше или ниже q3+1.5(q3-q1) и q1-1.5(q3-q1), где q1 и q3 25-е и 75-е процентили, соответственно.

Найдите наблюдения с невязками, которые являются 2,5 стандартными отклонениями выше и ниже среднего значения.

find(r > nanmean(r) + 2.5*nanstd(r))
ans = 7×1

    62
   252
   255
   330
   337
   341
   396

find(r < nanmean(r) - 2.5*nanstd(r))
ans = 3×1

   119
   324
   375

Больше о

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

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

| | |

Для просмотра документации необходимо авторизоваться на сайте