Returns a closure computing P(K=j | theta) for all components j, marginalized over the system failure time T. By Theorem 5 of the foundational paper, this equals E_T[P(K=j | T, theta)].
Value
a function with signature function(par, ...) returning an m-vector
where element j gives P(K=j | theta)
Details
The default method uses Monte Carlo integration via rdata().