I am analyzing data from animal experiments. Each experiment includes a response variable (continuous), group {A, B}, condition(X, Y, Z}. The same experiment is replicated multiple times (5 times in this dataset) to ensure reproducibility. I want the model to account for the variability between replicates. The experiment-specific estimates are NOT my interest. I am using it as a blocking variable. I initially used it as a random effect to account for the fact that group effects can vary across different experiments. Then, I included the experiment as a fixed effect in addition to the random effect, which improves the model tremendously (the covariance parameter estimates were reduced from 20.6 to 5.8, And the AIC for the model was reduced from 703 to 676. Consequently, the confidence intervals for LS-means estimates were tighter when the experiment was included as a fixed effect as well as a random effect. My question: Is it okay to include the 'experiment' variable in the model this way? Here is my code: proc mixed data=lumi maxiter=1000 covtest; class group(ref="flox") experiment condition (ref="Basal"); model value = group|condition experiment/solution outpred=check; random group/ subject=experiment g v vcorr; lsmeans group*condition/cl slice=condition pdiff=control('flox' 'Basal') adjust=dunnett; run; Any insight would be much appreciated!
... View more