What are the best way to get error estimates from a PDE solver?

There is not interface to EPSComputeError in the SLEPc plugin. You can use the options I gave in the earlier post in the sparams argument of your EPSSolve. If you want to store the errors in a file instead of being printed to screen, you should be able to do -eps_error_relative your_filename::ascii_info_detail.