The R codes and the Mlxtran model files can be found here


# Regular steps - version 1

In this first version, the value of the step function is defined by the amount of the input pulses. Tlag defines the length of the step. Using p=-1 is equivalent to define a negative amount

# Regular steps - version 2

In this second version, the value of the step function is defined as an input parameter A. The amount of the input pulse amtDose is ignored since it is multiplied by p=A/amtDose.

Remark that we define f_0=A as initial condition in this example while the signal is a sequence of negative and positive pulses.

# Irregular steps

In this example, the signal is a sequence of steps with irregular heights and lengths.

# Step function with continuous transitions

Using infusions with constant rates lead to replace the Heavyside functions (sudden jumps) with continuous transitions.