This is Tommy. I was trying to connect with you on Linkedin :) I know you rewrote Lisa and you have some script built-in to calculate the regulatory potential score (RP).
Given a bed file (peaks), a bigwig (intensity) and a gene model (e.g., gtf) file, I need to calculate the RP for every gene based on ATACseq data, do you have a script for that? I trust your implementation and do not want to re-invent the wheel.
Hi Allen,
This is Tommy. I was trying to connect with you on Linkedin :) I know you rewrote Lisa and you have some script built-in to calculate the regulatory potential score (RP).
Given a bed file (peaks), a bigwig (intensity) and a gene model (e.g., gtf) file, I need to calculate the RP for every gene based on ATACseq data, do you have a script for that? I trust your implementation and do not want to re-invent the wheel.
Thanks! Hope all is well with you PhD at Harvard!
Tommy