@@ -456,30 +456,36 @@ def _make_matlab_command(self, _):
456
456
contrasts [i ].weights = cont [3 ]
457
457
if len (cont ) >= 5 :
458
458
contrasts [i ].sessions = cont [4 ]
459
- script = ["""\
460
- % generated by nipype.interfaces.spm
459
+ script = [
460
+ """\
461
+ %% generated by nipype.interfaces.spm
461
462
spm_defaults;
462
463
jobs{1}.stats{1}.con.spmmat = {'%s'};
463
464
load(jobs{1}.stats{1}.con.spmmat{:});
464
465
SPM.swd = '%s';
465
466
save(jobs{1}.stats{1}.con.spmmat{:},'SPM');
466
467
names = SPM.xX.name;
467
- """ % (self .inputs .spm_mat_file , os .getcwd ())]
468
+ """
469
+ % (self .inputs .spm_mat_file , os .getcwd ())
470
+ ]
468
471
# get names for columns
469
472
if isdefined (self .inputs .group_contrast ) and self .inputs .group_contrast :
470
473
script += ["condnames=names;" ]
471
474
else :
472
475
if self .inputs .use_derivs :
473
476
script += [r"pat = 'Sn\([0-9]*\) (.*)';" ]
474
477
else :
475
- script += [r"""\
476
- pat = 'Sn\([0-9]*\) (.*)\*bf\(1\)|Sn\([0-9]*\) .*\*bf\([2-9]\)|Sn\([0-9]*\) (.*)';""" ]
478
+ script += [
479
+ r"""\
480
+ pat = 'Sn\([0-9]*\) (.*)\*bf\(1\)|Sn\([0-9]*\) .*\*bf\([2-9]\)|Sn\([0-9]*\) (.*)';"""
481
+ ]
477
482
478
483
script += ["t = regexp(names,pat,'tokens');" ]
479
484
# get sessidx for columns
480
485
script += [r"pat1 = 'Sn\(([0-9].*)\)\s.*';" ]
481
486
script += ["t1 = regexp(names,pat1,'tokens');" ]
482
- script += ["""\
487
+ script += [
488
+ """\
483
489
for i0=1:numel(t)
484
490
condnames{i0}='';
485
491
condsess(i0)=0;
@@ -488,7 +494,8 @@ def _make_matlab_command(self, _):
488
494
condsess(i0)=str2num(t1{i0}{1}{1});
489
495
end;
490
496
end;
491
- """ ]
497
+ """
498
+ ]
492
499
493
500
# BUILD CONTRAST SESSION STRUCTURE
494
501
for i , contrast in enumerate (contrasts ):
@@ -497,37 +504,44 @@ def _make_matlab_command(self, _):
497
504
script += ["consess{%d}.tcon.convec = zeros(1,numel(names));" % (i + 1 )]
498
505
for c0 , cond in enumerate (contrast .conditions ):
499
506
script += ["idx = strmatch('%s',condnames,'exact');" % cond ]
500
- script += ["""\
507
+ script += [
508
+ """\
501
509
if isempty(idx)
502
510
throw(MException('CondName:Chk', sprintf('Condition %%s not found in design','%s')));
503
511
end;
504
- """ % cond ]
512
+ """
513
+ % cond
514
+ ]
505
515
if contrast .sessions :
506
516
for sno , sw in enumerate (contrast .sessions ):
507
517
script += ["sidx = find(condsess(idx)==%d);" % (sno + 1 )]
508
- script += ["consess{%d}.tcon.convec(idx(sidx)) = %f;" % (
509
- i + 1 ,
510
- sw * contrast .weights [c0 ],
511
- ) ]
518
+ script += [
519
+ "consess{%d}.tcon.convec(idx(sidx)) = %f;"
520
+ % ( i + 1 , sw * contrast .weights [c0 ],)
521
+ ]
512
522
else :
513
- script += ["consess{%d}.tcon.convec(idx) = %f;" % (
514
- i + 1 ,
515
- contrast .weights [c0 ],
516
- ) ]
523
+ script += [
524
+ "consess{%d}.tcon.convec(idx) = %f;"
525
+ % ( i + 1 , contrast .weights [c0 ],)
526
+ ]
517
527
for i , contrast in enumerate (contrasts ):
518
528
if contrast .stat == "F" :
519
529
script += ["consess{%d}.fcon.name = '%s';" % (i + 1 , contrast .name )]
520
530
for cl0 , fcont in enumerate (contrast .conditions ):
521
531
tidx = cname .index (fcont [0 ])
522
- script += ["consess{%d}.fcon.convec{%d} = consess{%d}.tcon.convec;" %
523
- (i + 1 , cl0 + 1 , tidx + 1 )]
532
+ script += [
533
+ "consess{%d}.fcon.convec{%d} = consess{%d}.tcon.convec;"
534
+ % (i + 1 , cl0 + 1 , tidx + 1 )
535
+ ]
524
536
script += ["jobs{1}.stats{1}.con.consess = consess;" ]
525
- script += ["""\
537
+ script += [
538
+ """\
526
539
if strcmp(spm('ver'),'SPM8')
527
540
spm_jobman('initcfg');
528
541
jobs=spm_jobman('spm5tospm8',{jobs});
529
542
end;
530
- """ ]
543
+ """
544
+ ]
531
545
script += ["spm_jobman('run',jobs);" ]
532
546
return "\n " .join (script )
533
547
0 commit comments