kaplan_meier
Calculate the state survival function using the Kaplan-Meier estimator.
Given that a state transition occurred at time \(t_0\), calculate the probability that a compound is still in the new state at time \(t_0 + \Delta t\).
Options
- -f
Name of the file containing the discrete trajectory. The discrete trajectory must be stored as
numpy.ndarray
either in a binary NumPy .npy file or in a (compressed) NumPy .npz archive. Seemdtools.file_handler.load_dtrj()
for more information about the requirements for the input file.- -o
Output filename.
- -b
First frame to read from the discrete trajectory. Frame numbering starts at zero. Default:
0
.- -e
Last frame to read from the discrete trajectory. This is exclusive, i.e. the last frame read is actually
END - 1
. A value of-1
means to read the very last frame. Default:-1
.- --every
Read every n-th frame from the discrete trajectory. Default:
1
.- --intermittency
Maximum number of frames a compound is allowed to leave its state whilst still being considered to be in this state provided that it returns to this state after the given number of frames. In other words, a compound is only considered to have left its state if it has left it for at least the given number of frames.
- --discard-neg-start
If provided, discard state leavings starting from negative states. Transitions from positive to negative states are regarded as proper state leaving.
- --discard-all-neg
If provided, discard all state leavings starting from or ending in a negative state. This is different to --discard-neg-start in the sense that transitions from positive to negative states are treated as censored.
See also
mdtools.dtrj.kaplan_meier()
The underlying function that calculates the Kaplan-Meier estimate of the survival function
scripts.discretization.kaplan_meier_discrete
Calculate the state survival function using the Kaplan-Meier estimator resolved with respect to the states in the second discrete trajectory
Notes
For more information about the survival function and the Kaplan-Meier
estimator refer to mdtools.dtrj.kaplan_meier()
.