https://github.com/philjrdb/HCP
Tip revision: 6df52a9f08fe8150b87b53f004c461ea768bd60f authored by philjrdb on 05 April 2021, 01:56:00 UTC
Update README.md
Update README.md
Tip revision: 6df52a9
Phil_HCPques.m
%% HumanCondPun questionnaires
ques_names = {'dass_anx' 'dass_dep' 'bis_sr' 'bis_imp' 'aor_rno' 'aor_rpo'...
'bisbas_drv' 'bisbas_rew' 'bisbas_fun' 'bisbas_bis'...
'ipip_e' 'ipip_a' 'ipip_c' 'ipip_n' 'ipip_i'}; % questionnaire names
ques{1} = [1 0 1 0 1 1 0 0 1 0 0 1 1];
ques{2} = [0 1 0 1 0 0 1 1 0 1 1 0 0];
ques{3} = [1 0 0 1 1 1 0 0];
ques{4} = [0 1 1 0 0 0 1 1];
ques{5} = [1 0 1 0 1 0 1 0 1 0 1 0];
ques{6} = [0 1 0 1 0 -1 0 -1 0 -1 0 1];
ques{7} = [1 0 0 0 1 0 1 0 0 0 0 0 0];
ques{8} = [0 1 1 0 0 0 0 0 0 0 0 1 0];
ques{9} = [0 0 0 0 0 1 0 0 1 0 1 0 0];
ques{10} = [0 0 0 1 0 0 0 1 0 1 0 0 1];
ques{11} = [1 0 0 0 0 -1 0 0 0 0 1 0 0 0 0 -1 0 0 0 0];
ques{12} = [0 1 0 0 0 0 -1 0 0 0 0 1 0 0 0 0 -1 0 0 0];
ques{13} = [0 0 1 0 0 0 0 -1 0 0 0 0 1 0 0 0 0 -1 0 0];
ques{14} = [0 0 0 1 0 0 0 0 -1 0 0 0 0 1 0 0 0 0 -1 0];
ques{15} = [0 0 0 0 1 0 0 0 0 -1 0 0 0 0 -1 0 0 0 0 -1];
a_leng = length(HCP_aggr);
%%
for q = 1:length(ques_names)
n_idx = regexp(ques_names{q},'_');
ques_field = ['ques_' ques_names{q}(1:n_idx-1)];
q_idx = find(ques{q} ~= 0);
q_scoring = ques{q}(q_idx);
for r = 1:a_leng
HCP_aggr(r).(ques_names{q}) = mean(HCP_aggr(r).(ques_field)(q_idx).*q_scoring);
end
end
clearvars -except HCP* work_folder *keep;