[coeff,SE,PVAL,INMODEL,STATS,NEXTSTEP,HISTORY]= stepwisefit(XMatrix(:,2:end),Y,'maxiter',200,'display','off');
    coeff=[STATS.intercept; coeff.*INMODEL'];
    STATS;
    Y_predicted =[];