Save coefficients from fitlme function in a table

Question

alphabetagamma 2022-7-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1770925-save-coefficients-from-fitlme-function-in-a-table

评论： Abderrahim. B 2022-7-30

I am running a fixed effects model using the lme function and would like to save the coefficients in a table. I tried the following

coefficient_table = cell2table(cell(0,4), 'VariableNames', {'Estimate', 'SE', 'tStat', 'pValue'});
fe = fitlme(data_table, strcat(var1, '~', var2));
value = fe.Coefficients(2:end, :);

However, class(value) is 'classreg.regr.lmeutils.titleddataset', and not a table. The same problem does not arise when I fit a linear regreression model with fitlm function. Then, the class(value) = table

Ultimately, I want to save all the coefficients in

coefficient_table = [coefficient_table; value];

How can I convert the class of value to a table? Or is there another way to save the coefficients into a table? I'd really appreciate your help.

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

请先登录，再进行评论。

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

Answer 1

Abderrahim. B 2022-7-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1770925-save-coefficients-from-fitlme-function-in-a-table#answer_1018040

在 MATLAB Online 中打开

Hi

You need to use dataset2table.

Demo:

load imports-85
tbl = table(X(:,12),X(:,14),X(:,24),'VariableNames',{'Horsepower','CityMPG','EngineType'});
head(tbl, 5)
ans = 5×3 table
    Horsepower    CityMPG    EngineType
    __________    _______    __________

       111          21           13    
       111          21           13    
       154          19           37    
       102          24           35    
       115          18           35    

% 
lme = fitlme(tbl,'CityMPG~Horsepower+(1|EngineType)+(Horsepower-1|EngineType)');
% Retrieve and convert to table fitlme coefficients
coef_tbl = dataset2table(lme.Coefficients)
coef_tbl = 2×8 table
         Name          Estimate      SE       tStat     DF       pValue       Lower        Upper  
    _______________    ________    _______    ______    ___    __________    ________    _________

    {'(Intercept)'}      37.276     2.8556    13.054    201    1.3147e-28      31.645       42.906
    {'Horsepower' }    -0.12631    0.02284     -5.53    201    9.8848e-08    -0.17134    -0.081269
% In case you want to convert Name from cell to categorical. 
coef_tbl.Name = categorical(coef_tbl.Name)
coef_tbl = 2×8 table
       Name        Estimate      SE       tStat     DF       pValue       Lower        Upper  
    ___________    ________    _______    ______    ___    __________    ________    _________

    (Intercept)      37.276     2.8556    13.054    201    1.3147e-28      31.645       42.906
    Horsepower     -0.12631    0.02284     -5.53    201    9.8848e-08    -0.17134    -0.081269