Bootstrap multi linear model

Question

Fabrice Lambert 2021-10-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1466394-bootstrap-multi-linear-model

评论： Fabrice Lambert 2021-10-5

Hi. I'm trying to bootstrap a multilinear model using bootstrp. I can get it to work with regress, but not with fitlm. Unfortunatly, the help doesn't provide an example for this. Here's some model code:

x = (1:1000)';

y1 = 2*x+1 + random('normal',0,2,1000,1);

y2 = 0.5*x - 2 + random('normal',0,2,1000,1);

z = 3 * y1 + 2 * y2 + 1 + random('normal',0,2,1000,1);

b1 = bootstrp(1000, @regress, z, [y1 y2 ones(1000,1)] ); % <-- This works

b2 = bootstrp(1000, @fitlm, [y1 y2], z); % < -- This doesn't work

Could someone please explain how to use fitlm withn bootstrp?

Thanks a lot,

Fabrice

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

请先登录，再进行评论。

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

Answer 1

Jeff Miller 2021-10-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1466394-bootstrap-multi-linear-model#answer_801511

The problem is that fitlm produces a complicated output structure and bootstrp doesn't know what to do with it.

x = (1:1000)';
y1 = 2*x+1 + random('normal',0,2,1000,1);
y2 = 0.5*x - 2 + random('normal',0,2,1000,1);
z = 3 * y1 + 2 * y2 + 1 + random('normal',0,2,1000,1);
% look at the very different outputs of regress vs fitlm:
b11 = regress(z,[y1 y2 ones(1000,1)])
b21 = fitlm([y1 y2],z)
% One way around the problem is to define & bootstrap your own version
% of the function fitlm that returns only a vector of values
b2 = bootstrp(1000, @myfitlm, [y1 y2], z); % < -- This works
function myout = myfitlm(preds,tobepred)
% select out the coefficients that are to be bootstrapped.
b21 = fitlm(preds,tobepred);
myout = b21.Coefficients.Estimate;
end

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

Jeff Miller 2021-10-5

In the myfitlm function, you could select out different or additional information from the b21 structure if you wanted to bootstrap something other than the coefficients. I think you can select whatever / as much as / you want from that structure as long as you assemble all of it into a numerical vector myout.

Fabrice Lambert 2021-10-5