Regression Analysis Confusion!

%%Regression analysis for ice drift anomalies    
for i=1:361,
        for j=1:361,
           if ((~isnan(anom_wint_r(i,j,:))))
               stats = regstats(squeeze(anom_wint_r(i,j,:)),years','linear',{'beta' 'rsquare' 'fstat'});        
              idanom_trend_per_winter(i,j) = stats.beta(2); 

             else
                  idanom_trend_per_winter(i,j) = NaN;
              end
          end
  end

anom_wint_r is in cm/s, but i am plotting idanom_trend_per_winter i it unitless or is it cm/s/year .. Am i plotting the p-value here

This is the part that confuses me, any explanation would be helpful idanom_trend_per_winter(i,j) = stats.beta(2);