Type:	Package
Title:	Personalised Synthetic Controls
Version:	1.3.0
Maintainer:	Richard Jackson <richJ23@liverpool.ac.uk>
Description:	Allows the comparison of data cohorts (DC) against a Counter Factual Model (CFM) and measures the difference in terms of an efficacy parameter. Allows the application of Personalised Synthetic Controls.
License:	GPL-3
Encoding:	UTF-8
LazyData:	true
RoxygenNote:	7.3.2
URL:	https://github.com/richjjackson/psc/, https://github.com/richJJackson/psc
BugReports:	https://github.com/richJJackson/psc/issues
Depends:	R (≥ 3.5.0)
Imports:	mvtnorm, survival, enrichwith, stats, flexsurv, ggplot2, survminer, gtsummary, RColorBrewer, waffle
Suggests:	knitr, rmarkdown, devtools, ggpubr, testthat (≥ 3.0.0)
VignetteBuilder:	knitr
Config/testthat/edition:	3
NeedsCompilation:	no
Packaged:	2025-05-29 14:49:07 UTC; richardjackson
Author:	Richard Jackson [cre, aut, cph]
Repository:	CRAN
Date/Publication:	2025-05-29 15:10:02 UTC

psc: Personalised Synthetic Controls

Description

Allows the comparison of data cohorts (DC) against a Counter Factual Model (CFM) and measures the difference in terms of an efficacy parameter. Allows the application of Personalised Synthetic Controls.

Author(s)

Maintainer: Richard Jackson richJ23@liverpool.ac.uk (ORCID) [copyright holder]

See Also

Useful links:

• https://github.com/richjjackson/psc/

• https://github.com/richJJackson/psc

• Report bugs at https://github.com/richJJackson/psc/issues

acc

Description

Function to accept (or not) a proposed solution used as part of the MCMC procedure

Usage

acc(old, new)

Arguments

Details

A function for the evaluation of two likelihoods as part of the MCMC procedure

Value

returns the an evaluation of old/new > U where U is a draw from the uniform distribution

Example model for a survival outcome

Description

A generated model with a binary endpoint and a logistic link function. Data for the model were synthetically generated and are based on a dataset to evaulate the use of Sorafenib in HCC akin to the PROSASH model (see ?psc::surv.mod for more details)

Usage

bin.mod

Format

A model of class 'glm':

vi

vascular invasion

ecog

ECOG performance Status

logafp

AFP - log scale

alb

albumin

logcreat

Creatinine - log scale

allmets

metastesis

Source

simulated

Summarising data within a CFM

Description

The pscCFM creates a model object which is stripped of identifiable information. The cfmDataSumm fucntion supplies a tabulated form of the dataset for summary information

Usage

cfmDataSumm(cfm)

Arguments

Value

a summary table

Visualising data within a CFM

Description

The pscCFM creates a model object which is stripped of identifiable information. The cfmDataVis fucntion supplies a visualised form of the dataset for summary information

Usage

cfmDataVis(cfm)

Arguments

Value

a list of grobs for each model covariate

Returns the coefficeint estimate of a psc object.

Description

Returns basic measures of the posterior distribution obtained from the pscfit.R function

Usage

## S3 method for class 'psc'
coef(object, ..., level = 0.05)

Arguments

Value

The summary of the posterior distribution for the efficacy parameter in terms of the median and 95

Example model for a survival outcome

Description

A generated model with a continuous data endpoint and a identity link function. Data for the model were synthetically generated and are based on a dataset to evaulate the use of Sorafenib in HCC akin to the PROSASH model (see ?psc::surv.mod for more details)

Usage

cont.mod

Format

A model of class 'glm':

ecog

ECOG performance Status

logafp

AFP - log scale

alb

albumin

logcreat

Creatinine - log scale

Source

simulated

Example model for a survival outcome

Description

A generated model with a count data endpoint and a log link function. Data for the model were synthetically generated and are based on a dataset to evaulate the use of Sorafenib in HCC akin to the PROSASH model (see ?psc::surv.mod for more details)

Usage

count.mod

Format

A model of class 'glm':

ecog

ECOG performance Status

logafp

