dat <- data.frame(AGE = c(20,40,65),
TBM = c(50,70,90),
HGT = c(150,170,200),
MALE = c(TRUE,FALSE,TRUE))
schnider_poppk(dat, rand = FALSE, rate = FALSE)
schnider_poppk(dat, rand = TRUE, rate = TRUE)
Run the code above in your browser using DataLab