How to fit a general-linear mixed-effects model with categorical variables?

7 次查看（过去 30 天）

Falk Lieder 2015-10-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/246950-how-to-fit-a-general-linear-mixed-effects-model-with-categorical-variables

编辑： Falk Lieder 2015-10-6

Hi,

I am using the function fitglme from the statistics toolbox to fit a mixed-effects model with repeated measurements and categorical predictor variables as follows:

data_nr_acquisitions=table(nr_acquisitions,problem_type,block,subject);
%mixed-effects GLM that allows the effects of the problem type, the offset, and the
%block to vary randomly between subjects.
glme = fitglme(data_nr_acquisitions,'nr_acquisitions ~ problem_type + block + (problem_type| subject) + (block| subject) + (1|subject)');
anova(glme)

The variable problem_type is categorical, but when I run anova it says that problem_type has only one degree of freedom even though it has four possible values. This suggests that Matlab is treating it as a continuous regressors rather than as a categorical variable. Hence, something went wrong.

I tried to instruct fitglme to treat problem_type as a categorical variable with the argument "CategoricalVars" but unlike fitglm the function fitglme does not accept this argument. Can fitglme handle categorical variables and how can I get it to treat a variables as categorical?