The problem is: I have bivariate, normally distributed observations (x,y) to a model. But as I can't observe the original variables, I will observe them as (x+p,y+q) where p and q are some error terms which are assumed to be normal. Will the new data follow bivariate normal, or something else?