Abstract
 The objective of this paper is to describe general approaches of diagnostic test accuracy (DTA) that are available for the quantitative synthesis of data using R software. We conduct a DTA that summarizes statistics for univariate analysis and bivariate analysis. The package commands of R software were “metaprop” and “metabin” for sensitivity, specificity, and diagnostic odds ratio; forest for forest plot; reitsma of “mada” for a summarized receiveroperating characteristic (ROC) curve; and “metareg” for metaregression analysis. The estimated total effect sizes, test for heterogeneity and moderator effect, and a summarized ROC curve are reported using R software. In particular, we focus on how to calculate the effect sizes of target studies in DTA. This study focuses on the practical methods of DTA rather than theoretical concepts for researchers whose fields of study were nonstatistics related. By performing this study, we hope that many researchers will use R software to determine the DTA more easily, and that there will be greater interest in related research.

Keywords: Metaanalysis; Diagnostic test accuracy; Receiveroperating characteristic curve; Likelihood ratios; Mada; Reitsma
INTRODUCTION
 General pairwise metaanalysis calculates the effect size, such as relative risk and odds ratio (OR) for binary data and the mean difference for continuous data. By contrast, the diagnostic test accuracy (DTA) simultaneously combines two effect sizes, such as the sensitivity and specificity or positive predictive value and negative predictive value [13].
 Therefore, DTA is more complex than pairwise metaanalysis, which has one result value. The expansion to multivariate analysis with more than two results inevitably leads to the introduction of the multilayer concept, which requires some degree of mathematical understanding as well as an ability to use statistical programs.
 This study focuses on the procedures involved in running the R software (Supplementary Material 1) as well as the concepts of producing summary statistics, which need to be understood for DTA.
 In this study, the previous metaanalysis studies performed by the authors [13] are reviewed using R software. Furthermore, this study requires prior knowledge about the metaanalysis of diagnostic tests because it first deals with the types and changes of the effect size to calculate the summary statistics for DTA.
UNDERSTANDING DIAGNOSTIC TEST ACCURACY
 The data for DTA assumes a 2×2 table form in which the row cells are distinguished by the presence or absence of a test, and the column cells are distinguished by the presence or absence of a disease (Figure 1).
 Summary statistics for diagnostic test accuracy
 The DTA is represented by the summary statistics and summary line from four sets of basic data, namely true positive (TP), false positive (FP), false negative (FN), and true negative (TN). Representative summary statistics are the sensitivity, specificity, diagnostic odds ratio (DOR), and forest plot, and an example summary curve is the summary receiver operating characteristic (SROC) curve (Table 1).
 Diagnostic test accuracy model
 To calculate the summary statistics for the DTA, an appropriate model should be selected, as with pairwise metaanalysis. Models that simultaneously consider the sensitivity and specificity include the Moses–Littenberg SROC model [4,5], bivariate model [6], and hierarchical SROC (HSROC) model [7].
 The Moses–Littenberg model is a simple model that was created in the early stage to determine the DTA, and it estimates the SROC using simple linear regression. This is similar to the fixedeffect model in pairwise metaanalysis, and cannot estimate the heterogeneity between studies. Furthermore, this model cannot distinguish between withinstudy and betweenstudy variations in all variations, and can perform limited analysis because it only provides the SORC curve without parameter estimates, standard deviation, or confidence intervals (CIs).
 To overcome the disadvantages of the Moses–Littenberg model, the bivariate model and HSROC model were developed based on the hierarchical model. These two models provide the same value mathematically when there is no covariate [8,9]. This is similar to the randomeffect model in pairwise metaanalysis. Both models can estimate the withinstudy and betweenstudy variation of studies, that is, the heterogeneity.
 The bivariate model assumes a binominal distribution that directly models the sensitivity and specificity for withinstudy variations, while assuming a bivariate normal distribution for betweenstudy variation. However, the HSROC model assumes a binominal distribution for withinstudy variations, while assuming a hierarchical distribution for parameters included in the logistic model by applying the logistic regression model to determine the probability of a binominal distribution for betweenstudy variation.
 The R “mada” package reitsma model, which we will practice in this book, calculates the summary statistics and estimates the SROC curve using the bivariate model by default.
 Calculation of effect size
 Examine the sensitivity and specificity in Table 1.
 The sensitivity is TP/(TP+FN), and the specificity is TN/(TN+FP), which are proportiontype data.
 For these proportiontype data, logittransformed data are used more often than raw data. The logit transformation is a method of adjusting the data distribution according to statistical assumptions. The proportiontype data are limited between the lower and upper limits of 0 and 1, respectively. To convert these data to make them appropriate for the assumptions of statistics, their upper and lower limits should be released by performing multiplication and log transformations, respectively. This is called logit transformation.
 Upon completion of the calculation of the summary statistics for DTA, they are reverted to their original values for interpretation. In the practice using R, we will calculate the logittransformed sensitivity and specificity using the “metaprop” function of the “meta” package and then revert them to their original values to interpret them. Thus, we should understand why the effect size is transformed.
DIAGNOSTIC TEST ACCURACY USING THE “mada” AND “meta” PACKAGES OF R

