These data are from: http://www.ncdc.noaa.gov/cag/time-series/global
The global component of Climate at a Glance (GCAG) provides global-scale temperature information using data from the Global Historical Climatology Network-Monthly (GHCN-M) data set and International Comprehensive Ocean-Atmosphere Data Set (ICOADS). The tool provides near real-time analysis of monthly and annual temperatures for the globe and is intended for the study of climate variability and change. Data is provided globally, by hemisphere, and by land and ocean surface components. The interactive mapping tool allows analysis of the spatial patterns of global temperature anomalies and regional temperature analysis. Recent data are preliminary and may be modified after appropriate quality control has been performed.
gt <- matrix(c(1990, 0.28, 1991, 0.41, 1992, 0.38, 1993, 0.23, 1994, 0.23, 1995,
0.4, 1996, 0.41, 1997, 0.33, 1998, 0.57, 1999, 0.61, 2000, 0.43, 2001, 0.43,
2002, 0.6, 2003, 0.59, 2004, 0.63, 2005, 0.56, 2006, 0.65, 2007, 0.63, 2008,
0.52, 2009, 0.55, 2010, 0.62, 2011, 0.62, 2012, 0.53, 2013, 0.61, 2014,
0.62), byrow = TRUE, ncol = 2, dimnames = list(c(), c("year", "global.temp")))
The regression line below only uses data from 1998–2013, inclusive.
(gt.lm <- lm(global.temp ~ year, data.frame(gt), subset = (year >= 1998 & year <=
2013))) # 2014 data must be incomplete
##
## Call:
## lm(formula = global.temp ~ year, data = data.frame(gt), subset = (year >=
## 1998 & year <= 2013))
##
## Coefficients:
## (Intercept) year
## -7.53860 0.00404
plot(global.temp ~ year, gt, type = "b", xlab = "Year", ylab = "Global Temperature Anomolay (Degrees Celsius)",
ylim = c(0, 1))
abline(gt.lm, col = "red")
