N.I.M.R.O.D.  
Functions/Subroutines

observationModel.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine SystemResolution (Y, iout, simul, deriv, ifail)
 Observation model definition.
subroutine putObservationModel (data)
subroutine removeObservationModel (data)

Function Documentation

subroutine putObservationModel ( double precision,dimension(nbpatienta,tdef2,nbcomp),intent(inout)  data)

Definition at line 64 of file observationModel.f90.

Referenced by readPK().

Here is the caller graph for this function:

subroutine removeObservationModel ( double precision,dimension(nbpatienta,tdef2,nbcomp),intent(inout)  data)

Definition at line 74 of file observationModel.f90.

Referenced by EBCurves().

Here is the caller graph for this function:

subroutine SystemResolution ( double precision,dimension(2*nbcomp),intent(inout)  Y,
integer,intent(in)  iout,
double precision,dimension(tdef2,npmcomp),intent(inout)  simul,
double precision,dimension(tdef2,npmcomp),intent(inout)  deriv,
integer,intent(out)  ifail 
)

Observation model definition.

AUTHOR : Melanie Prague Daniel Commenges Julia Drylewicz Jeremy guedj Rodolphe Thiebaut

DESCRIPTION :

In practice, we do not observe $X^i(t)$ directly, but we have discrete-time observations $Y^i(t_{ij})$ of some functions of $X^i(t)$. We assume that there are known link functions $g_m(.),$ $m=1,\dots , M$, leading to an additive measurement error model. For $i=1, \dots ,n$, $m=1,\dots , M$ and $j=1, \dots T_i$ (the number of observation time for subject $i$), we observe:

\[ Y^i_{m}(t_{ij})=g_m(X^i(t_{ij}))+\epsilon_{ijm}, \qquad \epsilon_{ijm}\sim \mathcal{N}\left(0,\sigma^2_m\right). \nonumber \]

This subroutine can also catch solver error or inconsistency regarding plausible values.

CAUTION : Must "certainly/maybe" be modified when the model changes

MODIFICATION:

01/09/2012 - Prague - Refactoring

INFORMATIONS:

Parameters:
[in,out]YBiomarker trajectories through ODE system value (untransformed)
[in]ioutTime of computation
[in,out]simulODE system transformed trajectories values
[in,out]derivDerivatives of ODE system transformed trajectories values
[out]ifailFailure and inconsistency indicator

Definition at line 42 of file observationModel.f90.

References WorkingSharedValues::systeme.

Referenced by solution().

Here is the caller graph for this function: