         offset = 0;
        end

        % Transform data
        tr = (sel_data-offset)/m(1);
        
        % Apply anscrombe transform
        antr = 2*sqrt(tr+3/8);
        
        % Calculate mean and std
        s = std(antr,0,2);
        mn = mean(antr,2);
        
        % Determine sum of square residuals
        r = norm(s-1).^2;
                
        if display_progress && mod(c,10) == 0
            hold(ax,'off');
            plot(ax,mn,ones(size(s)),'k'); 
           