function [avg,avgr]=spe_avg(yo,th)
if nargin<=1, th=0.0498; end

[r,c]=size(yo);
avg=  (sum(yo')/c)';

validc=fix(c*(1-th));  %% exp(-5)= 6.74e-3
yos=sort(yo');
avgr=sum(yos(1:validc,:))'/validc; 


