When you're fitting a model $y=f(x)$ to data (${x_i, y_i}$) with errorbars on both the independent ($x$) and response ($y$) variables, it's standard that you can define an 'effective variance' when calculating the chi-squared
$$(\delta y_i)^2 + \left(\delta x_i \left(\frac{df}{dx}\right)_{x=x_i} \right)^2$$
Here the differential is evaluated at the value of $x_i$ of interest.
My question is if this can be extended to the case where there are two independent variables $(x, z)$, so that you're fitting a function $y=g(x, z)$. Something like:
$$(\delta y_i)^2 + \left(\delta x_i \left(\frac{dg}{dx}\right)_{x=x_i} \right)^2 + \left(\delta z_i \left(\frac{dg}{dz}\right)_{z=z_i} \right)^2$$