4  ANOVAs

Authors
Affiliations

John T West

The Pennsylvania State University

Rebecca Wagner

The Pennsylvania State University

Published

March 29, 2024

cat.drm.data = cat.drm.data %>%
  pivot_longer(cols = c(cat.drm.FA.Related.Rates, cat.drm.FA.Unrelated.Rates), names_to = "Relatedness", values_to = "FA.Rates") %>%
  mutate(Relatedness = case_when(
    Relatedness == "cat.drm.FA.Related.Rates" ~ "Related", TRUE ~ "Unrelated"))

object.sim.data = object.sim.data %>%
  pivot_longer(cols = c(object.sim.FA.Related.Rates, object.sim.FA.Unrelated.Rates), names_to = "Relatedness", values_to = "FA.Rates") %>%
  mutate(Relatedness = case_when(
    Relatedness == "object.sim.FA.Related.Rates" ~ "Related", TRUE ~ "Unrelated"))

schematic.data = schematic.data %>%
  pivot_longer(cols = c(schematic.FA.Related.Rates, schematic.FA.Unrelated.Rates), names_to = "Relatedness", values_to = "FA.Rates") %>%
  mutate(Relatedness = case_when(
    Relatedness == "schematic.FA.Related.Rates" ~ "Related", TRUE ~ "Unrelated"))

abstract.img.data = abstract.img.data %>%
  pivot_longer(cols = c(abstract.img.FA.Related.Rates, abstract.img.FA.Unrelated.Rates), names_to = "Relatedness", values_to = "FA.Rates") %>%
  mutate(Relatedness = case_when(
    Relatedness == "abstract.img.FA.Related.Rates" ~ "Related", TRUE ~ "Unrelated"))
#factor data
cat.drm.data$Group <- as.factor(cat.drm.data$Group)
cat.drm.data$Relatedness <- as.factor(cat.drm.data$Relatedness)

#anova & anova table output
cat.drm.anova = cat.drm.data %>%
  anova_test(dv = FA.Rates, wid = numeric.id, between = Group, within = Relatedness, effect.size = "pes", detailed = TRUE, type = 2)
Warning: NA detected in rows: 159,160,213,214,269,270,289,290,361,362,393,394,457,458,603,604,611,612,619,620,655,656,661,662.
Removing this rows before the analysis.
get_anova_table(cat.drm.anova)
ANOVA Table (type II tests)

             Effect DFn DFd    SSn   SSd       F        p p<.05   pes
1       (Intercept)   1 349 14.444 8.569 588.237 7.34e-77     * 0.628
2             Group   1 349  0.200 8.569   8.133 5.00e-03     * 0.023
3       Relatedness   1 349  8.997 4.023 780.488 5.07e-91     * 0.691
4 Group:Relatedness   1 349  0.011 4.023   0.914 3.40e-01       0.003
#means & sd
cat.drm.data %>%
  group_by(Group, Relatedness) %>%
  get_summary_stats(FA.Rates, type = "mean_sd")
# A tibble: 4 × 6
  Group Relatedness variable     n  mean    sd
  <fct> <fct>       <fct>    <dbl> <dbl> <dbl>
1 OA    Related     FA.Rates   148 0.232 0.159
2 OA    Unrelated   FA.Rates   148 0.015 0.047
3 YA    Related     FA.Rates   203 0.274 0.186
4 YA    Unrelated   FA.Rates   203 0.041 0.09 
#bar plot for quick visualization
cat.drm.bar.plot = ggplot(cat.drm.data, aes(x = Relatedness, y = FA.Rates, fill = Group)) +
  stat_summary(fun = "mean", geom = "bar", position = "dodge") +
  stat_summary(fun.data = mean_cl_normal, geom = "errorbar", fun.args = list(mult = 1), 
               position = position_dodge(width = 0.90), width = 0.2)
cat.drm.bar.plot
Warning: Removed 24 rows containing non-finite outside the scale range
(`stat_summary()`).
Warning: Removed 24 rows containing non-finite outside the scale range
(`stat_summary()`).

object.sim.data$Group <- as.factor(object.sim.data$Group)
object.sim.data$Relatedness <- as.factor(object.sim.data$Relatedness)

object.sim.anova = object.sim.data %>%
  anova_test(dv = FA.Rates, wid = numeric.id, between = Group, within = Relatedness, effect.size = "pes", detailed = TRUE)
Warning: NA detected in rows: 23,24,151,152,159,160,189,190,191,192,405,406,431,432,603,604,619,620,635,636,661,662.
Removing this rows before the analysis.
get_anova_table(object.sim.anova)
ANOVA Table (type III tests)

             Effect DFn DFd    SSn   SSd        F         p p<.05   pes
