# standard normal homogeneity test for trends(SNHT) for Turin 1961-2003
testTurin<- SNHTabsoluteTrend(yearly.average.temperature.Turin.Milan[,2])$SNHTdata
# standard normal homogeneity test for trends(SNHT) for Milan 1961-2003
testMilan<- SNHTabsoluteTrend(yearly.average.temperature.Turin.Milan[,3])$SNHTdata
plot(1961:2002,testTurin,type='l',col='cyan',ylim=c(0,30),
main='standard normal homogeneity test (SNHT) for trends, 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