An angular analysis of the rare decay B0s → φμ+μ− is presented, using proton-proton collision data collected by the LHCb experiment at centre-of-mass energies of 7, 8 and 13 TeV, corresponding to an integrated luminosity of 8.4 fb−1. The observables describing the angular distributions of the decay B0s → φμ+μ− are determined in regions of q2, the square of the dimuon invariant mass. The results are consistent with Standard Model predictions.