# Craddock test for testing the homogeneity of yearly average temperatures of Milan vs Turin
# example based on:
# Training session on homogenisation methods (Bologna 17-18 May, 2005)
# Historical Climatology Group
# Istituto di Scienze dell'Atmosfera e del Clima (Institute of Atmospheric Sciences and Climate)

craddock.result <- CraddockTest(yearly.average.temperature.Turin.Milan,2,3)
# Example of inhomogeneity
tmp <- yearly.average.temperature.Turin.Milan
tmp[20:43,3] <- tmp[20:43,3]+1 # adding an error of +1 degree to Milan's yearly average temperatures from 1980
craddock.result <- CraddockTest(tmp,2,3)

Try the climtrends package in your browser

Any scripts or data that you put into this service are public.

climtrends documentation built on May 29, 2017, 11:58 p.m.