function data=trim(data,outval)
%data=trim(data,outval)
% 祛除壞數(shù)據(jù)，包括NaN, Inf, 和異常大小數(shù)據(jù)
% data: 列狀數(shù)據(jù)，每列來自一個(gè)總體,返回凈化后的數(shù)據(jù)
% outval: 系數(shù)因子，離均值超過outval倍標(biāo)準(zhǔn)差被判為異常大小，默認(rèn)為4

% L.J.HU 8-17-1999

if nargin<2,outval=4;end
outliers=(isnan(data)|abs(data)==inf);
[n,m]=size(data);mu=mean(data);sigma=std(data);
outliers=outliers+(abs(data-ones(n,1)*mu)>outval*ones(n,1)*sigma);
if m>1,
   data(any(outliers'),:)=[];
else
  data(find(outliers'),:)=[];
end