Question: Running prepare.py on remote server
0
gravatar for suneetamodekurty77
4 weeks ago by
suneetamodekurty770 wrote:

Hi
Greetings
I am trying to run prepDE.py on remote server of my university but
There is an error

python: can't open file 'prepDE.py': [Errno 2] No such file or directory

Please advise an alternative to get the count matrix with stringtie.
I used the -A option with stringtie and got individual count matrices for each file ( controlled and treated separately) I want to do DeSeq2
How do I proceed?
Please give inputs
Thank you!!

prepde.py rna-seq stringtie • 163 views
ADD COMMENTlink modified 4 weeks ago by kristoffer.vittingseerup2.0k • written 4 weeks ago by suneetamodekurty770
  1. Use absolute path for everything to ensure you're not making path based mistakes
  2. Talk to your university sysadmin
  3. Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.
    code_formatting
ADD REPLYlink written 4 weeks ago by RamRS22k

You can use featurecounts (http://bioinf.wehi.edu.au/featureCounts/) to generate read count.

ADD REPLYlink written 4 weeks ago by waqaskhokhar99950
1

OP says

I want to do DeSeq2

Your advice is to

use featurecounts

I don't see how this answers OP's question.

ADD REPLYlink modified 4 weeks ago • written 4 weeks ago by RamRS22k

Thank you for the advise. I also read the manual for stringtie and learnt the -A option but my question is,how to get a single gene count matrix for DESeq 2 . I am unable to do that.

ADD REPLYlink written 4 weeks ago by suneetamodekurty770

Do not add answers unless you're answering the top level question. This content belongs as a reply to waqaskhokhar999 's answer. I'm moving it to a comment on the top-level post now, but please be more careful in the future.

ADD REPLYlink written 4 weeks ago by RamRS22k

thank you will remember the advise.

ADD REPLYlink written 4 weeks ago by suneetamodekurty770
2
gravatar for kristoffer.vittingseerup
4 weeks ago by
European Union
kristoffer.vittingseerup2.0k wrote:

You can use the R package tximport specifically look at this section of the vignette (which also explain how to continue with DESeq2 analysis a bit further down). An alternative is the importIsoformExpression() function from my R package IsoformSwitchAnalyzeR which will enable you to find and analyze isoform switches (both in individual genes and genome wide systematic changes).

ADD COMMENTlink written 4 weeks ago by kristoffer.vittingseerup2.0k

Got it! Thank you!

ADD REPLYlink written 4 weeks ago by suneetamodekurty770

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one answer if they all work.

Upvote|Bookmark|Accept

ADD REPLYlink written 4 weeks ago by RamRS22k
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: 1463 users visited in the last hour