fitlme "Index exceeds matrix dimensions."

Question

Bruno Bianchi 2015-11-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/256510-fitlme-index-exceeds-matrix-dimensions

评论： Bruno Bianchi 2015-11-20

Hi, I'm trying to use the new fitlme bult-in function (R2015a) and it keeps going on an error:

Index exceeds matrix dimensions.
Error in classreg.regr.LinearLikeMixedModel.makeSparseZ (line 1312)
                idxr     = accumarray(subs,vals(I),[],@(x) {x});
Error in LinearMixedModel/fitStandardLMEModel (line 1225)
        Zs = LinearMixedModel.makeSparseZ(Z,q,lev,Gid,N);
Error in LinearMixedModel/fitter (line 822)
        model.slme = fitStandardLMEModel(model);
Error in classreg.regr.FitObject/doFit (line 220)
            model = fitter(model);
Error in LinearMixedModel.fit (line 2395)
        model = doFit(model);
Error in fitlme (line 224)
    lme = LinearMixedModel.fit(ds,formula,varargin{:});

Also, running the fnuction takes too long (the error appears after several minutes of wating) and consume a lot of RAM. This is the line I'm running:

lmm1 = fitlme(datos_tabla, 'E1_T1 ~ pred + (1|SUJ)' );

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Gautam Pendse 2015-11-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/256510-fitlme-index-exceeds-matrix-dimensions#answer_200542

Hi Bruno,

This doesn't look like the intended behavior but it is difficult to say why you are getting that error. Could you post your data that reproduces the problem?

Gautam

3 个评论
显示 1更早的评论隐藏 1更早的评论

Gautam Pendse 2015-11-20

在 MATLAB Online 中打开

Hi Bruno,

Your table is of size 1-by-3 with 3 predictors. Each variable in the table is of size 1-by-N:

The error message can be improved but fitlme expects your table to be of size N-by-3 where each row represents 1 observation. In the code below, I modify your table by transposing variables and then fit the model:

%%Display the original table
test         = load('tabla_export.mat');
tabla_export = test.tabla_export;
size(tabla_export)
tabla_export
% Variables in your table are of size 1-by-N. The expectation is that your
% table will have N rows and each row would represent one observation. You
% can do this by just transposing the variables in your table as shown
% below.
%%Transpose 1-by-N variables into N-by-1 variables
E1_T1 = tabla_export.E1_T1;
pred  = tabla_export.pred;
SUJ   = tabla_export.SUJ;
E1_T1        = E1_T1';
pred         = pred';
SUJ          = SUJ';
newtbl       = table();
newtbl.E1_T1 = E1_T1;
newtbl.pred  = pred;
newtbl.SUJ   = SUJ;
%%Fit the model
rng(0,'twister')
lme = fitlme(newtbl,'E1_T1 ~ pred + (1|SUJ)','verbose',1,'startmethod','random')