Error while using EdgeR
1
0
Entering edit mode
6.7 years ago
KVC_bioinfo ▴ 600

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 • 2.6k views
ADD COMMENT
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 REPLY
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 REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode
ADD COMMENT

Login before adding your answer.

Traffic: 1352 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6