User: Jakub

gravatar for Jakub
Jakub30
Reputation:
30
Status:
New User
Location:
United Kingdom
Last seen:
1 month, 3 weeks ago
Joined:
4 years, 1 month ago
Email:
j***********@ndcn.ox.ac.uk

Posts by Jakub

<prev • 17 results • page 1 of 2 • next >
0
votes
1
answer
249
views
1
answer
merging RUV-corrected datasets in DESEq
... Dear Bioconductor Community, I have a workflow for multiple plates single-cell RNASeq that works as follows QC Use RUVSeq to eliminate library complexity bias from data (this works really well, with pre and post) Account for Zero Inflation using Zinbwave Model differential expression using DE ...
deseq2 ruvseq glm zinbwave written 11 months ago by Jakub30 • updated 11 months ago by Michael Love26k
1
vote
2
answers
1.5k
views
2
answers
Answer: A: Question: GoSeq: opposite pwf plot for up- and downregulated genes
... I am just going through some old posts. Just for completeness (since it is still being looked at by people who might look for an answer), the answer to this question was that the data was subject to 5'/3' bias in one of the experimental groups in particular, yielding to an overexpression of short g ...
written 13 months ago by Jakub30
0
votes
2
answers
771
views
2
answers
Comment: C: is it necessary to check batch effect in this case? and how to?
... Glad to help. I assume you used default trimmomatic settings, which are reasonable. The result looks fine. You will have to decide how much QC you do on 'lanes', and whether you want to know per lane GC content etc... Trimming is only one part of QC. There are many ways to explore biases from seque ...
written 3.1 years ago by Jakub30
0
votes
1
answer
549
views
1
answers
Answer: A: Change the FDR in DEXSeq
... You can subset the results class with any cuttoff you like, as the results object has a "padj" column, which contains the the adjusted p value (also known as the FDR adjusted p-value), by some multiple testing adjustment, here using the Benjamini-Hochberg method. At the end of chapter 4 of the vign ...
written 3.1 years ago by Jakub30
1
vote
2
answers
771
views
2
answers
Answer: A: is it necessary to check batch effect in this case? and how to?
... If I understand correctly, you have run the same 88 samples in each of the flow cell lanes? In that case that's an optimal experimental design and goes a long way in preventing the sort of batch effects you might be worried about. I would argue that an optimal block design is the main answer to your ...
written 3.1 years ago by Jakub30
0
votes
1
answer
742
views
1
answers
Comment: A: DESeq2 normFactors method: Error in fitBetaWrapper
... It does indeed, sorry for not being explicit about this! Thanks again. ...
written 3.6 years ago by Jakub30
0
votes
2
answers
1.3k
views
2
answers
Comment: C: Accounting for 5'/3' Bias in DESeq 2
... Thanks! I've done exactly this and computed TINs for each gene, and performed the loess regression. I now have the raw logcounts and corrected logcounts. I guess I am not clear in my head which value is best to use in a normFactor offset matrix, before normalising each row to a geometric mean of 1 ...
written 3.6 years ago by Jakub30
0
votes
1
answer
742
views
1
answers
Comment: C: DESeq2 normFactors method: Error in fitBetaWrapper
... Many thanks! Indeed, not all values were finite (lots of dividing by zero). I have added the line: normFactors[!is.finite(normFactors)] <- 1  A typical column now looks like this: Min. :0.1656 1st Qu.:0.9380 Median :1.0000 Mean :1.0208 3rd Qu.:1.1189 Max. :4.2724 ...
written 3.6 years ago by Jakub30
2
votes
1
answer
742
views
1
answer
DESeq2 normFactors method: Error in fitBetaWrapper
... Dear all, this is a follow up from my last question. When my run my DESeq2 analysis using the usual estimateSizeFactors, it runs fine, however if I provide a normalisation matrix, I get the following error on estimateDispersions: Error in fitBetaWrapper(ySEXP = counts(object), xSEXP = modelMatrix ...
normalization deseq2 written 3.6 years ago by Jakub30
2
votes
1
answer
283
views
1
answers
Answer: A: Issue with SampleSheet.csv during illumina methylatiom microarrays analysis by R
... Hi, I think there are indeed two issues. I would rename the columns to the standard column names (see below). I have added the Tissue column, but you can add any number of columns to the right the Sentrix_Position column. Sample_Name,Sample_Well,Sample_Plate,Sample_Group,Pool_ID,Sentrix_ID,Sentri ...
written 3.7 years ago by Jakub30

Latest awards to Jakub

Popular Question 3.1 years ago, created a question with more than 1,000 views. For Question: GoSeq: opposite pwf plot for up- and downregulated genes
Popular Question 3.1 years ago, created a question with more than 1,000 views. For Accounting for 5'/3' Bias in DESeq 2
Scholar 3.1 years ago, created an answer that has been accepted. For A: Issue with SampleSheet.csv during illumina methylatiom microarrays analysis by R
Scholar 3.7 years ago, created an answer that has been accepted. For A: Issue with SampleSheet.csv during illumina methylatiom microarrays analysis by R

Help
Access

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