Hi,
I want to report what I think it is a bug in the Table.write
API.
For testing purposes, I am writing a slice of gnomad’s variants dataset like so:
# hail 0.2.126
import hail as hl
ht = hl.read_table(
"gs://gcp-public-data--gnomad/release/3.1.2/ht/genomes/gnomad.genomes.v3.1.2.sites.ht",
_load_refs=False,
)
subset = ht.select_globals().head(20).to_spark(flatten=False)
subset.write(MY_OUTPUT_PATH)
Originally, I see that gnomad’s table does contain a list of globals. However, when I inspect my exported table, no globals are set:
my_subset = hl.read_table(MY_OUTPUT_PATH)
my_subset.describe()
>>> ----------------------------------------
Global fields:
None
----------------------------------------
...
Am I doing something wrong at the moment of reading/writing the data?
Thank you very much!
Irene