One natural choice is to use absolute reconstruction error, like anomaly detection with deterministic auto-encoders [1]}. The anomaly score \(e_{m,t}\)  for observation \(x_{m,t}\)  is computed as:
\(e_{m,t} = \Vert x_{m,t} - \mu _{m,t} \Vert \) 
