Question: HISAT2, Stringtie,DEseq2 differential gene expression
0
gravatar for sanchi.andrea
13 months ago by
sanchi.andrea0 wrote:

Hi there, I have got the .gtf file from stringtie and now I would like to run DEseq2. Unfortunately I should have a count table to use as input for DEseq2 and I actually havn't it. I have tried to us the python script prepDE.py as follows:

python prepDE.py -i myfile.gtf

I have got this error:

Error: Text file with sample ID and path invalid (1    havana    transcript    868240    870201    .    -    .    gene_id "ENSG00000230368"; transcript_id "ENST00000432963"; ref_gene_name "FAM41C"; cov "0.0"; FPKM "0.000000"; TPM "0.000000";)

The script can be checked here: http://ccb.jhu.edu/software/stringtie/dl/prepDE.py

Anyone know how to fix it? Thanks in advance

rna-seq • 1.7k views
ADD COMMENTlink modified 13 months ago by Devon Ryan91k • written 13 months ago by sanchi.andrea0

Which python version are you using?

ADD REPLYlink written 13 months ago by h.mon27k

I think is python 2 i think

ADD REPLYlink written 13 months ago by sanchi.andrea0

If you issue python --version, you will be sure.

ADD REPLYlink written 13 months ago by h.mon27k
1
gravatar for Devon Ryan
13 months ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k wrote:

Read the usage carefully. The -i option is described as follows:

the parent directory of the sample sub-directories or a textfile listing the paths to GTF files

So you need a text file with the name of the GTF file in it or the directory name that has the stringTie results.

As an aside, since you have human data there is unlikely to be any benefit to using stringTie. Either use something like featureCounts or salmon/kallisto.

ADD COMMENTlink written 13 months ago by Devon Ryan91k

I have tried putting directly the directory as follows:

python prepDE.py -i /media/andrea/Andrea/myfile.gtf

The error is the same:

Error: Text file with sample ID and path invalid (1 havana transcript 868240 870201 . - . gene_id "ENSG00000230368"; transcript_id "ENST00000432963"; ref_gene_name "FAM41C"; cov "0.0"; FPKM "0.000000"; TPM "0.000000";)

My python version is: Python 2.7.15rc1

thanks for any suggestions

ADD REPLYlink written 13 months ago by sanchi.andrea0

Try actually reading my reply.

ADD REPLYlink written 13 months ago by Devon Ryan91k

UPDATE: i fix the problem creating a txt file with tthe name of my file and then the path. Thanks

ADD REPLYlink written 13 months ago by sanchi.andrea0

You can accept @Devon's answer so everyone knows that your issue is solved

ADD REPLYlink written 13 months ago by grant.hovhannisyan1.7k
Please log in to add an answer.

Help
Access

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