Skip to content

Commit

Permalink
added m*MMD correction for quadratic time MMD statistic
Browse files Browse the repository at this point in the history
  • Loading branch information
karlnapf committed Aug 11, 2012
1 parent d3c42f6 commit 42992b0
Showing 1 changed file with 2 additions and 2 deletions.
Expand Up @@ -58,14 +58,14 @@
if "sample_null_spectrum" in dir(QuadraticTimeMMD):
mmd.set_null_approximation_method(MMD2_SPECTRUM)
mmd.set_statistic_type(BIASED)
null_samples_spectrum=mmd.sample_null_spectrum(num_bootstrap, m-10)/m
null_samples_spectrum=mmd.sample_null_spectrum(num_bootstrap, m-10)

# fit gamma distribution, biased statistic
mmd.set_null_approximation_method(MMD2_GAMMA)
mmd.set_statistic_type(BIASED)
gamma_params=mmd.fit_null_gamma()
# sample gamma with parameters
null_samples_gamma=array([gamma(gamma_params[0], gamma_params[1]) for _ in range(num_bootstrap)])/m
null_samples_gamma=array([gamma(gamma_params[0], gamma_params[1]) for _ in range(num_bootstrap)])


# plot
Expand Down

0 comments on commit 42992b0

Please sign in to comment.