print.copas: Print results of Copas selection model

In metasens: Statistical Methods for Sensitivity Analysis in Meta-Analysis

Print results of Copas selection model

Description

Print method for objects of class copas.

This function prints the main results of a Copas analysis, performed using the function copas. It complements the graphical summary of the results, generated using plot.copas.

Specifically it prints a table where the:

first column corresponds to the x-axis in plots 3 & 4 from plot.copas;

second column corresponds to the treatment effect displayed in plot 3 from plot.copas;

third and fourth columns give the confidence intervals for this treatment effect,

fifth colum gives the p-value for an overall treatment effect,

sixth column gives the p-value for residual publication bias (the y-axis of plot 4 from plot.copas (see plot.copas under plot 4 for a further explanation of this p-value))

seventh column gives an approximate estimate of the number of studies the model suggests remain unpublished if the probability of publishing the study with the largest SE is as in column 1.

Below this is displayed the results of the Copas analysis (Adjusted estimate) for the smallest degree of selection for which the p-value for evidence of residual selection bias exceeds sign.rsb (default: 0.1). This is simply extracted from the corresponding row in the table above.

Lastly, the unadjusted random effects estimate and 95% confidence interval is printed.

Usage

## S3 method for class 'copas'
print(
  x,
  backtransf = x$backtransf,
  digits = gs("digits"),
  digits.pval = max(gs("digits.pval"), 2),
  digits.prop = gs("digits.prop"),
  digits.tau2 = gs("digits.tau2"),
  digits.tau = gs("digits.tau"),
  scientific.pval = gs("scientific.pval"),
  big.mark = gs("big.mark"),
  header = TRUE,
  legend = TRUE,
  text.tau2 = gs("text.tau2"),
  text.tau = gs("text.tau"),
  ...
)

Arguments

x	An object of class copas.
backtransf	A logical indicating whether results should be back transformed in printouts and plots. If backtransf = TRUE (default), results are printed as odds ratios rather than log odds ratio, for example.
digits	Minimal number of significant digits, see print.default.
digits.pval	Minimal number of significant digits for p-value of overall treatment effect, see print.default.
digits.prop	Minimal number of significant digits for proportions, see print.default.
digits.tau2	Minimal number of significant digits for between-study variance τ^2, see print.default.
digits.tau	Minimal number of significant digits for τ, the square root of the between-study variance τ^2.
scientific.pval	A logical specifying whether p-values should be printed in scientific notation, e.g., 1.2345e-01 instead of 0.12345.
big.mark	A character used as thousands separator.
header	A logical indicating whether information on title of meta-analysis, comparison and outcome should be printed at the beginning of the printout.
legend	A logical indicating whether a legend should be printed.
text.tau2	Text printed to identify between-study variance τ^2.
text.tau	Text printed to identify τ, the square root of the between-study variance τ^2.
...	Additional arguments (ignored).

Author(s)

James Carpenter [e-mail], Guido Schwarzer [e-mail]

See Also

copas, plot.copas, summary.copas

Examples

data(Fleiss1993bin, package = "meta")

# Perform meta analysis, effect measure is odds ratio (OR)
#
m1 <- metabin(d.asp, n.asp, d.plac, n.plac, data = Fleiss1993bin, sm = "OR")

# Perform Copas analysis
#
cop1 <- copas(m1)
cop1

metasens documentation built on March 7, 2023, 7:51 p.m.