r3fang / SnapATAC

Analysis Pipeline for Single Cell ATAC-seq
GNU General Public License v3.0
301 stars 125 forks source link

Pull out feature topic probability matrix #97

Open StatsAreForLosers opened 5 years ago

StatsAreForLosers commented 5 years ago

Is there a place that snapATAC stores the feature topic probabilities matrix? I understand that LDA outputs a matrix of the probabiltiy a selected document (cell) has a specific topic (set of features) and from your source code that is stored in obj@smat@dmat. However, I am also curious if there is an output of the probability each future (genomic region) is used in the topic (set of features)? This can be useful in trying trying to find what genomic regions are distinguishing each topic and seeing which topics are enriched in each cluster.

r3fang commented 5 years ago

Hi there,

SnapATAC does not perform topic modeling, instead, it uses other type of dimensionality reduction technique that does not report the feature importance. Instead, you can use the function called findDAR to identify the differential peaks between clusters.

On Sep 13, 2019, at 10:40 AM, StatsAreForLosers notifications@github.com wrote:

Is there a place that snapATAC stores the feature topic probabilities matrix? I understand that LDA outputs a matrix of the probabiltiy a selected document (cell) has a specific topic (set of features) and from your source code that is stored in obj@smat@dmat. However, I am also curious if there is an output of the probability each future (genomic region) is used in the topic (set of features)? This can be useful in trying trying to find what genomic regions are distinguishing each topic and seeing which topics are enriched in each cluster.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/r3fang/SnapATAC/issues/97?email_source=notifications&email_token=ABT6GG33TXVYMQZB5CBHNIDQJPGB5A5CNFSM4IWSOMNKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4HLJKVKA, or mute the thread https://github.com/notifications/unsubscribe-auth/ABT6GG7GFKMYPO2UHF4GNR3QJPGB5ANCNFSM4IWSOMNA.