broadinstitute / epi-SHARE-seq-pipeline

Epigenomics Program pipeline to analyze SHARE-seq data.
MIT License
17 stars 3 forks source link

update insert size plot script #137

Closed mei-knudson closed 1 year ago

sidwekhande commented 1 year ago

@mei-knudson it looks like the branch is failing at the joint_qc task as well, perhaps because the script hasn't been updated?

error here:

File "/usr/local/bin/joint_cell_plotting.py", line 179, in <module>
main()
File "/usr/local/bin/joint_cell_plotting.py", line 170, in main
plot_cells(metrics_df, pkr, min_umis, min_genes, min_tss, min_frags, plot_file)
File "/usr/local/bin/joint_cell_plotting.py", line 115, in plot_cells
max_x = max(df["frags"])
ValueError: max() arg is an empty sequence

Looks like the input atac_barcode_metadata file has different columns now, with the following columns: barcode fragments_promoter reads_tss reads_promoter tss_enrichment reads_tss_total reads_unique reads_duplicate pct_duplicates reads_peaks fragment_peaks raw_reads_nonmito raw_reads_mito pct_reads_promoter pct_reads_peaks pct_mito_reads

could you take a look and update the script?

mei-knudson commented 1 year ago

@mei-knudson it looks like the branch is failing at the joint_qc task as well, perhaps because the script hasn't been updated?

error here:

File "/usr/local/bin/joint_cell_plotting.py", line 179, in <module>
main()
File "/usr/local/bin/joint_cell_plotting.py", line 170, in main
plot_cells(metrics_df, pkr, min_umis, min_genes, min_tss, min_frags, plot_file)
File "/usr/local/bin/joint_cell_plotting.py", line 115, in plot_cells
max_x = max(df["frags"])
ValueError: max() arg is an empty sequence

Looks like the input atac_barcode_metadata file has different columns now, with the following columns: barcode fragments_promoter reads_tss reads_promoter tss_enrichment reads_tss_total reads_unique reads_duplicate pct_duplicates reads_peaks fragment_peaks raw_reads_nonmito raw_reads_mito pct_reads_promoter pct_reads_peaks pct_mito_reads

could you take a look and update the script?

fixed