search

gwastro / pycbc-inference-paper

Data Release associated with the PyCBC Inference paper. This includes posterior samples for GW150914, LVT151012, and GW151226 along with configuration files and run scripts,
7 stars 5 forks source link

How to thinned the posterior samples? #1

Closed yanyuechuixue closed 6 years ago

yanyuechuixue commented 6 years ago

I use the script in run_files and get "gw150914_inference.hdf" file, but it is not same with "gw150914_posteriors_thinned.hdf".

Could you introduce the method of how to thinned the posterior samples?

cdcapano commented 6 years ago

You need to run pycbc_inference_extract_samples on gw150914_inference.hdf file, but I'm not sure what thin interval settings @soumide1102 used. @soumide1102 can you comment?

It would be good to add the extract samples command to the repo, either as separate scripts, or as an addendum to the run scripts.

soumide1102 commented 6 years ago

Hi @yanyuechuixue The "gw150914_inference_thinned.hdf" contains the thinned chains or independent samples extracted from the chains stored in the output file from the run "gw150914_inference.hdf". I've now added the commands for thinning the chains in scripts calledextract_independent_samples.sh in the events folders under run_files. You'll have to run the script for thinning with a version of PyCBC that is post v1.9.4 ( https://github.com/gwastro/pycbc/releases ).

yanyuechuixue commented 6 years ago

Thank you!