two-way Reapted mesure ANCOVA in R

Question

I am trying to run a two-way repeated measure ANCOVA (mix design) in R. My data set looks somethign like this:

Sub <- c(1,1,1,1,2,2,2,2,3,3,3,3,...)
Medici