Biostar Beta. Not for public use.
Error while using EdgeR
0
Entering edit mode
15 months ago
KVC_bioinfo • 390
Boston

Hello All,

I am trying to do differential expression using voom package. I constantly get the following error:

> Error in .isAllZero(counts) :    count matrix must be integer or
> double-precision

Below is the code I am trying:

library(limma)
library(preprocessCore)
library(edgeR)
library(tidyverse)

raw.data <- read.table(file="input.txt", header=TRUE)

tail(raw.data)
sample_info <- data.frame(condition = gsub( "_[0-93]+", "",names(raw.data)
),row.names = names(raw.data))
#sample_info

sample_info.edger <- factor(c( rep("sample", 84) , rep("Control", 8)))
sample_info.edger <- relevel( sample_info.edger , ref = "Control")

edgeR.DGElist <- DGEList(counts = raw.data , group = sample_info.edger )

Could someone help me here?

TIA

edgeR • 866 views
ADD COMMENTlink
0
Entering edit mode

Which step exactly causes the error? Is it edgeR.DGElist <- DGEList(counts = raw.data.m , group = sample_info.edger )?

What is the output of head(raw.data.m)?

ADD REPLYlink
0
Entering edit mode
edgeR.DGElist <- DGEList(counts = raw.data.m , group = sample_info.edger )

causes error.

raw.data.m was a typo sorry about that

ADD REPLYlink
0
Entering edit mode

What is the output of head(raw.data.m)?

ADD REPLYlink
0
Entering edit mode

Do you have NAs or Infs in your matrix? How did you normalize your data?

ADD REPLYlink
0
Entering edit mode

I do not have NA. I used quantile normalization using Bioconductor.

ADD REPLYlink
0
Entering edit mode
ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3