Apply formula to a set array of each group
3 views (last 30 days)
Show older comments
Hello everyone,
I have a set of 6672 observations that I was able to split into groups. Now I am trying to calculate the auto-correlation of a variable present in each group by using the code as follows:
Work = readtable('ARcal.xlsx');
Work.Groups = findgroups(Work.Numero);
ACor = splitapply(@(a,b) {corrcoef(a,b)},Work.MarketReturn(1:length(Work.MarketReturn)-1),Work.MarketReturn(2:length(Work.MarketReturn)),Work.Groups);
However, I get the following error message:
Error using splitapply (line 99)
The data variables must have the same number of rows as the vector of group numbers. The group number vector has 6672 row(s), and data variable 1 has 6671 row(s).
Error in ARcal (line 56)
ACor = splitapply(@(a,b) {corrcoef(a,b)},Work.MarketReturn(1:length(Work.MarketReturn)-1),Work.MarketReturn(2:length(Work.MarketReturn)),Work.Groups);
Any help would be much appreciated.
0 Comments
Answers (1)
Cris LaPierre
on 3 Sep 2020
It looks like you are using a table. If so, every variable in the table has the same number of rows. It looks like your table has 6672 rows. However, you inputs a and b only have 6671. This is because you are indexing Work.MarketReturn to remove 1 row (last for a, first for b), but you pass in the grouping variable, Work.Groups, without indexing.
Try this
ACor = splitapply(@(a,b) {corrcoef(a,b)},Work.MarketReturn(1:length(Work.MarketReturn)-1),Work.MarketReturn(2:length(Work.MarketReturn)),Work.Groups(1:end-1));
2 Comments
Cris LaPierre
on 4 Sep 2020
Edited: Cris LaPierre
on 4 Sep 2020
I'm afraid you're in a better position to answer that question. You need to specify 6671 groups. The obvious choices are
- Work.Groups(1:end-1)
- Work.Groups(2:end)
Basically, do you want to group your data by the first input to corrcoef or the second input?
See Also
Categories
Find more on Language Fundamentals in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!