Estep {frailtyEM} | R Documentation |
Perform the E step calculations
Description
This is an inner wrapper for the C++ functions which perform the E step and is not intended to be used directly.
This function does not check the input.
For a data set with K
clusters,
Usage
Estep(c, c_lt, delta, alpha, bbeta, pvfm, dist)
Arguments
c |
Vector of length |
c_lt |
Vector of length |
delta |
Vector of integers of length |
alpha , bbeta |
Parameters of the frailty distribution |
pvfm |
Parameter for the PVF distribution, only matters in that case |
dist |
One of 0 (for gamma), 1 (for stable) or 2 (for PVF) |
Value
A K x 3
matrix where the first column and the second column are the numerators
and the denominators of the frailty fraction (without the Laplace transform) and the
last column is the log(denominator) + log-Laplace transform, i.e. the log-likelihood contribution