17
$\begingroup$

I have been trying to brush up on my stats knowledge, especially in relation to Sample size determination and Statistical Power Analysis. But it seems that the more I read the more I need to read.

Anyway I found a tool called G*Power which seems to do everything I need but I am having an issue understanding Noncentrality Parameter, what is it, what does it do, what would be a suggested value etc?

The information on wikipedia etc is either incomplete or I'm not doing a very good job understanding it.

I am conducting a series of two tailed z-tests if that is any help.

p.s. Can anyone add better tags to this question?

$\endgroup$

1 Answer 1

15
$\begingroup$

In power calculations, we calibrate tests using knowledge of what the sampling distribution of the test statistic would be under the null hypothesis. Usually, it follows a $\chi^2$ or normal distribution. This allows you to calculate "critical values" for which, values exceeding this are deemed to be too highly inconsistent with what would be expected if the null were true.

The power of a statistical test is calculated by specifying the probability model for the data generating process under an alternative hypothesis, and calculating the sampling distribution for the same test statistic. This now takes on a different distribution.

For test statistics having a $\chi^2$ distribution under the null, they take a non-central $\chi^2$ distribution under the alternative that you create. These are very complicated distributions, but standard software can calculate density, distribution, and quantiles for them easily. The trick is that they are a convolution of standard $\chi^2$ densities and Poisson densities. In R, the dchisq, pchisq, and rchisq functions all have an optional ncp argument which is, by default, 0.

If the test statistic has a standard normal distribution under the null hypothesis, it will have a nonzero mean normal distribution under the alternative. Here that mean is the noncentrality parameter. For a t-test under an equal variance assumption, the mean is given by:

$$\delta = \frac{\mu_1 -\mu_2}{\sigma_{pooled}/\sqrt{n}}$$

enter image description here

In either case, data generated according to an alternative hypothesis will have test statistics following some noncentral distribution with noncentrality parameter ($\delta$). The $\delta$ is a sometimes unknown, often complicated function of the other data generating parameters.

$\endgroup$
4
  • 1
    $\begingroup$ I get why random sampling would lead to a normally distributed mean if the null hypothesis were true (your black line). But the web has given me conflicting descriptions of the distribution under the alternative (i.e., when $\mu_2$ is assumed to be different to $\mu_1$) - in yours it is also normal (red line) but e.g., real-statistics.com has shown it to be skewed (see image half-way down the page). Surely, I've missed a trick. Can you help clear things up for me? $\endgroup$
    – Ben
    Commented Apr 25, 2018 at 12:14
  • 1
    $\begingroup$ @ben I didn't draw a non-central T, I drew the power of a statistical test (red area, shaded). The non-central Chi-sq distribution describes that area when doing power calculation. $\endgroup$
    – AdamO
    Commented May 5, 2020 at 3:54
  • $\begingroup$ @AdamO Neither of those distributions look like chi-squared. Why not draw a t-distribution with an ncp. It's not that difficult with R which I guessing you used. The only reason I "here" is that there's an SO question that cites this Q. $\endgroup$
    – DWin
    Commented Nov 14, 2023 at 1:50
  • $\begingroup$ @AdamO, is mu_1 always zero, or can it take the value of the starting mean? I'm trying to calculate the ncp in R for a power analysis on detecting a % change in a mean value (density). Is my starting mean value = mu_1 in this case? $\endgroup$
    – Nate
    Commented Nov 16, 2023 at 0:48

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.