This is a test version of Biostars. For the public version, visit https://www.biostars.org.
edgeR glmRT code error message

Hi everyone, I did edgeR to analyze no replicated sample with the following code. Unfortunately, I met some of the problems that I can't fix.

count <- read.table("D:/NDM-12-4/test/Vero.CH.CH-C100.normalisedReadCounts.txt", sep = "\t", header = T, row.names = 1)
view(count)
d1 <- estimateGLMCommonDisp(count, design=NULL, verbose=TRUE)
d2 <- estimateCommonDisp(count, verbose=TRUE)
fit <- glmFit(count,dispersion = 0.26664)
results <- glmLRT(fit, coef=2)

The result shows that

Error in glmLRT(fit, coef = 2) : 
Need at least two columns for design, usually the first is the intercept column

But actually, I put two samples inside the count file. So what does that mean??

r gmlrt edger
• 2,000 views • link
updated by Mark Mark • written by jimmy0958073736 • jimmy0958073736

Please do not post links to Google Drive files and not even explain why you are doing this.

• 0 views • link
written by Kevin Blighe Kevin Blighe

1 answer

You need to create a design matrix. Why did you set it to NULL?

integer or character vector indicating which coefficients of the linear model are to be tested equal to zero. Values must be columns or column names of design

I highly suggest you read the edger manual it is fantastic.

• 0 views • link
written by Mark Mark

Dear Mark,

I have the same issue even providing the "designmatrix". I have only one sample for each condition and use the approach below:

designMat <- matrix(1,2,1)
dgList <- estimateGLMCommonDisp(dgList, design=designMat, method="deviance", robust = TRUE, subset=NULL)

fit <- glmFit(dgList)
lrt <- glmLRT(fit)

edgeR_result <- topTags(lrt, n=Inf)

I get "Need at least two columns for design, usually the first is the intercept column" error in lrt <- glmLRT(fit) step. I have tried a lot and do not know how to fix it - could you please help me with this? Thank you.

• 0 views • link
written by Apex92 • Apex92

Like I said above please read the manual. It's impossible to do DE with 1 sample in each condition. You need at minimum 2 samples per condition.

Please read the manual.

• 0 views • link
written by Mark Mark

Log in to answer this question.

More posts like this
  • RUVg correction
    written by Shaimaa Gamal • Shaimaa Gamal
    ![RUVg correction][1] [1]: /media/images/32d6c138-f261-4b94-aff7-34ca87e6 How can I improve RUVg correction? The RLE plot isn't still perfect ```r design = model.matrix(~factor, data = pData(set1)) design=as.matrix(design) y …
  • How to improve batch correction
    written by Shaimaa Gamal • Shaimaa Gamal
    RUVG correction: ![image: RUVG correction][1] [1]: /media/images/cb6d0e9c-8b35-4871-9080-46bb4fe4 I am not satisfied with the RLE plot and I wonder if there are anything to improve it? …
  • edgeR error
    written by rupali • rupali
    finding error to run edgeR , please check my code to be helpful for finding error and solving it. difficulty in finding the next steps …
  • Error when using edgeR with only two samples
    written by Apex92 • Apex92
    Dear all, I have only two samples (control vs KD) from an RNA-Seq experiment and want to see which genes are differentially expressed in the …
  • How to model a pre and post treatment in limma while adjusting for control
    written by AST • AST
    I have RNA-seq data from the GEO database that has three groups - Control (C), Disease1 (D1), and Disease2 (D2). The patients of D1 and …
  • Large sample size in edger
    written by Peter • Peter
    Hello, I have a spreadsheet that I intend to analyze using the edgeR package. Usually the number of samples is small, so I use a …
  • Differential expression of RNA seq data using EdgeR
    written by GuriGuri1565 • GuriGuri1565
    Hi everyone, I'm using edgeR for differential expression genes analysis. I make 4 groups that get differential expression genes each and then compare the number …
  • RNA-seq, EdgeR and batch effects
    written by Ekarl2 Ekarl2
    My experimental design consists of two experimental factors, four conditions (+/+, -/+, +/-, -/-) and four biological replicate per condition. The batch effect in this …
  • edgeR / deSeq parallel usage
    written by russhh russhh
    Hi. I'm using [edgeR][1] for differential expression analysis of RNASeq data (counts downloaded from TCGA). The workflow I'm using follows the GLM functionality of edgeR, …
  • EdgeR-GLMs correcting for batch effects
    written by Matina Matina
    Hi all, I have RNA-seq data for 3 conditions (control, breast, endometrial). Some of the samples where collected and sequenced in two different batches. I …