Calculation of dynamical and many-body observables for spin-fermion models using the polynomial expansion method