AFP - log scale

alb

albumin

logcreat

Creatinine - log scale

Source

simulated

Example Dataset of patients with aHCC receiving Lenvetanib

Description

A dataset containing 100 simulated patients. Data are based on the data used to generate PROSASH survival model -see ?psc::surv.mod for more detials.

Usage

data

Format

A model of class 'flezsurvreg':

gamma

cumulative baseline hazard parameters

vi

vascular invasion

age60

patient age (centred at 60)

ecog

ECOG performance Status

logafp

AFP - log scale

alb

albumin

logcreat

Creatinine - log scale

allmets

metastesis

ageVasInv

centred age nested within vascular invasion

time

survival time

cen

censoring indicator

os

survival time

count

exapmple outcome for count data

trt

exapmple identifier for mulitple treatment comparisons

aet

Aetiology

Source

simulated

A generic function for cleaning data ready for analysis

Description

A generic function for cleaning data ready for analysis

Usage

dataComb(CFM, DC, id = NULL, trt = NULL, cfmOb = FALSE)

Arguments

Value

datComb returns a list containing objects which detial the components of both the Counter Factual Model (CFM) and the Data Cohort (DC) the required exported components of the model and a cleaned data cohort.

Examples

bin.mod <- psc::bin.mod
data <- psc::data
dc <- dataComb(bin.mod,data)

Function for cleaning the data of a model with class 'flexsurvreg'

Description

The purpose of this function is to prepare the dataset and the counter-factual model for estimation and is the first step pf the pscfit.R process. The output is a complete-case dataset where the data names match the variables used in the CFM.

Usage

## S3 method for class 'flexsurvreg'
dataComb(CFM, DC, id = NULL, trt = NULL, cfmOb = FALSE)

Arguments

Value

a list containing objects which specifiy the required exported components of the model and a cleaned data cohort. Attributes include

• 'model.type' specifying the class of model to be used as the CFM

• 'model_extract' sepcifying the model componets required for estimation

• 'cov' a cleaned dataset of covariates

• 'outcome' a cleaned dataset containing the outcomes

Function for cleaning the data of a model with class 'flexsurvreg'

Description

The purpose of this function is to prepare the dataset and the counter-factual model for estimation and is the first step pf the pscfit.R process. The output is a complete-case dataset where the data names match the variables used in the CFM.

Usage

## S3 method for class 'glm'
dataComb(CFM, DC, id = NULL, trt = NULL, cfmOb = FALSE)

Arguments

Value

a list containing objects which specifiy the required exported components of the model and a cleaned data cohort. Attirbutes of the 'cleaned' object include: Attributes include

• 'model.type' specifying the class of model to be used as the CFM

• 'model_extract' sepcifying the model componets required for estimation

• 'cov' a cleaned dataset of covariates

• 'outcome' a cleaned dataset containing the outcomes

Function for cleaning the data of a model with class

Description

The purpose of this function is to prepare the dataset and the counter-factual model for estimation and is the first step pf the pscfit.R process. The output is a complete-case dataset where the data names match the variables used in the CFM.

Usage

## S3 method for class 'pscCFM'
dataComb(CFM, DC, id = NULL, trt = NULL, cfmOb = TRUE)

Arguments

Value

a list containing objects which specifiy the required exported components of the model and a cleaned data cohort. Attirbutes of the 'cleaned' object include: Attributes include

• 'model.type' specifying the class of model to be used as the CFM

• 'model_extract' sepcifying the model componets required for estimation

• 'cov' a cleaned dataset of covariates

• 'outcome' a cleaned dataset containing the outcomes

A function which provides basic summary information of a matrix x

Description

A function which provides basic summaries of data provided within a data frame of covariets. Summaries are in the form of frequencies of counts and associated percentages for categorical data and median (IQR) for continuous data. Intended for use with the pscCFM.R function to define the setting in which a model has been generated

Usage

dataSumm(x)

Arguments

Details

Categorical/Character data are summarised by a table and Continuous data are summarised as median (IQR)

A function to ensure that data from the cfm and data cohort are compatible

Description

The purpose of this function is to run a series of checks to ensure that the data included in the data cohort is comparable to the counter-factual model

