# area below curve from -infinity to -2 lower <- pt(-2, df = 1) # prop of area below curve within t = +/-2 1 - (2 * lower) [1] 0.7048 # critical values for inner 95% qt(c(.025, .975), df = 1) [1] -12.7062 12.7062