I have read a couple of places that it is possible to do a $2D$ (or $3D$) maximum likelihood fit, but I can't seem to understand how this would work. Suppose I'm considering a probability distribution function depending on $2$ observable variables, $x$ and $y$, given a multi-dimensional set of variable parameters, $\alpha, \beta, \gamma$, ..., PDF($x_i,y_i$|$\alpha, \beta, \gamma, \ldots$).
I would say that I have a good understanding of the ideas behind a maximum likelihood fit, but for some reason I cannot wrap my head around a multi-dimensional maximum likelihood fit.
How do I understand how a $2D$ maximum likelihood fit to determine the best estimate of the parameters $\alpha, \beta, \gamma$ works?