function(regressionmatrix)
{
	if(ncol(regressionmatrix) == 38) {
		n <- nrow(regressionmatrix)
		patternHS <- rep(0, n + 37)
		patternIHS <- rep(0, n + 37)
		patternBT <- rep(0, n + 37)
		patternBB <- rep(0, n + 37)
		patternTT <- rep(0, n + 37)
		patternTB <- rep(0, n + 37)
		patternRT <- rep(0, n + 37)
		patternRB <- rep(0, n + 37)
		patternDT <- rep(0, n + 37)
		patternDB <- rep(0, n + 37)
		for(t in 1:n) {
			print(n - t)
			patternHS <- HS(regressionmatrix[t,  ], t, patternHS)
			patternIHS <- IHS(regressionmatrix[t,  ], t, patternIHS)
			patternBT <- BT(regressionmatrix[t,  ], t, patternBT)
			patternBB <- BB(regressionmatrix[t,  ], t, patternBB)
			patternTT <- TT(regressionmatrix[t,  ], t, patternTT)
			patternTB <- TB(regressionmatrix[t,  ], t, patternTB)
			patternRT <- RT(regressionmatrix[t,  ], t, patternRT)
			patternRB <- RB(regressionmatrix[t,  ], t, patternRB)
			patternDT <- DT(regressionmatrix[t,  ], t, patternDT)
			patternDB <- DB(regressionmatrix[t,  ], t, patternDB)
		}
		return(matrix(c(patternHS, patternIHS, patternBT, patternBB, patternTT, patternTB, patternRT, patternRB, 
			patternDT, patternDB), 10, n + 37, byrow = T))
	}
	else print("Invalid regression matrix.")
}
