Класс: LinearMixedModel
Вектор отклика линейной модели смешанных эффектов
Загрузите выборочные данные.
load('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')
