Calculate intercepts and associated p-value
calc_p_apply(dat, t2, t1)
dat
data.frame
t2
upper intercept age
t1
lower intercept age