Класс: LinearMixedModel
Вектор отклика модели линейных смешанных эффектов
Загрузите выборочные данные.
load('weight.mat');weight содержит данные продольного исследования, где 20 субъектов случайным образом назначены для 4 программ упражнений, и их потеря веса регистрируется в течение двухнедельных периодов времени. Это моделируемые данные.
Сохраните данные в таблице. Определите 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')
