Moderator: i.sudbery

gravatar for i.sudbery
i.sudbery4.0k
Reputation:
4,050
Status:
Trusted
Location:
Sheffield, UK
Twitter:
IanSudbery
Last seen:
2 hours ago
Joined:
4 years ago
Email:
i********@sheffield.ac.uk

Posts by i.sudbery

<prev • 371 results • page 1 of 38 • next >
0
votes
4
answers
341
views
4
answers
Comment: C: Finding the best model for expression data (RNA-Seq counts)
... Depending on which normalisation you've used (VST or rlog), I was under the impression that there should be no zeros, and thus no need for +1. In fact rlog should is already adding be a pseudocount and taking the log. 1,000,000 is chosen because it it the value that puts the numbers and a "human sc ...
written 2 hours ago by i.sudbery4.0k
0
votes
4
answers
341
views
4
answers
Comment: C: Finding the best model for expression data (RNA-Seq counts)
... I would either fit a negative binomial, using gene length as an offset (not a co-variate, you need to force the coefficient to 1). Unfortunately getting the gene length is not entirely trival as you need to know which transcripts are expressed. If you can run the original data through salmon, that w ...
written 3 days ago by i.sudbery4.0k
0
votes
1
answer
85
views
1
answers
Answer: A: UMI Tools Dedup
... If you are doing drop-seq or 10x chromium, I highly recommend [`alevin`][1]. Other tools that can handle UMI deduplication are STARsolo, `umis` and `picard MarkDuplicatesWithCigar`. Not that the last two do not do error correction. Two things might cause excessive memory usage: 1. Many reads whos ...
written 3 days ago by i.sudbery4.0k
2
votes
4
answers
341
views
4
answers
Answer: A: Finding the best model for expression data (RNA-Seq counts)
... If I understand your question correctly you are trying to test whether genes that are more highly conserved tend to be more highly (or lowly) expressed. Thus would want to do what is conceptually similar to a linear correlation with conservation on one access and expression in another? Comparing ...
written 4 days ago by i.sudbery4.0k
0
votes
1
answer
106
views
1
answers
Comment: C: Htseq-count result assistance
... And where did you get the genome reference (fasta) from. The one you used to align the reads? ...
written 5 days ago by i.sudbery4.0k
1
vote
0
answers
143
views
0
answers
Comment: C: col as names
... I'm not sure what you mean by "the colnames of "my.names" are "3, 2, 3, 3, 3, 3...", with a row with the actual names I want" but read.table will not accept colnames that start with a number. Only column with names that are valid R variable names are allowed. As you have seen though, you can force ...
written 6 days ago by i.sudbery4.0k
0
votes
0
answers
71
views
0
answers
Comment: C: Read In Function for BAM file processing
... can you show us an outline of what this might look like for a single file? ...
written 6 days ago by i.sudbery4.0k
1
vote
1
answer
106
views
1
answers
Comment: C: Htseq-count result assistance
... Its possible nothing is wrong. You really need the full numbers. Depending on your GTF, most genes will have a count of zero, even in a good sample.... but how many can still be a mesaure of how good the sample is. If you are using the full ensembl GTF, which, I believe has about 50,000 genes in it, ...
written 6 days ago by i.sudbery4.0k
0
votes
1
answer
106
views
1
answers
Comment: C: Htseq-count result assistance
... Hi @citynsukka I edit your post to make the list of counts format correctly. You can do this by highlighting the part of the post to be formatted as-is and pressing the code button (labelled 101010). I've also changed the category - the tool category is used to announce/promote new tools. ...
written 6 days ago by i.sudbery4.0k
1
vote
2
answers
80
views
2
answers
Comment: C: help with MACS2 installation errors
... I think `pip --user` puts the binaries under `~/.local/bin` rather than `$HOME/bin`, so try the following in your `.bashrc` export PATH=$HOME/.local/bin:$PATH ...
written 6 days ago by i.sudbery4.0k

Latest awards to i.sudbery

Teacher 14 days ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Scholar 14 days ago, created an answer that has been accepted. For A: How to only change/substitute every QNAME(Read ID) from a Bam file?
Good Answer 20 days ago, created an answer that was upvoted at least 5 times. For A: Alternative to Cufflinks ?
Teacher 20 days ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Scholar 27 days ago, created an answer that has been accepted. For A: How to only change/substitute every QNAME(Read ID) from a Bam file?
Teacher 27 days ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Appreciated 27 days ago, created a post with more than 5 votes. For A: Alternative to Cufflinks ?
Commentator 4 weeks ago, created a comment with at least 3 up-votes. For C: if I am right about miRNA DE workflow
Commentator 4 weeks ago, created a comment with at least 3 up-votes. For C: if I am right about miRNA DE workflow
Teacher 5 weeks ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Commentator 6 weeks ago, created a comment with at least 3 up-votes. For C: if I am right about miRNA DE workflow
Scholar 6 weeks ago, created an answer that has been accepted. For A: How to only change/substitute every QNAME(Read ID) from a Bam file?
Appreciated 6 weeks ago, created a post with more than 5 votes. For A: Alternative to Cufflinks ?
Teacher 6 weeks ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Teacher 6 weeks ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Teacher 7 weeks ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Scholar 7 weeks ago, created an answer that has been accepted. For A: How to only change/substitute every QNAME(Read ID) from a Bam file?
Scholar 7 weeks ago, created an answer that has been accepted. For A: How to only change/substitute every QNAME(Read ID) from a Bam file?
Good Answer 8 weeks ago, created an answer that was upvoted at least 5 times. For A: Alternative to Cufflinks ?
Appreciated 8 weeks ago, created a post with more than 5 votes. For A: Alternative to Cufflinks ?
Teacher 8 weeks ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Teacher 8 weeks ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Teacher 8 weeks ago, created an answer with at least 3 up-votes. For A: Filtering during RNA Seq analysis?
Scholar 8 weeks ago, created an answer that has been accepted. For A: How to only change/substitute every QNAME(Read ID) from a Bam file?
Scholar 9 weeks ago, created an answer that has been accepted. For A: How to only change/substitute every QNAME(Read ID) from a Bam file?

Help
Access

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