Hello Romy,
You should have a look at your data, maybe they have different number of volumes, ID are different, etc .
I ran the MIGP option with 180 subjects, with and without automatic estimation of components and it worked very good in both cases. I should say that I look carefully at my raw data and after preprocessing. I ran it on a MacOs 16 GB RAM and 4 cores, therefore, your memory won´t be a problem at all...
I don´t think this message helps you a lot, but at least you can be sure that the MIGP works perfectly.
Regards
--
Yacila Deza Araujo, M.Sc. Neuropsych.
PhD Student
Technische Universität Dresden
Faculty of Medicine Carl Gustav Carus
Department of Psychiatry & Psychotherapy
Section of Systems Neuroscience
Würzburger Straße 35
01187 Dresden
Germany
Phone: +49 (0) 351 463 - 42217
Fax: +49 (0) 351 463 - 42202
mailto: [log in to unmask]
-----Ursprüngliche Nachricht-----
Von: FSL - FMRIB's Software Library [mailto:[log in to unmask]] Im Auftrag von Romy Lorenz
Gesendet: Freitag, 20. Februar 2015 13:06
An: [log in to unmask]
Betreff: [FSL] probelms with new group-PCA step in medodic - MIGP
Hi everyone,
as suggested by Steve in his post from 26.01.2015 I was running the new group-PCA step in Melodic (MIGP) for a large dataset (HCP data).
I used the following command and tested it on a small (n=20) dataset - not HCP data though:
melodic -i /group/HCP/HCP_MIGP_50random_input.txt --migp --sep_vn -o /home/rlorenz/scratch/HCP_Analysis/REST1_LR_hp2000_clean_20dim.mipg.gica -v --dim=20 --nomask --nobet --report --bgimage=/share/apps/fsl/data/standard/MNI152_T1_2mm_brain.nii.gz --Oall
However, when trying it to run on 50 subjects (HCP resting state data), I am getting the following error message:
terminate called after throwing an instance of 'std::bad_alloc'
I am running this command on a head node of a cluster with 32 cores and half a TB of RAM. So RAM should actually not be an issue, right?
I used /usr/bin/time in order to get a detailed error message (see below).
Could you please help me figuring out what is going wrong?
Thanks in advance and best wishes,
Romy
------------------------------------------------------------------------------
detailed error message:
terminate called after throwing an instance of 'std::bad_alloc'
what(): std::bad_alloc
Command terminated by signal 6
Command being timed: "melodic -i /group/HCP/HCP_MIGP_50random_input.txt --migp --sep_vn -o /home/rlorenz/scratch/HCP_Analysis/REST1_LR_hp2000_clean_20dim.mipg.gica -v --dim=20 --nomask --nobet --report --bgimage=/share/apps/fsl/data/standard/MNI152_T1_2mm_brain.nii.gz --Oall"
User time (seconds): 14078.95
System time (seconds): 36.19
Percent of CPU this job got: 99%
Elapsed (wall clock) time (h:mm:ss or m:ss): 3:55:47
Average shared text size (kbytes): 0
Average unshared data size (kbytes): 0
Average stack size (kbytes): 0
Average total size (kbytes): 0
Maximum resident set size (kbytes): 137120224
Average resident set size (kbytes): 0
Major (requiring I/O) page faults: 0
Minor (reclaiming a frame) page faults: 2392164
Voluntary context switches: 148
Involuntary context switches: 1427916
Swaps: 0
File system inputs: 3864824
File system outputs: 88
Socket messages sent: 0
Socket messages received: 0
Signals delivered: 0
Page size (bytes): 4096
Exit status: 0
|