This function computes the PIPs of all potential predictors
An object of class
an named vector with all PIPs
1 2 3 4 5 6 7 8 9 10 11 12 13 14
# extract the data: data("VAP_data") # the definition of the full model with three potential predictors: FULL <- outcome ~ ns(day, df = 4) + gender + type + SOFA # here we define time as a spline with 3 knots # computation of the posterior model probabilities: test <- PMP(fullModel = FULL, data = VAP_data, discreteSurv = TRUE, maxit = 150) class(test) #computation of the posterior inclusion probabilities: postInclusionProb(test)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.