N.I.M.R.O.D.  

moduleModelSpe.f90

Go to the documentation of this file.
00001 !------------------------------------------------------------------------------
00002 ! N.I.M.R.O.D. - Normal approximation Inference in Models with Random
00003 !                effects based on Ordinary Differential equations
00004 !------------------------------------------------------------------------------
00005 !
00006 ! VERSION : 1.0
00007 !
00008 ! MODULE: moduleModelSpe.f90
00024 
00025 !------------------------------------------------------------------------------
00028 module VectorLength
00029     use debugage
00030     implicit none
00031     integer, parameter:: npmbio=3 ! Number of biological parameters
00032     integer, parameter:: npmexpl=0 ! Number of explicative variables (estimated regression coefficient)
00033     integer, parameter:: ndim=3 ! Number of random effects
00034     integer,parameter :: npmcomp=1 ! Number of observed compartments
00035     integer,parameter :: nbcomp=2 ! Total Number of compartments in ODE systems
00036     integer,parameter :: tdef2=650 ! Maximum observation time
00037     integer,parameter :: nbpatienta=100 ! Maximum number of patient
00038     integer,parameter :: npm=npmbio+npmexpl+ndim+npmcomp
00039     integer,parameter :: nf=1
00040 end module VectorLength
00041 
00042 !------------------------------------------------------------------------------
00045 module systemODE
00046     use VectorLength
00047     implicit none
00048     ! Parametres variables Systeme ODE
00049      double precision,save ::cl0,V00,Vt0,ktp0,kpt0,ka0 ! Name if parameter id considered as constant (cf transfo)
00050      double precision,save ::cl,V0,Vt,ktp,kpt,dose0,ka! Name used in ODE
00051      double precision, save::dcl,dV0,dVt,dktp,dkpt,dka! Name of derivative (cf transfo)
00052 contains
00053     subroutine ODEinit()
00054         implicit none
00055         ! Parametres variables Systeme ODE
00056         cl=0.0D0
00057         V0=0.0D0
00058         Vt=0.0D0
00059         ktp=0.0D0
00060         kpt=0.0D0
00061         dose0=0.0D0
00062         cl0=0.0D0
00063         V00=0.0D0
00064         Vt0=0.0D0
00065         ktp0=0.0D0
00066         kpt0=0.0D0
00067         dcl=0.0D0
00068         dV0=0.0D0
00069         dVt=0.0D0
00070         dktp=0.0D0
00071         dkpt=0.0D0
00072         dka=0.0D0
00073         ka=0.0D0
00074         ka0=0.0D0
00075         return
00076     end subroutine ODEinit
00077 end module systemODE
00078 
00079 
00080 !------------------------------------------------------------------------------
00083 module Datareading
00084     use VectorLength
00085     implicit none
00086           ! Input Files
00087             character(len=30)::biomarkersFile,treatmentFile
00088 
00089           ! variables ReadData
00090             integer,dimension(nbpatienta,tdef2,nbcomp),save::schedule
00091             double precision,dimension(nbpatienta,tdef2,nbcomp),save::donnees
00092             integer,save::nbpatient
00093             integer,dimension(nbpatienta),save::idpat
00094             double precision,dimension(nbpatienta),save::dose
00095             integer,dimension(nbpatienta,nbcomp),save::nbobs
00096 end module Datareading
00097