# Loading required libraries
library(survival)
library(sandwich)
# Reading in a dataset
data(mri)
# Linear regression of LDL on age (with robust SE by default)
testReg <- regress ("mean", ldl~age, data = mri)
# 95% Prediction Interval for age 50
predict(testReg)
Run the code above in your browser using DataLab