function(coefs, tmat)
{
	K <- ncol(tmat) - 1
	mysum <- 0
	for(i in 1:K) {
		mysum <- mysum + coefs[i] * tmat[i, K + 1]
	}
	mysum <- mysum/tmat[K + 1, K + 1]
	mysum
}
