# Chapter 6 section 3.1 details the creation of compiled series. This appendix includes the code for the creation of such series.


## A single daily time series for the region


michdday=as.data.frame(matrix(nrow=15340,ncol=22))

names(michdday)=c('Time','Nstats','Meancl','Sdcl','Meanrl','Sdrl','Meanrt','Sdrt','RMSErl','RMSErt')

michdday$Time=1:15340

for (i in 1:15340){

	time=dataran1[dataran1$Time==i,] # Get all values from all stations where Time = i
	
	michdday$Nstats[i]=nrow(time)	
	
	michdday$Meancl[i]=mean(time$Cleanm)
	
	michdday$Sdcl[i]=sd(time$Cleanm)
	
	michdday$Meanrl[i]=mean(time$Released)

	michdday$Sdrl[i]=sd(time$Released)

	michdday$Meanrt[i]=mean(time$Returned)

	michdday$Sdrt[i]=sd(time$Returned)
	
        michdday$RMSErl[i]=rmse(time$Cleanm,time$Released)
	
	michdday$RMSErt[i]=rmse(time$Cleanm,time$Returned)}


## A single monthly time series for the region


michdmonth=as.data.frame(matrix(nrow=504,ncol=24))

names(michdmonth)=c('Time','Year','Month','Nstats','Meancl','Meancld','Sdcl','Sdcld','Meanrl','Meanrld','Sdrl','Sdrld','Meanrt','Meanrtd','Sdrt','Sdrtd','MAErl','MAEdrl','MAErt','MAEdrt','RMSErl','RMSEdrl','RMSErt','RMSEdrt')

michdmonth$Year=sort(rep(1970:2011,12))

michdmonth$Month=rep(1:12,42)

michdmonth$Time=1:504

for (i in 1970:2011){

	Year=dataran1[dataran1\$Year==i,]

	for (j in 1:12){

		Month=Year[Year$Month==j,]

		michdmonth$Nstats[(i-1970)*12+j]=nrow(Month)
		
		michdmonth$Meancl[(i-1970)*12+j]=mean(Month$Cleanm)	
		
		michdmonth$Sdcl[(i-1970)*12+j]=sd(Month$Cleanm)

		michdmonth$Meanrl[(i-1970)*12+j]=mean(Month$Released)

		michdmonth$Sdrl[(i-1970)*12+j]=sd(Month$Released)
		
		michdmonth$Meanrt[(i-1970)*12+j]=mean(Month$Returned)
		
		michdmonth$Sdrt[(i-1970)*12+j]=sd(Month$Returned)
				
		michdmonth$RMSErl[(i-1970)*12+j]=rmse(Month$Cleanm,Month$Released)
		
		michdmonth$RMSErt[(i-1970)*12+j]=rmse(Month$Cleanm,Month$Returned)

	}

}


## A time series of a single year for the region


michdday1=as.data.frame(matrix(nrow=366,ncol=22))

names(michdday1)=c('Dyear','Nstats','Mencl','Meancl','Sdcl','Sdcld','Menrl','Meanrl'
,'Sdrl','Sdrld', 'Menrt','Meanrt','Sdrt','Sdrtd','MAErl','MAEdrl','MAErt','MAEdrt','RMSErl','RMSEdrl','RMSErt','RMSEdrt')

michdday1$Dyear=1:366 

for (i in 1:366){

	time=dataran1[dataran1$Dyear==i,]
	
	michdday1$Nstats[i]=nrow(time)
	
	michdday1$Mencl[i]=mean(time$Cleanm)
	
	michdday1$Sdcl[i]=sd(time$Cleanm)
	
	michdday1$Menrl[i]=mean(time$Released)
	
	michdday1$Sdrl[i]=sd(time$Released)
	
	michdday1$Menrt[i]=mean(time$Returned)
	
	michdday1$Sdrt[i]=sd(time$Returned)
	
	michdday1$RMSErl[i]=rmse(time$Cleanm,time$Released)
	
	michdday1$RMSErt[i]=rmse(time$Cleanm,time$Returned)

}