How to average epochs_psd for multiple subjects?

Hi, everyone!
How to average the epochs_psd for multiple subjects? For evoked or TFR object, there is a way to average by mne.grand_average(). But what about the epochs_psd?