Title
metasim -- Tool to simulate a new two-arm clinical trial(s) or diagnostic accuracy study(ies) using results from current meta-analysis
Syntax
metasim varlist , n(#) es(numlist) var(numlist) type(string) [options]
options Description ------------------------------------------------------------------------- Main n(#) number of patients in group 1 of the new study (see below) es(numlist) pooled estimate/s from current meta-analysis var(numlist) variance/s for es type(string) type of study being simulated
Optional measure(string) outcome measure used in meta-analysis p(#) event rate or probability of being diseased dependent on (type) r(#) ratio of patients in two groups; treatment and control or diseased and healthy studies(#) number of new studies to be simulated model(string) meta-analysis model used on pre-existing data tausq(numlist) between study variance in current meta-analysis dist(string) distribution of effect sizes used to simulate the new study from corr(#) correlation between sensitivity and specificity (diagnostic study only) -------------------------------------------------------------------------
Description
metasim simulates a specified number of new studies based on the estimate/s obtained from a pre-existing meta-analysis. The user should input a maximum of 6 variables. The program will then generate a new set of variables and save it in a file called temppow within the working directory.
Options
+------+ ----+ Main +-------------------------------------------------------------
n(integer) relates to the number of patients in the new study. If simulating a new clinical trial, then this specifies the number of patients in the control group. If simulating a new diagnostic accuracy study then n is the number of diseased patients when using sensitivity and specificity as the outcome measure of accuracy, or the number of positive test results if using the diagnostic odds ratio (DOR).
es(numlist) specifies the pooled estimate/s from the meta-analysis of existing studies. If using the odds ratio, the diagnostic odds ratio or relative risk then ln(OR), ln(DOR) or ln(RR) estimates need to be specified. If using sensitivity and specificity then logit(sensitivity) and logit(specificity) estimates need to be specified in that order. var(numlist) specifies the variance/s for es. Need to specify two values if using sensitivity and specificity.
type(clinical/diagnostic) specifies the type of new study that the user would like to simulate; either a 2-arm clinical trial or a diagnostic test accuracy study.
+----------+ ----+ Optional +---------------------------------------------------------
measure(or/rr/rd/nostandard/dor/ss) specifies the outcome measure used in the meta-analysis to pool the results. The odds ratio (or), relative risk (rr), risk difference (rd) and unstandardised mean difference (nostandard) can only be used when simulating a new clinical study. The diagnostic odds ratio (dor) and sensitivity and specificity (ss) can only be used when simulating a new diagnostic accuracy study. The default for a clinical type study with 4 variables entered into the varlist is relative risk (rr, the default for a clinical type study with 6 variables entered into the varlist is unstandardised mean difference (nostandard and the default for a diagnostic type study is sensitivity and specificity (ss).
p(real) if simulating a new clinical study then this is the estimated event rate in the control group in the new study. When simulating a new diagnostic study this is the estimated probability of being diseased given a positive result in the new study. When this option is not specified by the user, the program will calculate this value by averaging the probabilities across the studies included in the dataset memory. Note that p is only relevant in the diagnostic framework when using the diagnostic odds ratio ( dor) as the option in measure.
r(real) is the ratio of patients in the control group to the treatment group when simulating a new clinical study. When simulating a new diagnostic accuracy study this is the ratio of diseased to healthy people if using sensitivity and specificity and the ratio of positive to negative results if using the DOR (default 1).
studies(integer) specifies the number of new studies to be simulated (default 1). When more than one are specified they are all assumed to have the same sample size.
model(fixed/fixedi/random/randomi/bivariate) defines the type of model used to meta-analyse the pre-existing data. The default is the fixed effect Mantel-Haenszel method (fixed) unless the outcome measure is the nonstandardised mean difference in which case the default is the inverse variance method (fixedi). The (fixedi) option specifies a fixed effect model using the inverse variance method. The (random) option uses the random effect DerSimonian & Laird method, taking the estimate for heterogeneity from the Mantel-Haenszel method. The (randomi) option specifies a random effects model using the method of DerSimonian and Laird, with the estimate of heterogeneity being taken from the inverse-variance fixed-effect model. All of the above options call on the metan command within the program. The final option is the bivariate random effects model (bivariate). This method calls on a combination of the metandi and midas commands. It may only be specified when simulating a new diagnostic accuracy study. tausq(numlist) is the measure of between study variance taken from the pre-existing meta-analysis (default 0). If the user specifies sensitivity and specificity (ss) as the measure then two values must be entered for tausq.
dist(normal/t) specifies the distribution of effect sizes used to sample a value from in order to simulate a new study(ies). The default for the (random) and (randomi) is a predictive distribution based on the t-distribution (t) allowing for heterogeneity between studies (and the uncertainty in the heterogeneity). The default for all other models is the (normal) distribution based on the mean and variance entered in es and var.
corr(real) is the correlation between the sensitivity and specificity (default 0). This option is only needed if the user chooses the bivariate model.
Examples
Simulating a new clinical trial with 100 patients in each group. The current fixed effect meta-analytic estimate of the pre-existing clinical trial has a pooled odds ratio of 2 (variance 0.5).
. metasim e_trt ne_trt e_ctrl ne_ctrl, model(fixed) measure(or) es(2) var(0.5) n(100) type(clinical)
Simulating a new clinical trial with 100 patients in the control group and 50 patients in the experimental treatment group. The current random effects meta-analytic estimate of the pre-existing clinical trail has a pooled relative risk of 1.5 (variance 0.2).
. metasim e_trt ne_trt e_ctrl ne_ctrl, model(random) measure(rr) es(1.5) var(0.2) tausq(0.1) n(100) type(clinical) r(0.5)
Simulating a new diagnostic test accuracy study assuming 80 patients in both the diseased group and the healthy group. A fixed effect meta-analysis of the pre-existing diagnostic test studies has a pooled logit sensitivity of 1.2 (variance 0.1) and a pooled logit specificity of 1.7 (variance 0.15).
. metasim TP FP FN TN, model(fixed) measure(ss) es(1.2 1.7) var(0.1 0.15) n(100) type(diagnostic)
Simulating a new diagnostic test accuracy study assuming 80 patients in both the diseased and the healthy group of the new study. A bivariate random effects meta-analysis of the pre-existing diagnostic test studies has a pooled logit sensitivity of 1.5 (variance 0.4) and a pooled logit specificity of 1.3 (variance 0.45).
. metasim TP FP FN TN, model(bivariate) measure(ss) es(1.5 1.3) var(0.4 0.45) n(80) type(diagnostic)
Authors
Michael J. Crowther, University of Leicester, United Kingdom. Email: michael.crowther@le.ac.uk.
Sally R. Hinchliffe, University of Leicester, United Kingdom. Email: srh20@le.ac.uk.
Alison Donald, University of Leicester, United Kingdom.
Alex J. Sutton, University of Leicester, United Kingdom. Email: ajs22@le.ac.uk.
References
Hinchliffe S, Crowther MJ, Phillips RS, Sutton AJ. Using meta-analysis to inform the design of subsequent studies of diagnostic test accuracy (Submitted)
Ferreira ML, Herbert RD, Crowther MJ, Verhagen A, Sutton AJ. When is another clinical trial justified? (Submitted)
Sutton AJ, Cooper NJ, Jones DR, Lambert PC, Thompson JR, Abrams KR. Evidence-based sample size calculations based upon meta-analysis. Statistics in Medicine 2007; 26:2479-2500.
Also see
Online: metapow, metapowplot