Skip to main content
Notice removed Authoritative reference needed by Numbermind
Bounty Ended with R Carnell's answer chosen by Numbermind
added 168 characters in body
Source Link
Numbermind
  • 217
  • 1
  • 7
  • 21

EDIT: I calculated the eta-square effect size in R and got a value of 0.39, which is a lot smaller than the previous one (1.0096). This would give a minimum sample size of 19 per group, i.e, 76 in total. This sample size sounds more correct to me. Should I use eta square and why? However, can I even use the 0.39 value in the "f" variable input for the pwr.anova.test? I think the "f" argument there is only for cohen's f and not for eta square.

EDIT: I calculated the eta-square effect size in R and got a value of 0.39, which is a lot smaller than the previous one (1.0096). This would give a minimum sample size of 19 per group, i.e, 76 in total. This sample size sounds more correct to me. Should I use eta square and why?

EDIT: I calculated the eta-square effect size in R and got a value of 0.39, which is a lot smaller than the previous one (1.0096). This would give a minimum sample size of 19 per group, i.e, 76 in total. This sample size sounds more correct to me. Should I use eta square and why? However, can I even use the 0.39 value in the "f" variable input for the pwr.anova.test? I think the "f" argument there is only for cohen's f and not for eta square.

Notice added Authoritative reference needed by Numbermind
Bounty Started worth 100 reputation by Numbermind
Improved information
Source Link
Numbermind
  • 217
  • 1
  • 7
  • 21

EDIT: I calculated the eta-square effect size in R and got a value of 0.39, which is a lot smaller than the previous one (1.0096). This would give a minimum sample size of 19 per group, i.e, 76 in total. This sample size sounds more correct to me. Should I use eta square and why?

> outcome <- c(6.3,2.8,7.8,7.9,4.9,9.9,4.1,3.9,6.3,6.9,5.1,2.9,3.6,5.7,4.5,1.0,2.8,4.8,3.9,1.6)  
# data
> treatment1 <- factor( c( 1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4 ))        
# grouping variable
> anova1 <- aov( outcome ~ treatment1 )                # run the ANOVA
> summary( anova1 )                                    # print the ANOVA 
table
            Df Sum Sq Mean Sq F value Pr(>F)  
treatment1   3  37.13  12.375    3.46 0.0414 *
Residuals   16  57.22   3.576                 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> etaSquared( anova1 ) 
             eta.sq eta.sq.part
 treatment1 0.3935058   0.3935058

EDIT: I calculated the eta-square effect size in R and got a value of 0.39, which is a lot smaller than the previous one (1.0096). This would give a minimum sample size of 19 per group, i.e, 76 in total. This sample size sounds more correct to me. Should I use eta square and why?

> outcome <- c(6.3,2.8,7.8,7.9,4.9,9.9,4.1,3.9,6.3,6.9,5.1,2.9,3.6,5.7,4.5,1.0,2.8,4.8,3.9,1.6)  
# data
> treatment1 <- factor( c( 1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4 ))        
# grouping variable
> anova1 <- aov( outcome ~ treatment1 )                # run the ANOVA
> summary( anova1 )                                    # print the ANOVA 
table
            Df Sum Sq Mean Sq F value Pr(>F)  
treatment1   3  37.13  12.375    3.46 0.0414 *
Residuals   16  57.22   3.576                 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> etaSquared( anova1 ) 
             eta.sq eta.sq.part
 treatment1 0.3935058   0.3935058
added 300 characters in body
Source Link
Numbermind
  • 217
  • 1
  • 7
  • 21

I'm calculating the minimum sample size to conduct a One-way ANOVA test that I will follow with a Tukey-HSD post hoc analysis. I'm using software R to do my calculations.

From the 'pwr' package, I need to give as input the number of groups (k), the effect size (f) the significance level, and the power to get the minimum sample size. However, to get the effect size, one needs a pilot study or a guess. From a small pilot study let's say I got the following trial data:

enter image description here

However, how should I compute the effect size from the pilot study? I've seen several "effect sizes": Cohen d, Eta-Squared, Partial Eta Squared, Omega Squared, Epsilon Squared, etc.. I also read in the following book: Practical Statistical Power Analysis using WebPower and R, that they have an online calculator: https://webpower.psychstat.org/models/means03/effectsize.php which gave me the following:

enter image description here

