Difference between revisions of "2009 Summer Project Week Statistical Toolbox"

From NAMIC Wiki
Jump to: navigation, search
Line 24: Line 24:
  
 
<h3>Approach, Plan</h3>
 
<h3>Approach, Plan</h3>
Due to the complexity of this toolbox the plan is to break it down in several releases or milestones.
+
Due to the complexity of this toolbox the plan is to break it down in several releases or milestones.  
  
<h4>First Milestone</h4>
+
In general terms I plan to use:
 +
* TinyXML for the XML processing
 +
* VNL (the one available in ITK) for matrices processing
 +
* Boost shared pointers for smart pointers
 +
* Some undetermined thread library for asynchronous processing (linear regression stuff...)
 +
 
 +
Any suggestion to change these libraries is more than welcome.
 +
 
 +
{{First Milestone}}
 
* The images should be spatially normalized and registered to an atlas
 
* The images should be spatially normalized and registered to an atlas
 
* The experiments are read/written using XML
 
* The experiments are read/written using XML
Line 36: Line 44:
 
* The integration with Slicer is resolved
 
* The integration with Slicer is resolved
  
<h4>Second Milestone</h4>
+
{{Second Milestone}}
 
* Summary statistics are available (Global/local)
 
* Summary statistics are available (Global/local)
 
* Extent threshold is available
 
* Extent threshold is available
Line 43: Line 51:
 
* Likelihood Ratio Test is available
 
* Likelihood Ratio Test is available
  
<h4>Third Milestone</h4>
+
{{Third Milestone}}
 
* Precalculated T-maps and F-maps can be compared
 
* Precalculated T-maps and F-maps can be compared
 
* Global and local statistic plots are available
 
* Global and local statistic plots are available
 
* Confidence intervals for the mean are available
 
* Confidence intervals for the mean are available
  
<h4>Fourth Milestone</h4>
+
{{Fourth Milestone}}
 
* Principal Component Analysis is available
 
* Principal Component Analysis is available
 
</div>
 
</div>
Line 62: Line 70:
  
 
==References==
 
==References==
[Statistical Parametrical Mapping]
+
* http://www.fil.ion.ucl.ac.uk/spm/
 
+
* http://www.math.mcgill.ca/keith/surfstat/
 
</div>
 
</div>

Revision as of 17:19, 4 June 2009

Home < 2009 Summer Project Week Statistical Toolbox


Key Investigators

  • Robarts Research/UWO: Diego Cantor


Objective

The objective of this project is to build a statistical toolbox for brain images, and integrate this toolbox into Slicer. Current tools such as SurfStat and SPM allow to run statistical experiments in two levels. The first level corresponds to a single subject study where the progression of a disease, the effect of a drug or the functional activation (fMRI) is analyzed over time. The second level corresponds to cross-sectional studies where multiple subjects are compared using designs such as ANOVA and T-tests in all its varieties.

The proposal is then to construct a third level where single-subject or multiple-subject experiments can be carried out using information from multiple modalities simultaneously. In such environment, every subject has a set of images(anatomical MR, functional MR, PET data, etc...) and the statistics are taken on vectors. Each vector is composed by the respective intensities in the image set.


Approach, Plan

Due to the complexity of this toolbox the plan is to break it down in several releases or milestones.

In general terms I plan to use:

  • TinyXML for the XML processing
  • VNL (the one available in ITK) for matrices processing
  • Boost shared pointers for smart pointers
  • Some undetermined thread library for asynchronous processing (linear regression stuff...)

Any suggestion to change these libraries is more than welcome.

Template:First Milestone

  • The images should be spatially normalized and registered to an atlas
  • The experiments are read/written using XML
  • Trilinear interpolation is used when the images in the set have different resolutions
  • There are not covariates
  • The only test available is a T^2 Hotelling test comparing all subjects against the mean
  • Uncorrected p-values are available
  • When the dimension is 1. The T^2 test is reduced to a regular t-test and this is verifiable
  • The integration with Slicer is resolved

Template:Second Milestone

  • Summary statistics are available (Global/local)
  • Extent threshold is available
  • Family-wise p-values are available
  • Two-sample T^2 test is available
  • Likelihood Ratio Test is available

Template:Third Milestone

  • Precalculated T-maps and F-maps can be compared
  • Global and local statistic plots are available
  • Confidence intervals for the mean are available

Template:Fourth Milestone

  • Principal Component Analysis is available

Progress

This is a long-term project. Right now I am working on the first milestone. I expect to learn about the Slicer architecture (modules, gui, event handling, processing threads, etc..) during the programming week so the posterior releases are developed entirely into Slicer.