Added power_spectrum member function to ar1_posterior.py