Entering edit mode
z.biranvand1983
•
0
@zbiranvand1983-7126
Last seen 9.0 years ago
Hi so I'm trying to do limma fitting on my data but i keep getting the following error
> library(limma)
> data<-read.table("C:/Users/test/Desktop/book1.txt",header = T, row.names = 1)
> dim(data)
[1] 720 14
> pairs(data)
> pairs(data)# it is better to use log2
> design <- model.matrix(~ 0+factor(c(rep(1,6),rep(2,8))))
> colnames(design) <- c("treat1","control1")
> #fit a linear model
> fit <- lmFit(data, design)
Error in rowMeans(y$exprs, na.rm = TRUE) : 'x' must be numeric
