Dear David,


For one-dimensional contrasts (e.g. t-tests) SPM asks you for two parameters for Bayesian inference at the second level (i) Effect Size Threshold (Default 0.1) and (ii) Log Odds Threshold (Default 10).


Other reasonable choices would be 0 and 3.


The effect size threshold, T, tells SPM that you are only interested in voxels with contrast values C^beta > T. ie. that your experimental effect is bigger than T.


The Log Odds Threshold, L, tells SPM that you are only interested in voxels where SPM is sure (with posterior probability 1/(1+exp(-L)) )

that this is the case.


Note that L=3 gives you p=0.95.

L=10 is much, much more stringent giving p=0.99995.


I would advise you use the most recent version of SPM when doing this.


Also, you don't have to do a first level Bayesian analysis if you want to a second-level one.


All the best,


Will.



From: SPM (Statistical Parametric Mapping) <[log in to unmask]> on behalf of David Hofmann <[log in to unmask]>
Sent: 31 January 2017 10:52
To: [log in to unmask]
Subject: [SPM] How to run a (1st + 2nd level) Bayesian analysis in SPM
 
Hi all,

I have an fMRI event-related design in which subjects viewed fearful and neutral faces. I want to run a 1st level and a second level Bayesian analysis in SPM. For this, I did the following steps:

1. 1st level Bayesian analysis with standard settings as described in the manual
2. Contrast fear faces > neutral faces
3. For the 2nd level analysis, I smoothed the con-files and ran a one-sample t-test (estimated the model first with the classical and then with the 2nd level Bayesian option) 
4. I specified a t-contrast (i.e. [1]) for the one-sample t-test of the subjects
5. I chose apply masking - none

Now SPM is asking me for the Effect size threshold for PPM at the 2nd level and suggests 0.99. Whereas the meaning of the effect size threshold was clearly explained in the manual for the 1st level analysis, I not sure what value to choose for the 2nd level analysis and what this value means.When I select the suggested value (0.99) and choose as Log Odds Threshold 10, which should correspond to 95 % certainty, then there is no effect. There are also no effects for a value as low as 0.2. This is very strange since in the classical analysis there are very strong effects (fusiform gyrus) which survive an FWE correction at 0.01.

The questions are as follows:

1. Are the analysis steps I did correct or is there a better way to test for group effect by means of Bayesian analysis (e.g. Bayesian model comparison, Rosa, M.J. et al., 2010)

2. What does the effect size threshold at the 2nd level mean and what are reasonable values?


Here is an overview of posts with topic Bayesian analysis, which did not help me answering my questions:
1.https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=spm;888fe64.1503
2.https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=SPM;41144d5.1403
3.https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=SPM;5f9a54e5.1405
4.https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=ind1603&L=spm&F=&S=&P=639757
5.https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=spm;2e6e6dca.1405
6.https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=spm;a5ab6e97.1603
7.https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=SPM;377114fa.0909


greetings

David