convergence.rate: Numerical estimation of the convergence rate of Rosenthal's...

View source: R/convergence.rate.R

Numerical estimation of the convergence rate of Rosenthal's fail-safe number NrR Documentation

Numerical estimation of the convergence rate of Rosenthal's fail-safe number Nr

Description

Numerical estimation of the convergence rate of Rosenthal's fail-safe number Nr.

Usage

convergence.rate(k = seq( 10, 5000, by = 10), R = 1000, alpha = 0.05)

Arguments

k

A grid of number of studies to consider.

R

The number of repeats for each number of studies.

alpha

The significance level, set to 0.05 by default.

Details

This function replicates the Figures 6 and 5 in Fragkos, Tsagris and Frangos (2017).

Value

Two plots, the absolute relative error of Nr agains the number of studies and the logarithm of absolute relative error against the logarithm of the number of studies and the coefficients of the regression model of the second plot. The second coefficient is the numerically estimated convergence rate of Nr.

Author(s)

Michail Tsagris and Constantinos Frangos

R implementation and documentation: Michail Tsagris mtsagris@uoc.gr and Constantinos Frangos kfragkos@outlook.com.

References

Konstantinos C. Fragkos, Michail Tsagris and Christos C. Frangos (2017). Exploring the distribution for the estimator of Rosenthal's "fail-safe" number of unpublished studies in meta-analysis. Communications in Statistics-Theory and Methods, 46(11):5672–5684.

See Also

den.plot, truncnorm.nr.density, rosenthal

Examples

convergence.rate(k = seq( 50, 500, by = 50), R = 1000, alpha = 0.05)

fsn documentation built on March 18, 2022, 6:02 p.m.