Problem with a rolling regression programme please help

Question

Andrew Wileman 2015-1-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/169308-problem-with-a-rolling-regression-programme-please-help

评论： Star Strider 2015-1-8

Hi,

I'm new to Matlab and was wondering if anyone may be able to help with a program I'm writing as I have seem to have hit a stumbling block. I'm trying to write a program that will will move along a data set using a window, at each increment, the program calls the function exp1fit, which calculates the regression coefficients for the exponential A and B within the window. Therefore a exponential function that fits the data in the window will be produced. The problem is getting the moving window to work (this is adapted from wreg code by Léon Bueckins), it keeps coming up with a 'Index exceeds matrix dimensions' warning and I can't seem to resolve it. I was wondering please, if anyone had any ideas.

Thanks, Andy

2 个评论
显示无隐藏无

Star Strider 2015-1-8

Please attach your .m file, and post the line throwing the error.

Andrew Wileman 2015-1-8

The problem seems to lie with ,: in the code

[A(:,i),B(:,i)] = exp1fit(x(s:s+window_length-1,:),y(s:s+window_length-1,:));

if this is removed it works, however it is not calculating the regression just in the window, but the whole function.

请先登录，再进行评论。

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

Answer 1

Star Strider 2015-1-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/169308-problem-with-a-rolling-regression-programme-please-help#answer_164344

在 MATLAB Online 中打开

One way to deal with that sort of problem is to use the minimum of s+window_length-1 and length(x).

Change the references to:

x(s:min(s+window_length-1, length(x)),:),y(s:min(s+window_length-1, length(y)),:)

and see if that makes a difference. (I didn’t test this specifically, but I’ve used it successfully in the past.) The idea is that it will use the shortest of the two options, so if all goes well, you will not read beyond the ends of your variable vectors. You can make that a bit more efficient by creating a separate index range:

idxrng = s:min(s+window_length-1, length(x));

and then just referring to it in your ‘x’ and ‘y’ subscripts. It looks like the author was attempting to do that by calculating ‘stoppingpoint’, but for some reason, it’s not working as it should.