Figure 2 shows the flow of the DTA. First, when coding the data, we must change the variable name appropriately for the corresponding function. After selecting a metaanalysis model (fixed or random), the total effect size is presented, the heterogeneity is verified, and the publication bias is then verified and reported.
 The “mada” package is required to analyze the DTA in R. After “mada” is installed, you will be promoted to install “mvtnorm,” “ellipse,” and “mvmeta.” Thus, you should install them in advance as follows:
 ·install.packages(“mada”)
 ·install.packages(“mvtnorm”)
 ·install.packages(“ellipse”)
 ·install.packages(“mvmeta”)
 In addition, you should install the “meta,” “metafor,” and “rmeta” packages for general pairwise intervention metaanalysis in R as follows:
 ·install.packages(“meta”)
 ·install.packages(“metafor”)
 ·install.packaqes(“rmeta”)
 The main explanations are applicable to the “mada” and “meta” packages. For detailed explanations about the “mada” package, refer to detailed codes, documents, and references for the package [10].
 We mark R commands with a dot (‘·’) in front of them, to distinguish them from the main text. When long commands are extended to the next line, there is no dot at the beginning of the next line. Thus, when you enter the command in the R software, you must type them without the dot (‘·’) in front of them.
 Data coding and loading
 As an example of the DTA, the urine sample measuring the albumin concentration method (Table 2) was selected from among the test methods for microalbuminuria in diabetes patients [2,3,11]. Subgroup (g) 1 consists of Western European counties, and 0 consists of countries other than Western European countries.
 Load the example file from the working folder with the following command. Note that R prefers commaseparated values (csv) file format. Thus, you should save Table 2 as “dta_shim.csv” in the specified working folder.
 ·dta_shim < read.csv (“dta_shim.csv”, header=TRUE)
 read.csv is a function for loading a csv file. The above command means to load the file “dta_shim.csv” and use the first variable name (header=TRUE). This loaded file is saved as “dta_shim” in the R memory. To confirm this, enter the specified data in the View() function.
 Summary statistics
 The “mada” package, which is a bivariate model for calculating the summary statistics for the DTA, does not provide the total effect sizes of summary statistics (sensitivity, specificity, and DOR) and only provides the effect size of individual studies as a forest plot, which is inconvenient.
 Therefore, it is more natural to check the value of each summary statistic by performing univariate analysis using the “meta” package first, and then to present an SROC curve using the “mada” package.
 Univariate analysis
 Calculate the sensitivity, specificity, and DOR and plot them using the univariate analysis model.
 Load the meta package to perform metaanalysis:
 ·library(meta)
Sensitivity
 The “meta” package includes many functions. Among them, the “metaprop” function calculates the total effect size using the number of events (event) and the number of samples (n) from proportiontype data.
 · sensitivity_logit < metaprop(dta_shim$TP, dta_shim$TP+ dta_shim$FN, comb.fixed=FALSE, comb.random=TRUE, sm=“PLOGIT”, method.ci=“CP”, studlab=dta_shim$id, byvar=dta_shim$g)
 ·print(sensitivity_logit, digits=3)
 In sensitivity analysis, the number of events is TP and the number of samples is TP+FN. The variables of these data in R can be indicated by using the symbol ‘$’ (for example, write “dta_shim$TP” to indicate the TP variable of the dta_shim data). After sequentially entering the number of events (dta_shim$TP) and the number of samples (dta_shim$TP+dta_shim$FN) in the metaprop function, input other optional arguments at the end.
 To calculate the effect size from proportiontype data, the method of reverting after logit transformation was used. Besides, you can enter sm=“PRAW” to use the raw data without transformation, or sm=“PLN” to find the reverted value after log transformation.
 For consistency with the assumptions of the statistic model, and to consider the symmetricity and distribution of data, it is desirable to transform proportiontype data (log transformation or logit transformation) as they produce conservative results. However, many previous studies and statistical models have broadened the operation scope of researchers. Thus, it is necessary to find and use an appropriate method for the research results.
 Even if data transformation is performed, the “metaprop” function automatically reverts and shows the total effect size that can be interpreted.
 Furthermore, there are several methods for calculating the confidence interval, but the default Clopper–Pearson method is recommended as it is not too complex (method.ci=“CP”).
 The random effect model was used, and comb.fixed=FALSE and comb.random=TRUE are also entered. The desired model can be selected by using FALSE or TRUE.
 Studlab=study indicates the name of individual studies. To show the result by subgroup, enter “byvar=g” where g is the variable name representing the subgroup. The results obtained when using the “metaprop” function are assigned to sensitivity_logit, and the result is shown in Figure 3.
 The results from sensitivity_logit in Figure 3 are examined below onebyone.
 ① Shows the total effect size of all nine studies. The proportion of the random effect model was 0.841 (95% CI, 0.788 to 0.882).
 ② Shows the result corresponding to the subgroup. The random model shows slight differences in sensitivity according to the subgroup (0 vs. 1). Based on the random effect model, the proportion is 0.816 for Western Europe countries and 0.855 for other countries. These values need to be tested using meta regression analysis according to country group later.
 ③ Shows the heterogeneity of all studies. The Higgins’ I^{2} of the heterogeneity is determined by subtracting the number of degrees of freedom from the Cochrane Q statistics, and then again dividing the resulting value by the Cochrane Q statistics. Thus, it quantifies the heterogeneity in a consistent manner. Values between 0% and 40% indicate that the heterogeneity may not be important; values between 30% and 60% indicate moderate heterogeneity; values between 50% and 90% indicate substantial heterogeneity; and values between 75% and 100% indicate considerable heterogeneity. The pvalue of the Cochrane Q statistics is 0.1, which is a somewhat wide range of significance [3].
 In this sensitivity analysis, the Higgins’ I^{2} is 32.5% and the Cochrane Q statistics pvalue is 0.158, which suggest weak heterogeneity.
 In addition, the calculation process for the results is revealed at the bottom of Figure 3. The inverse variance method is a basic metaanalysis method, and uses the inverse variance of the relevant study when calculating the weights of individual studies. The DerSimonianLaird estimator indicates that the tau value was used when calculating the betweenstudy variance.
 Furthermore, logit transformation and Clopper–Pearson method were used.
