So I stored a run of about 1m samples (takes up about 8gb) using the hdfbackend. I've used it several times previously. However, today I went to open the data again with the following commands:
"
reader = emcee.backends.HDFBackend("QHO_3_sample1.h5")
tau_3 = reader.get_autocorr_time()
burnin_3 = int(2 * np.max(tau_3))
thin_3 = int(0.5 * np.min(tau_3))
samples_3 = reader.get_chain()
"
And it tells me "
You must run the sampler with 'store == True' before accessing the results." Further, when I try to run reader.iteration, it tells me that there are no iterations. However, I very clearly see the 8gb file in my documents... not sure what is going on. Any ideas?