Класс: 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')