■ Forest plot
· forest(sensitivity_logit, digits=3, rightcols=c(“effect”, “ci”), xlab=“Sensitivity”)
 Enter the corresponding metaanalysis model (sensitivity_logit) in the forest function. Then, various options can be entered to facilitate identification. “digits=3” indicates that it shows only down to three decimal places, and “rightcols=c(“effect,” “ci”))” indicates that it shows the effect size and CI while omitting only the weight at the right side of the forest plot.
 In addition, the addition of colors or the addition/removal of certain information is possible at one’s discretion. You can learn more details by practicing the meta package.
 The forest plot provides the same information as the abovementioned total effect size. Furthermore, withinstudy and betweenstudy variation can be easily identified by the graphic representation of the effect size of individual studies.
 For example, it can be seen that Gansevoort, Ng, Wiegmann, and Ahn have large withinstudy variations, and Wiegmann and Incerti have large betweenstudy variations.
Specificity
· specificity_logit < metaprop(dta_shim$TN, dta_shim$TN+ dta_shim$FP, comb.fixed=FALSE, comb.random=TRUE, sm=“PLOGIT”, method.ci=“CP”, studlab=dta_shim$id, byvar=dta_shim$g)
·print(specificity_logit, digits=3)
 In the specificity analysis, the number of events is TN and the number of samples is TN+FP. After sequentially entering the number of events (dta_shim$TN) and the number of samples (dta_shim$TN+dta_shim$FP) in the metaprop function, input other optional arguments, respectively. The explanation after this is identical to that of the sensitivity analysis.
 We will examine the results of specificity_logit onebyone.
 The total effect size of all nine studies is shown. The proportion of the random effect model was 0.861 (95% CI, 0.794 to 0.909).
 The random model shows almost no difference in terms of the effect size between the subgroup (0 vs. 1). The Higgins’ I^{2} in this specificity analysis is 78.3%, and the pvalue of Cochrane Q statistics is <0.0001, which indicates the existence of heterogeneity.
■ Forest plot
· forest(specificity_logit, digits=3, rightcols=c(“effect”, “ci”), xlab=“Specificity”)
 The explanation for this command is the same as that for the sensitivity analysis.
Diagnostic odds ratio
 The “meta” package includes several functions. Among them, the “metabin” function calculates the total effect size from binary data when there exist all of the raw data. The respective sensitivity and specificity are proportiontype data, but the DOR of the 2×2 format is binary data.
 · DOR_model < metabin(TP,TP+FP,FN,FN+TN, sm=”OR”, comb.fixed=FALSE,comb.random=TRUE, method=“Inverse,” id, byvar=g, data=dta_shim)
 ·print(DOR_model)
 For binary data, enter TP, TP+FP, FN, and FN+TN in this order.
 Write OR for effect size (sm=“OR”) and use the general inverse variance method for weights of individual studies (method=“Inverse”).
 To set the random effect model considering the betweenstudy variations, additionally enter “comb.fixed=FALSE” and “comb.random=TRUE.”
 “id” indicates the name of the individual study, and “data=dta_shim” specifies the data “dta_shim” loaded to the R memory. To show the result for each g, enter “byvar=g,” where g is the variable name for the g. The results of the metabin function are assigned to the DOR model.
■ Forest plot
· forest(DOR_model, digits=3, rightcols=c(“effect”, “ci”), xlab =“Diagnostic Odds Ratio”)
 We will examine the results of the DOR_model in Figure 4 onebyone.
 The total effect sizes of all nine studies are shown. The OR of the random effect model is 37.935 (95% CI, 18.186 to 79.132) and pvalue <0.0001. In this diagnosis test, the OR for the positive result among persons with a disease was approximately 38 times higher than the OR for positive results among persons with no disease.
 It appears that the random model has almost no difference according to subgroup (0 vs. 1).
 The Higgins’ I^{2} of all studies is 72.7%, and the pvalue of the Cochrane Q statistics is 0.0003, indicating that there is heterogeneity.
 Bivariate analysis
 The “mada” package for bivariate analysis does not directly present the sensitivity, specificity, and DOR as in MetaDiSc or STATA, which are other DTA applications. Thus, to show the combined overall statistics with “mada” package you should check the source code and calculate it manually.
 Therefore, in this study, the summary statistics were analyzed separately for sensitivity, specificity, and DOR by performing univariate analysis. In the following bivariate analysis, only the SROC curve is estimated using the “mada” package.
 Before loading the “mada” package, the “meta” package that was used before should be unloaded, because “mada” and “meta” both use the “forest” function, which may not be executed if it is called simultaneously by both packages.
 ·detach(package:meta)
