How to compute a p value and extract a critical value in R