This work deals with the local and global sensitivity of PSD for random responses analytically. The work is conducted utilizing the Pseudo-Excitation Method (PEM). Analytical formulae are derived on the basis of a direct algebraic method. The analytical formulae are compact, devoid of modal truncation problem and efficient to program. Global sensitivity analysis is implemented analytically by incorporating Gaussian process model.