Thankyou for getting back to me on this . So, I wanted to put the full piece of code there because there is a sense of continuity.If I were to show you what the vsd_adjusted_log2 graph is , then it goes way back and the code would be longer.But I will update with the smallest amount of code so that you get a better understanding.
But this is what the vsd_adjusted_log2 contains.
vsd_adjusted_log2 <- removeBatchEffect(assay(vsd), batch=vsd$Time_point, covariates = vsd$Age)