I have some single-cell data I want to analyze in DESeq2 with a model of interaction terms of treatment and time variables. See the image for the metadata of samples I have. The data comes from several hashed captures, each with 2 samples. I wanted to model the capture to account for batch effects, but I've tried several permutations and keep getting linear combinations of terms. Is there a way to model both the individual and the capture it came from? Note that the same 2 individuals are always paired in a capture.
My thought for the model was something like this:
~ 0 + cluster_labels + condition:ind.condition.ID + timepoint:cap.timepoint.ID + condition * timepoint
The model works if I don't include the timepoint:cap.timepoint.ID term. Cluster labels isn't shown in the metadata table, but it's the celltypes (CM, TFH, etc.). One capture failed, so there are 6 captures for the 'pre' timepoint and 5 for the 'post'.
0 answers
No answers yet.
Log in to answer this question.
-
Finding differences between individual patients in a time-series RNA-seq?
written by txema.heredia •
-
Bulk RNA-seq, seeing strong differences between groups at a timepoint post-infection, but no strong…
written by Derrik •
-
DESeq2 for differential microbial abundances
written by ymj •
-
DESeq2: "model matrix is not full rank"
written by bart •
-
What is the appropriate Limma design matrix for the interaction between sex and age?
written by lintek08 •
-
Design in DESeq2: measruing effects of treatment while adjusting for sample type
written by Ridha •
-
Use DESeqDataSetFromMatrix from a matrix of normalized read counts
written by celia.escher •
-
Design for paired Deseq2 analysis
written by Pin.Bioinf
-
GOexpress analysis, removing genes with a score of 0?
written by draper •
-
How to identify human samples with a specific mouse gene expression signature?
written by Avro
