Open maxlarosa opened 5 years ago
When I try the following code
library(GDCRNATools) library(limma) data(lncTarget) data("pcTarget") data("DEGAll") data("mirCounts") data("rnaCounts")
rnaExpr <- gdcVoomNormalization(counts = rnaCounts, filter = FALSE) ####### Normalization of miRNAs data ####### mirExpr <- gdcVoomNormalization(counts = mirCounts, filter = FALSE)
deLNC <- gdcDEReport(deg = DEGAll, gene.type = 'long_non_coding') dePC <- gdcDEReport(deg = DEGAll, gene.type = 'protein_coding') ceOutput <- gdcCEAnalysis(lnc = rownames(deLNC), pc = rownames(dePC), lnc.targets = lncTarget, pc.targets = pcTarget, rna.expr = rnaExpr, mir.expr = mirExpr)
I get this error: Error in rna.expr[pc, ] : subscript out of bounds
Hello! Have you solved this problem now ?I encountered the same error with Error in mir.expr[mir, ] : subscript out of bounds. If you already know why and how to avoid such an error, I appreciate your reply!
When I try the following code
library(GDCRNATools) library(limma) data(lncTarget) data("pcTarget") data("DEGAll") data("mirCounts") data("rnaCounts")
Normalization of RNAseq data
rnaExpr <- gdcVoomNormalization(counts = rnaCounts, filter = FALSE) ####### Normalization of miRNAs data ####### mirExpr <- gdcVoomNormalization(counts = mirCounts, filter = FALSE)
deLNC <- gdcDEReport(deg = DEGAll, gene.type = 'long_non_coding') dePC <- gdcDEReport(deg = DEGAll, gene.type = 'protein_coding') ceOutput <- gdcCEAnalysis(lnc = rownames(deLNC), pc = rownames(dePC), lnc.targets = lncTarget, pc.targets = pcTarget, rna.expr = rnaExpr, mir.expr = mirExpr)
I get this error: Error in rna.expr[pc, ] : subscript out of bounds