# standard normal homogeneity test, SD different from 1(SNHT) for Turin 1961-2003
testTurin<- SNHTabsoluteII(yearly.average.temperature.Turin.Milan[,2])
# standard normal homogeneity test, SD different from 1(SNHT) for Milan 1961-2003
testMilan<- SNHTabsoluteII(yearly.average.temperature.Turin.Milan[,3])
plot(1961:2002,testTurin,type='l',col='cyan',ylim=c(0,30),
main='standard normal homogeneity test (SNHT), SD different from 1 for 1961-2003',
xlab='Years',ylab='test value')
par(new=TRUE)
plot(1961:2002,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