Thank you for taking my question into consideration.
I have RNA-seq dataset, a time series experiment, 3 time points with 8 replicates each. I wanna check the variation in gene expression across time. But I want to take into consideration the surgeon effect (2 surgeons were performing the surgery on mice).
I will run LRT test but I'm not really sure of what the reduced model should be: reduced~surgeon or reduced~1, the variable of interest should be removed which is time.