The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

Type: Package
Title: Conditional Aalen-Johansen Estimation
Version: 1.0
Maintainer: Martin Bladt <martinbladt@math.ku.dk>
Description: Provides the conditional Nelson-Aalen and Aalen-Johansen estimators. The methods are based on Bladt & Furrer (2023), in preparation.
License: GPL-2 | GPL-3 [expanded from: GPL (≥ 2)]
Encoding: UTF-8
RoxygenNote: 7.2.1
Suggests: knitr, rmarkdown
VignetteBuilder: knitr
NeedsCompilation: no
Packaged: 2023-02-28 18:01:12 UTC; martinbladt
Author: Martin Bladt [aut, cre], Christian Furrer [aut]
Repository: CRAN
Date/Publication: 2023-03-01 10:42:09 UTC

Compute the conditional Aalen-Johansen estimator.

Description

Compute the conditional Aalen-Johansen estimator.

Usage

aalen_johansen(
  data,
  x = NULL,
  a = NULL,
  p = NULL,
  alpha = 0.05,
  collapse = FALSE
)

Arguments

data

A list of trajectory data for each individual.

x

A numeric value for conditioning.

a

A bandwidth. Default uses an asymmetric version using alpha.

p

An integer representing the number of states. The absorbing state is last.

alpha

A probability around the point x, for asymmetric sub-sampling.

collapse

Logical, whether to collapse the last state of the model.

Value

A list containing the Aalen-Johansen estimator, the Nelson-Aalen estimator, and related quantities.


Calculate the product integral of a matrix function

Description

Calculate the product integral of a matrix function

Usage

prodint(start, end, step_size, lambda)

Arguments

start

Start time.

end

End time.

step_size

Step size of the grid.

lambda

A given matrix function.

Value

The product integral of the given matrix function.


Simulate the path of a time-inhomogeneous (semi-)Markov process until a maximal time

Description

Simulate the path of a time-inhomogeneous (semi-)Markov process until a maximal time

Usage

sim_path(i, rates, dists, t = 0, u = 0, tn = Inf, abs = numeric(0), bs = NA)

Arguments

i

The initial state, integer.

rates

The total transition rates out of states, a function with arguments state (integer), time (numeric), and duration (numeric) returning a rate (numeric).

dists

The distribution of marks, a function with arguments state (integer), time (numeric), and duration (numeric) returning a probability vector.

t

The initial time, numeric.

u

The initial duration (since the last transition), numeric. By default equal to zero

tn

The maximal time, numeric. By default equal to inifinity

abs

Vector indicating which states are absorbing. By default the last state is absorbing.

bs

Vector of upper bounds on the total transition rates. By default the bounds are determined using optimize, which might only identify a local maximum.

Value

A list concerning jump times and states, with the first time being the initial time t and state and the last time being tn (if not absorbed)

Examples


jump_rate <- function(i, t, u){if(i == 1){3*t} else if(i == 2){5*t} else{0}}
mark_dist <- function(i, s, v){if(i == 1){c(0, 1/3, 2/3)} else if(i == 2){c(1/5, 0, 4/5)} else{0}}
sim <- sim_path(sample(1:2, 1), t = 0, tn = 2, rates = jump_rate, dists = mark_dist)
sim

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.