Entering edit mode
6.8 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
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)
?causes error.
raw.data.m was a typo sorry about that
What is the output of
head(raw.data.m)
?Do you have NAs or Infs in your matrix? How did you normalize your data?
I do not have NA. I used quantile normalization using Bioconductor.