y = response(lme)
Navigate to a folder containing sample data.
Load the sample data.
weight contains data from a longitudinal
study, where 20 subjects are randomly assigned to 4 exercise programs,
and their weight loss is recorded over two-week time periods. This
is simulated data.
Store the data in a table. Define
tbl = table(InitialWeight,Program,Subject,Week,y); tbl.Subject = nominal(tbl.Subject); tbl.Program = nominal(tbl.Program);
Fit a linear mixed-effects model where the initial weight, type of program, week, and the interaction between the week and type of program are the fixed effects. The intercept and week vary by subject.
lme = fitlme(tbl,'y ~ InitialWeight + Program*Week + (Week|Subject)');
Compute the fitted values and the response.
F = fitted(lme); y = response(lme);
Plot the response versus the fitted values.
plot(F,y,'bs') xlabel('Fitted Values') ylabel('Response')