This work addresses the challenge of analysing the quality of human movements from visual information which has use in a broad range of applications, from diagnosis and rehabilitation to movement optimisation in sports science. Traditionally, such assessment is performed as a binary classification between normal and abnormal by comparison against normal and abnormal movement models, e.g. [5]. Since a single model of abnormal movement cannot encompass the variety of abnormalities, another class of methods only compares against one model of normal movement, e.g. [4]. We adopt this latter strategy and propose a continuous assessment of movement quality, rather than a binary classification, by quantifying the deviation from a normal model. In addition, while most methods can only analyse a movement after its completion e.g. [6], this assessment is performed on a frame-by-frame basis in order to allow fast system response in case of an emergency, such as a fall.Methods such as [4,6] are specific to one type of movement, mostly due to the features used. In this work, we aim to represent a large variety of movements by exploiting full body information. We use a depth camera and a skeleton tracker [3] to obtain the position of the main joints of the body, as seen in Fig. 1. We normalise this skeleton for global position and orientation of the camera, and for the varying height of the subjects, e.g. using Procrustes analysis.The normalised skeletons have high dimensionality and tend to contain outliers. Thus, the dimensionality is reduced using Diffusion Maps [1] which is modified by including the extension that Gerber et al. [2] presented to deal with outliers in Laplacian Eigenmaps. The resulting high level feature vector Y, obtained from the normalised skeleton at one frame, represents an individual pose and is used to build a statistical model of normal movement.Our statistical model is made up of two components that describe the normal poses and the normal dynamics of the movement. The pose model is in the form of the probability density function (pdf) f Y (y) of a random variable Y that takes as value y = Y our pose feature vector Y. The pdf is learnt from all the frames of training sequences that contain normal instances of the movement, using a Parzen window estimator. The quality of a new pose y t at frame t is then assessed as the log-likelihood of being described by the pose model, i.e. llh pose = log f Y (y t ) .(The dynamics model is represented as the pdf f Y t (y t |y 1 , . . . , y t−1 ) which describes the likelihood of a pose y t at a new frame t given the poses at the previous frames. In order to compute it, we introduce X t with value x t ∈ [0, 1], which is the stage of the (periodic or non-periodic) movement at frame t. Note, in the case of periodic movements, this movement stage can also be seen as the phase of the movement's cycle. Based on Markovian assumptions, we find thatwithx t an approximation of x t that minimises f {X 0 ,...,X t } (x 0 , . . . , x t |y 1 , . . . , y t ). f Y t (y t |x t ) is...