I have some clinical diagnostic data from three sites that I have been looking at performance against a reference standard using the McNemar test.
How would I set up the analysis to investigate if there are any differences between the three sites?
My data (in R) look like:
myData <- read.table(header=TRUE, text=' Site ReferenceStandard NewDiagnostic Count A R R 47 A S R 2 A R S 0 A S S 387 B R R 3 B S R 1 B R S 7 B S S 161 C R R 13 C S R 0 C R S 0 C S S 108') ct <- xtabs(Count ~ ReferenceStandard + NewDiagnostic + Site, data=myData) ftable(ct) Site A B C ReferenceStandard NewDiagnostic R R 47 3 13 S 0 7 0 S R 2 1 0 S 387 161 108
Do I use the Cochran–Mantel–Haenszel?