Calculation of Dynamical and Many-Body Observables with the Polynomial Expansion Method for Spin-Fermion Models