Deleted:How to sort out “subscript out of bounds” error in GDCRNATools R package
0
0
Entering edit mode
Bhavya • 0
@d707106c
Last seen 4.4 years ago
India

I want to use the GDCRNATools R package. I am trying to use its gdcCEAnalysis function. In the example of this function, the following code is mentioned.

library(GDCRNATools)
 deLNC <- c('ENSG00000260920','ENSG00000242125','ENSG00000261211')
 dePC <- c('ENSG00000043355','ENSG00000109586','ENSG00000144355')
 genes <- c(deLNC, dePC)
 samples <- c('TCGA-2F-A9KO-01', 'TCGA-2F-A9KP-01', 
             'TCGA-2F-A9KQ-01', 'TCGA-2F-A9KR-01', 
             'TCGA-2F-A9KT-01', 'TCGA-2F-A9KW-01')
 rnaExpr <- data.frame(matrix(c(2.7,7.0,4.9,6.9,4.6,2.5,
                     0.5,2.5,5.7,6.5,4.9,3.8,
                     2.1,2.9,5.9,5.7,4.5,3.5,
                     2.7,5.9,4.5,5.8,5.2,3.0,
                     2.5,2.2,5.3,4.4,4.4,2.9,
                     2.4,3.8,6.2,3.8,3.8,4.2),6,6), 
                     stringsAsFactors=FALSE)
 rownames(rnaExpr) <- genes
 colnames(rnaExpr) <- samples

 mirExpr <- data.frame(matrix(c(7.7,7.4,7.9,8.9,8.6,9.5,
                     5.1,4.4,5.5,8.5,4.4,3.5,
                     4.9,5.5,6.9,6.1,5.5,4.1,
                     12.4,13.5,15.1,15.4,13.0,12.8,
                     2.5,2.2,5.3,4.4,4.4,2.9,
                     2.4,2.7,6.2,1.5,4.4,4.2),6,6),
                     stringsAsFactors=FALSE)
 colnames(mirExpr) <- samples
 rownames(mirExpr) <- c('hsa-miR-340-5p','hsa-miR-181b-5p',
                     'hsa-miR-181a-5p', 'hsa-miR-181c-5p',
                     'hsa-miR-199b-5p','hsa-miR-182-5p')

 ceOutput <- gdcCEAnalysis(lnc       = deLNC, 
                         pc          = dePC, 
                         lnc.targets = 'starBase', 
                         pc.targets  = 'starBase', 
                         rna.expr    = rnaExpr, 
                         mir.expr    = mirExpr)

The above mentioned code gives the following output:

Step 1/3: Hypergenometric test done !
Step 2/3: Correlation analysis done !
Step 3/3: Regulation pattern analysis done !

However, if the values in deLNC are changed, the gdcCEAnalys function shows the following error:

deLNC <- c('ENSG00000000003','ENSG00000099869','ENSG00000261211')
#rest of the code is same.

Error in `colnames<-`(`*tmp*`, value = c("lncRNAs", "Genes", "Counts",  : 
  attempt to set 'colnames' on an object with less than two dimensions

Please help me in understanding the cause of this error. The details of the package and function can be found here: https://rdrr.io/bioc/GDCRNATools/src/R/gdcCEAnalysis.R https://rdrr.io/bioc/GDCRNATools/man/gdcCEAnalysis.html

mir.expr gdcCEAnalysis subscript_out_of_bounds GDCRNATools • 751 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 897 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6