:orphan: 





.. _blq_pk_tut:



Depot + One compartment PK with BLQ
###################################

[Generated automatically as a Tutorial summary]

Inputs
******



Description
===========

:Name: blq_pk

:Title: Depot + One compartment PK with BLQ

:Author: J.R. Hartley

:Abstract: 

| Depot One Comp PK model, with BLQ (below level of quantification) observations.

:Keywords: tutorial; pk; advan4; dep_two_cmp; blq

:Input Script: :download:`blq_pk_tut.pyml <blq_pk_tut.pyml>`

:Diagram: 

Failed to create compartment diagram

True f[X] values
================

.. code-block:: pyml

    f[KA] = 0.2000
    f[CL] = 2.0000
    f[V1] = 50.0000
    f[KA_isv,CL_isv,V1_isv] = [
        [ 0.1000, 0.0200, 0.0100 ],
        [ 0.0200, 0.0300, 0.0200 ],
        [ 0.0100, 0.0200, 0.0900 ],
    ]
    f[PNOISE] = 0.1500
    f[ANOISE] = 0.0100



Starting f[X] values
====================

.. code-block:: pyml

    f[KA] = 1.0000
    f[CL] = 1.0000
    f[V1] = 20.0000
    f[KA_isv,CL_isv,V1_isv] = [
        [ 0.0500, 0.0100, 0.0100 ],
        [ 0.0100, 0.0500, 0.0100 ],
        [ 0.0100, 0.0100, 0.0500 ],
    ]
    f[PNOISE] = 0.1000
    f[ANOISE] = 0.0100



Outputs
*******



Generating and Fitting Summaries
================================

* Gen: :ref:`blq_pk_gen` (gen)
* Fit: :ref:`blq_pk_fit` (fit)

Fitted f[X] values
==================

.. code-block:: pyml

    f[KA] = 0.1957
    f[CL] = 2.0245
    f[V1] = 47.0079
    f[KA_isv,CL_isv,V1_isv] = [
        [ 0.1365, 0.0159, 0.0617 ],
        [ 0.0159, 0.0390, 0.0200 ],
        [ 0.0617, 0.0200, 0.1256 ],
    ]
    f[PNOISE] = 0.1389
    f[ANOISE] = 0.0100



Plots
*****



Dense comp plots
================



.. thumbnail:: images/comp_dense/000001.svg
    :width: 200px


.. thumbnail:: images/comp_dense/000002.svg
    :width: 200px


.. thumbnail:: images/comp_dense/000003.svg
    :width: 200px


Alternatively see :ref:`blq_pk_dense_comp_plots`

Comparison
**********



True objective value
====================


.. code-block:: pyml

    -756.5645



Final fitted objective value
============================


.. code-block:: pyml

    -763.9910



Compare Main f[X]
=================



.. csv-table:: 
    :file: fx_comp_main.csv
    :header-rows: 1


Compare Noise f[X]
==================



.. csv-table:: 
    :file: fx_comp_noise.csv
    :header-rows: 1


Compare Variance f[X]
=====================



.. csv-table:: 
    :file: fx_comp_variance.csv
    :header-rows: 1
