# Buishand range test for Turin 1961-2003
BuishandTurin<- BuishandRangeTest(yearly.average.temperature.Turin.Milan[,2])
# Buishand range test for Milan 1961-2003
BuishandMilan<- BuishandRangeTest(yearly.average.temperature.Turin.Milan[,3])
plot(1961:2003,BuishandTurin,type='l',col='cyan',ylim=c(-2.5,0.5),
main='Buishand range test 1961-2003',xlab='Years',ylab='Buishand value')
par(new=TRUE)
plot(1961:2003,BuishandMilan,type='l',col='magenta',ylim=c(-2.5,0.5),xlab='',ylab='')
legend("top",c('Turin','Milan'),text.col=c('cyan','magenta'))
Run the code above in your browser using DataLab