compare_psm_likes {psm3mkv}

Compare likelihoods of PSMs

Description

Compare the total log-likelihood values for the patient-level dataset after fitting PSM-simple and PSM-complex models to each combination of endpoint distributions

Usage

compare_psm_likes(ptdata, fitslist, cuttime = 0)

Arguments

ptdata

Dataset of patient level data. Must be a tibble with columns named:

• ptid: patient identifier

• pfs.durn: duration of PFS from baseline

• pfs.flag: event flag for PFS (=1 if progression or death occurred, 0 for censoring)

• os.durn: duration of OS from baseline

• os.flag: event flag for OS (=1 if death occurred, 0 for censoring)

• ttp.durn: duration of TTP from baseline (usually should be equal to pfs.durn)

• ttp.flag: event flag for TTP (=1 if progression occurred, 0 for censoring).

fitslist

List of distribution fits to relevant endpoints, after calling fit_ends_mods_par() or fit_ends_mods_spl()

cuttime

Time cutoff - this is nonzero for two-piece models.

Value

List containing

• results: Dataset of calculation results for each model

• bests: Tibble indicating which is the best fitting model individually or jointly, to each endpoint, according to AIC or BIC

Examples

# Fit parametric distributions to a dataset
bosonc <- create_dummydata("flexbosms")
parfits <- fit_ends_mods_par(bosonc)

splfits <- fit_ends_mods_spl(bosonc)
# Present comparison of likelihood calculations
compare_psm_likes(bosonc, parfits)
compare_psm_likes(bosonc, splfits)