The study of global event shape variables can provide sensitive tests of predictions for multijet production in proton-proton collisions. This paper presents a study of several event shape variables calculated using jet four momenta in proton-proton collisions at a centre-of-mass energy of 13 TeV and uses data recorded with the CMS detector at the LHC corresponding to an integrated luminosity of 2.2 fb$^{−1}$. After correcting for detector effects, the resulting distributions are compared with several theoretical predictions. The agreement generally improves as the energy, represented by the average transverse momentum of the two leading jets, increases.