## User: Jakub

Jakub30
Reputation:
30
Status:
New User
Location:
United Kingdom
Last seen:
1 year, 9 months ago
Joined:
2 years, 9 months ago
Email:
j***********@ndcn.ox.ac.uk

#### Posts by Jakub

<prev • 15 results • page 1 of 2 • next >
2
551
views
2
... 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 22 months ago by Jakub30
1
337
views
1
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 22 months ago by Jakub30
2
551
views
2
... 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 22 months ago by Jakub30
1
520
views
1
written 2.3 years ago by Jakub30
2
889
views
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 2.3 years ago by Jakub30
1
520
views
1
... 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 2.3 years ago by Jakub30
1
520
views
1
... 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 ...
written 2.3 years ago by Jakub30
1
283
views
1
... 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 2.4 years ago by Jakub30
2
889
views
2
... I didn't find an answer to this searching the forums. I have RNASeq samples with 5'/3' biases that are unevenly distributed amongst the samples. Some of my conditions have more samples with the bias, some less - the reason for these biases is almost certainly different levels of RNA sample fragmenta ...
written 2.4 years ago by Jakub30 • updated 2.4 years ago by Ryan C. Thompson6.8k
1
973
views
1
... I have thought about this very hard, and have run goseq in a variety of ways. Similar to what goseq package suggests: genes=as.integer(p.adjust(res$padj[(res$log2FoldChange!=0) & (!is.na(res$log2FoldChange))], method="BH")<.05) names(genes)= row.names(res[(res$log2FoldChange!=0) & ( ...
written 2.8 years ago by Jakub30

#### Latest awards to Jakub

Scholar 2.4 years ago, created an answer that has been accepted. For A: Issue with SampleSheet.csv during illumina methylatiom microarrays analysis by R

Content
Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.