Usage

data_match(cls, lev, dc.data)

Arguments

Value

a dataset which is checked and compatible with the CFM

Example Dataset of patients treated with GemCap in the ESPAC-4 trial

Description

A dataset containing 346 simulated patients. Data are based on the patietns randomised to revceive GemCap in the ESPAC-4 trial

Usage

e4_data

Format

A model of class 'flezsurvreg':

time

survival time

cen

censoring indicator

nodes

negative (n=1) or positive (n=2) lymph nodes

grade

tumour grade (1,2 or 3)

lca199

log transformed ca19.9

t

T-stage (1,2 or 3)

Source

simulated

Visualising Categorical Data

Description

A function which summarises categorical data using a waffle plot

Usage

facVis(x, nm)

Arguments

Value

a ggplot object

Visualising Categorical Data

Description

A function which compares visually a new categorical covariate against equivalent data from a CFM

Usage

facVisComp(p, x)

Arguments

Value

a ggplot object

Model for a survival outcome based on Gemcitbine patietns from ESPAC-3

Description

A generated model with a survival endpoint and a cuymulative hazard function estimated using flexible parametric splines. Data for the model were obtained from the ESPAC-3 trials

Usage

gemCFM

Format

A model of class 'pscCFM' containg a 'flexsurvreg' model:

gamma

cumulative baseline hazard parameters

nodes

negative (n=1) or positive (n=2) lymph nodes

grade

tumour grade (1,2 or 3)

lca199

log transformed ca19.9

t

T-stage (1,2 or 3)

Source

simulated

Fucntion for estimating initial parameter values 'flexsurvreg'

Description

A generic function for estimating the initial parameters for estimation as part of the pscfti.R function. Parameter estimates are obtained using standard optomisation methods provided by the 'optim' set of functions. For a single parameter a Brent method is applied. For mutliple treatment comparisons the 'BFGS' method is applied.

Usage

initParm(CFM, DC_clean, trt)

Arguments

Details

This function takes the liklihood for a 'flexsurvreg' model and uses 'optim' to fit the likelihood.

Value

an 'optim' output giving the parameter values to be supplied as a starting value for the mcmc routine.

Fucntion for estimating initial parameter values 'flexsurvreg'

Description

Fucntion for estimating initial parameter values 'flexsurvreg'

Usage

## S3 method for class 'flexsurvreg'
initParm(CFM, DC_clean, trt = NULL)

Arguments

Details

This function takes the liklihood for a 'flexsurvreg' model and uses 'optim' to fit the likelihood.

Value

an 'optim' output giving the parameter values to be supplied as a starting value for the mcmc routine.

Fucntion for estimating initial parameter values

Description

Fucntion for estimating initial parameter values

Usage

## S3 method for class 'glm'
initParm(CFM, DC_clean, trt = trt)

Arguments

Details

This function takes the liklihood for a 'flexsurvreg' model and uses 'optim' to fit the likelihood.

Value

an 'optim' output giving the parameter values to be supplied as a starting value for the mcmc routine.

Fucntion for estimating initial parameter values

Description

Fucntion for estimating initial parameter values

Usage

## S3 method for class 'pscCFM'
initParm(CFM, DC_clean, trt = NULL)

Arguments

Details

This function takes the liklihood for a 'flexsurvreg' model and uses 'optim' to fit the likelihood.

Value

an 'optim' output giving the parameter values to be supplied as a starting value for the mcmc routine.

Likelihood function for a psc model of class 'flexsurvreg'

Description

A function which defines the likelihood for a PSC model where the Counter Factual Model (CFM) takes the form of a 'flexsurvreg' object and an efficacy parameter (\beta) is being estimated. For more details on fitting please see ?pscfit and ?pscEst

Usage

lik.flexsurvreg(beta, DC_clean)

Arguments

Details

A likelihood function for use by pscfit for a model of class 'flexsurvreg'

Likelihood function for a psc model of class 'flexsurv' with multiple treatment comparisons

Description

A function which defines the likelihood for a PSC model where the Counter Factual Model (CFM) takes the form of a 'flexsurvreg' object and a mulitple efficacy parameters (\beta) is being estimated. For more details on fitting please see ?pscfit and ?pscEst

