search for: simulweib

...roduced Has anybody else come across this issue, or could suggest why parfm struggles with low event rates? Or could someone please run my code to see if they get the same issue? Full reproducible code is presented below. Many thanks for any help, Alex CODE: ### Create function to generate data simulWeib <- function(N, lambda, rho, beta1, beta2, beta3, beta4, rateC, sigma) { # covariate --> N Bernoulli trials x1 <- sample(x=c(0, 1), size=N, replace=TRUE, prob=c(0.5, 0.5)) # Now create random effect stuff # Create one vector of length N, all drawn from same normal distribution ra...

...the observation level? Many thanks for any solutions/or references which may be helpful. Reproducible example here (using coxme): setwd("/mnt/ja01-home01/mbrxsap3/phd_risk/R/p4_run_analysis/") library(coxme) library(survival) ### Create data with a group level random effect simulWeib.group <- function(N, lambda, rho, beta1, beta2, beta3, beta4, rateC, sigma, M) { # covariate --> N Bernoulli trials x1 <- sample(x=c(0, 1), size=N, replace=TRUE, prob=c(0.5, 0.5)) x2 <- sample(x=c(0, 1), size=N, replace=TRUE, prob=c(0.5, 0.5)) x3 <- sample(x=c(0, 1)...