********************
> B <- c(60,45,42,47,55) #Jan, Sept-Dec 1998 GISS land/sea
> mean(B)
[1] 49.8
> C <- c(50,40,40,52,61,60,51,66,59,62,49,59,66,54,57,59) #1998-2013 annual GISS land/sea values with 1998 editing out the 7 high months
> fit <- lm(C~Years)
> summary(fit)

Call:
lm(formula = C ~ Years)

Residuals:
Min 1Q Median 3Q Max
-10.5162 -4.5129 0.6235 5.7254 11.1235

Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1693.6015 758.4956 -2.233 0.0424 *
Years 0.8721 0.3782 2.306 0.0369 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 6.974 on 14 degrees of freedom
Multiple R-squared: 0.2752, Adjusted R-squared: 0.2235
F-statistic: 5.317 on 1 and 14 DF, p-value: 0.03694

************************