For multimodal likelihood functions of fixed dimension the unimodal algorithm could be extended to build only contiguous (simply connected) regions using order statistics. We have not pursued this as we have devised the variable dimension posterior solution, which is more general (see next section). However, a multimodal solution based on order statistics could be investigated in the future.