Skip to content
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

inferRegulations function example #18

Open
Ying-Lab opened this issue Sep 23, 2024 · 0 comments
Open

inferRegulations function example #18

Ying-Lab opened this issue Sep 23, 2024 · 0 comments

Comments

@Ying-Lab
Copy link

Thank you for your great work!
I wonder if there is any examples in using the inferRegulations() function. I noticed that the A549 dataset was used to evaluated the reconstructed GRN. However, there exists no detailed usages in the example Walkthrough.

I tried to use this function after the scAI object is created.
My pipeline is as follows:
load("/home/zwh/program/scAI-master/data/data_A549.rda")
X <- data_A549$data # List of data matrix
labels <- data_A549$labels # the collected time of cells, which is used for validation
scAI_outs <- create_scAIobject(raw.data = X)
scAI_outs <- preprocessing(scAI_outs, assay = NULL, minFeatures = 200, minCells = 1,
libararyflag = F, logNormalize = F)
scAI_outs <- addpData(scAI_outs, pdata = labels, pdata.name = "Time")
scAI_outs <- run_scAI(scAI_outs, K = 5, nrun = 32)
markers_ATAC <- identifyFactorMarkers(scAI_outs, assay = 'ATAC', n.top = 5)
gene.use = c('GPR153','ACOT7')
scAI_GRN <- inferRegulations(scAI_outs,gene.use = gene.use ,candinate_loci = markers_ATAC)
image

However, the returned scAI_GRN has Null values. In addition, how to directly get the TF-GENE regulatory inference results. In scAI_GRN, I noticed that the results are presented in the TF-factor-loci format. But in the published paper, the reconstructed GRN on A549 dataset was shown in the TF-gene level.
Thank you in advance.
Scotty Zhang

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant