Question: count matrices for genes and transcripts for all samples
gravatar for fatimarasool135
10 months ago by
fatimarasool13560 wrote:

Hi According to this protocol

I have generated the gtf files for ballgown but now i want to use the Deseq. for the Deseq tool i have to generate the text file of gene and transcripts files. I have run the python script for one samples as python ./ballgown/sample1/sample1.gtf

I gave me 2 output files. I run this command for one sample. My number of samples is 3. Can i have to run the same command for other 2 samples?

ADD COMMENTlink modified 10 months ago by Haci370 • written 10 months ago by fatimarasool13560

I would use tximport to aggregate the abundances to the gene level and then feed this into DESeq2:

ADD REPLYlink written 10 months ago by ATpoint36k

Hi AT point, I have run this command for all samples python ./ballgown/sample1/sample1.gtf python ./ballgown/sample2/sample2.gtf python ./ballgown/sample3/sample3.gtf

Two output files are generated in each run. Then i renamed files according to sample in command. and now i have 6 files. for next step , I have to run this command for all 3 samples separately? Load gene(/transcript) count matrix and labels

countData <- as.matrix(read.csv("gene_count_matrix.csv", row.names="gene_id"))

I do not know how to work on Deseq2 . Can you please guide me in detail.

ADD REPLYlink written 10 months ago by fatimarasool13560
gravatar for Haci
10 months ago by
Haci370 wrote:

You do not have to run prepDE as many times as the number of your samples. The -i option (i meaning input) enables you to provide a list of samples, generating a single count file with genes/transcripts as rows and samples as columns, exactly what Deseq2 uses as input. Here is a

ADD COMMENTlink written 10 months ago by Haci370

Hi Haci, I have use this method. which is written in manual of stringtie. I did not use the i option with list of samples because in this method i got error of incorrect path so i use the other method. that is list below...I have run this command for 3 samples . now i have 6 files for gene and transcripts. These six files used in DEseq?




I am new to Deseq and R laguage Can you please gave me script of Deseq , I have ead the manual but cant understand it how to start it. I have just install it in my system

ADD REPLYlink modified 10 months ago • written 10 months ago by fatimarasool13560

Here is the Deseq2 tutorial:

The tutorial is very comprehensive and covers a lot scenarios for how to go from counts to the Deseq2 object. One way is the tximport, as suggested by ATPoint above. If you repeatedly encounter problems with prepDE you might consider using tximport instead.

Speaking of encountering problems with prepDE, please pay attention to the error messages, for example an error of "incorrect path" tells you to provide the "correct" path of the gtf files in the list of samples file. You can go to the directory where the gtf files sit and type "pwd" to print the full path of the directory. You should then append something like "/fileX.gtf" for each file.

ADD REPLYlink written 10 months ago by Haci370
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2016 users visited in the last hour