Usage

lik.flexsurvreg.mtc(beta, DC_clean)

Arguments

Details

A likelihood function for use by pscfit for a model of class 'flexsurvreg' where multiple treatment comparisons are required

Likelihood function for a psc model of class 'glm'

Description

A function which defines the likelihood for a PSC model where the Counter Factual Model (CFM) takes the form of a 'glm' object and an efficacy parameter (\beta) is being estimated. For more details on fitting please see ?pscfit and ?pscEst

Usage

lik.glm(beta, DC_clean)

Arguments

Details

A likelihood function for use by pscfit for a model of class 'glm'

Likelihood function for a psc model of class 'glm' with multiple treatment comparisons

Description

A function which defines the likelihood for a PSC model where the Counter Factual Model (CFM) takes the form of a 'glm' object and a mulitple efficacy parameters (\beta) is being estimated. For more details on fitting please see ?pscfit and ?pscEst

Usage

lik.glm.mtc(beta, DC_clean)

Arguments

Details

A likelihood function for use by pscfit for a model of class 'flexsurvreg' where multiple treatment comparisons are required

Estimates the linear predictor of a psc object

Description

A function which created a linear predictor from a model and an external dataset. If required, linear predictors will be provided on the response scale. For a survival outcome, this will take the form of median survival estimates.

Usage

linPred(DC_clean, resp = FALSE)

Arguments

Details

A function which combines the data from the data cohort against the model parameters of the PSC

Value

Extracts the linear predictor from a object containing both a counter factual model and a data cohort which is created using the dataComb() fucntion.

A generic function for extracting model information

Description

A generic function for extracting model information

Usage

modelExtract(CFM)

Arguments

Details

A function for extracting the model information required for using pscfit

Value

a list of extracted model components

A generic function for extracting model information

Description

This function extracts model information for use with the pscfit.R function.

Usage

## S3 method for class 'flexsurvreg'
modelExtract(CFM)

Arguments

Details

A function for extracting the model information required for using pscfit

Value

a list of extracted model components. Included are

• mod_class: The class of the model object

• terms: The terms included in the model

• cov_cls: The classes of the model covariates

• cov_lev: The levels of categorical variables

• cov_co: covariate coefficients

• sig: variance-covariance matrix

• haz_co: hazard parameter coefficients

• k: number of knots

• knots: knot position

• lam: lambda parameter

• form: model formula

A generic function for extracting model information

Description

A generic function for extracting model information

Usage

## S3 method for class 'glm'
modelExtract(CFM)

Arguments

Details

A function for extracting the model information required for using pscfit

Value

a list of extracted model components. Included are

• mod_class: The class of the model object

• terms: The terms included in the model

• cov_cls: The classes of the model covariates

• cov_lev: The levels of categorical variables

• cov_co: covariate coefficients

• sig: variance-covariance matrix

• form: model formula

• family: model family

• out.nm: outcome covariates names

modp

Description

A function which rrturns either the input value (if positive) or zero (if negative)

Usage

modp(x)

Arguments

Details

A fucntion which returns a version of x with negative values replacd with 0

Value

a numeric vector with negative values replaced with 0

Visualising Categorical Data

Description

A function which summarises categorical data using a waffle plot

Usage

numVis(x, nm)

Arguments

Value

a ggplot object

Visualising Categorical Data

Description

A function which summarises numerical data using a density plots

Usage

numVisComp(p, x)

Arguments

Value

a ggplot object

Function for Plotting PSC objects

Description

A function which illsutrates the predicted response under the Counter Factual Model (CFM) and the observed response under the experimental treatment(s). Form of the output will depend on the form of the CFM used

Usage

## S3 method for class 'psc'
plot(x, ...)

Arguments

Details

This function plots the expected response of the control treatment along with the observe response rates of the experimental arms

Value

a survival plot corresponding to the psc fit

Examples

bin.mod <- psc::bin.mod
data <- psc::data
bin.psc <- pscfit(bin.mod,data,nsim=3000)
plot(bin.psc)

Function for Plotting PSC objects

Description

