Regression analysis

[Yuko]

Dear all,

I would like to perform a regression analysis using brain measures and behavior assessments. So, I’m wondering if the BRAPH can extract each brain measure from individual participants.

Any comments are highly appreciated.

Best wishes,
Yuko

[Mite Mijalkov]

Dear Yuko,

Yes, you can extract the brain measures from the individual participants by running a script in Matlab. Only thing that you would need for this script to work is to run Braph in Matlab (so that the relevant functions can be loaded) and save your graph analysis file in .fga format (in this case the file is named trial.fga). Then, you can check the following sample code that extracts the modularity for all patients in a given group:

clear all; close all; clc;
% this code will extract the results for all subjects in a given group at
% all densities specfied as a range. The code currently extracts modularity
% measure.
% the values are extracted in the array meas_extracted
% each row is a subject in the group, while each column represents
% different density

%% load the fga file with already calculated measures
load(‘trial.fga’,’-mat’)

%% specify groups for which measures can be extracted
% the groups can be checked by typing ga.cohort on the command line
group = 1;

%% specify density range
densities = 5:1:15;

%% specify the measure number
% the measure number is found by typing Graph.NAME on the command line
% 34 – modularity
measure_numbers = 34;

%% calculate the number of subjects in this group
sub_num = numel(find(ga.cohort.getGroup(group).getProp(Group.DATA)));

%% initialize the array that holds the measure
meas_extracted = zeros(sub_num,length(densities));

%% extract the measures
for d = 1:1:length(densities)

%% get all measures from the particular group
[ms_gr1,mi_gr1] = ga.getMeasures(measure_numbers,group);
% get only the values of the measure at the given density
value_gr1 = ms_gr1{densities(d)}.getProp(fMRIMeasureBUD.VALUES1);

%% put the measures at the corresponding place in the array
meas_extracted(:,d) = value_gr1;
end

Hope this helps. If you need some other specific analyses, I can try to help with that.

Best,
Mite

[Author]

Posts