1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
function.MeanPrev <-
function(data, marker, status, tag.healthy = 0, direction = c("<", ">"), control = control.cutpoints(), pop.prev, ci.fit = FALSE, conf.level = 0.95, measures.acc){
direction <- match.arg(direction)
if (measures.acc$cutoffs < 0 || measures.acc$cutoffs > 1) {
warning("Diagnostic marker values are not between 0 and 1 for this \n criterion. A data transformation has been performed.", call. = FALSE, immediate. = TRUE)
tmarker <- (measures.acc$cutoffs - min(measures.acc$cutoffs))/(max(measures.acc$cutoffs)-min(measures.acc$cutoffs))
difference <- abs(tmarker-mean(tmarker))
} else {
difference <- abs(measures.acc$cutoffs-mean(measures.acc$cutoffs))
}
cMeanPrev <- measures.acc$cutoffs[which(round(difference,10) == round(min(difference,na.rm=TRUE),10))]
optimal.cutoff <- obtain.optimal.measures(cMeanPrev, measures.acc)
res <- list(measures.acc = measures.acc, optimal.cutoff = optimal.cutoff)
res
}
|
