I have an experimental design where we change the genotype and the treatments. The replicates were collected at different times, so I am assmuing there is a batch effect. I have two questions here -
To build the DESeq object, I am using the ~Batch + treatment + genotype + treatment:genotype model. DESeq2 automatically converts the treatment and genotype to factor variables, but not the batch effect. So, how would it matter if Batch was considered as a factor or as a continuous variable?
When I do the PCA after removing batch effect using limma's removeBatcheffect, the samples still cluster into distinct groups as we would expect it to. So, do we still need to model for the batch effect? Is there anything else that determines if the batch effect would change the further downstream analysis drastically?
Thank you so much!