
%find points outside 2 sd and remove them - do this by SELECTING AND
%KEEPING points within 2 sd.
[row col val]=find(abs(BadNoisy-modlBadNoisy)<SigScale*stdBadNoisy);
linearidx = sub2ind(size(BadNoisy),row,col);
cleanedBadNoisy=BadNoisy(linearidx);
%need to make new x vector (and y vector for plot - only use once 
%below) with just the cleaned points (x and y vectors have to be same sizes
%polyfit, polyvag, and plotting)
xclean=x(linearidx);
yclean=y(linearidx);
