# standard normal homogeneity test (SNHT) for Turin 1961-2003
testTurin<- SNHTabsolute2SDs(yearly.average.temperature.Turin.Milan[,2])
# standard normal homogeneity test (SNHT) for Milan 1961-2003
testMilan<- SNHTabsolute2SDs(yearly.average.temperature.Turin.Milan[,3])
plot(1961:2000,testTurin,type='l',col='cyan',ylim=c(0,30),
main='standard normal homogeneity test (SNHT) 2 SDs 1961-2003',xlab='Years',ylab='test value')
par(new=TRUE)
plot(1961:2000,testMilan,type='l',col='magenta',ylim=c(0,30),xlab='',ylab='')
legend("top",c('Turin','Milan'),text.col=c('cyan','magenta'))
Run the code above in your browser using DataLab