# Create sample data
library(dplyr)
df <- data.frame(
id = c(1, 2, 3, 4, 4, 4, 4, 5, 5),
time = c(8, 1, 5, 2, 6, 7, 8, 3, 3),
cause = c(0, 0, 2, 1, 1, 1, 0, 1, 2)
) |>
arrange(id, time)
# Calculate MCC
mcc_result <- mcc(df, "id", "time", "cause")
# Extract MCC estimates
estimates <- mcc_estimates(mcc_result)
print(estimates)
# For grouped analysis
df_grouped <- df |>
mutate(group = c("A", "A", "B", "B", "B", "B", "B", "A", "A"))
mcc_grouped <- mcc(df_grouped, "id", "time", "cause", by = "group")
estimates_grouped <- mcc_estimates(mcc_grouped)
print(estimates_grouped)
# Clean up
rm(df, df_grouped, mcc_result, mcc_grouped, estimates, estimates_grouped)
Run the code above in your browser using DataLab