calculating effect of condition taking into account patient origin and clone_origin
1
0
Entering edit mode
idetoma • 0
@646f29c0
Last seen 10 months ago
Spain

I have a complex design where I have 3 karyotypes. Each genotype has samples from different patients (patient_name in the metadata), and from each patient they derived different clones. I am interested in the effect of the karyotype (47vs46, and 49vs46), taking into account the patient origin and clone origin.

This is my metadata:

enter image description here

I try to create my model with a custom matrix by doing like this:

mod.matrix=model.matrix(~   karyotype:pat.n:ind.n + karyotype, metadata_NSC)
mod.matrix=mod.matrix[,which(colMeans(mod.matrix)!=0)]

dds_NSC <- DESeq(dds_NSC, betaPrior = FALSE, full=mod.matrix)

However I get a not full rank model error:

Error in checkFullRank(full) : 
  the model matrix is not full rank, so the model cannot be fit as specified.
  One or more variables or interaction terms in the design formula are linear
  combinations of the others and must be removed.

Am I using the nesting wrongly?

PS. The clones are unfrozen independently and therefore are not proper technical replicates (even if they derive from the same patient) Thanks

DESeq2 • 376 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 10 hours ago
United States

For questions about setting up the statistical design for your study, I recommend speaking with a local statistician or someone familiar with linear models in R. I have to restrict my time on the support site for software related questions.

ADD COMMENT

Login before adding your answer.

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