2012年2月13日 星期一

SPM8 分析流程備忘

1. Import function and T1 DICOM file, respectory, in order to convert to .img and .hdr format

2.Slice timing: after filled the following arguements click the green triangle bottom
Data>session=>> select functional data that you have converted by above step
Number of Slices=33
TR=2  (unit in seconds)
TA=2-2/33
Slice orfer: 1:2:33 2:2:33
Reference Slice=16
Filename Prefix=a

3.Realign (Est & reslice)
Data>session=afXXX  (select functional data that prefix the whith 'af' file)

4.Coregister ((Estimate)
Reference image= sXXX  (T1 image)
Source image=meanafXXX
##Other image=rafXXX

5.Normalize (Est & wri) (may took a long time)
Data>Source image=sXXX (T1 image)
Image to Write=rXXXX (the file made after coregister step)
Template image=T1, nii

6.Smooth
Image to smooth=wXXX

7.Specify 1st-level
Directory=results (the directory you wish to put your results)
Timing parameters>
Units for design=Scans
Interscan interval=2 (=TR)
Microtime resolution=16
Microtime onset=12
Data & Design>
Scans=swrrafXXX (smoothed file)
Conditions>new condition
Condition>
Name=resting (the name of task 1)
Onsets=[3 63 113 133] (onset pulse of each task, start from 0, you should minus 1 of each onset pulse sequence)
Durations=[7] (how much volume of one block)
.
.
.
add new condition...

8.Estimate
Select SPM.mat in results directionary

9. Results
Select SPM.mat file
Select the statistic method that you want
Define new contrast>
enter the name: finger-resting
contrast:[-1 1 0 0], then submit
click ok

Continous define new contrast...
enter the name: toe-resting
contrast:[-1 0 1 0]...
               [task1 task2 task3 task4]
... completed all contrast

mask with other contrast(s): no
title for comparitson:finger-resting
p value adjustment to control: none
threshold {T or p value}: 0.001
&extent thrshold {voxels}: 0

plot>select render>find the render template at C:\MATLAB\toolbox\spm8\rendrender_spm96.mat
style: old

沒有留言:

張貼留言