Problem with a rolling regression programme please help

1 次查看（过去 30 天）
Andrew Wileman
Andrew Wileman 2015-1-8
评论： 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
Star Strider 2015-1-8
Please attach your .m file, and post the line throwing the error.
Andrew Wileman
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.

请先登录，再进行评论。

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

采纳的回答

Star Strider
Star Strider 2015-1-8
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.

更多回答（1 个）

Andrew Wileman
Andrew Wileman 2015-1-8
Thanks for getting back to me.
The problem seems to be with lines 66 & 77 in the window_regression mfile.

类别

在 Help Center 和 File Exchange 中查找有关 Linear Regression 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by