The aim of the snSMART R package is to consolidate data simulation, sample size calculation and analysis functions for several snSMART (small sample sequential, multiple assignment, randomized trial) designs under one library.
An snSMART is a multi-stage trial design where for a two-stage design, randomization in the second stage depends on the outcome to first stage treatment. snSMART designs require that the same outcome is measured at the end of the first stage and at the end of the second stage. Additionally, the length of the first stage of the trial must be the same amount of time as the length for the second stage. snSMARTs are motivated by obtaining more information from a small sample of individuals with the primary goal to identify the superior first stage treatment or dosage level using both stages of data. Data are shared across the two stages of the snSMART design to more precisely estimate the effect of the treatments given in the first stage.
Before using this package, please install the JAGS library on your device.
You can install the snSMART package from CRAN:
install.packages("snSMART", repos = "http://cran.us.r-project.org")##
## The downloaded binary packages are in
## /var/folders/5_/sp2r7r_s5snf4hq634xmk_0r0000gn/T//RtmpUpvbHg/downloaded_packages
library(snSMART)## Loading required package: EnvStats
##
## Attaching package: 'EnvStats'
## The following objects are masked from 'package:stats':
##
## predict, predict.lm
Or get the development version from GitHub:
# Install devtools first if you haven't done so
library(devtools)
# install snSMART
devtools::install_github("sidiwang/snSMART")
library(snSMART)- snSMART with binary outcome (3 active treatments or placebo and 2 dose
level, non-responders re-randomized)
- Bayesian Joint Stage Model (BJSM) analysis function -
BJSM_binary - Log Poisson Joint Stage Model (JSRM) analysis function -
JSRM_binary - sample size calculation function -
sample_size - Group Sequential snSMART BJSM analysis function -
group_seq
- Bayesian Joint Stage Model (BJSM) analysis function -
- snSMART with mapping function (3 active treatments, re-randomization
depends on continuous outcome at stage 1; continuous outcomes)
- BJSM analysis function -
BJSM_c
- BJSM analysis function -
BJSM_binary: We call the BJSM_binary function using data from an
snSMART with 30 total individuals. We assumed a six beta model with the
priors for
mydata <- data_binary
BJSM_result <- BJSM_binary(
data = mydata, prior_dist = c("beta", "beta", "pareto"),
pi_prior = c(0.4, 1.6, 0.4, 1.6, 0.4, 1.6), beta_prior = c(1.6, 0.4, 3, 1),
n_MCMC_chain = 1, n.adapt = 1000, MCMC_SAMPLE = 20, ci = 0.95,
six = TRUE, DTR = TRUE, verbose = FALSE
)
summary(BJSM_result)##
## Treatment Effects Estimate:
## Estimate Std. Error C.I. CI low CI high
## trtA 0.4230557 0.10984677 0.95 0.2132406 0.6239648
## trtB 0.4621068 0.09676214 0.95 0.2859294 0.6378900
## trtC 0.5949919 0.04873058 0.95 0.5357735 0.6742898
##
## Differences between Treatments:
## Estimate Std.Error C.I. CI low CI high
## diffAB -0.03905106 0.1733779 0.95 -0.3686263 0.22405873
## diffBC -0.13288514 0.1246969 0.95 -0.3755066 0.07647332
## diffAC -0.17193620 0.1133842 0.95 -0.3321949 0.04989511
##
## Linkage Parameter Estimate:
## Estimate Std. Error C.I. CI low CI high
## beta0A 0.9377270 0.07330169 0.95 0.7715834 0.9999638
## beta0B 0.8822049 0.08253566 0.95 0.7357460 0.9998799
## beta0C 0.7042976 0.26361834 0.95 0.1427238 0.9771399
## beta1A 1.3582895 0.27307278 0.95 1.0028447 1.9003715
## beta1B 1.3271310 0.29322491 0.95 1.0072391 2.2666544
## beta1C 1.3889213 0.19685318 0.95 1.0015164 1.7286500
##
## Expected Response Rate of Dynamic Treatment Regimens (DTR):
## Estimate Std. Error C.I. CI low CI high
## rep_AB 0.5003874 0.07678636 0.95 0.3585359 0.6151457
## rep_AC 0.5671207 0.07019227 0.95 0.4567067 0.6971029
## rep_BA 0.4983502 0.09710728 0.95 0.3015048 0.6378734
## rep_BC 0.5780166 0.09550426 0.95 0.4319781 0.7034776
## rep_CA 0.6077879 0.10240957 0.95 0.4230910 0.7724611
## rep_CB 0.6205831 0.08853430 0.95 0.4310405 0.8091465
LPJSM_binary: Here, we call the LPJSM_binary mirroring our example for
the BJSM_binary above.
LPJSM_result <- LPJSM_binary(data = data_binary, six = TRUE, DTR = TRUE)
summary(LPJSM_result)##
## GEE output:
##
## Call:
## geepack::geeglm(formula = Y ~ alphaA + alphaB + alphaC + gamma1A +
## gamma2A + gamma1B + gamma2B + gamma1C + gamma2C - 1, family = poisson(link = "log"),
## data = geedata, id = ptid, corstr = "independence")
##
## Coefficients:
## Estimate Std.err Wald Pr(>|W|)
## alphaA -1.2155 0.4030 9.096 0.00256 **
## alphaB -0.9845 0.3403 8.370 0.00381 **
## alphaC -0.8444 0.3036 7.737 0.00541 **
## gamma1A 1.2155 0.4030 9.096 0.00256 **
## gamma2A 0.7029 0.3664 3.680 0.05508 .
## gamma1B 0.9845 0.3403 8.370 0.00381 **
## gamma2B 0.8271 0.3497 5.594 0.01803 *
## gamma1C 0.8444 0.3036 7.737 0.00541 **
## gamma2C -42.2449 0.6085 4820.409 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation structure = independence
## Estimated Scale Parameters:
##
## Estimate Std.err
## (Intercept) 0.4096 0.1259
## Number of clusters: 30 Maximum cluster size: 2
##
## Treatment Effect Estimate:
## Estimate Std. Error
## trtA 0.2966 0.1195
## trtB 0.3736 0.1271
## trtC 0.4298 0.1305
##
## Expected Response Rate of Dynamic Treatment Regimens (DTR):
## Estimate Std. Error
## rep_AB 0.8274 0.1328
## rep_AC 0.9072 0.2341
## rep_BA 0.7984 0.1403
## rep_BC 0.9892 0.1427
## rep_CA 0.4298 0.1305
## rep_CB 0.4298 0.1305
sample_size: In this example, we call the function to request the
sample size needed per arm with the following assumptions: the response
rates for treatments A, B, and C are 0.7, 0.5 and 0.25, respectively;
sampleSize <- sample_size(
pi = c(0.7, 0.5, 0.25), beta1 = 1.4, beta0 = 0.5,
coverage = 0.9, power = 0.8, mu = c(0.65, 0.55, 0.25),
n = c(4, 2, 3)
)## With given settings, the estimated sample size per arm for an snSMART is: 34
## This implies that for an snSMART with sample size of 34 per arm (102 in total for three treatments):
## The probability of successfully identifying the best treatment is 0.8 when the difference of response rates between the best and second best treatment is at least 0.2, and the response rate of the best treatment is 0.7
group_seq: This function either outputs which treatment arm should be
dropped if or provides a full BJSM analysis based on the complete
dataset if . The dataset and are provided in the package for
illustration purpose.
result1 <- group_seq(
data = groupseqDATA_look1, interim = TRUE, drop_threshold_pair = c(0.5, 0.4),
prior_dist = c("beta", "beta", "pareto"), pi_prior = c(0.4, 1.6, 0.4, 1.6, 0.4, 1.6),
beta_prior = c(1.6, 0.4, 3, 1), MCMC_SAMPLE = 6000, n.adapt = 1000, n_MCMC_chain = 1
)##
## Interim Analysis Outcome:
## Threshold tau_l is set to:
## 0.5
##
## Threshold psi_l is set to:
## 0.4
##
## Step 1: Arm C's interim posterior probability of having the greatest response is bigger than threshold
## 0.5
##
## Step 2: Arm A's interim posterior probability of having the lowest response is higher
## Arm A is dropped
##
result2 <- group_seq(
data = groupseqDATA_full, interim = FALSE, prior_dist = c("beta", "beta", "pareto"),
pi_prior = c(0.4, 1.6, 0.4, 1.6, 0.4, 1.6), beta_prior = c(1.6, 0.4, 3, 1),
MCMC_SAMPLE = 60000, BURN.IN = 10000, n_MCMC_chain = 1, ci = 0.95, DTR = TRUE
)
summary(result2)##
## Treatment Effects Estimate:
## Estimate Std. Error C.I. CI low CI high
## trtA 0.3010 0.04901 0.95 0.2064 0.3973
## trtB 0.4734 0.03969 0.95 0.3967 0.5518
## trtC 0.6829 0.04032 0.95 0.6039 0.7587
##
## Differences between Treatments:
## Estimate Std.Error C.I. CI low CI high
## diffAB -0.1724 0.06288 0.95 -0.2934 -0.04701
## diffBC -0.2095 0.05548 0.95 -0.3176 -0.10193
## diffAC -0.3819 0.06203 0.95 -0.5016 -0.25930
##
## Linkage Parameter Estimate:
## Estimate Std. Error C.I. CI low CI high
## beta0A 0.8687 0.11044 0.95 0.6615 1.0000
## beta0B 0.7141 0.11250 0.95 0.4944 0.9419
## beta0C 0.8718 0.12633 0.95 0.6187 1.0000
## beta1A 1.4751 0.36568 0.95 1.0000 2.1793
## beta1B 1.3613 0.16992 0.95 1.0266 1.6774
## beta1C 1.4494 0.08872 0.95 1.2860 1.6265
BJSM_c: Below, we call the function assuming the mean and standard
deviation of the normal prior being 50 and 50 for all three treatments,
and the standard deviation of the prior distribution of
BJSM_result <- BJSM_c(
data = trialDataMF, xi_prior.mean = c(50, 50, 50), xi_prior.sd = c(50, 50, 50),
phi3_prior.sd = 20, n_MCMC_chain = 1, n.adapt = 1000, MCMC_SAMPLE = 5000, BURN.IN = 1000,
ci = 0.95, n.digits = 5
)
summary(BJSM_result)##
## Parameter Estimation:
## Estimate CI CI_low CI_high
## V1[1,1] 111.4730 0.95 6.127e+01 170.443
## V1[2,1] 85.5731 0.95 4.313e+01 131.885
## V1[1,2] 85.5731 0.95 4.313e+01 131.885
## V1[2,2] 77.7046 0.95 4.287e+01 119.640
## V2[1,1] 123.2901 0.95 6.477e+01 194.723
## V2[2,1] 29.2337 0.95 -1.714e+01 78.510
## V2[1,2] 29.2337 0.95 -1.714e+01 78.510
## V2[2,2] 107.6312 0.95 5.585e+01 171.342
## phi1 0.1817 0.95 2.822e-04 0.387
## phi3 4.0143 0.95 2.807e+00 5.400
## xi_[A] 51.1199 0.95 4.716e+01 54.890
## xi_[B] 62.0682 0.95 5.819e+01 65.932
## xi_[C] 68.9650 0.95 6.520e+01 72.650
This R package will continue to be updated as more snSMART designs and methods are developed. We hope that this package translates snSMART design and methods into finding more effective treatments for rare disease.
Chao, Y.C., Trachtman, H., Gipson, D.S., Spino, C., Braun, T.M. and Kidwell, K.M., 2020. Dynamic treatment regimens in small n, sequential, multiple assignment, randomized trials: An application in focal segmental glomerulosclerosis. Contemporary clinical trials, 92, p.105989.
Chao, Y.C., Braun, T.M., Tamura, R.N. and Kidwell, K.M., 2020. A Bayesian group sequential small n sequential multiple‐assignment randomized trial. Journal of the Royal Statistical Society: Series C (Applied Statistics), 69(3), pp.663-680.
Fang, F., Hochstedler, K.A., Tamura, R.N., Braun, T.M. and Kidwell, K.M., 2021. Bayesian methods to compare dose levels with placebo in a small n, sequential, multiple assignment, randomized trial. Statistics in Medicine, 40(4), pp.963-977.
Hartman, H., Tamura, R.N., Schipper, M.J. and Kidwell, K.M., 2021. Design and analysis considerations for utilizing a mapping function in a small sample, sequential, multiple assignment, randomized trials with continuous outcomes. Statistics in Medicine, 40(2), pp.312-326.
Wei, B., Braun, T.M., Tamura, R.N. and Kidwell, K., 2020. Sample size determination for Bayesian analysis of small n sequential, multiple assignment, randomized trials (snSMARTs) with three agents. Journal of Biopharmaceutical Statistics, 30(6), pp.1109-1120.
Wei, B., Braun, T.M., Tamura, R.N. and Kidwell, K.M., 2018. A Bayesian analysis of small n sequential multiple assignment randomized trials (snSMARTs). Statistics in medicine, 37(26), pp.3723-3732.