-
Notifications
You must be signed in to change notification settings - Fork 77
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
gene_peak matrix in RP_model #93
Comments
Hi, how did you get the gene activity score? Did you use the enhanced model or simple model? if a peak is at the promoter, it will have a weight of 1, and then adding up other more distal peaks' regulatory potential, so it could be great than 1. If you used the I am doing a benchmark for those different ways to calculate gene activity score and will update here. Thanks! |
thanks, but i mean the gene_peak without multiplication by count matrix, the element in this gene_peak matrix represent the peak regulatory the gene potential. Should max(gene_peak) is equal to 1? But i got some values greater than 1. |
how did you get the gene_peak matrix? you are right, for every peak, the max weight should be 1. |
The value range of gene_peak matrix in RP_model should between 0 and 1. Why does "genes_peaks_score_csr" in scATAC_Genescore.py have a value greater than 1, with a maximum of 1600+.
The text was updated successfully, but these errors were encountered: