We consider maximal operators acting on vector valued functions, that is, functions taking values on \(\mathbb {C}^d,\) that incorporate matrix weights in their definitions. We show vector valued estimates, in the sense of Fefferman–Stein inequalities, for such operators. These are proven using an extrapolation result for convex body valued functions due to Bownik and Cruz-Uribe. Finally, we show an \(\textrm{H}^1\) - \(\textrm{BMO}\) duality for matrix valued functions and we apply the previous vector valued estimates to show upper bounds for biparameter paraproducts. For the reader’s convenience, we include an appendix explaining how to adapt the extrapolation for real convex body valued functions of Bownik and Cruz-Uribe to the setting of complex convex body valued functions that we treat.
