AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

|
May 19, 2009 10:12AM
> Thanks for the response, Rick. I have 12 subjects, but I am looking
> at the distributions across voxels. 60,000 voxels should be plenty to
> get an even distribution that would allow 5% of the voxels to exceed
> the p < .05 threshold by random chance, would it not?

This is not quite correct. The t-distribution approaches the normal
distribution as the degrees of freedom approaches infinity (subjects
minus 1 in your case), but you have only 12 subjects. 60000 voxels is
multiple comparisons, but of the same non-normal test.

Keep in mind that these results depend on the variance of your data,
which comes not only from BOLD responses to your stimuli, but from how
the experiment was designed and how the data was processed.

---

To get a feel for the importance of the number of subjects, try this:

First, compute your test using 3dttest instead of 3dANOVA*, as it is
easier to manipulate for a single test. Just grab all the subject
volumes that went into the ANOVA and list them as '-set2' in a 3dttest
command. Run the 3dttest command. It should give the same results.

Now for comparison, input all of your subject datasets 4 times in the
otherwise same 3dttest command, so it is as if you had 4 times as many
subjects. Clearly the mean beta will be the same, but now the t-stats
should be basically doubled (as expected by the definition, the stat
will grow as sqrt(N)).

For a similar test, drop down to 3 subjects.

- rick

Subject Author Posted

2% of voxels reaching p < .05

cdm May 14, 2009 09:14AM

Re: 2% of voxels reaching p < .05

Gang Chen May 14, 2009 05:44PM

Re: 2% of voxels reaching p < .05

rick reynolds May 14, 2009 06:26PM

Re: 2% of voxels reaching p < .05

cdm May 19, 2009 09:04AM

Re: 2% of voxels reaching p < .05

rick reynolds May 19, 2009 10:12AM

Re: 2% of voxels reaching p < .05

cdm May 20, 2009 11:04AM

Re: 2% of voxels reaching p < .05

Bob Cox May 15, 2009 09:14AM

Re: 2% of voxels reaching p < .05

cdm May 15, 2009 09:32AM