Fantastic!
I really appreciate your quality of life tips! I’m really excited to try this now!
Last night I was thinking about how if I dotted b_fourier_monthly @ fourier.T (without indexing into b_fourier_monthly first), I would have shape [2, 732]. I would then want to index it with something in order to match the 732 long observed. I saw this comment of yours about multiindex. I wonder if that could have helped me here. If the multiindex had been ordered groups, date, then I think I might have been able to index with that directly. But I’m on thin ice here.