thyroids {clust.bin.pair} | R Documentation |
Following surgery which confirmed the absence of hyperparathyroidism two diagnostic tests, PET and SPECT, were performed. Their measures of true negatives and false positives are reported. Data reported in Obuchowki 1998.
data(thyroids)
A data frame with 21 rows and 6 variables:
ID of the patient
number of glands tested from the patient
number of true negatives from the PET test
individual results per gland from the PET test
number of true negatives from the SPECT test
individual results per gland from the SPECT test
Obuchowski, N. A. (1998). On the comparison of correlated proportions for clustered data. Statistics in medicine, 17(13), 1495-1507.
data(thyroids)
thyroids$n.glands == sapply(thyroids$x.pet, length)
thyroids$n.glands == sapply(thyroids$x.spect, length)
thyroids$n.pet == sapply(thyroids$x.pet, function(x) length(which(x == 1)))
thyroids$n.spect == sapply(thyroids$x.spect, function(x) length(which(x == 1)))
tc <- nested.to.contingency(thyroids$x.pet, thyroids$x.spect)
clust.bin.pair(tc[,'ak'], tc[,'bk'], tc[,'ck'], tc[,'dk'], method="obuchowski")
do.call(clust.bin.pair, data.frame(tc))