Problem with readDGE
2
0
Entering edit mode
iroda_55 • 0
@iroda_55-22861
Last seen 4.6 years ago

I need to analyze differentially expressed genes, for this, I try to use readDGE. Here is all commands what I did

if (!requireNamespace("BiocManager", quietly = TRUE))
  install.packages("BiocManager")
BiocManager::install(version = "3.11")
BiocManager::install("limma")
BiocManager::install("Glimma")
a
BiocManager::install("edgeR")
a
BiocManager::install("RNAseq123")
a
BiocManager::install((c("gplots", "RColorBrewer", "R.utils")))
a
files <- c("GSM1545535_10_6_5_11.txt", "GSM1545536_9_6_5_11.txt",
           "GSM1545538_purep53.txt", "GSM1545539_JMS8-2.txt",
           "GSM1545540_JMS8-3.txt", "GSM1545541_JMS8-4.txt",
           "GSM1545542_JMS8-5.txt", "GSM1545544_JMS9-P7c.txt",
           "GSM1545545_JMS9-P8c.txt")
read.delim(files[1], nrow = 5)
x <- readDGE(files, columns = c(1, 3))
BiocManager::install("readDGE")
a

But I got this message:

Warning message:
package ‘readDGE’ is not available (for R Under development) 

how can I handle this problem?

edger • 1.4k views
ADD COMMENT
1
Entering edit mode
Mike Smith ★ 6.5k
@mike-smith
Last seen 8 hours ago
EMBL Heidelberg

readDGE() is a function within the package edgeR.

You use BiocManager::install() to install new packages, which it looks like you've done in your example. This step only needs to be done once, they will still be installed next time you start R.

However the error you're seeing is because you're trying to install a package that doesn't exist. readDGE() is a function defined by the edgeR package.

Instead you need to load the package in your R session and then run the function. To load a package you need to use the library() command, and then you should be able to execute functions provided by that package. For example, the code below will load the edgeR package, and then run the readDGE() function.

library(edgeR)
x <- readDGE(files, columns = c(1, 3))
ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 5 hours ago
WEHI, Melbourne, Australia

It seems that you are following the limma RNA-seq workflow package:

limmaWorkflow.html

but you have omitted the "Setup" commands listed in Section 3.

In a workflow like this, the later commands depend on earlier steps so you can't generally run later commands without running the preceding code first.

ADD COMMENT

Login before adding your answer.

Traffic: 592 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