Gene-peak association fails with Cellranger GTF file
jordanc17 opened this issue · 0 comments
jordanc17 commented
I'm trying to use the SEACells.genescores.get_gene_peak_correlations
function using the GTF file included in the Cell Ranger GRCh38-2020-A reference. When I try to use this GTF file, all of the gene-peak correlations are 0, but I can successfully get the gene-peak correlations if I use the hg38 GTF included in the tutorial, but not changing anything else about the function call. I'm not exactly sure why this is happening, since it seems like this function should work with other GTF files than the one in the tutorial. This is an example of my code and output
Code:
gene_set = rna_meta_ad.var_names[:1000]
gene_peak_cors = SEACells.genescores.get_gene_peak_correlations(atac_meta_ad, rna_meta_ad,
path_to_gtf=path_to_repository + '/data/genes.gtf',
span=100000,
n_jobs=1,
gene_set=gene_set)
gene_peak_cors[gene_peak_cors == 0]
Output:
LINC00115 0
NOC2L 0
KLHL17 0
ISG15 0
..
FCER1G 0
TOMM40L 0
NR1I3 0
MPZ 0
SDHC 0
Length: 1000, dtype: int64