Diagnostic test accuracy summary line (summary receiver operating characteristic curve)
 Load the “mada” package for bivariate analysis:
 ·library(mada)
 To see the forest plots of univariate analysis for sensitivity, specificity, and DOR using the “mada” package, enter the following commands:
 · forest(madad(dta_shim), type=“sens”, xlab=“Sensitivity”, snames=dta_shim$id)
 · forest(madad(dta_shim), type=“spec”, xlab=“Specificity”, snames=dta_shim$id)
 ·forest(madauni(dta_shim))
 These plots are the same as those obtained in the univariate analysis, and are not recommended because they do not show the overall effect size of the summary statistics.
 In the “mada” package, use the reitsma function, which is appropriate for a bivariate model.
 ·fit < reitsma(dta_shim, correction.control=“single”)
 ·summary(fit)
 Enter the dta_shim data in the reitsma function. It becomes impossible to calculate if there is ‘0’ in a data cell. To prevent this, you can enter 0.5 in all cells of every study (correction.control= “all”), or correct only the cell of the corresponding study (horizontal) (correction.control=“single”). In the options, you can adjust it to a random value such as ‘correction=0.5,’ where 0.5 is the default value. For models using the reitsma function, ‘fit’ is assigned.
 In addition, you can refer to the area under the curve (AUC), which is 0.906, in the middle of the console window and the values corresponding to the HSROC model.
 Now, we will draw the SROC curve (Figure 5). The graphs will be drawn in the order of commands by overlapping because the first SROC curve remains in the memory.
 · plot(fit, sroclwd=2, xlim=c(0,1), ylim=c(0,1), main=“SROC curve (bivariate model) for Diagnostic Test Accuracy”)
 “plot” is a graph drawing function. Enter the set model fit. “sroclwd=2” indicates the thickness of the SROC curve. Adjust the units of the x and y axes by adjusting xlim and ylim, respectively. The current graph shows the range from a minimum of 0 to a maximum of 1.
 ·points(fpr(dta_shim), sens(dta_shim), pch=2)
 Enter the individual study in points. fpr() and sens() respectively indicate the false positive rate and sensitivity of individual studies in the corresponding data. pch=2 indicates a triangle shape. You can choose from among various shapes: rectangle (0), circle (1), triangle (2), cross (3), scissors (4), rhombus (5), inverted triangle (6), star (8), and black dot (20). The black dot (20) appears to have the best discrimination (Figure 5).
 ·legend(“bottomleft,” c(“SROC,” “95% CI region”), lwd=c(2,1))
 There is an annotation for each curve at the left bottom of the SROC curve.
 Heterogeneity review
 Once the summary statistics and the SROC summary line are presented, we have the major components of the DTA. Then, if there is any significant heterogeneity of study, researchers should verify it and report the heterogeneity factors. The basic assumption of the SROC curve is that the shape of the ROC curve is identical in all studies. However, this basic assumption is not met if there is heterogeneity between studies. There are many causes of this heterogeneity such as chance, difference in cutoff value, difference in study design, prevalence, research environment, and the demographic factors of the sample population [3].
 The DTA presents various methods for diagnosing the heterogeneity [3].
 First, the asymmetry of the SROC curve may be a cause of heterogeneity.
 Second, heterogeneity may be suspected if the degree of scattering or variation of individual studies in the SROC curve is large.
 Third, heterogeneity may be suspected if the betweenstudy variation is greater than the withinstudy variation in the forest plot (sensitivity, specificity, DOR).
 Fourth, heterogeneity may be suspected if the correlation coefficient of sensitivity and specificity is larger than zero.
 The first to third factors only depend on visual distinction, so only the overall outline can be seen.
 The symmetry of the SROC curve indicates the agreement of the models of the divided SROC curves when the SROC curve is divided by a random line from the top of the yaxis to the right bottom of the xaxis. In other words, when the SROC curve is symmetrical and the inflection point is drawn to the top left corner and sharply turned, the area AUC of the SROC curve increases and the Youden’s J index (J=sensitivity+specificity1) becomes high, which indicate a good DTA.
 In visual verification, the SROC curve in this example does not appear to have a high symmetry, and the degree of scattering of individual studies also does not appear to be large.
 According to the withinstudy and betweenstudy variation in the forest plot (Figure 4), the betweenstudy variation does not appear to be large.
■ Sensitivity and specificity correlation coefficient
Finally, to examine the correlation coefficient of sensitivity and specificity, additional variables are created for the current data as follows:
· dta_shim$sn < dta_shim$TP/(dta_shim$TP+dta_shim$FN)
·dta_shim$sp < dta_shim$TN/(dta_shim$FP+dta_shim$TN)
·dta_shim$logitsn < log(dta_shim$sn/(1dta_shim$sn))
·dta_shim$logitsp < log(dta_shim$sp/(1dta_shim$sp))
 First, the sensitivity (dta_shim$sn) and specificity of each study are determined using the equations. Then, the sensitivity and specificity, which are proportion data, are logittransformed to meet the distribution assumption. Then, the variables are checked to determine whether they have been created properly.
 ·View(dta_shim)
 Once the variables are logittransformed, the correlation coefficient of the sensitivity and specificity is obtained as follows:
 ·cor(dta_shim$logitsn, dta_shim$logitsp)
 The correlation coefficient function is “cor”. When the logittransformed sensitivity and specificity are entered in this function, a correlation coefficient of 0.227 is obtained.
 If the sensitivity and specificity are mutually equal and have a normal symmetric distribution, they show a tradeoff relationship. The two are balanced against each other, and when one of them is lowered, the other one is raised. Therefore, the sizes of these two measurements differ in opposite directions depending on the cutoff value in the diagnostic test, and hence, these two values inevitably have a negative correlation.
 The correlation coefficient in this example is a negative value, indicating a low heterogeneity.