A function which illsutrates the predicted response under the counter factual model and the observed response under the experimental treatment(s).

Usage

## S3 method for class 'psc.binary'
plot(x, ...)

Arguments

Details

This function plots the expected response of the control treatment along with the observe response rates of the experimental arms

Value

a survival plot corresponding to the psc fit

Examples

bin.mod <- psc::bin.mod
data <- psc::data
bin.psc <- pscfit(bin.mod,data,nsim=3000)
plot(bin.psc)

Function for Plotting PSC objects

Description

A function which illsutrates the predicted response under the counter factual model and the observed response under the experimental treatment(s).

Usage

## S3 method for class 'psc.cont'
plot(x, ...)

Arguments

Details

This function plots the expected response of the control treatment along with the observe response rates of the experimental arms

Value

a survival plot corresponding to the psc fit

Examples

cont.mod <- psc::cont.mod
data <- psc::data
cont.psc <- pscfit(cont.mod,data,nsim=3000)
plot(cont.psc)

Function for Plotting PSC objects #' A function which illsutrates the predicted response under the counter factual model and the observed response under the experimental treatment(s).

Description

Function for Plotting PSC objects #' A function which illsutrates the predicted response under the counter factual model and the observed response under the experimental treatment(s).

Usage

## S3 method for class 'psc.count'
plot(x, ...)

Arguments

Details

This function plots the expected response of the control treatment along with the observe response rates of the experimental arms

Value

a survival plot corresponding to the psc fit

Examples

count.mod <- psc::count.mod
data <- psc::data
count.psc <- pscfit(count.mod,data,nsim=3000)
plot(count.psc)

Function for Plotting PSC objects

Description

Function for Plotting PSC objects

Usage

## S3 method for class 'psc.flexsurvreg'
plot(x, ...)

Arguments

Details

making use of 'ggsurvplot' in the survminer package, this function plots the expected survival funtion for the 'control' treatment estimated from the CFM along with the Kaplan Meier estimates of the observed events

Value

a survival plot corresponding to the psc fit

Personalised Synthetic Controls - print

Description

Personalised Synthetic Controls - print

Usage

## S3 method for class 'psc'
print(x, ...)

Arguments

Value

printing psc results

quiet_gglist

Description

Ensuring a quiet list of the grobs data are supplied to cfmDataVis

Usage

## S3 method for class 'quiet_gglist'
print(x, ...)

Arguments

Value

A quiet list

quiet_gtsumm

Description

Ensuring a quiet list of the grobs data are supplied to cfmDataVis

Usage

## S3 method for class 'quiet_gtsumm'
print(x, ...)

Arguments

Value

A quiet list

quiet_gtsumm

Description

Ensuring a quiet list of the grobs data are supplied to cfmDataVis

Usage

## S3 method for class 'quiet_list'
print(x, ...)

Arguments

Value

A quiet list

Fitted psc object

Description

An object returned by the pscfit function, inheriting from class psc and representing a fitted personlised synthetic control model.

Usage

psc.object

Format

An object of class NULL of length 0.

Author(s)

Richard Jasckson (richj23@liverpool.ac.uk)

Creating a CFM model which can be shared

Description

Standard R model objects contain within them the datasets used to create the model and as such care is needed when sharing these objects for research. The psc.cfm function creates an object with all identifiable information retracted and includes only the information required to use the models within the psc package

Usage

pscCFM(cfm, dataSumm = T, dataVis = T)

Arguments

Value

a list containing objects which specifiy the required exported components of the model.

Function for performing estimation procedures in 'pscfit'

Description

Function for performing estimation procedures in 'pscfit'

Usage

pscEst(CFM, DC_clean, nsim, start, start.se, trt)

Arguments

Details

Define the set of model parameters B to contain \Gamma which summarize the parameters of the CFM. Prior distributions are defined for B using a multivariate normal distribution \pi (B) \sim MVN(\mu ,\Sigma) where \mu| is the vector of coefficient estimates from the validated model and \Sigma is the variance-covariance matrix. This information is taken directly from the outputs of the parametric model and no further elicitation is required. The prior distirbution for the efficacy parameter (\pi{(\beta)}) is set as an uniformative N(0,1000).

