How to create expression set from .txt files using R
2
0
Entering edit mode
Isadeghi66 • 0
@isadeghi66-12255
Last seen 7.2 years ago

Hi

I have some .txt files of RNA seq data form an experiment with 3 groups with 3 replicates each of which has the same features including genes ID and read counts.  How can I make expression set from these files using R regarding for example following files names and groups:

a.txt, b.txt, c.txt, d.txt, e.txt, f.txt, g.txt, h.txt, I.txt

of which a,b, c are in the same group, d, e, f in one group and g, h and I in one group. 

Thanks

Iman

RNA seq anaylsis R • 1.8k views
ADD COMMENT
1
Entering edit mode

You almost surely don't want an ExpressionSet. In addition, telling us the names of the files is not really useful. Instead it would be helpful if you were to say what the files contain, and what you want to do with those data once you have read into R. It would also help to say what you have tried, and any error messages you might have got.

ADD REPLY
0
Entering edit mode

All files contain  rows of  "ENTREZ  ID" and columns including Gene length and counts. I want to make an expression set in order to do differential expression analysis(DE) using DEseq2. I`m so new to R and even coding and I try to learn a bit about it. I have done DE analysis with edgeR and limma with these files, but was wondering how I can make expression sets too. I just don`t know what codes I have to start with. I have read the command for making Exp. set but the files features and contents were different. So, any suggestion?

Tnx

ADD REPLY
0
Entering edit mode
@james-w-macdonald-5106
Last seen 7 hours ago
United States

If you want to use DESeq2, then you should read the help page for DESeqDataSetFromMatrix, which is a function for reading in a matrix of counts.

But as I already mentioned, you don't want an ExpressionSet. I don't think DESeq2 uses anything but a RangedSummarizedExperiment which is not the same thing.

ADD COMMENT
0
Entering edit mode

OK. Thats cool. I will give it a shot. Thank you very much for being patient.
 

ADD REPLY
0
Entering edit mode
@gordon-smyth
Last seen 1 hour ago
WEHI, Melbourne, Australia

Given that you've done DE analyses with edgeR, can you not simply use the readDGEList() function? That is specifically designed to collate files containing counts when you have one file for each sample, which is what you seem to have.

ADD COMMENT

Login before adding your answer.

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