How to use OutputFcn with fitnlm

Question

Adrian Batista Planas 2024-8-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2142791-how-to-use-outputfcn-with-fitnlm

评论： Matt J 2024-8-5

Hi all,

How do I execute an OutputFcn function in each iteration inside fitlm? I have tried to run the following code but it does not print anything or throw an exception.

  load carbig
    tbl = table(Horsepower,Weight,MPG);
    modelfun = @(b,x)b(1) + b(2)*x(:,1).^b(3) + ...
                        b(4)*x(:,2).^b(5);
    
    beta0 = [-50 500 -1 500 -1];
    opts = statset("fitnlm");
    opts.OutputFcn = @(x)disp(x);
    mdl = fitnlm(tbl,modelfun,beta0,"Options",opts);

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

dpb 2024-8-4

在 MATLAB Online 中打开

fitnlm doesn't document the 'OutputFcn' field in the options list; I'm guessing it's being ignored.

Try nlinfit instead.

This penchant for multiple nearly identical functions is maddening, indeed...

which -all fitnlm
/MATLAB/toolbox/stats/classreg/fitnlm.m
which -all nlinfit
/MATLAB/toolbox/stats/stats/nlinfit.m

请先登录，再进行评论。

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

Answer 1

Matt J 2024-8-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2142791-how-to-use-outputfcn-with-fitnlm#answer_1494431

编辑：Matt J 2024-8-5

在 MATLAB Online 中打开

Only Optimization Toolbox solvers have an OutputFcn option. Perhaps you are thinking of lsqcurvefit,

    load carbig
    tbl = table(Horsepower,Weight,MPG);
    modelfun = @(b,x)b(1) + b(2)*x(:,1).^b(3) + ...
                        b(4)*x(:,2).^b(5);
    
    beta0 = [-50 500 -1 500 -1];
    
    opts = optimoptions('lsqcurvefit',OutputFcn=@outfun);
    beta1 = lsqcurvefit(modelfun,tbl{:,1:2},tbl{:,3},beta0,"Options",opts);
function stop = outfun(x,optimValues,state)
    stop = false;
    % Check whether directional derivative norm is less than .01.
    if strcmp(state,'iter')
        disp(x)
    end
end

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

Torsten 2024-8-5

If statistics for the fitted parameter were accessible when using "lsqcurvefit", it would be mentionned on the documentation page as possible output from the solver:

https://uk.mathworks.com/help/optim/ug/lsqcurvefit.html

But it isn't.

Matt J 2024-8-5

在 MATLAB Online 中打开

@Adrian Batista Planas,

once you have solved with lsqcurvefit, you can take the solution beta1 and use it for beta0 with fitnlm:

    load carbig
    tbl = table(Horsepower,Weight,MPG);
    modelfun = @(b,x)b(1) + b(2)*x(:,1).^b(3) + ...
                        b(4)*x(:,2).^b(5);
    
    beta0 = [-50 500 -1 500 -1];
    
    opts = optimoptions('lsqcurvefit',OutputFcn=@outfun);
    beta1 = lsqcurvefit(modelfun,tbl{:,1:2},tbl{:,3},beta0,"Options",opts);
 
    mdl = fitnlm(tbl,modelfun,beta1,"Options",statset("fitnlm"));