I am searching for correlations between symptoms severity and the parameters of neural activity in a disease model. As predictors, I have two (slightly collinear) measures of disease severity on each side of the body that may differ (sometimes substantially) on the right and left body sides. Also, in each brain hemisphere, I have several recordings (4-16) from the structure involved in pathology – and for each recording, we have several estimated parameters of activity. And I want to study how these parameters correlate with symptoms severity.
So, I am wondering what is the best way to analyze the data that would account for all these conditions. Now I am thinking about using mixed linear model function for my task (nlme package):
fit <- lme(parameterA ~ measure1 + measure2 + random=~1|animal/side)
Here I wish to account for a subset of parameters derived from the given hemisphere. All the parameters in the hemisphere may be seen as repeated measures. But two subsets of recordings from the right and left hemisphere may also be dependent as they come from one animal.
Is my approach correct? Do I need to account for 'measure1 * measure2' instead of 'measure1 + measure2'? Do the variables for random effect estimates (animal and side) appear in the right order? Some of the neural parameters I want to study also have collinearity – what is the best way to account for it? I am going to analyze them in separate models and then adjust the estimates for multiple comparisons. Am I wrong?
If you know a more appropriate way to analyze such data please let me know.