Greetings,
I'm new to BRMS but interested in tweaking the code underlying MBA to run some type of Connectivity-behavior analysis. The idea would be to look at how for example changes in connectivity between two conditions relates to changes in some behavioral measure.
Assuming the following default code:
brm(deltaY ~ 1 + (1|Subj) + (1|ROI1:ROI2) + (1|mm(ROI1, ROI2, weights = cbind(w, w), scale=FALSE)) + (1|mm(ROI1:Subj, ROI2:Subj, weights = cbind(w, w), scale=FALSE)) ...
one can control for covariates by running I believe:
brm(deltaY ~ 1 + Covar + (1|Subj) + (1|ROI1:ROI2) + (1|mm(ROI1, ROI2, weights = cbind(w, w), scale=FALSE)) + (1|mm(ROI1:Subj, ROI2:Subj, weights = cbind(w, w), scale=FALSE)) ...
however if purposefully want to investigate behavioral correlates should I run:
brm(deltaY ~ 1 + (1|Subj) + (1|ROI1:ROI2:**deltaBehavior**) + (1|mm(ROI1, ROI2, weights = cbind(w, w), scale=FALSE)) + (1|mm(ROI1:Subj, ROI2:Subj, weights = cbind(w, w), scale=FALSE)) ...
if so, would I then just calculate the Average Marginal Effect of deltaBehavior predicting deltaY at each ROI1:RO2 pairing?