library(smooth)

## Warning: package 'smooth' was built under R version 4.0.2

## Loading required package: greybox

## Warning: package 'greybox' was built under R version 4.0.2

## Registered S3 method overwritten by 'quantmod':
##   method            from
##   as.zoo.data.frame zoo

## Package "greybox", v0.6.0 loaded.

## This is package "smooth", v2.6.0

globaltemp=read.csv(file = "C://temp//GlobalTemperatures.csv", header = TRUE)
attach(globaltemp)
names(globaltemp)

## [1] "dt"                                       
## [2] "LandAverageTemperature"                   
## [3] "LandAverageTemperatureUncertainty"        
## [4] "LandMaxTemperature"                       
## [5] "LandMaxTemperatureUncertainty"            
## [6] "LandMinTemperature"                       
## [7] "LandMinTemperatureUncertainty"            
## [8] "LandAndOceanAverageTemperature"           
## [9] "LandAndOceanAverageTemperatureUncertainty"

Creating time-series objects:

landavg1=ts(data = LandAverageTemperature, start = c(1850, 1), end = c(2015,12), frequency = 12)
landmax1=ts(data = LandMaxTemperature, start = c(1850, 1), end = c(2015,12), frequency = 12)
landoceanavg1=ts(data = LandAndOceanAverageTemperature, start = c(1850, 1), end = c(2015,12), frequency = 12)

Plotting timeseries objects:

plot(decompose(landavg1))

plot(decompose(landmax1))

plot(decompose(landoceanavg1))

Positive trend as well as seasonality present. Some random variation is also present.

10-year moving average on Land and Ocean Average Temperature:

sma(landoceanavg1, order = 10)

## Time elapsed: 0.02 seconds
## Model estimated: SMA(10)
## Initial values were produced using backcasting.
## 
## Loss function type: MSE; Loss function value: 2.1596
## Error standard deviation: 1.4703
## Sample size: 1992
## Number of estimated parameters: 2
## Number of degrees of freedom: 1990
## Information criteria:
##      AIC     AICc      BIC     BICc 
## 7190.693 7190.699 7201.887 7201.910