Dataset for "Inference of Neutron Star Mass Distributions and the Equation of State from Multi-messenger Observations"

This dataset contains the sampled data produced using the affine-invariant Markov Chain Monte Carlo (MCMC) method. The samples represent the posterior distributions in a Bayesian inference, where the likelihoods were computed using BAMR and the MCMC was performed using its dependency O2SCL. Both codes are written in C++ and the simulations were run on local Linux and HPC systems. All outputs are stored in HDF5 (Hierarchical Data Format, version 5), a structured and portable format that organizes data in a tree-like hierarchy of groups and datasets. Each file includes a table-like group named markov_chain_0, with the thinned samples stored in the object markov_chain_0/data. In this structure, rows correspond to individual samples and columns correspond to the parameters and derived output quantities. This table object can be accessed as follows: import h5pyf = h5py.File("filename", "r")data = f["markov_chain_0/data"] # full table of samplescol_names = list(data.keys()) # get column namessamples = data[col_names[i]] # access specific columns These files provide the complete set of samples needed to reproduce the posterior distributions, statistical summaries, and figures associated with the corresponding analysis.