Apply formula to a set array of each group
3 次查看(过去 30 天)
显示 更早的评论
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 个评论
回答(1 个)
Cris LaPierre
2020-9-3
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 个评论
Cris LaPierre
2020-9-4
编辑:Cris LaPierre
2020-9-4
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?
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!