feComputeEvidenceFromTracts.m

[francopestilli]

write a function that takes:

• fg of a set of tracts
• fg of whole brain tractography
• DWI nifti, bvecs, bvals
• coregistered anatomical t1 (for plotting)
• a whole brain white matter segmentation

The function should:

1. compute the minmax(X,Y,Z) of the coordinates of all tracts passed it.
2. restrict the white matter mask to the cube contained by those coordinates.
3. build endpoints roi from the tracts
4. clip the whole brain fiber group to this new white-matter volume.
5. find tracts similar to the ones passed in present inside the whole brain fiber group by looking at fibers that share both endpoints with the tracts fibers.
6. add the found tracts in 5. to the tracts passed in.
   1. combine the whole brain fiber group and the tracts keeping track of the indices of each individual tract inside the new fiber group.
   2. build a life model from the new fiber group.
   3. fit the model.
   4. test the evidence for each tract.

This can be thought as a general mechanisms for any test when building a whole brian connectome is not possible due to computational limits.