Biostar Beta. Not for public use.
Question: Can't make sva comparison despite very distinct batch effects
0
Entering edit mode

I continuously receive this error when I try to run sva with my count data

R code for sva:

# Loading all variables 
library(sva)
dat  <- dds_matrix
idx  <- rowMeans(dat) > 1
dat  <- dat[idx,]
mod  <- model.matrix(~batch, colData(dds))
mod0 <- model.matrix(~1, colData(dds))
svseq <- svaseq(dat, mod, mod0)

I receive this error message.

Error in solve.default(t(mod) %*% mod) : Lapack routine dgesv: system is exactly singular: U[4,4] = 0

I have very distinct batches which can be seen in this plot.

PCA plot of my data set, the colors indicate my batches

I have clearly dictated the batches as factors and this exact process has worked before. This is with a new, but very similar data set.

If it is simply not possible to remove these batch effects, how would I go about bringing these two batches closer together, they are definitely the same condition and the variance is likely due to batch effects, not biological significance.

ADD COMMENTlink 2.0 years ago jjrin • 10
Entering edit mode
0

what is your contrast of interest? if you model sva on your known batch, it will look for covariates that are independet of it.

ADD REPLYlink 2.0 years ago
Martombo
♦ 2.4k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0