1       (Intercept)   1 350 28.204 7.802 1265.310 2.81e-118     * 0.783
2             Group   1 350  0.011 7.802    0.503  4.79e-01       0.001
3       Relatedness   1 350 23.785 5.761 1445.061 2.66e-126     * 0.805
4 Group:Relatedness   1 350  0.048 5.761    2.918  8.90e-02       0.008
object.sim.data %>%
  group_by(Group, Relatedness) %>%
  get_summary_stats(FA.Rates, type = "mean_sd")
# A tibble: 4 × 6
  Group Relatedness variable     n  mean    sd
  <fct> <fct>       <fct>    <dbl> <dbl> <dbl>
1 OA    Related     FA.Rates   147 0.394 0.166
2 OA    Unrelated   FA.Rates   147 0.004 0.018
3 YA    Related     FA.Rates   205 0.385 0.204
4 YA    Unrelated   FA.Rates   205 0.029 0.07 
object.sim.bar.plot = ggplot(object.sim.data, aes(x = Relatedness, y = FA.Rates, fill = Group)) +
  stat_summary(fun = "mean", geom = "bar", position = "dodge") +
  stat_summary(fun.data = mean_cl_normal, geom = "errorbar", fun.args = list(mult = 1), 
               position = position_dodge(width = 0.90), width = 0.2)
object.sim.bar.plot
Warning: Removed 22 rows containing non-finite outside the scale range
(`stat_summary()`).
Warning: Removed 22 rows containing non-finite outside the scale range
(`stat_summary()`).

schematic.data$Group <- as.factor(schematic.data$Group)
schematic.data$Relatedness <- as.factor(schematic.data$Relatedness)

cat.anova = schematic.data %>%
  anova_test(dv = FA.Rates, wid = numeric.id, between = Group, within = Relatedness, effect.size = "pes", detailed = TRUE)
Warning: NA detected in rows: 1,2,3,4,7,8,27,28,33,34,43,44,59,60,63,64,65,66,79,80,81,82,91,92,93,94,99,100,103,104,109,110,111,112,117,118,119,120,121,122,123,124,137,138,161,162,165,166,179,180,183,184,185,186,193,194,195,196,197,198,199,200,203,204,213,214,225,226,233,234,237,238,239,240,241,242,247,248,251,252,253,254,263,264,269,270,273,274,275,276,277,278,279,280,287,288,291,292,293,294,297,298,319,320,331,332,337,338,339,340,347,348,351,352,361,362,365,366,367,368,375,376,379,380,399,400,401,402,407,408,409,410,413,414,415,416,419,420,437,438,455,456,457,458,459,460,463,464,471,472,475,476,477,478,483,484,493,494,495,496,497,498,499,500,505,506,507,508,513,514,523,524,537,538,553,554,557,558,559,560,567,568,575,576,581,582,589,590,591,592,593,594,603,604,619,620,639,640,641,642,649,650,655,656,669,670,677,678,679,680,693,694,701,702,711,712,721,722,725,726.
Removing this rows before the analysis.
get_anova_table(cat.anova)
ANOVA Table (type III tests)

             Effect DFn DFd    SSn   SSd        F         p p<.05   pes
1       (Intercept)   1 251 73.735 9.409 1966.947 9.33e-121     * 0.887
2             Group   1 251  0.094 9.409    2.506  1.15e-01       0.010
3       Relatedness   1 251  5.213 2.878  454.592  2.90e-58     * 0.644
4 Group:Relatedness   1 251  0.286 2.878   24.926  1.12e-06     * 0.090
schematic.data %>%
  group_by(Group, Relatedness) %>%
  get_summary_stats(FA.Rates, type = "mean_sd")
# A tibble: 4 × 6
  Group Relatedness variable     n  mean    sd
  <fct> <fct>       <fct>    <dbl> <dbl> <dbl>
1 OA    Related     FA.Rates   101 0.504 0.159
2 OA    Unrelated   FA.Rates   101 0.248 0.151
3 YA    Related     FA.Rates   152 0.483 0.146
4 YA    Unrelated   FA.Rates   152 0.324 0.168
cat.bar.plot = ggplot(schematic.data, aes(x = Relatedness, y = FA.Rates, fill = Group)) +
  stat_summary(fun = "mean", geom = "bar", position = "dodge") +
  stat_summary(fun.data = mean_cl_normal, geom = "errorbar", fun.args = list(mult = 1), 
               position = position_dodge(width = 0.90), width = 0.2)
cat.bar.plot
Warning: Removed 220 rows containing non-finite outside the scale range
(`stat_summary()`).
Warning: Removed 220 rows containing non-finite outside the scale range
(`stat_summary()`).

