data=TemporalData[];

model=EstimatedProcess[data,MAProcess[3]]

residuals=Rest[data["PathStates"]-KalmanFilter[model, data]["PathStates"]];
Histogram[residuals]

DistributionFitTest[residuals,Automatic,{"TestDataTable",All}]