The influence of the vessel shape, the initial conditions, and the vertical temperature gradient on dynamics and amount of disorder in convective patterns evolving in Bénard-Marangoni instability have been analyzed by using statistical tools, namely the density of defects, a disorder function, the order-disorder (m, σ) diagram introduced from the minimal spanning tree approach by Dussert et al.,[Phys. Rev. B 34, 3528 (1986)] and the entropy function recently defined by Loeffler (unpublished). Pattern disorder is studied for transient and steady states. Experimental results show that the disorder in the hexagonal patterns of Bénard-Marangoni convection (i) is minimized in a hexagonal vessel and (ii) can be described as a Gaussian noise superimposed on a perfect array of hexagonal cells. Starting from imposed arrays, both hexagonal and nonhexagonal, with a wavelength different from the one that is naturally selected, the final state is independent of initial conditions. Disorder increases with the distance from the threshold. Depending on the Prandtl number, different behaviors of the patterns are observed.© 1996 The American Physical Society.