Ultimately the aim is to estimate the posterior distribution for \beta conditional on the distribution of B and the observed data. A full form for the posterior distribution is then given as

P(\beta \vert B,D) \propto L(D \vert B,\beta) \pi(B) \pi(\beta)

Please see 'pscfit' for more details on liklihood formation.

For each iteration of the MCMC procedure, the following algorithm is performed

1. Set and indicator s=1, and define an initial state based on prior hyperparameters for \pi(B) and \pi(\beta) such that b_s = \mu and \tau_s=0

2. Update s = s+1 and draw model parameters b_s from \pi(B) and an draw a proposal estimate of \beta from some target distribution

3. Estimate \Gamma_(i,S)=\nu^T x_i where \nu is the subset of parameters from b_s which relate to the model covariates and define 2 new likelihood functions \Theta_(s,1)=L(D \vert B=b_s,\beta=\tau_(s-1) ) & \Theta_(s,2)= L(D \vert B=b_s,\beta=\tau_s)

4. Draw a single value \psi from a Uniform (0,1) distribution and estimate the condition \omega= \Theta_(s,1)/\Theta_(s,2). If \omega > \psi then accept \tau_s as belonging to the posterior distribution P(\beta \vert B,D) otherwise retain \tau_(s-1)

5. Repeat steps 2 – 4 for the required number of iterations

The result of the algorithm is a posterior distribution for the log hazard ratio, \beta, captures the variability in B through the defined priors \pi{(\beta)}.

Value

A matrix containing the draws form the posterior distribution

Function for estimating initial parameter values 'flexsurvreg'

Description

A function which performs the Bayesian MCMC estimation procedure for estimating the efficacy parameter (\beta) using personalised sunthetic controls methodology.

Usage

## S3 method for class 'flexsurvreg'
pscEst(CFM, DC_clean, nsim, start, start.se, trt = trt)

Arguments

Details

An MCMC routine for fitting a psc model

Value

A matrix containing the draws form the posterior distribution

Function for estimating initial parameter values 'flexsurvreg'

Description

A function which performs the Bayesian MCMC estimation procedure for estimating the efficacy parameter (\beta) using personalised sunthetic controls methodology.

Usage

## S3 method for class 'glm'
pscEst(CFM, DC_clean, nsim, start, start.se, trt = trt)

Arguments

Details

An MCMC routine for fitting a psc model

Value

a matrix containing the draws form the posterior distribution

Function for estimating initial parameter values 'flexsurvreg'

Description

A function which performs the Bayesian MCMC estimation procedure for estimating the efficacy parameter (\beta) using personalised sunthetic controls methodology.

Usage

## S3 method for class 'pscCFM'
pscEst(CFM, DC_clean, nsim, start, start.se, trt = NULL)

Arguments

Details

An MCMC routine for fitting a psc model

Value

A matrix containing the draws form the posterior distribution

Personalised Synthetic Controls - summary

Description

A function which provides a summary of a DC_clean object. To be used either in conjunction with dataComb.R or summary.psc.R

Usage

pscSumm(DC_clean)

Arguments

Value

psc summary results including an estimate of the linear predictor combing the data and the model, an estimate of patient level response and summary statistics of the average responses for the sythenthic and observed populations

Personalised Synthetic Controls model fit

Description

Function which allows comparison of a data cohort against a parametric Counter Factual Model (CFM). The function allows models of the type 'flexsurvreg' and 'glm' to be supplied. The function performs by calculating the linear predictor as a combination of the CFM and the dataset supplied and then selects a likelihood based on the type of model specified. Likelihood is estimated using a Baysian MCMC procedure wherebey the parameters of the CFM acts as informative priors.

Usage

pscfit(CFM, DC, nsim = 5000, id = NULL, trt = NULL)

Arguments

Details

