This thesis provides a nonstationary statistical analysis of annual maximum temperature records to evaluate global warming. It demonstrates that a nonstationary extreme value Weibull model with a linear trend in the location parameter best explains the data among various parametric and nonparametric models. However, other modeling techniques using splines in a generalized additive model previously showed that the trend in annual maxima is not simultaneously significant over time. The thesis develops theoretical backgrounds on state-of-the-art extreme value analysis methods and presents their careful application in a reusable R package and Shiny application.