SVA paired samples
Entering edit mode
whyw6948 • 0
Last seen 3.5 years ago

Hello, I have EPIC methylation data. This data is paired: we have two samples for each individual subject (one at timepoint 1, the other at timepoint 2). Some of these patients have a condition, others are healthy. I'm interested in the change in methylation at timepoint 2 versus timepoint 1 comparing diseased vs healthy people. There are 82 paired subjects, 164 total samples. Each of the variables is recognized as a factor.

I am having issues with SVA. I have set up my mod and mod0 matrix to adjust for the effects of patient not attributed to time or condition. However, I am receiving an error and am not exactly sure what it means.

mod <- model.matrix(~0 + ID + Time:Condition,data = pd) #interested in average effect of condition within each timepoint
mod0 <- model.matrix(~0 + ID, data = pd) #adjust for effect of individual patient
SVA <- sva(dat,mod,mod0,vfilter=10000) ##vfilter for reducing time 
Error in solve.default(t(mod) %*% mod) : 
  Lapack routine dgesv: system is exactly singular: U[85,85] = 0
SVA EPIC methylation limma • 787 views
Entering edit mode
Last seen 10 hours ago
United States

The simple answer is that your mod design matrix is not full rank. You probably need to do something similar to section 3.5 in the edgeR User's Guide in order to make it full rank.


Login before adding your answer.

Traffic: 469 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6