■ Meta regression analysis
The “mada” package does not provide functions for the meta regression analysis of the DTA. Therefore, the statistical significance of the moderating variable subgroup (Western European countries vs. other countries) is verified by performing meta regression analysis with the DOR as the effect size.
·library(meta)
·metareg(DOR_model, g, method.tau=“REML,” digits=3)
 Load the “meta” package into the memory again because it was unloaded before the “mada” package was loaded.
 Then, enter the DOR metaanalysis model (DOR_model) and the moderating variable g into the meta regression analysis function metareg. Next, determine the betweenstudy variation of restricted maximumlikelihood estimator, and check the value to only three decimal places.
 The meta regression analysis result confirmed that the pvalue of the moderating variable g was 0.922, indicating statistical insignificance.
CONCLUSION
 This study summarized statistical theory and focused on the actual performance of metaanalysis so that it is easily understandable to general researchers who do not have majors in statistics. In other words, this study aimed to allow general researchers to adequately use already developed statistical methods in their respective fields of study to interpret the results.
 Performing an analysis to determine the DTA in R software can be a complex task because one needs to use various packages. Therefore, we recommend that researchers learn the analysis method using STATA and MetaDiSc applications as well, which can be operated as a single package.
 Researchers who desire to perform an analysis of the DTA should establish the concepts of summary statistics and summary line.
 We hope that this study will help domestic researchers perform metaanalysis more easily, and that it will encourage related research.
SUPPLEMENTARY MATERIALS
NOTES

^{} The authors have no conflicts of interest to declare for this study.
ACKNOWLEDGEMENTS
None.
Figure 1.Summary statistics for diagnostic test accuracy.
Figure 2.Flow chart of diagnostic test accuracy (DTA) using R “mada” & “meta” package. TP, true positive; FP, false positive; FN, false negative; TN, true negative; DOR, diagnostic odds ratio; SROC, summary receiver operating characteristic.
Figure 3.Univariate analysis: sensitivity. CI, confidence interval; g, subgroup.
Figure 4.Univariate analysis: diagnostic odds ratio. OR, odds ratio; CI, confidence interval; g, subgroup.
Figure 5.Summary receiver operating characteristic (SROC) curve (bivariate model) for diagnostic test accuracy. CI, confidence interval; AUC, area under the curve; DOR, diagnostic odds ratio.
Table 1.Diagnostic test accuracy summary statistics [2]
Summary statistics 
Equation 
Definition 
Sn 
TP/(TP+FN) 
Proportion of persons who have positive test results to those with disease 
Sp 
TN/(FP+TN) 
Proportion of persons who have negative test result to those without disease 
PPV 
TP/(TP+FP) 
Proportion of persons with disease to those who have positive test result 
NPV 
TN/(FN+TN) 
Proportion of persons without disease to those who have negative test result 
LR+ 
Sn/(1Sp) 
Ratio of the probability of a positive test result among those with disease to that of a positive test result among those without disease 
LR 
(1Sn)/Sp 
Ratio of the probability of a negative test result among those with disease to that of a negative test result among those without disease 
Accuracy of index test 
(TP+TN)/(TP+FP+FN+TN) 
The proportion of persons who are true positive and persons who are true negative among all subjects 
DOR 
(TP*TN)/(FP*FN) 
The ratio of the OR for a positive test result among persons with disease to that among persons without disease 
Table 2.Sample data for diagnostic test accuracy [2]
Id 
TP 
FP 
FN 
TN 
g 
Wiegmann 
21 
1 
9 
104 
1 
Bouhanick 
49 
21 
7 
110 
1 
Schwab 
24 
5 
3 
31 
1 
Zelmanovitz 
39 
6 
5 
48 
0 
Ahn 
23 
9 
7 
41 
0 
Ng 
12 
7 
2 
44 
0 
Gansevoort 
10 
13 
3 
40 
1 
Incerti 
82 
12 
7 
177 
0 
Sampaio 
99 
45 
21 
128 
0 
REFERENCES
 1. Hwang SD, Shim SR. Metaanalysis: from forest plot to network metaanalysis. Seoul: Hannarae; 2018. p 224246 (Korean).
 2. Shim SR. Diagnostic test accuracy using R & MetaDiSc software. Gwacheon: SDB Lab; 2019. (Korean).
 3. Shim SR, Shin IS, Bae JM. Metaanalysis of diagnostic tests accuracy using STATA software. J Health Info Stat 2015;40:190199 (Korean).
 4. Moses LE, Shapiro D, Littenberg B. Combining independent studies of a diagnostic test into a summary ROC curve: data analytic approaches and some additional considerations. Stat Med 1993;12:12931316.ArticlePubMed
 5. Littenberg B, Moses LE. Estimating diagnostic accuracy from multiple conflicting reports: a new meta analytic method. Med Decis Making 1993;13:313321.ArticlePubMed
 6. Reitsma JB, Glas AS, Rutjes AW, Scholten RJ, Bossuyt PM, Zwinderman AH. Bivariate analysis of sensitivity and specificity produces informative summary measures in diagnostic reviews. J Clin Epidemiol 2005;58:982990.ArticlePubMed
 7. Rutter CM, Gatsonis CA. A hierarchical regression approach to metaanalysis of diagnostic test accuracy evaluations. Stat Med 2001;20:28652884.ArticlePubMed
 8. Arends LR, Hamza TH, van Houwelingen JC, HeijenbrokKal MH, Hunink MG, Stijnen T. Bivariate random effects metaanalysis of ROC curves. Med Decis Making 2008;28:621638.ArticlePubMed
 9. Harbord RM, Deeks JJ, Egger M, Whiting P, Sterne JA. A unification of models for metaanalysis of diagnostic accuracy studies. Biostatistics 2007;8:239251.ArticlePubMedPDF
 10. Comprehensive R Archive Network. mada: metaanalysis of diagnostic accuracy. [cited 2019 May 8]. Available from: https://cran.rproject.org/web/packages/mada/index.html.
 11. Wu HY, Peng YS, Chiang CK, Huang JW, Hung KY, Wu KD, et al. Diagnostic performance of random urine samples using albumin concentration vs ratio of albumin to creatinine for microalbuminuria screening in patients with diabetes mellitus: a systematic review and metaanalysis. JAMA Intern Med 2014;174:11081115.ArticlePubMed