Model currently supports estimation of more than one treatment (using the 'trt') option and esitmation restricted to sub-groups of the data cohort (using the 'id' option.

the pscfit function compares a dataset ('DC') against a parametric model. This is done by selecting a likelihood which is identified by the type of CFM that is supplied. At present, two types of model are supported, a flexible parmaeteric survival model of type 'flexsurvreg' and a geleneralised linear model of type 'glm'.

Where the CFM is of type 'flexsurvreg' the likeihood supplied is of the form:

L(D \vert \Lambda, \Gamma_i) = \prod^{n}_{i=1} f(t_i \vert \Lambda, \Gamma_i)^{c_i} S(t_i|\Lambda, \Gamma_i)^{(1-c_i)}

Where \Lambda defines the cumulative baseline hazard function, \Gamma is the linear predictor and t and c are the event time and indicator variables.

Where the CFM is of the type 'glm' the likelihood supplied is of the form:

L(x \vert \Gamma_i) = \prod^{n}_{i=1} b(x \vert \Gamma_i) \exp{\{\Gamma_i^T t(x) - c(\Gamma_i)\} }

Where b(.), t(.) and c(.) represent the functions of the exponential family. In both cases, \Gamma is defined as:

\Gamma = \gamma x + \beta

Where \gamma are the model coefficients supplied by the CFM and \beta is the parameter set to measure the difference between the CFM and the DC.

Estimation is performed using a Bayesian MCMC procedure. Prior distributions for \Gamma (& \Lambda) are derived directly from the model coefficients (mean and variance covariance matrix) or the CFM. A bespoke MCMC routine is performed to estimate \beta. Please see '?mcmc' for more detials.

For the standard example where the DC contains information from only a single treatment, trt need not be specified. Where comparisons between the CFM and multiple treatments are require, a covariate of treamtne allocations must be specified sperately (using the 'trt' option).

Value

a object of class 'psc' with attributes model.type, the cleaned Dataset and the posterior distribution of the fitted model

Attributes include

• A 'cleaned' dataset including extracted components of the CFM and the cleaned DC included in the procedure

• An object defining the class of model (and therefore the procedure applied - see above)

• A matrix containing the draws of the posterior distributions

Examples

library(psc)
e4_data <-psc::e4_data
gemCFM <- psc::gemCFM
psc <- pscfit(gemCFM,e4_data)
summary(psc)

Personalised Synthetic Controls - summary

Description

A generic function to provide a summary of a 'psc' object obtained from pscfit.R

Usage

## S3 method for class 'psc'
summary(object, ...)

Arguments

Value

A summary of a psc object obtained using pscSumm and a copy of the pscfit object

Examples

library(psc)
e4_data <-psc::e4_data
gemCFM <- psc::gemCFM
psc <- pscfit(gemCFM,e4_data)
summary(psc)

Example model for a survival outcome

Description

A generated model with a survival endpoint and a cuymulative hazard function estimated using flexible parametric splines. Data for the model were synthetically generated and are based on a dataset to evaulate the use of Sorafenib in HCC akin to the PROSASH model

Usage

surv.mod

Format

A model of class 'flezsurvreg':

gamma

cumulative baseline hazard parameters

vi

vascular invasion

age60

patient age (centred at 60)

ecog

ECOG performance Status

logafp

AFP - log scale

alb

albumin

logcreat

Creatinine - log scale

allmets

metastesis

ageVasInv

centred age nested within vascular invasion

time

survival time

cen

censoring indicator

os

survival time

de

count

exapmple outcome for count data

trt

exapmple identifier for mulitple treatment comparisons

aet

Aetiology

Source

simulated

References

Using prognostic and predictive clinical features to make personalised survival prediction in advanced hepatocellular carcinoma patients undergoing sorafenib treatment. Berhane S, et al., Br J Cancer. 2019 Jul;121(2):117-124

A function to obtain survival estimates from a flexible parametric model

Description

This function provides basic survival estimates from a flexible parametric survival model

Usage

surv_fpm(DC_clean, beta = 0, s = NULL)

Arguments

Details

A fucntion which extracts survival probabilities from a flexsurvreg object

Value

a list of times and assoicated survival probabilities

Visualising Comparisons between a CFM and a DC

Description

The visComp function takes the data visualisations supplied as part of the CFM model and appends summaries of the equivalent datapoints from the Data Cohort.

Usage

visComp(CFM, DC, id = NULL)

Arguments

Value

a list of grobs for each model covariate