Hello,
I'm setting up an experiment where I measure changes in cell numbers in response to the concentration of 8 different cell culture media ingredients (i.e. continuous factors). A DSD with 8 extra runs gives me a 29 run experiment for this. Since the main time sink of this experiment is making up the different media mixes, I'm planning to run 3 exact replicates, preferrably blocked since the replicates will be in different 96-well plates, and differences in handling (e.g. time not spent in the incubator, etc) may lead to differential growth. All easy enough to set up by augmenting the DSD.
Run | Media Mix | Block | Response |
1 | 1 | 1 | ? |
... | ... | ... | ... |
28 | 28 | 1 | ? |
29 | 29 | 1 | ? |
30 | 1 | 2 | ? |
However, I'm also gathering data over several days (non-invasively measuring cell numbers) of these runs, up to 7 times. I'm expecting the difference between treatments to increase due to exponential growth of the cells up to a point, then decrease again as they run out of space. I could fit models to the data from the individual days and use the one with the best PValues, but it seems like a waste of data not to integrate it all into a single model.
However, the way I understand it, if I just used one blocking variable and put all replicate/day combinations into seperate blocks (see below), I would be aliasing the variability in response due to different days and different replicas.
Run | Media Mix | Day | Block | Response |
1 | 1 | 1 | 1 | ? |
... | ... | ... | ... | ... |
28 | 28 | 1 | 1 | ? |
29 | 29 | 1 | 1 | ? |
30 | 1 | 1 | 2 | ? |
... | ... | ... | ... | ... |
87 | 29 | 1 | 3 | ? |
88 | 1 | 2 | 4 | ? |
Is it possible to use two blocking variables to get around this problem?
Run | Media Mix | Day | Block A | Block B | Response |
1 | 1 | 1 | 1 | 1 | ? |
... | ... | ... | ... | ... | ... |
28 | 28 | 1 | 1 | 1 | ? |
29 | 29 | 1 | 1 | 1 | ? |
30 | 1 | 1 | 2 | 1 | ? |
... | ... | ... | ... | ... | ... |
87 | 29 | 1 | 3 | 1 | ? |
88 | 1 | 2 | 1 | 2 | ? |
I hope what I'm asking makes sense. I think I could alternatively put the day as a categorical factor, but I don't want any model to include it, so it feels unnecessary.