3
$\begingroup$

am writing a seminar on estimation of population ratio using delta method and am having problem on the literature review. i have written the introduction but need help on on how to use delta method to estimate the population ratio and the methodology

$\endgroup$
4
  • $\begingroup$ What ratio is this? A ratio of means? A population attributable fraction? $\endgroup$
    – AdamO
    Commented Jun 12, 2018 at 15:53
  • $\begingroup$ @gung and others - this question is terribly worded I agree, but Lucas's answer is the only one on the site with the full derivation and explanation - this question deserves to be better worded and re-opened $\endgroup$ Commented Feb 13, 2020 at 3:02
  • $\begingroup$ @XavierBourretSicotte maybe a conversation I wasn't privy to but if my answer would be better suited on another page/question so that this one may be closed please let me know. I agree question is not clearly stated. My answer was written assuming the OP is asking for a ratio of means. $\endgroup$ Commented Feb 14, 2020 at 20:16
  • $\begingroup$ Yes - this question (mine :) ) is a good place ! stats.stackexchange.com/questions/398436/… . So what happened with my question above is that I used a combination of bootstrap and delta method (using your exact formula bellow) to validate the distribution of the metric which was indeed normal but had a different variance to the Binomial $\endgroup$ Commented Feb 15, 2020 at 1:17

1 Answer 1

14
$\begingroup$

The multivariate delta method has a heuristic justification here: https://en.wikipedia.org/wiki/Delta_method#Multivariate_delta_method. For the multivariate delta method you have a specific function $f$ that takes a vector argument which is $p$ dimensional and maps this to a $k$ dimensional space. In the case of a ratio estimator $p=2$ and $k=1$. The function $f$ is

$$f\left(\begin{bmatrix} \bar{y} \\ \bar{x} \\ \end{bmatrix}\right) = \bar{y}/\bar{x}$$ Now what are needed are a few more quantities, the first is:

$$f(\vec{\mu})=f\left(\begin{bmatrix} \mu_{y} \\ \mu_{x} \\ \end{bmatrix}\right) = \mu_{y}/\mu_{x}$$

These are the $h(B)$ and $h(\beta)$ respectively in notation in the Wikipedia link.
Next you need the vector of partial derivatives of $f(\vec{\mu})$, this is:

$$\nabla f(\vec{\mu})=\begin{bmatrix} \frac{1}{\mu_{x}} \\ \frac{-\mu_{y}}{\mu_{x}^2} \\ \end{bmatrix}$$ Also we need the variance covariance matrix of the vector

$$\begin{bmatrix} \bar{y} \\ \bar{x} \\ \end{bmatrix}$$ which is

$$\begin{bmatrix} \sigma^2_{y}/n & \sigma_{yx} \\ \sigma_{yx} & \sigma^2_{x}/n\\ \end{bmatrix}.$$ Note this variance-covariance matrix is the $\Sigma/n$ in the Wikipedia notation. For a proof that $\mathbb{C}ov(\bar{y},\bar{x}) =\mathbb{C}ov(x,y)$ see Estimating the covariance of the means from two samples? Now the only thing left is to calculate the quadratic form:

$$\nabla f(\vec{\mu})^T\begin{bmatrix} \sigma^2_{y}/n & \sigma_{yx} \\ \sigma_{yx} & \sigma^2_{x}/n\\ \end{bmatrix}\nabla f(\vec{\mu}) = \begin{bmatrix} \frac{1}{\mu_{x}} \\ \frac{-\mu_{y}}{\mu_{x}^2} \\ \end{bmatrix}^T \begin{bmatrix} \sigma^2_{y}/n & \sigma_{yx} \\ \sigma_{yx} & \sigma^2_{x}/n\\ \end{bmatrix} \begin{bmatrix} \frac{1}{\mu_{x}} \\ \frac{-\mu_{y}}{\mu_{x}^2} \\ \end{bmatrix}.$$

Which when I worked this out gives you the equation:

$$\sigma^2_R=\frac{\sigma_y^2}{n\mu_x^2} - 2\frac{\mu_y\sigma_{xy}}{\mu_x^3}+\frac{\sigma^2_x\mu_y^2}{n\mu_x^4},$$

where this quantity is the variance of the delta method normal distribution. Putting this altogether gives us that

$$\left(\frac{\bar{y}}{\bar{x}}-\frac{\mu_y}{\mu_x}\right) \sim N(0,\sigma^2_R)$$

So you can estimate the ratio of the population means by the ratio of the sample means provided you can estimate variances and the covariance, or equivalently, the correlation $\rho = \frac{\sigma_{xy}}{\sigma_x\sigma_y}$, by substitution, $\sigma_x\sigma_y\rho = \sigma_{xy}$. This is how the delta method is most commonly used in the derivation of the ratio estimator distribution.

$\endgroup$
9
  • $\begingroup$ I believe that the evaluation of the quadratic form should give a slightly different result. $$\sigma^2_R=\frac{\sigma_y^2}{n\mu_x^2} - 2\frac{\sigma_{xy}\mu_y}{\mu_x^3}+\frac{\sigma^2_x\mu_y^2}{n\mu_x^4}$$ $\endgroup$ Commented Mar 28, 2019 at 19:19
  • $\begingroup$ @MartijnWeterings you are correct that is a typo going from the quadratic form to the equation for $\sigma^2_R$. Will fix that now, thanks for catching my error. $\endgroup$ Commented Mar 29, 2019 at 0:52
  • $\begingroup$ @LucasRoberts do you need the $sqrt(n)$ in front of the expression (y¯x¯−μyμx)∼N(0,σ2R)? Isn't $1/n$ already built into σ2R? $\endgroup$
    – Amazonian
    Commented Oct 26, 2019 at 20:05
  • $\begingroup$ @Amazonian No that is a copy paste error, I've removed it. Also, you can use $\LaTeX$ in comments just as in the regular posts if you'd like your math rendered more legibly. Also please keep in mind this is an approximation and if you decide to use your mileage may vary e.g. if you do not have finite first or second moments for $X$ or $Y$ the approximation is likely not a good one. $\endgroup$ Commented Oct 26, 2019 at 23:58
  • $\begingroup$ Should it be possible to get a negative number from this formulation? Not sure what I'm doing wrong.. this is what I have coded out in Python: (var_y / (n * (mean_x ** 2))) - ((2 * mean_y * cov_xy) / (mean_x ** 3)) + ((var_x * (mean_y ** 2)) / (n * (mean_x ** 4))) $\endgroup$
    – getup8
    Commented Sep 6, 2022 at 5:00