External Email - Use Caution
Dear experts,
We run the following analysis in Freesurfer:
*Group* (2 levels) x *X1* (continuous IV) x *X2 *(continuous IV) + *Sex *(categorical covariate) + *Age* (continuous covariate)
We exported the vertex with the maximum statistic after applying FWE correction to replicate the analysis in R. We want to make sure about the structure of the model to later do a Bayesian analysis using the same structure.
We got significant results for most of our analyses but not for all of them, likely because we are not doing exactly the same thing.
We have fitted several models in R, but the one that gave the closest results was the following:
model <- lm (vertex ~ NewFactor (4 levels: Group*Sex) + X1X2 (resulting from X1*X2) * Group + Age*NewFactor, data=mydata)
How should we introduce the interaction term in order to make the models fitted by Freesurfer and R comparable?
Thanks a lot in advance,
Marina