def calcR2(predTst, yTest, axis=<NUM_LIT:0>):
rss = np.sum((yTest - predTst) ** <NUM_LIT:2>, axis=axis)<EOL>tss = np.sum((yTest - yTest.mean()) ** <NUM_LIT:2>, axis=axis)<EOL>return <NUM_LIT:1> - rss/tss<EOL>
calculate coefficient of determination. Assumes that axis=0 is time

        Parameters
        ----------
        predTst : np.array, predicted reponse for yTest
        yTest : np.array, acxtually observed response for yTest
        Returns
        -------
        aryFunc : np.array
            R2