I was not sure how to modify the last line in your code in order to calculate adjusted values accounting for age, sex, and ICV
for each participant (this is a three group design).
Your code:
There is not an easy way to do it, but you can do it in
matlab
cd glmdir
X = load('Xg.dat');
y = load('perm.th30.abs.y.ocn.dat');
beta = inv(X'*X)*X'*y;
yhat = X(:,[i j k ...]*beta([i j k ...])
where i j k ... indicate the factors you want keep in your data. Each
factor gets a column in the design matrix (and a beta). If you need help
figuring this out, send the fsgd file.