Help with linear mixed effects model formula

16 次查看（过去 30 天）
Bryce Grier
Bryce Grier 2020-2-8
评论： Bryce Grier 2020-2-11
Hello,
I have a repeated measures data set with a response variable and two predictor variables (stimulus frequency and stimulus number). There are missing data, but it is only random for stimulus. Explained further, if I tested a given frequency in a subject, there is an observation for every possible stimulus number with that frequency, but not every frequency was tested in every subject. Due to the missing data, I can't run a repeated measures two-way ANOVA, so, as I understand it, I need to fit a linear mixed effects model.
My data are arranged in a table with 4 columns (Subject, Frequency, Stimulus, Response), with NaNs in the place of missing data. Additionally, I have reason to believe that there is an interaction between the predictors.
The formula that I'm currently using is:
lme = fitlme(mixedModelData,'Response ~ 1 + Stimulus*Frequency + (1|Subject)');
Does this model appropriately address the situation I have described?
Thanks!

请先登录，再回答此问题。

采纳的回答

Jeff Miller
Jeff Miller 2020-2-9
That model is appropriate if you want to assume these things (in addition to normality, independence, etc):
(a) subjects differ only in their overall response levels; that is, all subjects are affected identically by stimulus, frequency, and their interaction.
(b) the numbers used in the Stimulus and Frequency columns are real numbers to which the response is linearly related (i.e., they are not just arbitrary indicators of categories).
(c) the multiplication of the stimulus*frequency numbers gives a meaningful numerical interaction quantity to which the response is linearly related.
  3 个评论
显示 1更早的评论
Jeff Miller
Jeff Miller 2020-2-10
I think you probably want this, which allows Ss to be affected differently by stim, freq, & the interaction:
lme = fitlme(mixedModelData,'Response ~ 1 + Stimulus*Frequency + (Stimulus*Frequency|Subject)');
Bryce Grier
Bryce Grier 2020-2-11
Great, thanks again. The AIC and BIC suggest that your model is the most appropriate as well.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Analysis of Variance and Covariance 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by