tl;dr: I was confused about the precise expression for the quotient of two multivariate Gaussian densities, so I’m writing it up here.
tl;dr: I ported an R function to Python that helps avoid some numerical issues in polynomial regression.
Lorem ipsum blah blah blah. Let’s test some math: \[\sum_{n=1}^\infty 1/n^2 = \frac{\pi^2}{6}.\]
And inline: \(A^{-1}X=b\)
.