Statistics for Proper orthogonal decomposition for reduced order modeling: 2D heat flow