Home > brede > brede_mat_mrtm.m

brede_mat_mrtm

PURPOSE ^

brede_mat_mrtm - Ichise's MRTM analysis

SYNOPSIS ^

function Mout = brede_mat_mrtm(M, Mref, varargin)

DESCRIPTION ^

 brede_mat_mrtm       - Ichise's MRTM analysis

       BP = brede_mat_mrtm(M, Mref, varargin)

       Input:    M       'Mat' structure or numerical matrix
                 Mref    'Mat' strutcure or numerical matrix with
                         reference value/plasma curve

       Property: Info    [ {0} | value ] Amount of feedback/debug
                         information 
                 k2prime [ {[]} | Inf | positive value ] k'2 paramere
                 TimingType [ end | {middle] ] Determines where the
                         integration (cumulation) of the values in the
                         matrix should stop when the 'Times' are given.
                 Times   'Mat' structure with times: the first column 
                         should be the starting times of the scan
                         while the second column should be the end
                         time. There should be as many rows as there
                         are rows in the M and Mref matrices.

       Output:   Mout    Kinetic modeling parameters in a 'mat'
                         structure 

       Ichise's multilinear reference tissue model (MRTM/MRTM1/MRTM2)
       analyses are types of kinetic modelings for dynamic PET
       scans. This function will estimate the binding potential as well
       as other parameters of the kinetic model.

       The function expects an input where the activity curve for each
       region/voxel is in a single column and each row correspond to a
       time point/period. The outputted binding potential will be in a
       single row. Other rows are R1, k2 and k'2.

       The normal equation is used for estimation of the parameters,
       i.e., least squares with no correction for different amount of
       noise in the different scans. All time points are used for
       estimation. 

       The function may not estimate all voxels in a volume well. A
       spatial smoothing with a median filter as implemented in the
       brede_vol_smooth function may remove the badly estimated voxels. 

       Note that this function is not properly validated.

       Reference:
         Ichise et al., Linearized Reference Tissue Parametric Imaging
         Methods: Application to [11C]DASB Positron Emission Tomography
         Studies of the Serotonin Transporter in Human Brain, Journal of
         Cerebral Blood Flow & Metabolism, 23:1096+, 2003.

         PMOD, Ichise's Multilinear Reference Tissue Model (MRTM)
         http://www.pmod.com/technologies/doc/pkin/2316.htm

       Example:
         % "0.066 per minute" is from Kim (2006) J. Nucl. Med., 47:208+
         R = brede_mat_mrtm(M, Mref, 'Times', Mtimes, 'k2prime', 0.066/60);

       See also BREDE, BREDE_MAT, BREDE_MAT_CUMSUM, BREDE_MAT_LOGANPLOT,
                BREDE_VOL_SMOOTH.

 $Id: brede_mat_mrtm.m,v 1.5 2008/06/10 17:11:19 fn Exp $

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:
Generated on Fri 27-Nov-2009 18:11:22 by m2html © 2005