create likelihood function that takes v1, v2, and a visit as input

[davidwhogg]

Notes to selves: The likelihood is something like:

-0.5 * chi_squared - 0.5 * detATA

where detATAinv is the determinant of the design matrix outer producted with itself. With the covariance matrix in there. And/or maybe with a sign error