Класс: LinearMixedModel
Вектор отклика линейной модели смешанных эффектов
y = response(lme)
y = response(lme), y раньше соответствовал линейной модели lme смешанных эффектов.
Загрузите выборочные данные.
load(fullfile(matlabroot,'examples','stats','weight.mat'));
weight содержит данные из продольного исследования, где 20 предметов случайным образом присвоены 4 программам подготовки, и их потеря веса зарегистрирована по 2D недельным периодам времени. Это - моделируемые данные.
Храните данные в таблице. Задайте 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); y = response(lme);
Постройте ответ по сравнению с подходящими значениями.
plot(F,y,'bs') xlabel('Fitted Values') ylabel('Response')
