Once the weights of samples are obtained, it is possible to estimate a multi-dimensional DoS for e(x) and other variables by making multidimensional bins and summing up the weight of all samples falling into each bin. In our study, to obtain h i (e), we use PA where the resampling scheme is implemented when the inverse temperature is changed. MCMC calculations with the fixed inverse temperature are performed by Metropolis sampling. EPA is the method to estimate DoS by combining PA and MHM.