(1) cond_probs <- cond_logits_exp / rowSums(cond_logits_exp) cause <- apply(cond_probs, 1, \(prob) sample.int(n_causes + 1L, 1L, prob = prob)) # compute the observed time if needed obs_time <- mapply(\(cause, idx){ if(cause > n_causes) return(delta) # can likely be...