Dump demographics

Here's a matlab script that dumps the age & sex of a list of subjects.

% Dumps sexes and ages of a list of subjects, extracted from the DICOM files in /mridata/cbu
% Also produces mean +/- sd of ages across group

subjs={'CBU090515','CBU090516','CBU090518','CBU090519','CBU090520','CBU090829','CBU090830','CBU090837','CBU090838','CBU090842'};
rawpath='/mridata/cbu';

age=[];
for subjind=1:length(subjs)
    fn=dir(fullfile(rawpath,[subjs{subjind} '*']));
    try
        daydir=dir(fullfile(rawpath,fn(1).name,'20*'));
        seriesdir=dir(fullfile(rawpath,fn(1).name,daydir(1).name,'Series*'));
        dicomfn=dir(fullfile(rawpath,fn(1).name,daydir(1).name,seriesdir(1).name,'*.dcm'));
        dicomfn=fullfile(rawpath,fn(1).name,daydir(1).name,seriesdir(1).name,dicomfn(1).name);
        fprintf('%s\t',subjs{subjind});
        H=spm_dicom_headers(dicomfn);
        fprintf('%s\t%s\n',H{1}.PatientsSex,H{1}.PatientsAge);
        age=[age str2num(H{1}.PatientsAge(1:end-1))];
    catch
            fprintf('Cannot find %s\n',subjs{subjind});
    end;
end;
    
fprintf('Age %f +/- %f\n',mean(age),std(age));