Is it possible to add a custom parameter for the degrees of freedom in the PowerDivergenceTest function which will be an adjustment to the df for the p-value? This may be useful for the goodness of fit tests ($G^2, \chi^2$) as the current implementation does not account for the number of parameres of the reference distribution. I think that implementation proposed in scipy.stats.power_divergence seems be the right way to do so (ddof parameter). Now, the user needs to remember to correct p-value after the test.