Landmark-based geometric morphometric methods are probably the most widely used approaches for shape analysis. Much work has been done for static or cross-sectional shape analysis while considerably less research has focused on dynamic or longitudinal shapes. The question of analysing shape changes over time is a fundamental issue in many research fields. In this paper, as a motivating example, we consider the problem of describing the dynamics of facial expressions for which medical and sociological studies call for a proper differential analysis to distinguish their different characteristics. We address the problem from an inferential point of view testing whether landmark positions change over time, within each facial expression, and whether these changes are different between different expressions. As the shape changes over time completely depend on geometrical landmarks, part of the problem becomes finding the subset of landmarks which best describes the dynamics of the expressions. In this paper, we show by means of a motivating example related to the analysis of the FG-NET (Face and Gesture Recognition Research Network) database with facial expressions and emotions from the Technical University Munich [Wallhoff, F. (2006), Database with Facial Expressions and Emotions from Technical University of Munich (FEEDTUM)'], that NonParametric Combination (NPC) tests can be effective tools when testing whether there is a difference between dynamics of facial expressions or testing which of the landmarks are more informative in explaining their dynamics. In particular, we start analysing data by means of bivariate linear mixed-effects models and then we improve inferential results using the NPC methodology