It is also said in the book "Cohen defined the size of effect as: small 0.1, medium 0.25, and large 0.4", and here I got >1, which is really strange. I don't understand why are there so many different effect sizes as output and what is being used for the "overall effect size". Why are there interaction effect sizes, i.e, the effect size for group $x$ vs group $y$? I thought we only analyze the difference between each pair of groups in a Tukey-HSD test.

For example, if the effect size was f=1.0096, I would get the following:

> pwr.anova.test(k =4 , f =1.0096 , sig.level =0.05 , power =0.80)

 Balanced one-way analysis of variance power calculation 

          k = 4
          n = 3.806971
          f = 1.0096
  sig.level = 0.05
      power = 0.8

  NOTE: n is number in each group

This would give around 4 for each group, i.e, 4*4=16 in total. And The sample is so small because the effect size of 1.0096 is really big. This leads me to another question: is the effect size of 1.0096 correct?

And for a two-way ANOVA, is the effect size "formula" the same?

I'm calculating the minimum sample size to conduct a One-way ANOVA test that I will follow with a Tukey-HSD post hoc analysis. I'm using software R to do my calculations.

From the 'pwr' package, I need to give as input the number of groups (k), the effect size (f) the significance level, and the power to get the minimum sample size. However, to get the effect size, one needs a pilot study or a guess. From a small pilot study let's say I got the following trial data:

enter image description here

However, how should I compute the effect size from the pilot study? I've seen several "effect sizes": Cohen d, Eta-Squared, Partial Eta Squared, Omega Squared, Epsilon Squared, etc.. I also read in the following book: Practical Statistical Power Analysis using WebPower and R, that they have an online calculator: https://webpower.psychstat.org/models/means03/effectsize.php which gave me the following:

enter image description here

I don't understand why are there so many different effect sizes as output and what is being used for the "overall effect size". Why are there interaction effect sizes, i.e, the effect size for group $x$ vs group $y$? I thought we only analyze the difference between each pair of groups in a Tukey-HSD test.

For example, if the effect size was f=1.0096, I would get the following:

> pwr.anova.test(k =4 , f =1.0096 , sig.level =0.05 , power =0.80)

 Balanced one-way analysis of variance power calculation 

          k = 4
          n = 3.806971
          f = 1.0096
  sig.level = 0.05
      power = 0.8

  NOTE: n is number in each group

This would give around 4 for each group, i.e, 4*4=16 in total. And for a two-way ANOVA, is the effect size "formula" the same?

I'm calculating the minimum sample size to conduct a One-way ANOVA test that I will follow with a Tukey-HSD post hoc analysis. I'm using software R to do my calculations.

From the 'pwr' package, I need to give as input the number of groups (k), the effect size (f) the significance level, and the power to get the minimum sample size. However, to get the effect size, one needs a pilot study or a guess. From a small pilot study let's say I got the following trial data:

enter image description here

However, how should I compute the effect size from the pilot study? I've seen several "effect sizes": Cohen d, Eta-Squared, Partial Eta Squared, Omega Squared, Epsilon Squared, etc.. I also read in the following book: Practical Statistical Power Analysis using WebPower and R, that they have an online calculator: https://webpower.psychstat.org/models/means03/effectsize.php which gave me the following:

enter image description here

It is also said in the book "Cohen defined the size of effect as: small 0.1, medium 0.25, and large 0.4", and here I got >1, which is really strange. I don't understand why are there so many different effect sizes as output and what is being used for the "overall effect size". Why are there interaction effect sizes, i.e, the effect size for group $x$ vs group $y$? I thought we only analyze the difference between each pair of groups in a Tukey-HSD test.

For example, if the effect size was f=1.0096, I would get the following:

> pwr.anova.test(k =4 , f =1.0096 , sig.level =0.05 , power =0.80)

 Balanced one-way analysis of variance power calculation 

          k = 4
          n = 3.806971
          f = 1.0096
  sig.level = 0.05
      power = 0.8

  NOTE: n is number in each group

This would give around 4 for each group, i.e, 4*4=16 in total. The sample is so small because the effect size of 1.0096 is really big. This leads me to another question: is the effect size of 1.0096 correct?

And for a two-way ANOVA, is the effect size "formula" the same?

Source Link
Numbermind
  • 217
  • 1
  • 7
  • 21
Loading