# call data
data(laparotomy)
# define the function parameters
markers <- cbind(laparotomy$ddimer, laparotomy$log_leukocyte)
status <- factor(laparotomy$group, levels = c("not_needed", "needed"))
neg.set <- markers[status == levels(status)[1], ]
pos.set <- markers[status == levels(status)[2], ]
lambda <- 0.5
stat <- helper_PT(lambda, neg.set = neg.set, pos.set = pos.set)
Run the code above in your browser using DataLab