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 $

- brede_mat_cumsum brede_mat_cumsum - Cumulated sum
- brede_mat_mat2mat brede_mat_mat2mat - Convert or clean matrix
- brede_struct_guesstype brede_struct_guesstype - Try to guess type of structure

Generated on Fri 27-Nov-2009 18:11:22 by