meta.exact: Exact fixed effect meta analysis

In exactmeta: Exact fixed effect meta analysis

Description

Compute the exact confidence interval and p-value for risk difference, risk ratio, rate difference or rate ratio in meta analysis using fixed effect model.

Usage

meta.exact(data, type="risk difference", BB.grdnum=2000, B.sim=20000, cov.prob=0.95, 

midp=T, print=T, studyCI=T, ratio.upper=1000)

Arguments

data	The data matrix summarizing information from multiple studies needed in the meta analysis. The first column consists of the numbers of events in arm 1; the second column consists of the numbers of events in arm 2; the third column consists of the sample sizes or total follow up times in arm 1 and the fourth column consists of the sample sizes or total follow up times in arm 2. Each row contains data from a study.
type	The parameter of interest. It can be “risk difference" (the difference of two binomial proportions); "risk ratio" (the ratio of two binomial proportions); "rate difference" (the difference of two Poission rates) and "rate ratio" (the ratio of two Poisson rates).
BB.grdnum	The number of grid points used in searching the ends of the confidence interval. A bigger value increases the precision of the constructed confidence interval.
B.sim	The number of Monte-carlo simulations used to generate the null-distribution based on which the confidence interval and p-value are obtained.
cov.prob	The coverage probability of the confidence interval to be constructed.
midp	The binary indicator for if the mid-p values are used in calculating the exact p-values and constructing the exact confidence intervals. The results are more conservative wihout using mid-p values.
print	The binary indicator for if the computing progress is printed during the computation.
studyCI	The binary indicator for if the confidence interval and p-value for individual study are reported.
ratio.upper	A value indicating the limits in computing the confidence interval for either the risk ratio or rate ratio. The default value of 1000 means that the widest confidence interval for the ratio is (0.001, 1000).

Value

ci.fixed	A matrix reporting the point estimators, confidence intervals, and p-values from three methods: constant weighting, inverse-variance weighting and Fisher's weighting. Normally, the result from the inverse-variance weighting method is recommended. It also reports the corresponding result from the asymptotical MH method with or without continuity correction as specified by users. When type is either "risk ratio" or "rate ratio", the continuity correction are automatically made for zero event. The last column of the matrix provides the range used for searching the confidence limits. When a confidence limit is the same as the range limit, the true confidence interval goes beyond the specified range.
study.ci	When "studyCI=T", "study.ci" is a matrix reporting the point estimators (the first column), confidence intervals (the second and thrid columns), p values (the fourth column) and if the confidence interval goes beyond the pre-specified range: 1 if within the range and 0 otherwise (the fifth column).
precision	The precision used to search the limit of the confidence interval. It can be improved by increasing "BB.grdnum" in the arguments.

Author(s)

Yilei Yu

References

Tian, L., Cai, T., Pfeffer, M., Piankov, N., Cremieux, P., and Wei, LJ. (2008). Exact and effcient inference procedure for meta-analysis and its application to the analysis of independent two by two tables with all available data but without artificial continuity correction. Biostatistics, 10(2):275-281.

Examples

data=cbind(c(2, 1), c(1, 1), c(10, 10), c(12, 12))

meta.exact(data, type="risk difference", BB.grdnum=100, studyCI=FALSE)
meta.exact(data, type="rate difference", BB.grdnum=100, studyCI=FALSE)
meta.exact(data, type="risk ratio", BB.grdnum=100, studyCI=FALSE)
meta.exact(data, type="rate ratio", BB.grdnum=100)

Example output

study= 1 
study= 2 
$ci.fixed
            constant inverse-variance     fisher asymptotical-MH  range
est      -0.08163265      -0.08163265 -0.1020408     -0.06666667     NA
lower CI -0.26530612      -0.26530612 -0.2653061     -0.25702418     -1
upper CI  0.14285714       0.14285714  0.1428571      0.12369085      1
p         0.76210000       0.73260000  0.7036000      0.49245141     NA

$study.ci
NULL

$precision
[1] "+/- 0.02"

study= 1 
study= 2 
$ci.fixed
            constant inverse-variance     fisher asymptotical-MH     range
est      -0.08159091       -0.1087879 -0.1087879       0.9355070        NA
lower CI -0.29916668       -0.2991667 -0.2991667      -0.2719697 -1.359849
upper CI  0.16318183        0.1631818  0.1631818       0.1386364  1.359849
p         0.64370000        0.6222000  0.5953000       0.5244859        NA

$study.ci
NULL

$precision
[1] "+/- 0.0269276941163715"

study= 1 
study= 2 
$ci.fixed
          constant inverse-variance    fisher asymptotical-MH range
est      0.5754399        0.5754399 0.5754399       0.5555556    NA
lower CI 0.1258925        0.1258925 0.1258925       0.1033767 1e-03
upper CI 3.0199517        3.0199517 3.0199517       2.9856039 1e+03
p        0.5803000        0.5607000 0.5245000       0.4932867    NA

$study.ci
NULL

$precision
[1] "+/- 0.136787233247171"

study= 1 
study= 2 
$ci.fixed
         constant inverse-variance  fisher asymptotical-MH  range
est       0.57544          0.57544 0.57544    5.555600e-01     NA
lower CI  0.09550          0.09550 0.09550    0.000000e+00  1e-03
upper CI  3.46737          3.46737 3.46737    1.235196e+07  1e+03
p         0.51630          0.49790 0.46680    9.457100e-01     NA

$study.ci
               est   lower CI  upper CI         p
study  1 0.4166667 0.01412729  5.478235 0.5259204
study  2 0.8333333 0.02136752 32.500000 0.9090909

$precision
[1] "+/- 0.136787233247171"

exactmeta documentation built on May 2, 2019, 3:24 a.m.