Citations
Citations to this article as recorded by

Identifying the culprit artery via 12‐lead electrocardiogram in inferior wall
ST
‐segment elevation myocardial infarction: A meta‐analysis
Peng Zhou, Yingying Wu, Meng Wang, Yikai Zhao, Yangjie Yu, Maieryemu Waresi, Huiyang Li, Bo Jin, Xinping Luo, Jian Li
Annals of Noninvasive Electrocardiology.2023;[Epub] CrossRef  Clinical value of alarm features for colorectal cancer: a metaanalysis
Leonardo Frazzoni, Liboria Laterza, Marina La Marca, Rocco Maurizio Zagari, Franco Radaelli, Cesare Hassan, Alessandro Repici, Antonio Facciorusso, Paraskevas Gkolfakis, Cristiano Spada, Konstantinos Triantafyllou, Franco Bazzoli, Mario DinisRibeiro, Lor
Endoscopy.2023; 55(05): 458. CrossRef  CA125 for the Diagnosis of Advanced Urothelial Carcinoma of the Bladder: A Systematic Review and MetaAnalysis
HsuanJen Lin, RouhMei Hu, HungChih Chen, ChungChih Lin, ChiYu Lee, CheYi Chou
Cancers.2023; 15(3): 813. CrossRef  Performance of screening tools for cervical neoplasia among women in low and middleincome countries: A systematic review and metaanalysis
Sabrina K. Smith, Oguchi Nwosu, Alex Edwards, Meseret Zerihun, Michael H. Chung, Kara Suvada, Mohammed K. Ali, Nebiyu Dereje
PLOS Global Public Health.2023; 3(2): e0001598. CrossRef  Diagnostic accuracy of pointofcare lung ultrasound for COVID19: a systematic review and metaanalysis
Ashley Matthies, Michael Trauer, Karl Chopra, Robert David Jarman
Emergency Medicine Journal.2023; 40(6): 407. CrossRef  Influence of seasonal and operator variations on diagnostic accuracy of lateral flow devices during the COVID19 pandemic: a systematic review and metaanalysis
Ashwin Krishnamoorthy, Subashini Chandrapalan, Gohar JalayeriNia, Yaqza Hussain, Ayman Bannaga, Ian Io Lei, Ramesh Arasaradnam
Clinical Medicine.2023; 23(2): 144. CrossRef  Prognostic role of different findings at echocardiography in acute pulmonary embolism: a critical review and metaanalysis
Ludovica Anna Cimini, Matteo Candeloro, Magdalena Pływaczewska, Giorgio Maraziti, Marcello Di Nisio, Piotr Pruszczyk, Giancarlo Agnelli, Cecilia Becattini
ERJ Open Research.2023; 9(2): 006412022. CrossRef  Virotyping and genetic antimicrobial susceptibility testing of porcine ETEC/STEC strains and associated plasmid types
Nick Vereecke, Sander Van Hoorde, Daniel Sperling, Sebastiaan Theuns, Bert Devriendt, Eric Cox
Frontiers in Microbiology.2023;[Epub] CrossRef  The Diagnostic Power of Circulating miR1246 in Screening Cancer: An Updated Metaanalysis
Khanh Quang Huynh, Anh Tuan Le, Thang Thanh Phan, Toan Trong Ho, Suong Phuoc Pho, Hang Thuy Nguyen, Binh Thanh Le, Thuc Tri Nguyen, Son Truong Nguyen, Ihtisham Bukhari
Oxidative Medicine and Cellular Longevity.2023; 2023: 1. CrossRef  A Systematic Review and MetaAnalysis Comparing the Diagnostic Accuracy Tests of COVID19
Juan Jeferson VilcaAlosilla, Mayron Antonio CandiaPuma, Katiusca CoronelMonje, Luis Daniel GoyzuetaMamani, Alexsandro Sobreira Galdino, Ricardo Andrez MachadodeÁvila, Rodolfo Cordeiro Giunchetti, Eduardo Antonio Ferraz Coelho, Miguel Angel ChávezFu
Diagnostics.2023; 13(9): 1549. CrossRef  Prediction Models for Intrauterine Growth Restriction Using Artificial Intelligence and Machine Learning: A Systematic Review and MetaAnalysis
Riccardo Rescinito, Matteo Ratti, Anil Babu Payedimarri, Massimiliano Panella
Healthcare.2023; 11(11): 1617. CrossRef  Reliability of machine learning to diagnose pediatric obstructive sleep apnea: Systematic review and meta‐analysis
Gonzalo C. Gutiérrez‐Tobal, Daniel Álvarez, Leila Kheirandish‐Gozal, Félix del Campo, David Gozal, Roberto Hornero
Pediatric Pulmonology.2022; 57(8): 1931. CrossRef  How to Analyze the Diagnostic Performance of a New Test? Explained with Illustrations
Deepak Dhamnetiya, Ravi Prakash Jha, Shalini Shalini, Krittika Bhattacharyya
Journal of Laboratory Physicians.2022; 14(01): 090. CrossRef  Diagnostic accuracy of magnetic resonance imaging targeted biopsy techniques compared to transrectal ultrasound guided biopsy of the prostate: a systematic review and metaanalysis
E. J. Bass, A. Pantovic, M. J. Connor, S. Loeb, A. R. Rastinehad, M. Winkler, Rhian Gabe, H. U. Ahmed
Prostate Cancer and Prostatic Diseases.2022; 25(2): 174. CrossRef  Diagnostic performance and prognostic impact of coronary angiography‐based Index of Microcirculatory Resistance assessment: A systematic review and meta‐analysis
Weijia Li, Tatsunori Takahashi, Saul A. Rios, Azeem Latib, Joo Myung Lee, William F. Fearon, Yuhei Kobayashi
Catheterization and Cardiovascular Interventions.2022; 99(2): 286. CrossRef  Assessing the Knowledge of the Osteopathic Profession in New York City’s Eastern European Communities
Justin Chin, Lina Kleyn, Emily Dube, Mark Terrell, Christine M Lomiguen, Mikhail Volokitin
Cureus.2022;[Epub] CrossRef  The diagnostic accuracy of clinical tests for anterior cruciate ligament tears are comparable but the Lachman test has been previously overestimated: a systematic review and metaanalysis
Pawel A. Sokal, Richard Norris, Thomas W. Maddox, Rachel A. Oldershaw
Knee Surgery, Sports Traumatology, Arthroscopy.2022; 30(10): 3287. CrossRef  Comparison of Diagnostic Test Accuracy of ConeBeam Breast Computed Tomography and Digital Breast Tomosynthesis for Breast Cancer: A Systematic Review and MetaAnalysis Approach
Temitope Emmanuel Komolafe, Cheng Zhang, Oluwatosin Atinuke Olagbaju, Gang Yuan, Qiang Du, Ming Li, Jian Zheng, Xiaodong Yang
Sensors.2022; 22(9): 3594. CrossRef  Clinical Validity of 16α[18F]Fluoro17βEstradiol Positron Emission Tomography/Computed Tomography to Assess Estrogen Receptor Status in Newly Diagnosed Metastatic Breast Cancer
Jasper J.L. van Geel, Jorianne Boers, Sjoerd G. Elias, Andor W.J.M. Glaudemans, Erik F.J. de Vries, Geke A.P. Hospers, Michel van Kruchten, Evelien J.M. Kuip, Agnes Jager, Willemien C. Menkevan der Houven van Oordt, Bert van der Vegt, Elisabeth G.E. de V
Journal of Clinical Oncology.2022; 40(31): 3642. CrossRef  Diagnostic Performance of Antigen Rapid Diagnostic Tests, Chest Computed Tomography, and Lung PointofCareUltrasonography for SARSCoV2 Compared with RTPCR Testing: A Systematic Review and Network MetaAnalysis
Sung Ryul Shim, SeongJang Kim, Myunghee Hong, Jonghoo Lee, MinGyu Kang, Hyun Wook Han
Diagnostics.2022; 12(6): 1302. CrossRef  Diagnostic Accuracy of Machine Learning Models on Mammography in Breast Cancer Classification: A MetaAnalysis
Tengku Muhammad Hanis, Md Asiful Islam, Kamarul Imran Musa
Diagnostics.2022; 12(7): 1643. CrossRef  Association of No Evidence of Disease Activity With No Longterm Disability Progression in Multiple Sclerosis
Dalia Rotstein, Jacqueline Madeleine Solomon, Maria Pia Sormani, Xavier Montalban, Xiang Y. Ye, Dina Dababneh, Alexandra Muccilli, Prakesh Shah
Neurology.2022; 99(2): e209. CrossRef  Prospective study of Na[18F]F PET/CT for cancer staging in morbidly obese patients compared with [99mTc]TcMDP wholebody planar, SPECT and SPECT/CT
Sharjeel Usmani, Najeeb Ahmed, Gopinath Gnanasegaran, Fareeda Al kandari, Fahad Marafi, Ahmed BaniMustafa, Ahmed Musbah, Maryam Jassem Almashmoum, Tim Van den Wyngaert
Acta Oncologica.2022; 61(10): 1230. CrossRef  Structured data vs. unstructured data in machine learning prediction models for suicidal behaviors: A systematic review and metaanalysis
Danielle Hopkins, Debra J. Rickwood, David J. Hallford, Clare Watsford
Frontiers in Digital Health.2022;[Epub] CrossRef  Growth Differentiation Factor15 as a Candidate Biomarker in Gynecologic Malignancies: A Metaanalysis
Dipayan Roy, Anupama Modi, Purvi Purohit, Manoj Khokhar, Manu Goyal, Shailja Sharma, Puneet Setia, Antonio Facciorusso, Praveen Sharma
Cancer Investigation.2022; 40(10): 901. CrossRef  Diagnostic accuracy of MRI techniques for treatment response evaluation in patients with brain metastasis: A systematic review and metaanalysis
Wouter H.T. Teunissen, Chris W. Govaerts, Miranda C.A. Kramer, Jeremy A. Labrecque, Marion Smits, Linda Dirven, Anouk van der Hoorn
Radiotherapy and Oncology.2022; 177: 121. CrossRef  Accuracy of Diagnostic Tests for the Detection of Chagas Disease: A Systematic Review and MetaAnalysis
Mayron Antonio CandiaPuma, Laura Yesenia MachacaLuque, Brychs Milagros RoquePumahuanca, Alexsandro Sobreira Galdino, Rodolfo Cordeiro Giunchetti, Eduardo Antonio Ferraz Coelho, Miguel Angel ChávezFumagalli
Diagnostics.2022; 12(11): 2752. CrossRef  Assessment of the accuracy of 11 different diagnostic tests for the detection of Schistosomiasis mansoni in individuals from a Brazilian area of low endemicity using latent class analysis
Silvia Gonçalves Mesquita, Roberta Lima Caldeira, Tereza Cristina Favre, Cristiano Lara Massara, Lílian Christina Nóbrega Holsbach Beck, Taynãna César Simões, Gardênia Braz Figueiredo de Carvalho, Flória Gabriela dos Santos Neves, Gabriela de Oliveira, La
Frontiers in Microbiology.2022;[Epub] CrossRef  Diagnosis of pathological conditions through electronic nose analysis of urine samples: a systematic review and metaanalysis
Helga A.S. Afonso, Mariana V. Farraia, Mónica A. Vieira, João Cavaleiro Rufo
Porto Biomedical Journal.2022; 7(6): e188. CrossRef  Comparison of diagnostic test accuracy modeling of metaanalysis for data integration
Sung Ryul Shim
Epidemiology and Health.2022; : e2022083. CrossRef  A systematic review and metaanalysis of the diagnostic accuracy of biparametric prostate MRI for prostate cancer in men at risk
E. J. Bass, A. Pantovic, M. Connor, R. Gabe, A. R. Padhani, A. Rockall, H. Sokhi, H. Tam, M. Winkler, H. U. Ahmed
Prostate Cancer and Prostatic Diseases.2021; 24(3): 596. CrossRef  Diagnosis of Alzheimer’s Disease in Developed and Developing Countries: Systematic Review and MetaAnalysis of Diagnostic Test Accuracy
Miguel A. ChávezFumagalli, Pallavi Shrivastava, Jorge A. AguilarPineda, Rita NietoMontesinos, Gonzalo Davila DelCarpio, Antero PeraltaMestas, Claudia CaracelaZeballos, Guillermo ValdezLazo, Victor FernandezMacedo, Alejandro PinoFigueroa, Karin J.
Journal of Alzheimer's Disease Reports.2021; 5(1): 15. CrossRef  Perceptions of the Osteopathic Profession in New York City's Korean Communities
Justin Chin, Haeinn Woo, Diane Choi, Emily Dube, Mikhail Volokitin,, Christine Lomiguen
Osteopathic Family Physician.2021; : 12. CrossRef  Breath biopsy of breast cancer using sensor array signals and machine learning analysis
HsiaoYu Yang, YiChia Wang, HsinYi Peng, ChiHsiang Huang
Scientific Reports.2021;[Epub] CrossRef  MetaAnalysis and Systematic Review of the Application of Machine Learning Classifiers in Biomedical Applications of Infrared Thermography
Carolina Magalhaes, Joaquim Mendes, Ricardo Vardasca
Applied Sciences.2021; 11(2): 842. CrossRef  Predicting Clinical Outcomes in Acute Ischemic Stroke Patients Undergoing Endovascular Thrombectomy with Machine Learning
Yao Hao Teo, Isis Claire Z. Y. Lim, Fan Shuen Tseng, Yao Neng Teo, Cheryl Shumin Kow, Zi Hui Celeste Ng, Nyein Chan Ko Ko, ChingHui Sia, Aloysius S. T. Leow, Wesley Yeung, Wan Yee Kong, Bernard P. L. Chan, Vijay K. Sharma, Leonard L. L. Yeo, Benjamin Y.
Clinical Neuroradiology.2021; 31(4): 1121. CrossRef  Application of artificial intelligence in diagnosis of osteoporosis using medical images: a systematic review and metaanalysis
L. Gao, T. Jiao, Q. Feng, W. Wang
Osteoporosis International.2021; 32(7): 1279. CrossRef  The Accuracy of Visceral Adiposity Index for the Screening of Metabolic Syndrome: A Systematic Review and MetaAnalysis
Moniba Bijari, Sara Jangjoo, Nima Emami, Sara Raji, Mahdi Mottaghi, Roya Moallem, Ali Jangjoo, Amin Saberi, Pawel Grzmil
International Journal of Endocrinology.2021; 2021: 1. CrossRef  Immunofluorescence Targeting PBP2a Protein: A New Potential Methicillin Resistance Screening Test
Serenella Silvestri, Elisa Rampacci, Valentina Stefanetti, Michele Trotta, Caterina Fani, Lucia Levorato, Chiara Brachelente, Fabrizio Passamonti
Frontiers in Veterinary Science.2021;[Epub] CrossRef  Perceptions of the osteopathic profession in New York City’s Chinese Communities
Justin Chin, Sarah Li, Gregory Yim, YaQun Arlene Zhou, Peter Justin Wan, Emily R Dube, Mikhail Volokitin, Sonu Sahni, Mark A Terrell, Christine M Lomiguen
Family Medicine and Community Health.2020; 8(1): e000248. CrossRef  Development of a voiding diary using urination recognition technology in mobile environment
Gun Hyun Park, Su Jin Kim, Young Sam Cho
Journal of Exercise Rehabilitation.2020; 16(6): 529. CrossRef