You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
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)
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
The text was updated successfully, but these errors were encountered:
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)
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
The text was updated successfully, but these errors were encountered: