# Tag Info

### Understanding DESeq2 design, contrast and results

The simplest manner is to not use a wald test, but rather an LRT with a reduced model lacking the factor of interest: ...
• 19.2k
Accepted

### Removing PCR duplicates in RNA-seq Analysis

For normal RNA-seq PCR duplicates are normally kept in, but the duplication rate can be used as a quality control: The higher the duplication rate, the lower the quality. For expression analysis, it ...
• 3,211

### What does an FDR value of 1 in RNA-seq mean?

FDR stands for False Discovery Rate. It is a statistic tool used in multiple hypothesis testing. As you may know, when you use a p-value cutoff (usually 0.05) for your experiments, it means strictly ...
• 191

### error in heatmap using R

You have scale="col" in the code. What you are plotting is the z-score calculated based on the value distribution per column. Try changing it to ...
• 395
Accepted

### Understanding DESeq2 design, contrast and results

It seems that the "combining factors" trick described in part 3.3 of DESeq2 current "vignette" (as of may 2017) under the title "Interaction" is a way to access to the desired contrasts. It seems ...
• 2,960

### differential gene expression complex design no replicates

Mike Love is right. If the response you are looking for is linear in terms of change per minute, the most productive approach is likely to be fitting a linear model. You might get something out of ...
• 3,211
Accepted

### differential gene expression complex design no replicates

Although it is not recommended to use no replicates, in the edgeR manual they give some advice on how to go on with no replicates design. See page 21 of user guide. You can e.g. estimate a BCV value. ...
• 3,541
Accepted

### Missing genes and normalisation of RSEM output using EBSeq

Yes, that blog post does represent just one guy's opinion (hi!) and it does date all the way back to 2014, which is, like, decades in genomics years. :-) By the way, there is quite a bit of literature ...
• 4,938

### Getting a "system is computationally singular" error in sleuth

This happen when the variables (strain +batch) create a design matrix like this: batch strain 1 1 # 1 1 # 1 2 2 2 3 3 4 3 ... 16 72 Which means that some of the ...
• 4,602
Accepted

### Searching for gene expression data by cell line

Try the Gene Expression Omnibus - it looks like they have some datasets.
Accepted

### Running differential expression analyses on count matrices with many zeroes

I don't think that the issue is the low counts, but rather the number of features without any real variance (the black dots at the bottom). So what the heck is the dispersion plot and why does one ...
• 19.2k
Accepted

### Is removing samples based on clustering for downstream analysis a right choice?

I would be very hesitant to blindly exclude those samples based on the clustering. Check to see if the clusters actually denote some sort of batch effect, since it's not like all of the TCGA datasets ...
• 19.2k

### Running differential expression analyses on count matrices with many zeroes

Remove low-count features in advance. This is standard for most tools including DESeq2 and edgeR (see section 2.6). This will keep you from testing a lot of features that cannot be differentially ...
Accepted

### Differential Expression With 2 Treatments

If you're interested in looking at the differences between two treatments then you'll end up wanting to do both a direct contrast as well as the individual comparisons to baseline. The direct ...
• 19.2k

### Comparison of gene set enrichment statistics

I assume that you are talking about the implementation of these methods in the limma package. Otherwise this answer does not apply. I think that your questions can be answered with some simulations ...
• 4,602
Accepted

### Performing t.test to obtain p.values for each gene

You can simply apply() t.test() to your matrix. In general, though, I expect lmFit() from ...
• 19.2k
Accepted

### The biological meaning of the random variables and the responses in Seurat analysis

In the linked article the authors formalize microarray analysis as the study of the joint distributions of $\overrightarrow{X}_i$ and $Y_i$, where $\overrightarrow{X}_i$ is a vector of random ...
• 3,211

### How can I specify to DEseq2 to only perform comparisons on which I am interested with?

You can specify the exact comparisons you want in the results() function. So: ...
• 19.2k

### Removing PCR duplicates in RNA-seq Analysis

Generally you should just leave them as is. One does remove/mark duplicates in DNA seq. For further read check this Nature paper
• 51

### Missing genes and normalisation of RSEM output using EBSeq

Include all genes/transcripts in your analysis. A transcript that is not detected could be undetected through sampling error (i.e. the sequencer / library prep just happened to miss that transcript), ...
• 11.7k

### Differential Expression With 2 Treatments

Depending upon the type of treatment used the set of DEGs will change.If the treatments have similar kind of effect you will get a small list (less variable genes will have higher p-val) using a ...
• 584
Accepted

### Standard Cutoff for Moderated T-statistics

You're misinterpreting the moderated T-statistic, it's basically the fold-change divided by its variance. The p-value comes directly from that, so if you filter by moderated fold-change you're just ...
• 19.2k
Accepted

### Error in creating a volcano plot in MATLAB

I'm not sure if nowadays matlab is still a good option for informatics. I would rather go for a R script to make a volcano plot. If you don't want to have fancy code but functional one, you can even ...
• 56
Accepted

### Smallest group size for differential expression in limma (bulk RNA-Seq)

The problem with most of the methods is that use the gene's variance for each group, which can't be calculated (reliable) when the sample is <= 2. Also, statistically, it would have extremely low ...
• 4,602

### Smallest group size for differential expression in limma (bulk RNA-Seq)

While I share Llopis' concern about estimating variance from 2 samples, the statement you quoted is about avoiding false positives from genes that are only expressed in a few samples. It's fairly ...
• 1,806
Accepted

### Detecting differentially expressed genes with foldchange >= 2 and FDR < 0.05

Fold-change >= 2 is the same as logFC (log2(fold-change)) >= 1, so your example is doing exactly what you want. logFC is generally easier to think in and work with ...
• 19.2k

### Detecting differentially expressed genes with foldchange >= 2 and FDR < 0.05

The edgeR authors recommend that you use a relatively low logFC threshold for glmTreat such as lfc=log2(1.2). A ...
• 196
Accepted

### Selection of differential expressed genes

You can usually get away with FDR < 0.1, but that's as high as you can go. This all presumes you're doing follow-up experiments of some sort, of course. I guess you could increase the FDR more, but ...
• 19.2k

### How to normalise scRNASeq data for differential expression analysis

I would suggest using a likelihood ratio test for differential expression using logistic regression with batch as a latent variable. In Seurat you can do: ...
• 674