SVA paired samples
1
0
Entering edit mode
whyw6948 • 0
@whyw6948-23037
Last seen 6 months ago
United States

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 • 1.1k views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 2 days 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.

ADD COMMENT

Login before adding your answer.

Traffic: 365 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