1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
function.MaxSpSe <-
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)
m <- vector()
for(i in 1:length(measures.acc$cutoffs)) {
if (measures.acc$Se[i,1] <= measures.acc$Sp[i,1]) {
m[i] <- measures.acc$Se[i,1]
} else {
m[i] <- measures.acc$Sp[i,1]
}
}
M <- max(m,na.rm=TRUE)
optimal.index <- which(round(m,10) == round(M,10))
cMaxSpSe <- measures.acc$cutoffs[optimal.index]
optimal.cutoff <- obtain.optimal.measures(cMaxSpSe, measures.acc)
res <- list(measures.acc = measures.acc, optimal.cutoff = optimal.cutoff, criterion = m, optimal.criterion = M)
res
}
|
