# call data
data(laparotomy)
# define the function parameters
markers <- laparotomy[, -1]
status <- factor(laparotomy$group, levels = c("not_needed", "needed"))
event <- "needed"
score1 <- linComb(
markers = markers, status = status, event = event,
method = "logistic", resample = "none",
standardize = "none", direction = "<", cutoff.method = "Youden"
)
comb.score1 <- predict(score1, markers)
score2 <- nonlinComb(
markers = markers, status = status, event = "needed", include.interact = TRUE,
method = "polyreg", resample = "repeatedcv", nfolds = 5,
nrepeats = 10, cutoff.method = "Youden", direction = "auto"
)
comb.score2 <- predict(score2, markers)
score3 <- mathComb(
markers = markers, status = status, event = event,
method = "distance", distance = "euclidean", direction = "auto",
standardize = "tScore", cutoff.method = "Youden"
)
comb.score3 <- predict(score3, markers)
Run the code above in your browser using DataLab