5
$\begingroup$

I have two questions :

I am wondering if it is possible to know the size of the sample (n) for a multiple linear regression if I want the tests to be powerful enough ? To me, it seems very difficult in reality to use the calculation of power for the multiple linear model, because it would be necessary to enter valid values for the different parameters. Is it possible to do that with R ? If yes, How can I do this with R ? The R code ?

With R, there is this code for the multiple linear regression :

pwr.f2.test(u = NULL, v = NULL, f2 = NULL, sig.level = 0.05, power = NULL)

with :
u : degrees of freedom for numerator
v : degrees of freedom for denominator
f2 : effect size
sig.level : Significance level (Type I error probability)
power : Power of test (1 minus Type II error probability)

  1. Is it possible to improve this R code to take account of (consider) the size of the sample (n) and the different parameters of the regression ?
  2. What about the size of the sample (n) for a logistic regression (logit regression)? Is it possible to know the size of the sampling knowing some criterion like (sig.level=0.05, power=0.8…)?
$\endgroup$
3
$\begingroup$

Power analysis for multiple regression is quite complex as there are many moving parts and potentially several different tests of interest. The function pwr.f2.test is based on Cohen's book Statistical Power Analysis for the Behavioral Sciences and you can find detailed explanations and many examples there.

The most important insight is that the sample size is already captured by the coefficient v (degrees of freedom for the denominator). Exactly how depends on the details of the model. Consequently, the analysis already takes it into account.

Alternatively, another way to conduct a power analysis is to use simulation. It is particularly attractive for this sort of settings as you can play with each individual aspect of the design. For more on this, see Calculating statistical power (see also G. Jay Kern's post).

Once you get the hang of it, it's also quite easy to extend the simulation approach to numerous other tests. For logistic regression, I am not sure if there is any specific function in the pwr package but there is one in G*Power. I don't remember ever using it so I can't comment further on this part of the question.

$\endgroup$
  • $\begingroup$ Thanks for your very clear answer. What about the R script for a logistic regression (to identify the sample size) ? Thanks $\endgroup$ – varin sacha Jul 8 '13 at 17:43
  • $\begingroup$ I edited my answer with some remarks on that. $\endgroup$ – Gala Jul 8 '13 at 18:07
1
$\begingroup$

The best treatment of power for logistic regression I have seen was in Breslow and Day's volume 2 of "Statistical Methods In Cancer Research". The starting point is to realize that the simplest case is a binomial test and that the variance depends on N and the case proportion by way of var(t) = N*p*(1-p). Rather delightful is the fact that the publisher now makes the entirety of both volumes available online and the "Design Considerations" chapter is at the link I embedded. See section 5.6

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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