#For Related FA by Age
schematicFMRelated<- schematic.data %>% 
  filter(Relatedness == "Related") %>%
  aov(FA.Rates~Group, data = .)
Anova(schematicFMRelated, type = "III")
Anova Table (Type III tests)

Response: FA.Rates
             Sum Sq  Df  F value Pr(>F)    
(Intercept) 25.6264   1 1119.693 <2e-16 ***
Group        0.0260   1    1.137 0.2873    
Residuals    5.7446 251                    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
partial_eta_squared(schematicFMRelated)
      Group 
0.004509371 
#For Unrelated FA by Age
schematicFMUnrelated <- schematic.data %>% 
  filter(Relatedness == "Unrelated") %>%
  aov(FA.Rates~Group, data = .)
Anova(schematicFMUnrelated, type = "III")
Anova Table (Type III tests)

Response: FA.Rates
            Sum Sq  Df F value    Pr(>F)    
(Intercept) 6.2088   1  238.19 < 2.2e-16 ***
Group       0.3537   1   13.57 0.0002815 ***
Residuals   6.5427 251                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
partial_eta_squared(schematicFMUnrelated)
     Group 
0.05129078 
abstract.img.data$Group <- as.factor(abstract.img.data$Group)
abstract.img.data$Relatedness <- as.factor(abstract.img.data$Relatedness)

abstract.img.anova = abstract.img.data %>%
  anova_test(dv = FA.Rates, wid = numeric.id, between = Group, within = Relatedness, effect.size = "pes", detailed = TRUE)
Warning: NA detected in rows: 151,152,275,276,319,320,331,332,367,368,379,380,399,400,409,410,429,430,463,464,497,498,603,604.
Removing this rows before the analysis.
get_anova_table(abstract.img.anova)
ANOVA Table (type III tests)

             Effect DFn DFd    SSn    SSd        F         p p<.05   pes
1       (Intercept)   1 349 89.483 12.445 2509.462 1.93e-161     * 0.878
2             Group   1 349  0.136 12.445    3.804  5.20e-02       0.011
3       Relatedness   1 349 25.159  3.233 2715.462 1.02e-166     * 0.886
4 Group:Relatedness   1 349  0.261  3.233   28.180  1.97e-07     * 0.075
abstract.img.data %>%
  group_by(Group, Relatedness) %>%
  get_summary_stats(FA.Rates, type = "mean_sd")
# A tibble: 4 × 6
  Group Relatedness variable     n  mean    sd
  <fct> <fct>       <fct>    <dbl> <dbl> <dbl>
1 OA    Related     FA.Rates   150 0.586 0.136
2 OA    Unrelated   FA.Rates   150 0.164 0.147
3 YA    Related     FA.Rates   201 0.519 0.151
4 YA    Unrelated   FA.Rates   201 0.175 0.161
abstract.img.bar.plot = ggplot(abstract.img.data, aes(x = Relatedness, y = FA.Rates, fill = Group)) +
  stat_summary(fun = "mean", geom = "bar", position = "dodge") +
  stat_summary(fun.data = mean_cl_normal, geom = "errorbar", fun.args = list(mult = 1), 
               position = position_dodge(width = 0.90), width = 0.2)
abstract.img.bar.plot
Warning: Removed 24 rows containing non-finite outside the scale range
(`stat_summary()`).
Warning: Removed 24 rows containing non-finite outside the scale range
(`stat_summary()`).

#For Related FA by Age
nonsemFMRelated<- abstract.img.data %>% 
  filter(Relatedness == "Related") %>%
  aov(FA.Rates~Group, data = .)
Anova(nonsemFMRelated, type = "III")
Anova Table (Type III tests)

Response: FA.Rates
            Sum Sq  Df  F value    Pr(>F)    
(Intercept) 51.464   1 2457.952 < 2.2e-16 ***
Group        0.387   1   18.462  2.25e-05 ***
Residuals    7.307 349                       
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
partial_eta_squared(nonsemFMRelated)
     Group 
0.05024312 
#For Unrelated FA by Age
nonsemFMUnrelated <- abstract.img.data %>% 
  filter(Relatedness == "Unrelated") %>%
  aov(FA.Rates~Group, data = .)
Anova(nonsemFMUnrelated, type = "III")
Anova Table (Type III tests)

Response: FA.Rates
            Sum Sq  Df  F value Pr(>F)    
(Intercept) 4.0380   1 168.3529 <2e-16 ***
Group       0.0102   1   0.4243 0.5152    
Residuals   8.3710 349                    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
partial_eta_squared(nonsemFMUnrelated)
      Group 
0.001214202