modell.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpcextensions.h"
/*****************************************************************************/
#include "tpcmodels.h"
/*****************************************************************************/
 
/*****************************************************************************/
typedef struct TPC_MODEL {
  char *code;
  unsigned int id;
  unsigned int type; 
  unsigned int parNr;
  char *desc;
} TPC_MODEL;
/*****************************************************************************/
 
/*****************************************************************************/
static TPC_MODEL tpc_model[]={
  {"unknown",    0,    0, 0,  "unknown"},
  /* mathematical functions */
  {"level",       100,  2, 1,  "f(x)=A"},
  {"line",        101,  2, 2,  "line"},
  {"pol2",        102,  2, 3,  "2nd degree polynomial"},
  {"pol3",        103,  2, 4,  "3rd degree polynomial"},
  {"pol4",        104,  2, 5,  "4th degree polynomial"},
  {"pol5",        105,  2, 6,  "5th degree polynomial"},
  {"pol6",        106,  2, 7,  "6th degree polynomial"},
  {"pol7",        107,  2, 8,  "7th degree polynomial"},
  {"pol8",        108,  2, 9,  "8th degree polynomial"},
  {"pol9",        109,  2, 10, "9th degree polynomial"},
  {"ratf11",      211,  2, 0,  "1st/1st order rational function"},
  {"ratf21",      221,  2, 0,  "2nd/1st order rational function"},
  {"ratf22",      222,  2, 0,  "2nd/2nd order rational function"},
  {"ratf32",      232,  2, 0,  "3rd/2nd order rational function"},
  {"ratf33",      233,  2, 0,  "3rd/3rd order rational function"},
  {"1exp",        301,  2, 0,  "single exponential"},
  {"2exp",        302,  2, 0,  "sum of 2 exponentials"},
  {"3exp",        303,  2, 0,  "sum of 3 exponentials"},
  {"4exp",        304,  2, 0,  "sum of 4 exponentials"},
  {"5exp",        305,  2, 0,  "sum of 5 exponentials"},
  {"lundq1",      321,  2, 0,  "Lundqvist function"},
  {"lundq2",      322,  2, 0,  "sum of 2 Lundqvist functions"},
  {"lundq3",      323,  2, 0,  "sum of 3 Lundqvist functions"},
  {"ebolinf",     331,  2, 0,  "exponential bolus infusion function"},
  {"ebolinfrw",   332,  2, 0,  "Kudomi's exponential bolus infusion function for radiowater"},
  {"ebolinfaz",   334,  2, 0,  "exponential bolus function approaching zero"},
  {"ppfpk11195",  351,  2, 0,  "exp function for PK11195 plasma parent fraction"},
  {"igam",        401,  2, 2,  "Lower incomplete gamma function divided by gamma function"},
  {"igamc",       402,  2, 2,  "Upper incomplete gamma function divided by gamma function"},
  {"ppfigam",     403,  2, 5,  "Inverted gamma cdf for plasma parent fraction"},
  {"mpfhill",     841,  2, 0,  "Hill function"},
  {"ppfhill",     842,  2, 0,  "Hill function (1-f(x))"},
  {"ppfhille",    843,  2, 0,  "Hill function (1-f(x)) with ascending or descending end"},
  {"mpfb",        844,  2, 0,  "Hill function with background"},
  {"ppfb",        845,  2, 0,  "Hill function (A-f(x))"},
  {"ppfehill1",   846,  2, 0,  "Extended Hill function for plasma parent fraction"},
  {"mpfehill1",   847,  2, 0,  "Extended Hill function for plasma metabolite fraction"},
  {"ppfefill2",   848,  2, 0,  "Extended Hill function #2 for plasma parent fraction"},
  {"mpfehill2",   849,  2, 0,  "Extended Hill function #2 for plasma metabolite fraction"},
  {"mpfmamede",   851,  2, 0,  "Mamede function"},
  {"ppfmamede",   852,  2, 0,  "Mamede function (1-f(x)"},
  {"ppfmeyer",    861,  2, 0,  "Meyer function for plasma parent fraction"},
  {"mpfmeyer",    862,  2, 0,  "Meyer function for plasma metabolite fraction"},
  {"ppfemeyer",   863,  2, 0,  "extended Meyer function for plasma parent fraction"},
  {"mpfemeyer",   864,  2, 0,  "extended Meyer function for plasma metabolite fraction"},
  {"ppf3hill",    871,  2, 0,  "1-3 metabolite Hill function for plasma parent fraction"},
  {"mpf3hill1",   872,  2, 0,  "1-3 metabolite Hill function for plasma metab1 fraction"},
  {"mpf3hill2",   873,  2, 0,  "1-3 metabolite Hill function for plasma metab2 fraction"},
  {"mpf3hill3",   874,  2, 0,  "1-3 metabolite Hill function for plasma metab3 fraction"},
  {"ppf3pow",     881,  2, 0,  "1-3 metabolite power function for plasma parent fraction"},
  {"mpf3pow1",    882,  2, 0,  "1-3 metabolite power function for plasma metab1 fraction"},
  {"mpf3pow2",    883,  2, 0,  "1-3 metabolite power function for plasma metab2 fraction"},
  {"mpf3pow3",    884,  2, 0,  "1-3 metabolite power function for plasma metab3 fraction"},
  {"step",        1010, 2, 2,  "Step function with x values and y values after those"},
  {"ratf32d",     1232, 2, 0,  "3rd/2nd order rational function with delay"},
  {"fengm2s",     1312, 2, 5,  "PTAC function with two exponentials"},
  {"fengm2",      1313, 2, 7,  "PTAC function with three exponentials"},
  {"fengm2e",     1314, 2, 9,  "PTAC function with four exponentials"},
  {"gammav",      1401, 2, 4,  "Gamma variate"},
  {"gammavb",     1402, 2, 0,  "Gamma variate with background"},
  {"weibullcdfd", 1421, 2, 4,  "Weibull cdf with delay"},
  {"weibullcdfdd",1423, 2, 5,  "Weibull cdf and pdf sum with delay"},
  {"dmsurge",     1430, 2, 0,  "Sum of surge functions with delay"},
  {"surge",       1431, 2, 2,  "Surge function with AUC as parameter"},
  {"tradsurge",   1432, 2, 2,  "Traditional surge function"},
  {"surgerecirc", 1433, 2, 3,  "Surge function with recirculation"},
  {"p2bsrc",      1434, 2, 4,  "Surge function with recirculation for plasma-to-blood ratio"},
  {"surgefdgaif", 1435, 2, 5,  "Surge function for late FDG AIF with delay"},
  {"erlangpdf",   1441, 2, 3,  "Probability density function of Erlang distribution"},
  {"hilld",       1801, 2, 0,  "Hill function with delay"},
  {"hilldd",      1811, 2, 0,  "derivative of Hill function with delay"},
  {"hillssd",     1821, 2, 0,  "sum of Hill function and derivative with delay"},
  {"surgerecircd",1833, 2, 4,  "Surge function with recirculation and delay"},
  {"imgprofile",  2111, 2, 0,  "image profile"},
  {"p2brf",       2233, 2, 7,  "Rational function for plasma-to-blood ratio"},
  {"p2bfm2",      2313, 2, 7,  "Feng M2 function for plasma-to-blood ratio"},
  {"p2bhill",     2841, 2, 4,  "Hill function for plasma-to-blood ratio"},
  /* compartmental models */
  {"FRTM",                 0, 1, 4,  "full reference tissue model"},
  {"SRTM",                 0, 1, 3,  "simplified reference tissue model"},
  {"RRTM",                 0, 1, 3,  "reduced reference tissue model"},
  {"TRTM",                 0, 1, 3,  "transport-limited reference tissue model"},
  {"1TCM",                 0, 1, 5,  "one-tissue compartmental model"},
  {"2TCM",                 0, 1, 7,  "two-tissue compartmental model"},
  {"PAR3TCM",              0, 1, 9,  "parallel three-tissue compartmental model"},
  {"SER3TCM",              0, 1, 9,  "serial three-tissue compartmental model"},
  {"TTM",                  0, 1, 2,  "transit-time model"},
  {"radiowater",           0, 1, 4,  "radiowater perfusion model"},
  {"radiowater2",          0, 1, 6,  "radiowater perfusion model 2"},
  {"radiowater-lung",      0, 1, 5,  "radiowater perfusion model for lungs"},
  {"radiowater-liver",     0, 1, 7,  "radiowater perfusion model for liver"},
  {"radiowater-liver-TTM", 0, 1, 6,  "radiowater perfusion TTM model for liver"},
  {"oxygen-brain",         0, 1, 7,  "oxygen model for brain"},
  {"dispdelay",            0, 1, 2,  "dispersion and delay in radiowater studies"},
  /* pharmacokinetic compartment models */
  {"O1CM",        0,    3, 3,  "one-compartment model with first-order absorption and elimination"},
  {"O2CM",        0,    3, 5,  "two-compartment model with first-order absorption and elimination"},
  /* combinations of functions and models */
  {"ebolinfdd",   3331,12, 8,  "exponential bolus infusion function plus delay and dispersion"},
  {"",            0,    0, 0,  ""}
};
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int modelNr()
{
  unsigned int i=0;
  while(tpc_model[i].code[0]!='\0') i++;
  return(i);
}
/*****************************************************************************/
 
/*****************************************************************************/
char *modelCode(
  const unsigned int i
) {
  if(i>=modelNr()) return((char*)NULL);
  return(tpc_model[i].code);
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int modelOldId(
  const unsigned int i
) {
  return(tpc_model[i].id);
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int modelOld2New(
  const unsigned int i
) {
  unsigned int n=modelNr();
  if(i==0) return(0);
  for(unsigned int j=1; j<n; j++) if(tpc_model[j].id==i) return(j);
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/
char *modelDesc(
  const unsigned int i
) {
  if(i>=modelNr()) return((char*)NULL);
  return(tpc_model[i].desc);
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int modelCodeIndex(
  const char *s
) {
  if(s==NULL || *s=='\0') return(0);
  unsigned int i=0;
  while(tpc_model[i].code[0]!='\0') {
    if(!strcasecmp(tpc_model[i].code, s)) return(i);
    i++;
  }
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int modelParNr(
  const unsigned int code
) {
  if(code<1 || code>=modelNr()) return(0);
  return(tpc_model[code].parNr);
}
/*****************************************************************************/
 
/*****************************************************************************/