Alpha power estimated by FOOOF
Fondecyt 2017
Alvaro Rivera-Rei
05 septiembre, 2025
rm(list = ls()) # clean workspace
try(dev.off(), silent = TRUE) # close all plots
library(tidyverse)
library(afex)
library(lmerTest)
library(emmeans)
library(GGally)
library(easystats)theme_set(
theme_minimal()
)
a_posteriori <- function(afex_aov, sig_level = .05) {
factors <- as.list(rownames(afex_aov$anova_table))
for (j in 1:length(factors)) {
if (grepl(":", factors[[j]])) {
factors[[j]] <- unlist(strsplit(factors[[j]], ":"))
}
}
p_values <- afex_aov$anova_table$`Pr(>F)`
for (i in 1:length(p_values)) {
if (p_values[i] <= sig_level) {
cat(rep("_", 80), '\n', sep = "")
print(emmeans(afex_aov, factors[[i]], contr = "pairwise"))
} else {
print(emmeans(afex_aov, factors[[i]]))
}
}
}
a_posteriori_lmer <- function(lmer_obj, sig_level = .05) {
anova_lmer <- anova(lmer_obj)
factors <- as.list(row.names(anova_lmer))
for (j in 1:length(factors)) {
if (grepl(':', factors[[j]])) {
factors[[j]] <- unlist(strsplit(factors[[j]], ':'))
}
}
p_values <- anova_lmer$`Pr(>F)`
for (i in 1:length(p_values)) {
if (p_values[i] <= sig_level) {
cat(rep('_', 60), '\n', sep = '')
print(emmeans(lmer_obj, factors[[i]], contr = 'pairwise'))
}
}
}# master_dir <- '/home/alvaro/Insync/Drive/00LABWORKS_drive/Huepe/Fdcyt_2017/resting_huepe'
# data_dir <- paste(master_dir, 'FAA_results', sep = '/')
# alpha_power_data_name <- paste(data_dir, 'foof_data_2_to_48_Hz.csv', sep='/')
alpha_power_data <- read.table('foof_data_2_to_48_Hz.csv', header = TRUE, strip.white = TRUE, sep = ",")
alpha_power_data$Dataset <- factor(alpha_power_data$Dataset)
alpha_power_data$Electrode <- factor(alpha_power_data$Electrode)
alpha_power_data$Subject <- factor(alpha_power_data$Subject)
alpha_power_data$hemisphere[alpha_power_data$Electrode %in% c('E093-Fp1', 'E092-AF3a', 'E094-AF7', 'E089-F1a', 'E100-F3a', 'E101-F5a', 'E103-F7', 'E088-FC1a')] <- 'Left'
alpha_power_data$hemisphere[alpha_power_data$Electrode %in% c('E080-Fp2', 'E079-AF4a', 'E072-AF8', 'E076-F2a', 'E068-F4a', 'E069-F6a', 'E071-F8', 'E075-FC2a')] <- 'Right'
alpha_power_data$hemisphere <- factor(alpha_power_data$hemisphere)
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E093-Fp1' , 'E080-Fp2')] <- 'Fp1-Fp2'
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E092-AF3a', 'E079-AF4a')] <- 'AF3-AF4'
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E094-AF7' , 'E072-AF8')] <- 'AF7-AF8'
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E089-F1a' , 'E076-F2a')] <- 'F1-F2'
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E100-F3a' , 'E068-F4a')] <- 'F3-F4'
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E101-F5a' , 'E069-F6a')] <- 'F5-F6'
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E103-F7' , 'E071-F8')] <- 'F7-F8'
alpha_power_data$electrode_pair[alpha_power_data$Electrode %in% c('E088-FC1a', 'E075-FC2a')] <- 'FC1-FC2'
alpha_power_data$electrode_pair <- factor(alpha_power_data$electrode_pair, levels = c('Fp1-Fp2', 'AF3-AF4', 'AF7-AF8', 'F1-F2', 'F3-F4','F5-F6', 'F7-F8', 'FC1-FC2'))
group_id <- read.table("/home/alvaro/Insync/Drive/00EEG/Proyectos/Huepe/fdcyt_2017/Registro-Evaluaciones-FDCYT-DH-2017 - General ANONIMO.csv",
sep = ",", header = TRUE, col.names = c("full.id", "Subject", "Sex", "Group", "Stress"))
group_id$Sex <- factor(group_id$Sex)
group_id$Group <- factor(group_id$Group)
levels(group_id$Sex) <- list(female = "F", male = "M")
levels(group_id$Group) <- list(Control = "CN", "low-SES" = "EX")
group_id <- group_id[c('Subject', 'Group', 'Sex')]
alpha_power_data <- merge(alpha_power_data, group_id, by = 'Subject')
write.csv(alpha_power_data, 'alpha_power_data_clean.csv', row.names = FALSE)
asymmetry_Fp2_Fp1 <- c()
asymmetry_AF4_AF3 <- c()
asymmetry_AF8_AF7 <- c()
asymmetry_F2_F1 <- c()
asymmetry_F4_F3 <- c()
asymmetry_F6_F5 <- c()
asymmetry_F8_F7 <- c()
asymmetry_FC2_FC1 <- c()
Subject <- unique(alpha_power_data$Subject)
for (subj in Subject) {
subject_data <- subset(alpha_power_data, Subject == subj)
asymmetry_Fp2_Fp1 <- c(asymmetry_Fp2_Fp1, subject_data[which(subject_data$Electrode == 'E080-Fp2') , 5] - subject_data[which(subject_data$Electrode=='E093-Fp1') , 5])
asymmetry_AF4_AF3 <- c(asymmetry_AF4_AF3, subject_data[which(subject_data$Electrode == 'E079-AF4a'), 5] - subject_data[which(subject_data$Electrode=='E092-AF3a'), 5])
asymmetry_AF8_AF7 <- c(asymmetry_AF8_AF7, subject_data[which(subject_data$Electrode == 'E072-AF8') , 5] - subject_data[which(subject_data$Electrode=='E094-AF7') , 5])
asymmetry_F2_F1 <- c(asymmetry_F2_F1 , subject_data[which(subject_data$Electrode == 'E076-F2a') , 5] - subject_data[which(subject_data$Electrode=='E089-F1a') , 5])
asymmetry_F4_F3 <- c(asymmetry_F4_F3 , subject_data[which(subject_data$Electrode == 'E068-F4a') , 5] - subject_data[which(subject_data$Electrode=='E100-F3a') , 5])
asymmetry_F6_F5 <- c(asymmetry_F6_F5 , subject_data[which(subject_data$Electrode == 'E069-F6a') , 5] - subject_data[which(subject_data$Electrode=='E101-F5a') , 5])
asymmetry_F8_F7 <- c(asymmetry_F8_F7 , subject_data[which(subject_data$Electrode == 'E071-F8') , 5] - subject_data[which(subject_data$Electrode=='E103-F7') , 5])
asymmetry_FC2_FC1 <- c(asymmetry_FC2_FC1, subject_data[which(subject_data$Electrode == 'E075-FC2a'), 5] - subject_data[which(subject_data$Electrode=='E088-FC1a'), 5])
}
alpha_asymmetry_data <- data.frame(Subject, asymmetry_Fp2_Fp1, asymmetry_AF4_AF3, asymmetry_AF8_AF7, asymmetry_F2_F1, asymmetry_F4_F3, asymmetry_F6_F5, asymmetry_F8_F7, asymmetry_FC2_FC1)
alpha_asymmetry_data <- merge(alpha_asymmetry_data, group_id, by = 'Subject')
alpha_asymmetry_data_long <- pivot_longer(alpha_asymmetry_data,
cols = asymmetry_Fp2_Fp1:asymmetry_FC2_FC1,
names_to = 'Electrodes' ,
names_prefix = 'asymmetry_',
names_ptypes = factor(),
values_to = 'Asimmetry')
write.csv(alpha_asymmetry_data , 'alpha_asymmetry_data_clean.csv' , row.names = FALSE)
write.csv(alpha_asymmetry_data_long, 'alpha_asymmetry_data_clean_long.csv', row.names = FALSE)1 Spectral decomposition
options(width = 100)
mytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_data)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 19 20 39
low-SES 24 14 38
Sum 43 34 77- Infinity Reference or Reference Electrode Standardization
Technique (REST).
Dong, L., Li, F., Liu, Q., Wen, X., Lai, Y., Xu, P., & Yao, D. (2017). MATLAB Toolboxes for Reference Electrode Standardization Technique (REST) of Scalp EEG. Frontiers in Neuroscience, 11(OCT), 601. https://doi.org/10.3389/fnins.2017.00601 - 120 consecutive segments, 5 seconds each.
- PSD computed with Welch’s method.
1.1 Scalp Map, mean power 9-11 Hz
1.2 PSD topography, 1 to 55 Hz, grand average
1.3 PSD topography, 4 to 30 Hz, grand average
1.4 Frontal Electrodes, 4 to 30 Hz, grand average
- 1 standard error bands
2 Parameterization of neural power spectra
Each individual PSD is regarded as a combination of an aperiodic component, characterized by an offset and a slope, and putative periodic oscillatory peaks:
The process goes like this…
Power Spectrum P Fit: \[P = L + \sum_{n=0}^N G_n\]
Where: \[L = aperiodic\ signal\] \[G_n = Gaussian\ fit\ to\ a\ peak\] \[N = total\ peaks\ extracted\]
Aperiodic (background) Fit L: \[L = b - \log_{10}(k + F^\chi)\]
Where: \[b = broadband\ offset\] \[k = ``knee"\ parameter\] \[F = vector\ of\ frequencies\] \[\chi = aperiodic\ slope\]
Donoghue, T., Haller, M., Peterson, E. J., Varma, P., Sebastian, P.,
Gao, R., Noto, T., Lara, A. H., Wallis, J. D., Knight, R. T., Shestyuk,
A., & Voytek, B. (2020). Parameterizing neural power spectra into
periodic and aperiodic components. Nature Neuroscience,
23(12), 1655–1665. https://doi.org/10.1038/s41593-020-00744-x
2.1 Individual Electrodes
2.2 Entire Subject
3 General Asymmetry Description
Subject asymmetry_Fp2_Fp1 asymmetry_AF4_AF3 asymmetry_AF8_AF7 asymmetry_F2_F1
1 : 1 Min. :-0.626665 Min. :-0.412796 Min. :-0.437897 Min. :-0.21493
2 : 1 1st Qu.:-0.075688 1st Qu.:-0.093804 1st Qu.:-0.106656 1st Qu.:-0.02792
4 : 1 Median : 0.015353 Median : 0.006818 Median :-0.001352 Median : 0.02084
5 : 1 Mean :-0.005934 Mean :-0.009114 Mean :-0.006438 Mean : 0.02021
6 : 1 3rd Qu.: 0.068792 3rd Qu.: 0.069743 3rd Qu.: 0.096106 3rd Qu.: 0.04983
7 : 1 Max. : 0.347914 Max. : 0.380950 Max. : 0.498119 Max. : 0.20971
(Other):71 NA's :17 NA's :15 NA's :14 NA's :6
asymmetry_F4_F3 asymmetry_F6_F5 asymmetry_F8_F7 asymmetry_FC2_FC1 Group
Min. :-0.240622 Min. :-0.393810 Min. :-0.51034 Min. :-0.20251 Control:39
1st Qu.:-0.078740 1st Qu.:-0.095586 1st Qu.:-0.14622 1st Qu.:-0.02829 low-SES:38
Median : 0.005338 Median :-0.000607 Median :-0.02253 Median : 0.01477
Mean :-0.002026 Mean :-0.028220 Mean :-0.05082 Mean : 0.01972
3rd Qu.: 0.072646 3rd Qu.: 0.065410 3rd Qu.: 0.03538 3rd Qu.: 0.05133
Max. : 0.267287 Max. : 0.268483 Max. : 0.26139 Max. : 0.48270
NA's :10 NA's :12 NA's :17 NA's :5
Sex
female:43
male :34
asymmetry_pairs <- c('asymmetry_Fp2_Fp1', 'asymmetry_AF4_AF3', 'asymmetry_AF8_AF7', 'asymmetry_F2_F1', 'asymmetry_F4_F3', 'asymmetry_F6_F5', 'asymmetry_F8_F7', 'asymmetry_FC2_FC1')
asymmetry_pairs_pairs <- ggpairs(alpha_asymmetry_data,
columns = asymmetry_pairs,
aes(colour = Group, alpha = .25),
progress = FALSE,
lower = list(continuous = wrap("points")))
suppressWarnings(print(asymmetry_pairs_pairs))
4 Alpha Asymmetry = peak_power(right) - peak_power(left)
4.1 Mixed Model
Subject Group Sex Electrodes Asimmetry
1 : 8 Control:312 female:344 Fp2_Fp1: 77 Min. :-0.626665
2 : 8 low-SES:304 male :272 AF4_AF3: 77 1st Qu.:-0.072283
4 : 8 AF8_AF7: 77 Median : 0.006624
5 : 8 F2_F1 : 77 Mean :-0.006713
6 : 8 F4_F3 : 77 3rd Qu.: 0.065228
7 : 8 F6_F5 : 77 Max. : 0.498119
(Other):568 (Other):154 NA's :96 alpha_asymmetry_lmer = lmer(Asimmetry ~ Group*Sex + (1|Subject) + (1|Electrodes), alpha_asymmetry_data_long)
afex_plot(
alpha_asymmetry_lmer,
x = "Group",
trace = "Sex",
id = "Subject",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: Asimmetry ~ Group * Sex + (1 | Subject) + (1 | Electrodes)
Data: alpha_asymmetry_data_long
REML criterion at convergence: -690.4
Scaled residuals:
Min 1Q Median 3Q Max
-4.6138 -0.4514 -0.0109 0.5262 4.0657
Random effects:
Groups Name Variance Std.Dev.
Subject (Intercept) 0.0035418 0.05951
Electrodes (Intercept) 0.0003694 0.01922
Residual 0.0126293 0.11238
Number of obs: 520, groups: Subject, 73; Electrodes, 8
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 0.01270 0.01906 65.12076 0.666 0.5077
Grouplow-SES -0.06800 0.02368 75.21783 -2.871 0.0053 **
Sexmale 0.01345 0.02463 71.88706 0.546 0.5866
Grouplow-SES:Sexmale 0.04391 0.03493 70.53360 1.257 0.2129
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Correlation of Fixed Effects:
(Intr) Gr-SES Sexmal
Grouplw-SES -0.702
Sexmale -0.675 0.543
Grplw-SES:S 0.476 -0.678 -0.705# Indices of model performance
AIC | AICc | BIC | R2 (cond.) | R2 (marg.) | ICC | RMSE | Sigma
--------------------------------------------------------------------------
-676.4 | -676.2 | -646.6 | 0.284 | 0.062 | 0.236 | 0.106 | 0.112ANOVA-like table for random-effects: Single term deletions
Model:
Asimmetry ~ Group + Sex + (1 | Subject) + (1 | Electrodes) + Group:Sex
npar logLik AIC LRT Df Pr(>Chisq)
<none> 7 345.20 -676.40
(1 | Subject) 6 319.55 -627.10 51.301 1 7.924e-13 ***
(1 | Electrodes) 6 342.43 -672.85 5.552 1 0.01846 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
Group 0.087784 0.087784 1 70.516 6.9508 0.01030 *
Sex 0.051871 0.051871 1 70.568 4.1072 0.04648 *
Group:Sex 0.019954 0.019954 1 70.534 1.5799 0.21292
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1# Effect Size for ANOVA (Type III)
Parameter | Omega2 (partial) | 95% CI | Interpretation
------------------------------------------------------------
Group | 0.08 | [0.00, 0.22] | medium
Sex | 0.04 | [0.00, 0.16] | small
Group:Sex | 7.93e-03 | [0.00, 0.09] | very small
- Interpretation rule: field2013____________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Group emmean SE df lower.CL upper.CL
Control 0.0194 0.0141 37.1 -0.0091 0.04794
low-SES -0.0266 0.0141 37.1 -0.0553 0.00202
Results are averaged over the levels of: Sex
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Control - (low-SES) 0.046 0.0175 67 2.635 0.0104
Results are averaged over the levels of: Sex
Degrees-of-freedom method: kenward-roger
____________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Sex emmean SE df lower.CL upper.CL
female -0.0213 0.0137 35.0 -0.0491 0.00646
male 0.0141 0.0145 38.9 -0.0153 0.04349
Results are averaged over the levels of: Group
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
female - male -0.0354 0.0175 67 -2.026 0.0468
Results are averaged over the levels of: Group
Degrees-of-freedom method: kenward-roger 4.2 Fp2-Fp1 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_Fp2_Fp1", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 17 ID(s), which were removed before analysis:
4, 16, 17, 18, 20, 28, 34, 42, 43, 51, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 12 17 29
low-SES 19 12 31
Sum 31 29 60# alpha_asymmetry_rain <- ggplot(alpha_asymmetry_rep_anova$data$long, aes(y = asymmetry_Fp2_Fp1, x = Group, color = Sex, fill = Sex)) +
# # ggtitle("alpha_asymmetry") +
# ylab("power") +
# stat_halfeye(
# trim = FALSE,
# adjust = 1,
# .width = 0,
# justification = -.15,
# alpha = .5,
# point_colour = NA) +
# # theme(legend.position='none')
# # geom_boxplot(width = .15, alpha = .2, outlier.shape = NA) +
# geom_point(size = 2, alpha = .4, position = position_jitter(width = .05, height = 0))
# suppressWarnings(print(alpha_asymmetry_rain))
alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_Fp2_Fp1
Effect df MSE F ges p.value
1 Group 1, 56 0.02 2.64 .045 .110
2 Sex 1, 56 0.02 1.49 .026 .227
3 Group:Sex 1, 56 0.02 0.22 .004 .639
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control 0.0267 0.0257 56 -0.0248 0.0782
low-SES -0.0318 0.0252 56 -0.0822 0.0186
Results are averaged over the levels of: Sex
Confidence level used: 0.95 G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female -0.0245 0.0252 56 -0.0749 0.0259
male 0.0195 0.0257 56 -0.0321 0.0710
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female 0.01320 0.0394 56 -0.0657 0.092109
low-SES female -0.06224 0.0313 56 -0.1249 0.000472
Control male 0.04022 0.0331 56 -0.0261 0.106519
low-SES male -0.00131 0.0394 56 -0.0802 0.077603
Confidence level used: 0.95 4.3 AF4-AF3 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_AF4_AF3", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 15 ID(s), which were removed before analysis:
1, 4, 16, 17, 18, 20, 42, 43, 51, 52, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 13 18 31
low-SES 18 13 31
Sum 31 31 62alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_AF4_AF3
Effect df MSE F ges p.value
1 Group 1, 58 0.02 2.43 .040 .125
2 Sex 1, 58 0.02 2.22 .037 .141
3 Group:Sex 1, 58 0.02 1.00 .017 .322
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control 0.0216 0.0252 58 -0.0289 0.0721
low-SES -0.0340 0.0252 58 -0.0846 0.0165
Results are averaged over the levels of: Sex
Confidence level used: 0.95 G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female -0.0328 0.0252 58 -0.0834 0.0177
male 0.0204 0.0252 58 -0.0301 0.0709
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female 0.0128 0.0385 58 -0.0642 0.0898
low-SES female -0.0785 0.0327 58 -0.1439 -0.0131
Control male 0.0303 0.0327 58 -0.0351 0.0958
low-SES male 0.0104 0.0385 58 -0.0666 0.0874
Confidence level used: 0.95 4.4 AF8-AF7 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_AF8_AF7", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 14 ID(s), which were removed before analysis:
1, 4, 18, 20, 34, 42, 43, 51, 52, 64, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 13 17 30
low-SES 19 14 33
Sum 32 31 63alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_AF8_AF7
Effect df MSE F ges p.value
1 Group 1, 59 0.03 9.39 ** .137 .003
2 Sex 1, 59 0.03 1.30 .022 .258
3 Group:Sex 1, 59 0.03 2.66 .043 .108
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1________________________________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Group emmean SE df lower.CL upper.CL
Control 0.0671 0.0307 59 0.00565 0.12854
low-SES -0.0631 0.0294 59 -0.12184 -0.00436
Results are averaged over the levels of: Sex
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Control - (low-SES) 0.13 0.0425 59 3.065 0.0033
Results are averaged over the levels of: Sex G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female -0.0223 0.0300 59 -0.0823 0.0378
male 0.0263 0.0301 59 -0.0339 0.0864
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female 0.07751 0.0462 59 -0.0150 0.1700
low-SES female -0.12202 0.0382 59 -0.1985 -0.0455
Control male 0.05669 0.0404 59 -0.0242 0.1376
low-SES male -0.00418 0.0445 59 -0.0933 0.0850
Confidence level used: 0.95 4.5 F2-F1 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_F2_F1", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 6 ID(s), which were removed before analysis:
17, 20, 43, 65, 69, 78
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 17 18 35
low-SES 22 14 36
Sum 39 32 71alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_F2_F1
Effect df MSE F ges p.value
1 Group 1, 67 0.01 0.62 .009 .435
2 Sex 1, 67 0.01 0.19 .003 .666
3 Group:Sex 1, 67 0.01 0.01 <.001 .904
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control 0.0271 0.0125 67 0.00203 0.0521
low-SES 0.0131 0.0127 67 -0.01226 0.0384
Results are averaged over the levels of: Sex
Confidence level used: 0.95 G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female 0.0239 0.0120 67 1.66e-05 0.0478
male 0.0162 0.0132 67 -1.02e-02 0.0426
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female 0.0320 0.0180 67 -0.0039 0.0679
low-SES female 0.0158 0.0158 67 -0.0157 0.0474
Control male 0.0221 0.0175 67 -0.0128 0.0570
low-SES male 0.0103 0.0198 67 -0.0293 0.0499
Confidence level used: 0.95 4.6 F4-F3 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_F4_F3", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 10 ID(s), which were removed before analysis:
4, 18, 20, 34, 42, 43, 52, 65, 69, 78, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 16 17 33
low-SES 20 14 34
Sum 36 31 67alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_F4_F3
Effect df MSE F ges p.value
1 Group 1, 63 0.01 4.37 * .065 .041
2 Sex 1, 63 0.01 0.93 .015 .339
3 Group:Sex 1, 63 0.01 0.98 .015 .326
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1________________________________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Group emmean SE df lower.CL upper.CL
Control 0.0295 0.0194 63 -0.00927 0.0683
low-SES -0.0279 0.0194 63 -0.06666 0.0109
Results are averaged over the levels of: Sex
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Control - (low-SES) 0.0574 0.0274 63 2.091 0.0406
Results are averaged over the levels of: Sex G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female -0.0124 0.0187 63 -0.0498 0.0249
male 0.0140 0.0201 63 -0.0261 0.0542
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female 0.02984 0.0278 63 -0.0258 0.08548
low-SES female -0.05469 0.0249 63 -0.1045 -0.00492
Control male 0.02915 0.0270 63 -0.0248 0.08313
low-SES male -0.00107 0.0298 63 -0.0606 0.05842
Confidence level used: 0.95 4.7 F6-F5 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_F6_F5", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 12 ID(s), which were removed before analysis:
4, 7, 17, 18, 20, 34, 42, 43, 64, 65, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 16 17 33
low-SES 18 14 32
Sum 34 31 65alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_F6_F5
Effect df MSE F ges p.value
1 Group 1, 61 0.02 4.42 * .068 .040
2 Sex 1, 61 0.02 4.04 * .062 .049
3 Group:Sex 1, 61 0.02 0.40 .007 .528
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1________________________________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Group emmean SE df lower.CL upper.CL
Control 0.00688 0.0222 61 -0.0375 0.0513
low-SES -0.05985 0.0227 61 -0.1052 -0.0145
Results are averaged over the levels of: Sex
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Control - (low-SES) 0.0667 0.0317 61 2.102 0.0397
Results are averaged over the levels of: Sex
________________________________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Sex emmean SE df lower.CL upper.CL
female -0.05840 0.0219 61 -0.1022 -0.0146
male 0.00543 0.0230 61 -0.0405 0.0514
Results are averaged over the levels of: Group
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
female - male -0.0638 0.0317 61 -2.011 0.0487
Results are averaged over the levels of: Group
Group Sex emmean SE df lower.CL upper.CL
Control female -0.0150 0.0319 61 -0.0787 0.0487
low-SES female -0.1018 0.0300 61 -0.1619 -0.0418
Control male 0.0287 0.0309 61 -0.0331 0.0905
low-SES male -0.0179 0.0341 61 -0.0860 0.0502
Confidence level used: 0.95 4.8 F8-F7 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_F8_F7", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 17 ID(s), which were removed before analysis:
1, 4, 7, 17, 20, 34, 42, 43, 45, 51, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 11 17 28
low-SES 18 14 32
Sum 29 31 60alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_F8_F7
Effect df MSE F ges p.value
1 Group 1, 56 0.03 0.00 <.001 .983
2 Sex 1, 56 0.03 0.95 .017 .334
3 Group:Sex 1, 56 0.03 0.23 .004 .633
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control -0.0536 0.0306 56 -0.115 0.00771
low-SES -0.0528 0.0282 56 -0.109 0.00373
Results are averaged over the levels of: Sex
Confidence level used: 0.95 G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female -0.0734 0.0303 56 -0.1341 -0.0128
male -0.0329 0.0286 56 -0.0901 0.0243
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female -0.0839 0.0477 56 -0.179 0.0117
low-SES female -0.0630 0.0373 56 -0.138 0.0117
Control male -0.0234 0.0384 56 -0.100 0.0535
low-SES male -0.0425 0.0423 56 -0.127 0.0422
Confidence level used: 0.95 4.9 FC2-FC1 pair
options(width = 100)
alpha_asymmetry_rep_anova = aov_ez("Subject", "asymmetry_FC2_FC1", alpha_asymmetry_data, between = c("Group", "Sex"))G2;H2;Warningh: Missing values for 5 ID(s), which were removed before analysis:
17, 20, 43, 65, 69
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_asymmetry_rep_anova$data$long)
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 18 18 36
low-SES 22 14 36
Sum 40 32 72alpha_asymmetry_afex_plot <-
afex_plot(
alpha_asymmetry_rep_anova,
x = "Group",
trace = "Sex",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_asymmetry_afex_plot))
Anova Table (Type 3 tests)
Response: asymmetry_FC2_FC1
Effect df MSE F ges p.value
1 Group 1, 68 0.01 0.14 .002 .711
2 Sex 1, 68 0.01 4.44 * .061 .039
3 Group:Sex 1, 68 0.01 3.12 + .044 .082
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control 0.0283 0.0149 68 -0.00143 0.0580
low-SES 0.0203 0.0153 68 -0.01013 0.0508
Results are averaged over the levels of: Sex
Confidence level used: 0.95
________________________________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Sex emmean SE df lower.CL upper.CL
female 0.00184 0.0142 68 -0.0265 0.0302
male 0.04677 0.0159 68 0.0150 0.0785
Results are averaged over the levels of: Group
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
female - male -0.0449 0.0213 68 -2.107 0.0388
Results are averaged over the levels of: Group
Group Sex emmean SE df lower.CL upper.CL
Control female 0.0246 0.0211 68 -0.0174 0.0667
low-SES female -0.0210 0.0190 68 -0.0590 0.0170
Control male 0.0319 0.0211 68 -0.0101 0.0739
low-SES male 0.0616 0.0239 68 0.0140 0.1093
Confidence level used: 0.95 5 Alpha peak parameters and aperiodic activity
Subject Dataset Electrode Frequency Amplitude
1 : 16 rest_001_resting.csv: 16 E068-F4a : 77 Min. : 7.043 Min. :0.06018
2 : 16 rest_002_resting.csv: 16 E069-F6a : 77 1st Qu.: 9.152 1st Qu.:0.51268
4 : 16 rest_004_resting.csv: 16 E071-F8 : 77 Median :10.107 Median :0.76769
5 : 16 rest_005_resting.csv: 16 E072-AF8 : 77 Mean : 9.964 Mean :0.80731
6 : 16 rest_006_resting.csv: 16 E075-FC2a: 77 3rd Qu.:10.727 3rd Qu.:1.07556
7 : 16 rest_007_resting.csv: 16 E076-F2a : 77 Max. :12.046 Max. :1.94427
(Other):1136 (Other) :1136 (Other) :770 NA's :145 NA's :145
Width Offset Slope r_squared error
Min. :1.000 Min. :-0.7554 Min. :-0.4384 Min. :0.1708 Min. :0.01429
1st Qu.:1.606 1st Qu.: 0.7167 1st Qu.: 0.8150 1st Qu.:0.9509 1st Qu.:0.03191
Median :2.072 Median : 0.9958 Median : 1.0931 Median :0.9846 Median :0.04258
Mean :2.217 Mean : 1.0009 Mean : 1.0775 Mean :0.9474 Mean :0.04874
3rd Qu.:2.576 3rd Qu.: 1.3058 3rd Qu.: 1.3345 3rd Qu.:0.9928 3rd Qu.:0.05904
Max. :5.000 Max. : 3.1101 Max. : 2.6227 Max. :0.9988 Max. :0.19031
NA's :145
n_peaks f_inf f_sup alpha_inf alpha_sup hemisphere electrode_pair
Min. :0.000 Min. :2 Min. :48 Min. :7 Min. :13 Left :616 Fp1-Fp2:154
1st Qu.:2.000 1st Qu.:2 1st Qu.:48 1st Qu.:7 1st Qu.:13 Right:616 AF3-AF4:154
Median :3.000 Median :2 Median :48 Median :7 Median :13 AF7-AF8:154
Mean :2.962 Mean :2 Mean :48 Mean :7 Mean :13 F1-F2 :154
3rd Qu.:4.000 3rd Qu.:2 3rd Qu.:48 3rd Qu.:7 3rd Qu.:13 F3-F4 :154
Max. :4.000 Max. :2 Max. :48 Max. :7 Max. :13 F5-F6 :154
(Other):308
Group Sex
Control:624 female:688
low-SES:608 male :544
spec_params <- c('Amplitude', 'Frequency', 'Width', 'Offset', 'Slope')
spec_params_pairs <- ggpairs(alpha_power_data,
columns = spec_params,
aes(colour = Group, alpha = .25),
progress = FALSE,
lower = list(continuous = wrap("points")))
suppressWarnings(print(spec_params_pairs))
5.1 Alpha Power
options(width = 100)
alpha_pow_rep_anova = aov_ez("Subject", "Amplitude", alpha_power_data, between = c("Group", "Sex"), within = c("hemisphere", "electrode_pair"))G2;H2;Warningh: Missing values for 21 ID(s), which were removed before analysis:
1, 4, 7, 16, 17, 18, 20, 28, 34, 42, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_pow_rep_anova$data$long) / length(levels(alpha_pow_rep_anova$data$long$electrode_pair)) / length(levels(alpha_pow_rep_anova$data$long$hemisphere))
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 10 17 27
low-SES 17 12 29
Sum 27 29 56alpha_pow_afex_plot <-
afex_plot(
alpha_pow_rep_anova,
x = "Group",
trace = "Sex",
panel = "hemisphere",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_pow_afex_plot))
Anova Table (Type 3 tests)
Response: Amplitude
Effect df MSE F ges p.value
1 Group 1, 52 1.69 2.04 .034 .160
2 Sex 1, 52 1.69 0.49 .008 .487
3 Group:Sex 1, 52 1.69 0.67 .012 .416
4 hemisphere 1, 52 0.02 0.19 <.001 .665
5 Group:hemisphere 1, 52 0.02 6.33 * .002 .015
6 Sex:hemisphere 1, 52 0.02 3.54 + <.001 .065
7 Group:Sex:hemisphere 1, 52 0.02 1.78 <.001 .188
8 electrode_pair 2.91, 151.48 0.04 51.14 *** .056 <.001
9 Group:electrode_pair 2.91, 151.48 0.04 0.12 <.001 .946
10 Sex:electrode_pair 2.91, 151.48 0.04 0.87 .001 .455
11 Group:Sex:electrode_pair 2.91, 151.48 0.04 0.51 <.001 .670
12 hemisphere:electrode_pair 4.46, 231.96 0.01 2.79 * .001 .022
13 Group:hemisphere:electrode_pair 4.46, 231.96 0.01 2.13 + <.001 .071
14 Sex:hemisphere:electrode_pair 4.46, 231.96 0.01 0.52 <.001 .741
15 Group:Sex:hemisphere:electrode_pair 4.46, 231.96 0.01 1.04 <.001 .390
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1
Sphericity correction method: GG $emmeans
Group = Control:
hemisphere emmean SE df lower.CL upper.CL
Left 0.951 0.0659 52 0.819 1.084
Right 0.973 0.0644 52 0.844 1.102
Group = low-SES:
hemisphere emmean SE df lower.CL upper.CL
Left 0.851 0.0623 52 0.726 0.976
Right 0.820 0.0609 52 0.698 0.942
Results are averaged over the levels of: Sex, electrode_pair
Confidence level used: 0.95
$contrasts
Group = Control:
contrast estimate SE df t.ratio p.value
Left - Right -0.0218 0.0152 52 -1.432 0.1582
Group = low-SES:
contrast estimate SE df t.ratio p.value
Left - Right 0.0309 0.0144 52 2.146 0.0365
Results are averaged over the levels of: Sex, electrode_pair ________________________________________________________________________________$emmeans
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 0.856 0.0444 52 0.767 0.945
AF3.AF4 0.897 0.0460 52 0.804 0.989
AF7.AF8 0.811 0.0456 52 0.719 0.902
F1.F2 1.016 0.0458 52 0.924 1.108
F3.F4 0.931 0.0459 52 0.839 1.023
F5.F6 0.867 0.0471 52 0.772 0.961
F7.F8 0.786 0.0488 52 0.688 0.884
FC1.FC2 1.027 0.0442 52 0.938 1.115
Results are averaged over the levels of: Group, Sex, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 -0.0406 0.01110 52 -3.652 0.0131
Fp1.Fp2 - AF7.AF8 0.0452 0.01610 52 2.806 0.1158
Fp1.Fp2 - F1.F2 -0.1600 0.02140 52 -7.465 <.0001
Fp1.Fp2 - F3.F4 -0.0749 0.02330 52 -3.209 0.0439
Fp1.Fp2 - F5.F6 -0.0106 0.02190 52 -0.483 0.9997
Fp1.Fp2 - F7.F8 0.0698 0.02200 52 3.169 0.0485
Fp1.Fp2 - FC1.FC2 -0.1706 0.02540 52 -6.704 <.0001
AF3.AF4 - AF7.AF8 0.0858 0.01550 52 5.542 <.0001
AF3.AF4 - F1.F2 -0.1194 0.01610 52 -7.425 <.0001
AF3.AF4 - F3.F4 -0.0343 0.01850 52 -1.849 0.5910
AF3.AF4 - F5.F6 0.0300 0.01700 52 1.772 0.6414
AF3.AF4 - F7.F8 0.1104 0.01850 52 5.974 <.0001
AF3.AF4 - FC1.FC2 -0.1300 0.02100 52 -6.190 <.0001
AF7.AF8 - F1.F2 -0.2052 0.01920 52 -10.695 <.0001
AF7.AF8 - F3.F4 -0.1201 0.01920 52 -6.240 <.0001
AF7.AF8 - F5.F6 -0.0558 0.01630 52 -3.428 0.0245
AF7.AF8 - F7.F8 0.0246 0.01490 52 1.646 0.7204
AF7.AF8 - FC1.FC2 -0.2158 0.02250 52 -9.589 <.0001
F1.F2 - F3.F4 0.0851 0.01010 52 8.436 <.0001
F1.F2 - F5.F6 0.1494 0.01340 52 11.143 <.0001
F1.F2 - F7.F8 0.2298 0.01800 52 12.774 <.0001
F1.F2 - FC1.FC2 -0.0106 0.00855 52 -1.237 0.9170
F3.F4 - F5.F6 0.0643 0.01080 52 5.961 <.0001
F3.F4 - F7.F8 0.1447 0.01420 52 10.221 <.0001
F3.F4 - FC1.FC2 -0.0957 0.01200 52 -8.009 <.0001
F5.F6 - F7.F8 0.0804 0.01020 52 7.897 <.0001
F5.F6 - FC1.FC2 -0.1600 0.01530 52 -10.489 <.0001
F7.F8 - FC1.FC2 -0.2404 0.01980 52 -12.149 <.0001
Results are averaged over the levels of: Group, Sex, hemisphere
P value adjustment: tukey method for comparing a family of 8 estimates ________________________________________________________________________________$emmeans
hemisphere = Left:
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 0.858 0.0451 52 0.768 0.949
AF3.AF4 0.896 0.0475 52 0.801 0.991
AF7.AF8 0.812 0.0481 52 0.715 0.908
F1.F2 1.006 0.0463 52 0.913 1.099
F3.F4 0.929 0.0474 52 0.834 1.024
F5.F6 0.879 0.0487 52 0.781 0.977
F7.F8 0.815 0.0505 52 0.713 0.916
FC1.FC2 1.014 0.0443 52 0.925 1.103
hemisphere = Right:
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 0.854 0.0458 52 0.762 0.946
AF3.AF4 0.898 0.0464 52 0.804 0.991
AF7.AF8 0.810 0.0460 52 0.718 0.902
F1.F2 1.026 0.0460 52 0.934 1.119
F3.F4 0.933 0.0455 52 0.842 1.024
F5.F6 0.854 0.0470 52 0.760 0.949
F7.F8 0.758 0.0495 52 0.658 0.857
FC1.FC2 1.039 0.0451 52 0.949 1.130
Results are averaged over the levels of: Group, Sex
Confidence level used: 0.95
$contrasts
hemisphere = Left:
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 -0.037387 0.01380 52 -2.716 0.1409
Fp1.Fp2 - AF7.AF8 0.046677 0.01870 52 2.496 0.2199
Fp1.Fp2 - F1.F2 -0.147332 0.02220 52 -6.633 <.0001
Fp1.Fp2 - F3.F4 -0.070463 0.02610 52 -2.696 0.1471
Fp1.Fp2 - F5.F6 -0.020472 0.02390 52 -0.856 0.9886
Fp1.Fp2 - F7.F8 0.043758 0.02450 52 1.786 0.6320
Fp1.Fp2 - FC1.FC2 -0.155621 0.02650 52 -5.868 <.0001
AF3.AF4 - AF7.AF8 0.084064 0.02120 52 3.974 0.0050
AF3.AF4 - F1.F2 -0.109945 0.01710 52 -6.427 <.0001
AF3.AF4 - F3.F4 -0.033076 0.02120 52 -1.561 0.7704
AF3.AF4 - F5.F6 0.016915 0.02010 52 0.842 0.9897
AF3.AF4 - F7.F8 0.081145 0.02250 52 3.607 0.0149
AF3.AF4 - FC1.FC2 -0.118234 0.02240 52 -5.282 0.0001
AF7.AF8 - F1.F2 -0.194009 0.02340 52 -8.301 <.0001
AF7.AF8 - F3.F4 -0.117141 0.02530 52 -4.627 0.0006
AF7.AF8 - F5.F6 -0.067150 0.02080 52 -3.230 0.0415
AF7.AF8 - F7.F8 -0.002920 0.01960 52 -0.149 1.0000
AF7.AF8 - FC1.FC2 -0.202298 0.02600 52 -7.792 <.0001
F1.F2 - F3.F4 0.076869 0.01240 52 6.185 <.0001
F1.F2 - F5.F6 0.126859 0.01490 52 8.532 <.0001
F1.F2 - F7.F8 0.191089 0.01910 52 10.006 <.0001
F1.F2 - FC1.FC2 -0.008289 0.01090 52 -0.763 0.9943
F3.F4 - F5.F6 0.049991 0.01480 52 3.369 0.0288
F3.F4 - F7.F8 0.114221 0.01830 52 6.226 <.0001
F3.F4 - FC1.FC2 -0.085158 0.01510 52 -5.623 <.0001
F5.F6 - F7.F8 0.064230 0.01560 52 4.128 0.0031
F5.F6 - FC1.FC2 -0.135148 0.01740 52 -7.753 <.0001
F7.F8 - FC1.FC2 -0.199379 0.02080 52 -9.592 <.0001
hemisphere = Right:
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 -0.043811 0.01390 52 -3.163 0.0493
Fp1.Fp2 - AF7.AF8 0.043747 0.01770 52 2.469 0.2317
Fp1.Fp2 - F1.F2 -0.172678 0.02580 52 -6.703 <.0001
Fp1.Fp2 - F3.F4 -0.079262 0.02630 52 -3.010 0.0721
Fp1.Fp2 - F5.F6 -0.000652 0.02520 52 -0.026 1.0000
Fp1.Fp2 - F7.F8 0.095885 0.02690 52 3.569 0.0166
Fp1.Fp2 - FC1.FC2 -0.185532 0.02910 52 -6.370 <.0001
AF3.AF4 - AF7.AF8 0.087558 0.01720 52 5.091 0.0001
AF3.AF4 - F1.F2 -0.128867 0.02070 52 -6.211 <.0001
AF3.AF4 - F3.F4 -0.035451 0.02150 52 -1.646 0.7208
AF3.AF4 - F5.F6 0.043159 0.02140 52 2.019 0.4798
AF3.AF4 - F7.F8 0.139696 0.02370 52 5.900 <.0001
AF3.AF4 - FC1.FC2 -0.141721 0.02520 52 -5.615 <.0001
AF7.AF8 - F1.F2 -0.216425 0.02280 52 -9.473 <.0001
AF7.AF8 - F3.F4 -0.123009 0.02140 52 -5.752 <.0001
AF7.AF8 - F5.F6 -0.044398 0.01900 52 -2.340 0.2925
AF7.AF8 - F7.F8 0.052138 0.01900 52 2.743 0.1329
AF7.AF8 - FC1.FC2 -0.229278 0.02640 52 -8.692 <.0001
F1.F2 - F3.F4 0.093416 0.01240 52 7.563 <.0001
F1.F2 - F5.F6 0.172027 0.01760 52 9.752 <.0001
F1.F2 - F7.F8 0.268563 0.02350 52 11.429 <.0001
F1.F2 - FC1.FC2 -0.012854 0.00998 52 -1.287 0.8994
F3.F4 - F5.F6 0.078610 0.01400 52 5.608 <.0001
F3.F4 - F7.F8 0.175147 0.01840 52 9.543 <.0001
F3.F4 - FC1.FC2 -0.106270 0.01550 52 -6.874 <.0001
F5.F6 - F7.F8 0.096536 0.01280 52 7.516 <.0001
F5.F6 - FC1.FC2 -0.184880 0.01940 52 -9.517 <.0001
F7.F8 - FC1.FC2 -0.281417 0.02520 52 -11.187 <.0001
Results are averaged over the levels of: Group, Sex
P value adjustment: tukey method for comparing a family of 8 estimates 5.2 Alpha Frequency
options(width = 100)
alpha_freq_rep_anova = aov_ez("Subject", "Frequency", alpha_power_data, between = c("Group", "Sex"), within = c("hemisphere", "electrode_pair"))G2;H2;Warningh: Missing values for 21 ID(s), which were removed before analysis:
1, 4, 7, 16, 17, 18, 20, 28, 34, 42, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_freq_rep_anova$data$long) / length(levels(alpha_freq_rep_anova$data$long$electrode_pair)) / length(levels(alpha_freq_rep_anova$data$long$hemisphere))
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 10 17 27
low-SES 17 12 29
Sum 27 29 56alpha_freq_afex_plot <-
afex_plot(
alpha_freq_rep_anova,
x = "Group",
trace = "Sex",
panel = "hemisphere",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_freq_afex_plot))
Anova Table (Type 3 tests)
Response: Frequency
Effect df MSE F ges p.value
1 Group 1, 52 14.40 1.80 .031 .186
2 Sex 1, 52 14.40 0.77 .013 .384
3 Group:Sex 1, 52 14.40 0.48 .008 .492
4 hemisphere 1, 52 0.27 0.21 <.001 .646
5 Group:hemisphere 1, 52 0.27 0.06 <.001 .812
6 Sex:hemisphere 1, 52 0.27 0.09 <.001 .767
7 Group:Sex:hemisphere 1, 52 0.27 0.99 <.001 .325
8 electrode_pair 2.58, 134.26 0.25 3.38 * .003 .026
9 Group:electrode_pair 2.58, 134.26 0.25 0.78 <.001 .488
10 Sex:electrode_pair 2.58, 134.26 0.25 0.48 <.001 .671
11 Group:Sex:electrode_pair 2.58, 134.26 0.25 0.41 <.001 .717
12 hemisphere:electrode_pair 3.04, 157.88 0.12 0.74 <.001 .529
13 Group:hemisphere:electrode_pair 3.04, 157.88 0.12 0.73 <.001 .537
14 Sex:hemisphere:electrode_pair 3.04, 157.88 0.12 0.79 <.001 .504
15 Group:Sex:hemisphere:electrode_pair 3.04, 157.88 0.12 1.07 <.001 .366
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1
Sphericity correction method: GG G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control 9.79 0.189 52 9.41 10.2
low-SES 10.14 0.179 52 9.78 10.5
Results are averaged over the levels of: Sex, electrode_pair, hemisphere
Confidence level used: 0.95 G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female 9.85 0.189 52 9.47 10.2
male 10.08 0.179 52 9.72 10.4
Results are averaged over the levels of: Group, electrode_pair, hemisphere
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female 9.59 0.300 52 8.98 10.2
low-SES female 10.11 0.230 52 9.65 10.6
Control male 9.99 0.230 52 9.53 10.5
low-SES male 10.16 0.274 52 9.61 10.7
Results are averaged over the levels of: electrode_pair, hemisphere
Confidence level used: 0.95
hemisphere emmean SE df lower.CL upper.CL
Left 9.96 0.134 52 9.69 10.2
Right 9.97 0.129 52 9.71 10.2
Results are averaged over the levels of: Group, Sex, electrode_pair
Confidence level used: 0.95
Group hemisphere emmean SE df lower.CL upper.CL
Control Left 9.79 0.194 52 9.40 10.2
low-SES Left 10.13 0.184 52 9.76 10.5
Control Right 9.79 0.187 52 9.42 10.2
low-SES Right 10.15 0.177 52 9.80 10.5
Results are averaged over the levels of: Sex, electrode_pair
Confidence level used: 0.95
Sex hemisphere emmean SE df lower.CL upper.CL
female Left 9.84 0.194 52 9.45 10.2
male Left 10.08 0.184 52 9.71 10.4
female Right 9.86 0.187 52 9.49 10.2
male Right 10.08 0.177 52 9.73 10.4
Results are averaged over the levels of: Group, electrode_pair
Confidence level used: 0.95
Group Sex hemisphere emmean SE df lower.CL upper.CL
Control female Left 9.56 0.308 52 8.94 10.2
low-SES female Left 10.11 0.236 52 9.64 10.6
Control male Left 10.01 0.236 52 9.54 10.5
low-SES male Left 10.14 0.281 52 9.57 10.7
Control female Right 9.61 0.297 52 9.02 10.2
low-SES female Right 10.11 0.228 52 9.66 10.6
Control male Right 9.98 0.228 52 9.52 10.4
low-SES male Right 10.19 0.271 52 9.64 10.7
Results are averaged over the levels of: electrode_pair
Confidence level used: 0.95
________________________________________________________________________________
$emmeans
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 10.01 0.129 52 9.75 10.3
AF3.AF4 9.98 0.138 52 9.70 10.3
AF7.AF8 10.01 0.134 52 9.74 10.3
F1.F2 9.92 0.132 52 9.65 10.2
F3.F4 9.91 0.133 52 9.64 10.2
F5.F6 9.97 0.136 52 9.70 10.2
F7.F8 10.03 0.130 52 9.77 10.3
FC1.FC2 9.89 0.133 52 9.62 10.2
Results are averaged over the levels of: Group, Sex, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 0.03655 0.0299 52 1.223 0.9214
Fp1.Fp2 - AF7.AF8 0.00496 0.0246 52 0.202 1.0000
Fp1.Fp2 - F1.F2 0.09511 0.0485 52 1.961 0.5172
Fp1.Fp2 - F3.F4 0.10752 0.0458 52 2.350 0.2874
Fp1.Fp2 - F5.F6 0.03891 0.0466 52 0.836 0.9901
Fp1.Fp2 - F7.F8 -0.01798 0.0377 52 -0.477 0.9997
Fp1.Fp2 - FC1.FC2 0.12626 0.0513 52 2.459 0.2361
AF3.AF4 - AF7.AF8 -0.03159 0.0238 52 -1.326 0.8847
AF3.AF4 - F1.F2 0.05856 0.0476 52 1.230 0.9191
AF3.AF4 - F3.F4 0.07097 0.0449 52 1.581 0.7592
AF3.AF4 - F5.F6 0.00236 0.0379 52 0.062 1.0000
AF3.AF4 - F7.F8 -0.05453 0.0407 52 -1.341 0.8787
AF3.AF4 - FC1.FC2 0.08971 0.0495 52 1.814 0.6137
AF7.AF8 - F1.F2 0.09015 0.0436 52 2.065 0.4501
AF7.AF8 - F3.F4 0.10256 0.0420 52 2.442 0.2436
AF7.AF8 - F5.F6 0.03395 0.0349 52 0.974 0.9762
AF7.AF8 - F7.F8 -0.02294 0.0325 52 -0.706 0.9965
AF7.AF8 - FC1.FC2 0.12130 0.0448 52 2.708 0.1432
F1.F2 - F3.F4 0.01241 0.0198 52 0.627 0.9983
F1.F2 - F5.F6 -0.05619 0.0517 52 -1.086 0.9570
F1.F2 - F7.F8 -0.11309 0.0538 52 -2.103 0.4267
F1.F2 - FC1.FC2 0.03115 0.0138 52 2.259 0.3352
F3.F4 - F5.F6 -0.06861 0.0522 52 -1.314 0.8894
F3.F4 - F7.F8 -0.12550 0.0532 52 -2.359 0.2831
F3.F4 - FC1.FC2 0.01874 0.0208 52 0.899 0.9849
F5.F6 - F7.F8 -0.05689 0.0274 52 -2.075 0.4444
F5.F6 - FC1.FC2 0.08735 0.0495 52 1.764 0.6462
F7.F8 - FC1.FC2 0.14424 0.0545 52 2.644 0.1637
Results are averaged over the levels of: Group, Sex, hemisphere
P value adjustment: tukey method for comparing a family of 8 estimates
Group electrode_pair emmean SE df lower.CL upper.CL
Control Fp1.Fp2 9.84 0.188 52 9.46 10.2
low-SES Fp1.Fp2 10.19 0.178 52 9.83 10.5
Control AF3.AF4 9.81 0.200 52 9.40 10.2
low-SES AF3.AF4 10.15 0.189 52 9.77 10.5
Control AF7.AF8 9.81 0.194 52 9.42 10.2
low-SES AF7.AF8 10.20 0.184 52 9.83 10.6
Control F1.F2 9.73 0.192 52 9.34 10.1
low-SES F1.F2 10.11 0.181 52 9.75 10.5
Control F3.F4 9.72 0.193 52 9.33 10.1
low-SES F3.F4 10.09 0.183 52 9.73 10.5
Control F5.F6 9.83 0.198 52 9.44 10.2
low-SES F5.F6 10.12 0.187 52 9.74 10.5
Control F7.F8 9.90 0.188 52 9.53 10.3
low-SES F7.F8 10.16 0.178 52 9.80 10.5
Control FC1.FC2 9.69 0.193 52 9.30 10.1
low-SES FC1.FC2 10.09 0.182 52 9.72 10.5
Results are averaged over the levels of: Sex, hemisphere
Confidence level used: 0.95
Sex electrode_pair emmean SE df lower.CL upper.CL
female Fp1.Fp2 9.91 0.188 52 9.54 10.3
male Fp1.Fp2 10.11 0.178 52 9.76 10.5
female AF3.AF4 9.89 0.200 52 9.49 10.3
male AF3.AF4 10.07 0.189 52 9.69 10.4
female AF7.AF8 9.91 0.194 52 9.52 10.3
male AF7.AF8 10.11 0.184 52 9.74 10.5
female F1.F2 9.78 0.192 52 9.39 10.2
male F1.F2 10.06 0.181 52 9.70 10.4
female F3.F4 9.77 0.193 52 9.39 10.2
male F3.F4 10.04 0.183 52 9.67 10.4
female F5.F6 9.87 0.198 52 9.47 10.3
male F5.F6 10.08 0.187 52 9.71 10.5
female F7.F8 9.92 0.188 52 9.55 10.3
male F7.F8 10.14 0.178 52 9.78 10.5
female FC1.FC2 9.75 0.193 52 9.36 10.1
male FC1.FC2 10.03 0.182 52 9.66 10.4
Results are averaged over the levels of: Group, hemisphere
Confidence level used: 0.95
Group Sex electrode_pair emmean SE df lower.CL upper.CL
Control female Fp1.Fp2 9.62 0.298 52 9.02 10.2
low-SES female Fp1.Fp2 10.21 0.229 52 9.75 10.7
Control male Fp1.Fp2 10.06 0.229 52 9.60 10.5
low-SES male Fp1.Fp2 10.17 0.272 52 9.63 10.7
Control female AF3.AF4 9.63 0.317 52 9.00 10.3
low-SES female AF3.AF4 10.14 0.243 52 9.66 10.6
Control male AF3.AF4 9.98 0.243 52 9.49 10.5
low-SES male AF3.AF4 10.15 0.290 52 9.57 10.7
Control female AF7.AF8 9.61 0.308 52 8.99 10.2
low-SES female AF7.AF8 10.22 0.236 52 9.74 10.7
Control male AF7.AF8 10.02 0.236 52 9.55 10.5
low-SES male AF7.AF8 10.19 0.281 52 9.62 10.8
Control female F1.F2 9.49 0.304 52 8.88 10.1
low-SES female F1.F2 10.06 0.233 52 9.59 10.5
Control male F1.F2 9.96 0.233 52 9.50 10.4
low-SES male F1.F2 10.16 0.278 52 9.60 10.7
Control female F3.F4 9.49 0.307 52 8.87 10.1
low-SES female F3.F4 10.06 0.235 52 9.59 10.5
Control male F3.F4 9.95 0.235 52 9.47 10.4
low-SES male F3.F4 10.13 0.280 52 9.57 10.7
Control female F5.F6 9.66 0.314 52 9.03 10.3
low-SES female F5.F6 10.08 0.241 52 9.59 10.6
Control male F5.F6 10.01 0.241 52 9.53 10.5
low-SES male F5.F6 10.16 0.286 52 9.58 10.7
Control female F7.F8 9.71 0.299 52 9.12 10.3
low-SES female F7.F8 10.13 0.229 52 9.67 10.6
Control male F7.F8 10.09 0.229 52 9.63 10.6
low-SES male F7.F8 10.19 0.273 52 9.64 10.7
Control female FC1.FC2 9.48 0.306 52 8.87 10.1
low-SES female FC1.FC2 10.02 0.234 52 9.55 10.5
Control male FC1.FC2 9.89 0.234 52 9.42 10.4
low-SES male FC1.FC2 10.16 0.279 52 9.60 10.7
Results are averaged over the levels of: hemisphere
Confidence level used: 0.95
hemisphere electrode_pair emmean SE df lower.CL upper.CL
Left Fp1.Fp2 10.00 0.135 52 9.73 10.3
Right Fp1.Fp2 10.03 0.133 52 9.76 10.3
Left AF3.AF4 9.96 0.137 52 9.69 10.2
Right AF3.AF4 9.99 0.140 52 9.71 10.3
Left AF7.AF8 10.00 0.137 52 9.73 10.3
Right AF7.AF8 10.01 0.133 52 9.75 10.3
Left F1.F2 9.90 0.136 52 9.62 10.2
Right F1.F2 9.94 0.131 52 9.68 10.2
Left F3.F4 9.88 0.139 52 9.60 10.2
Right F3.F4 9.93 0.131 52 9.67 10.2
Left F5.F6 10.00 0.143 52 9.71 10.3
Right F5.F6 9.95 0.136 52 9.67 10.2
Left F7.F8 10.04 0.136 52 9.77 10.3
Right F7.F8 10.02 0.134 52 9.75 10.3
Left FC1.FC2 9.87 0.137 52 9.59 10.1
Right FC1.FC2 9.91 0.132 52 9.64 10.2
Results are averaged over the levels of: Group, Sex
Confidence level used: 0.95
Group hemisphere electrode_pair emmean SE df lower.CL upper.CL
Control Left Fp1.Fp2 9.83 0.197 52 9.44 10.2
low-SES Left Fp1.Fp2 10.17 0.186 52 9.79 10.5
Control Right Fp1.Fp2 9.84 0.193 52 9.46 10.2
low-SES Right Fp1.Fp2 10.22 0.182 52 9.85 10.6
Control Left AF3.AF4 9.78 0.199 52 9.38 10.2
low-SES Left AF3.AF4 10.14 0.188 52 9.76 10.5
Control Right AF3.AF4 9.83 0.204 52 9.42 10.2
low-SES Right AF3.AF4 10.16 0.193 52 9.77 10.5
Control Left AF7.AF8 9.82 0.200 52 9.42 10.2
low-SES Left AF7.AF8 10.19 0.189 52 9.81 10.6
Control Right AF7.AF8 9.81 0.193 52 9.42 10.2
low-SES Right AF7.AF8 10.22 0.183 52 9.85 10.6
Control Left F1.F2 9.74 0.198 52 9.34 10.1
low-SES Left F1.F2 10.05 0.187 52 9.68 10.4
Control Right F1.F2 9.72 0.191 52 9.33 10.1
low-SES Right F1.F2 10.16 0.181 52 9.80 10.5
Control Left F3.F4 9.69 0.202 52 9.29 10.1
low-SES Left F3.F4 10.06 0.191 52 9.68 10.4
Control Right F3.F4 9.74 0.191 52 9.36 10.1
low-SES Right F3.F4 10.13 0.180 52 9.76 10.5
Control Left F5.F6 9.83 0.208 52 9.42 10.3
low-SES Left F5.F6 10.17 0.197 52 9.77 10.6
Control Right F5.F6 9.83 0.198 52 9.43 10.2
low-SES Right F5.F6 10.06 0.187 52 9.69 10.4
Control Left F7.F8 9.90 0.198 52 9.51 10.3
low-SES Left F7.F8 10.18 0.187 52 9.81 10.6
Control Right F7.F8 9.90 0.195 52 9.51 10.3
low-SES Right F7.F8 10.14 0.184 52 9.77 10.5
Control Left FC1.FC2 9.69 0.199 52 9.29 10.1
low-SES Left FC1.FC2 10.05 0.188 52 9.67 10.4
Control Right FC1.FC2 9.69 0.191 52 9.30 10.1
low-SES Right FC1.FC2 10.13 0.181 52 9.77 10.5
Results are averaged over the levels of: Sex
Confidence level used: 0.95
Sex hemisphere electrode_pair emmean SE df lower.CL upper.CL
female Left Fp1.Fp2 9.87 0.197 52 9.47 10.3
male Left Fp1.Fp2 10.13 0.186 52 9.76 10.5
female Right Fp1.Fp2 9.96 0.193 52 9.57 10.3
male Right Fp1.Fp2 10.10 0.182 52 9.73 10.5
female Left AF3.AF4 9.87 0.199 52 9.48 10.3
male Left AF3.AF4 10.05 0.188 52 9.67 10.4
female Right AF3.AF4 9.90 0.204 52 9.49 10.3
male Right AF3.AF4 10.08 0.193 52 9.69 10.5
female Left AF7.AF8 9.90 0.200 52 9.50 10.3
male Left AF7.AF8 10.11 0.189 52 9.73 10.5
female Right AF7.AF8 9.92 0.193 52 9.54 10.3
male Right AF7.AF8 10.10 0.183 52 9.74 10.5
female Left F1.F2 9.73 0.198 52 9.33 10.1
male Left F1.F2 10.07 0.187 52 9.69 10.4
female Right F1.F2 9.83 0.191 52 9.44 10.2
male Right F1.F2 10.05 0.181 52 9.69 10.4
female Left F3.F4 9.74 0.202 52 9.33 10.1
male Left F3.F4 10.02 0.191 52 9.64 10.4
female Right F3.F4 9.81 0.191 52 9.43 10.2
male Right F3.F4 10.06 0.180 52 9.70 10.4
female Left F5.F6 9.92 0.208 52 9.50 10.3
male Left F5.F6 10.09 0.197 52 9.69 10.5
female Right F5.F6 9.82 0.198 52 9.42 10.2
male Right F5.F6 10.08 0.187 52 9.70 10.5
female Left F7.F8 9.95 0.198 52 9.56 10.3
male Left F7.F8 10.13 0.187 52 9.76 10.5
female Right F7.F8 9.89 0.195 52 9.50 10.3
male Right F7.F8 10.14 0.184 52 9.78 10.5
female Left FC1.FC2 9.72 0.199 52 9.32 10.1
male Left FC1.FC2 10.01 0.188 52 9.64 10.4
female Right FC1.FC2 9.78 0.191 52 9.39 10.2
male Right FC1.FC2 10.04 0.181 52 9.67 10.4
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex hemisphere electrode_pair emmean SE df lower.CL upper.CL
Control female Left Fp1.Fp2 9.54 0.312 52 8.91 10.2
low-SES female Left Fp1.Fp2 10.19 0.239 52 9.71 10.7
Control male Left Fp1.Fp2 10.12 0.239 52 9.64 10.6
low-SES male Left Fp1.Fp2 10.14 0.285 52 9.57 10.7
Control female Right Fp1.Fp2 9.70 0.306 52 9.08 10.3
low-SES female Right Fp1.Fp2 10.23 0.234 52 9.75 10.7
Control male Right Fp1.Fp2 9.99 0.234 52 9.52 10.5
low-SES male Right Fp1.Fp2 10.21 0.279 52 9.65 10.8
Control female Left AF3.AF4 9.60 0.315 52 8.96 10.2
low-SES female Left AF3.AF4 10.15 0.242 52 9.67 10.6
Control male Left AF3.AF4 9.97 0.242 52 9.48 10.5
low-SES male Left AF3.AF4 10.13 0.288 52 9.55 10.7
Control female Right AF3.AF4 9.67 0.323 52 9.02 10.3
low-SES female Right AF3.AF4 10.14 0.248 52 9.64 10.6
Control male Right AF3.AF4 9.99 0.248 52 9.49 10.5
low-SES male Right AF3.AF4 10.18 0.295 52 9.58 10.8
Control female Left AF7.AF8 9.60 0.317 52 8.96 10.2
low-SES female Left AF7.AF8 10.20 0.243 52 9.71 10.7
Control male Left AF7.AF8 10.04 0.243 52 9.56 10.5
low-SES male Left AF7.AF8 10.17 0.289 52 9.59 10.8
Control female Right AF7.AF8 9.61 0.306 52 9.00 10.2
low-SES female Right AF7.AF8 10.24 0.235 52 9.77 10.7
Control male Right AF7.AF8 10.00 0.235 52 9.53 10.5
low-SES male Right AF7.AF8 10.20 0.280 52 9.64 10.8
Control female Left F1.F2 9.49 0.314 52 8.86 10.1
low-SES female Left F1.F2 9.97 0.241 52 9.49 10.5
Control male Left F1.F2 9.99 0.241 52 9.51 10.5
low-SES male Left F1.F2 10.14 0.287 52 9.56 10.7
Control female Right F1.F2 9.50 0.303 52 8.89 10.1
low-SES female Right F1.F2 10.15 0.232 52 9.69 10.6
Control male Right F1.F2 9.93 0.232 52 9.47 10.4
low-SES male Right F1.F2 10.17 0.277 52 9.62 10.7
Control female Left F3.F4 9.46 0.321 52 8.81 10.1
low-SES female Left F3.F4 10.01 0.246 52 9.52 10.5
Control male Left F3.F4 9.93 0.246 52 9.43 10.4
low-SES male Left F3.F4 10.11 0.293 52 9.53 10.7
Control female Right F3.F4 9.52 0.302 52 8.91 10.1
low-SES female Right F3.F4 10.10 0.232 52 9.64 10.6
Control male Right F3.F4 9.97 0.232 52 9.50 10.4
low-SES male Right F3.F4 10.15 0.276 52 9.60 10.7
Control female Left F5.F6 9.64 0.330 52 8.97 10.3
low-SES female Left F5.F6 10.20 0.253 52 9.69 10.7
Control male Left F5.F6 10.03 0.253 52 9.53 10.5
low-SES male Left F5.F6 10.14 0.301 52 9.54 10.7
Control female Right F5.F6 9.68 0.314 52 9.05 10.3
low-SES female Right F5.F6 9.95 0.241 52 9.47 10.4
Control male Right F5.F6 9.98 0.241 52 9.50 10.5
low-SES male Right F5.F6 10.18 0.287 52 9.60 10.8
Control female Left F7.F8 9.69 0.313 52 9.06 10.3
low-SES female Left F7.F8 10.22 0.240 52 9.73 10.7
Control male Left F7.F8 10.12 0.240 52 9.63 10.6
low-SES male Left F7.F8 10.15 0.286 52 9.58 10.7
Control female Right F7.F8 9.74 0.309 52 9.12 10.4
low-SES female Right F7.F8 10.05 0.237 52 9.57 10.5
Control male Right F7.F8 10.06 0.237 52 9.59 10.5
low-SES male Right F7.F8 10.23 0.282 52 9.66 10.8
Control female Left FC1.FC2 9.47 0.316 52 8.83 10.1
low-SES female Left FC1.FC2 9.97 0.242 52 9.49 10.5
Control male Left FC1.FC2 9.90 0.242 52 9.42 10.4
low-SES male Left FC1.FC2 10.12 0.288 52 9.54 10.7
Control female Right FC1.FC2 9.49 0.304 52 8.88 10.1
low-SES female Right FC1.FC2 10.07 0.233 52 9.60 10.5
Control male Right FC1.FC2 9.88 0.233 52 9.41 10.3
low-SES male Right FC1.FC2 10.20 0.277 52 9.64 10.8
Confidence level used: 0.95 5.3 Alpha Bandwidth
options(width = 100)
alpha_bwd_rep_anova = aov_ez("Subject", "Width", alpha_power_data, between = c("Group", "Sex"), within = c("hemisphere", "electrode_pair"))G2;H2;Warningh: Missing values for 21 ID(s), which were removed before analysis:
1, 4, 7, 16, 17, 18, 20, 28, 34, 42, ... [showing first 10 only]
Below the first few rows (in wide format) of the removed cases with missing data.
g
G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = alpha_bwd_rep_anova$data$long) / length(levels(alpha_bwd_rep_anova$data$long$electrode_pair)) / length(levels(alpha_bwd_rep_anova$data$long$hemisphere))
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 10 17 27
low-SES 17 12 29
Sum 27 29 56alpha_bwd_afex_plot <-
afex_plot(
alpha_bwd_rep_anova,
x = "Group",
trace = "Sex",
panel = "hemisphere",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(alpha_bwd_afex_plot))
Anova Table (Type 3 tests)
Response: Width
Effect df MSE F ges p.value
1 Group 1, 52 6.04 0.01 <.001 .921
2 Sex 1, 52 6.04 0.47 .006 .496
3 Group:Sex 1, 52 6.04 0.48 .006 .493
4 hemisphere 1, 52 0.33 1.45 .001 .234
5 Group:hemisphere 1, 52 0.33 0.25 <.001 .620
6 Sex:hemisphere 1, 52 0.33 0.07 <.001 .786
7 Group:Sex:hemisphere 1, 52 0.33 0.11 <.001 .737
8 electrode_pair 2.74, 142.29 0.75 22.31 *** .088 <.001
9 Group:electrode_pair 2.74, 142.29 0.75 0.87 .004 .451
10 Sex:electrode_pair 2.74, 142.29 0.75 1.36 .006 .257
11 Group:Sex:electrode_pair 2.74, 142.29 0.75 1.01 .004 .386
12 hemisphere:electrode_pair 4.56, 236.88 0.15 0.94 .001 .450
13 Group:hemisphere:electrode_pair 4.56, 236.88 0.15 0.87 .001 .496
14 Sex:hemisphere:electrode_pair 4.56, 236.88 0.15 1.03 .002 .399
15 Group:Sex:hemisphere:electrode_pair 4.56, 236.88 0.15 0.95 .001 .441
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1
Sphericity correction method: GG G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control 2.25 0.122 52 2.01 2.5
low-SES 2.27 0.116 52 2.04 2.5
Results are averaged over the levels of: Sex, electrode_pair, hemisphere
Confidence level used: 0.95 G3;NOTE: Results may be misleading due to involvement in interactions
g Sex emmean SE df lower.CL upper.CL
female 2.32 0.122 52 2.07 2.56
male 2.20 0.116 52 1.97 2.44
Results are averaged over the levels of: Group, electrode_pair, hemisphere
Confidence level used: 0.95
Group Sex emmean SE df lower.CL upper.CL
Control female 2.25 0.194 52 1.86 2.64
low-SES female 2.39 0.149 52 2.09 2.68
Control male 2.25 0.149 52 1.95 2.55
low-SES male 2.15 0.177 52 1.80 2.51
Results are averaged over the levels of: electrode_pair, hemisphere
Confidence level used: 0.95
hemisphere emmean SE df lower.CL upper.CL
Left 2.24 0.0862 52 2.06 2.41
Right 2.28 0.0869 52 2.11 2.46
Results are averaged over the levels of: Group, Sex, electrode_pair
Confidence level used: 0.95
Group hemisphere emmean SE df lower.CL upper.CL
Control Left 2.22 0.125 52 1.97 2.47
low-SES Left 2.26 0.118 52 2.02 2.49
Control Right 2.29 0.126 52 2.03 2.54
low-SES Right 2.28 0.119 52 2.04 2.52
Results are averaged over the levels of: Sex, electrode_pair
Confidence level used: 0.95
Sex hemisphere emmean SE df lower.CL upper.CL
female Left 2.30 0.125 52 2.05 2.55
male Left 2.17 0.118 52 1.94 2.41
female Right 2.34 0.126 52 2.08 2.59
male Right 2.23 0.119 52 1.99 2.47
Results are averaged over the levels of: Group, electrode_pair
Confidence level used: 0.95
Group Sex hemisphere emmean SE df lower.CL upper.CL
Control female Left 2.22 0.199 52 1.82 2.62
low-SES female Left 2.38 0.152 52 2.08 2.69
Control male Left 2.22 0.152 52 1.91 2.53
low-SES male Left 2.13 0.181 52 1.76 2.49
Control female Right 2.29 0.200 52 1.88 2.69
low-SES female Right 2.39 0.154 52 2.08 2.70
Control male Right 2.29 0.154 52 1.98 2.59
low-SES male Right 2.18 0.183 52 1.81 2.55
Results are averaged over the levels of: electrode_pair
Confidence level used: 0.95
________________________________________________________________________________
$emmeans
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 2.04 0.0791 52 1.88 2.20
AF3.AF4 2.18 0.0935 52 1.99 2.37
AF7.AF8 2.04 0.0666 52 1.91 2.18
F1.F2 2.58 0.1170 52 2.35 2.82
F3.F4 2.38 0.1190 52 2.14 2.62
F5.F6 2.15 0.0839 52 1.98 2.32
F7.F8 2.06 0.0731 52 1.91 2.21
FC1.FC2 2.65 0.1280 52 2.40 2.91
Results are averaged over the levels of: Group, Sex, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 -0.13921 0.0628 52 -2.216 0.3594
Fp1.Fp2 - AF7.AF8 -0.00426 0.0385 52 -0.111 1.0000
Fp1.Fp2 - F1.F2 -0.54291 0.0863 52 -6.289 <.0001
Fp1.Fp2 - F3.F4 -0.34192 0.0756 52 -4.521 0.0009
Fp1.Fp2 - F5.F6 -0.11099 0.0555 52 -1.999 0.4927
Fp1.Fp2 - F7.F8 -0.02021 0.0583 52 -0.347 1.0000
Fp1.Fp2 - FC1.FC2 -0.61625 0.0928 52 -6.644 <.0001
AF3.AF4 - AF7.AF8 0.13495 0.0707 52 1.910 0.5505
AF3.AF4 - F1.F2 -0.40370 0.0714 52 -5.657 <.0001
AF3.AF4 - F3.F4 -0.20270 0.0572 52 -3.544 0.0177
AF3.AF4 - F5.F6 0.02823 0.0732 52 0.386 0.9999
AF3.AF4 - F7.F8 0.11900 0.0779 52 1.527 0.7895
AF3.AF4 - FC1.FC2 -0.47703 0.0766 52 -6.224 <.0001
AF7.AF8 - F1.F2 -0.53865 0.0955 52 -5.641 <.0001
AF7.AF8 - F3.F4 -0.33766 0.0878 52 -3.848 0.0074
AF7.AF8 - F5.F6 -0.10673 0.0588 52 -1.814 0.6138
AF7.AF8 - F7.F8 -0.01595 0.0555 52 -0.288 1.0000
AF7.AF8 - FC1.FC2 -0.61199 0.1040 52 -5.880 <.0001
F1.F2 - F3.F4 0.20100 0.0605 52 3.321 0.0327
F1.F2 - F5.F6 0.43193 0.0849 52 5.088 0.0001
F1.F2 - F7.F8 0.52270 0.0926 52 5.644 <.0001
F1.F2 - FC1.FC2 -0.07333 0.0380 52 -1.928 0.5386
F3.F4 - F5.F6 0.23093 0.0719 52 3.213 0.0434
F3.F4 - F7.F8 0.32171 0.0852 52 3.775 0.0091
F3.F4 - FC1.FC2 -0.27433 0.0553 52 -4.957 0.0002
F5.F6 - F7.F8 0.09078 0.0386 52 2.353 0.2861
F5.F6 - FC1.FC2 -0.50526 0.0866 52 -5.832 <.0001
F7.F8 - FC1.FC2 -0.59604 0.1000 52 -5.948 <.0001
Results are averaged over the levels of: Group, Sex, hemisphere
P value adjustment: tukey method for comparing a family of 8 estimates
Group electrode_pair emmean SE df lower.CL upper.CL
Control Fp1.Fp2 2.04 0.1150 52 1.81 2.27
low-SES Fp1.Fp2 2.03 0.1090 52 1.82 2.25
Control AF3.AF4 2.26 0.1360 52 1.99 2.54
low-SES AF3.AF4 2.09 0.1290 52 1.83 2.35
Control AF7.AF8 2.05 0.0968 52 1.85 2.24
low-SES AF7.AF8 2.04 0.0915 52 1.86 2.22
Control F1.F2 2.60 0.1710 52 2.26 2.94
low-SES F1.F2 2.56 0.1610 52 2.24 2.89
Control F3.F4 2.36 0.1720 52 2.02 2.71
low-SES F3.F4 2.40 0.1630 52 2.07 2.73
Control F5.F6 2.08 0.1220 52 1.83 2.32
low-SES F5.F6 2.22 0.1150 52 1.99 2.46
Control F7.F8 2.02 0.1060 52 1.80 2.23
low-SES F7.F8 2.10 0.1000 52 1.90 2.30
Control FC1.FC2 2.61 0.1860 52 2.24 2.98
low-SES FC1.FC2 2.70 0.1760 52 2.35 3.05
Results are averaged over the levels of: Sex, hemisphere
Confidence level used: 0.95
Sex electrode_pair emmean SE df lower.CL upper.CL
female Fp1.Fp2 2.05 0.1150 52 1.82 2.28
male Fp1.Fp2 2.03 0.1090 52 1.81 2.25
female AF3.AF4 2.23 0.1360 52 1.95 2.50
male AF3.AF4 2.13 0.1290 52 1.87 2.39
female AF7.AF8 2.07 0.0968 52 1.87 2.26
male AF7.AF8 2.02 0.0915 52 1.83 2.20
female F1.F2 2.68 0.1710 52 2.34 3.02
male F1.F2 2.48 0.1610 52 2.16 2.81
female F3.F4 2.50 0.1720 52 2.15 2.84
male F3.F4 2.26 0.1630 52 1.94 2.59
female F5.F6 2.18 0.1220 52 1.94 2.43
male F5.F6 2.12 0.1150 52 1.88 2.35
female F7.F8 2.03 0.1060 52 1.82 2.25
male F7.F8 2.09 0.1000 52 1.88 2.29
female FC1.FC2 2.81 0.1860 52 2.44 3.18
male FC1.FC2 2.50 0.1760 52 2.15 2.85
Results are averaged over the levels of: Group, hemisphere
Confidence level used: 0.95
Group Sex electrode_pair emmean SE df lower.CL upper.CL
Control female Fp1.Fp2 2.02 0.182 52 1.65 2.38
low-SES female Fp1.Fp2 2.08 0.140 52 1.80 2.36
Control male Fp1.Fp2 2.07 0.140 52 1.79 2.35
low-SES male Fp1.Fp2 1.99 0.167 52 1.66 2.32
Control female AF3.AF4 2.35 0.216 52 1.92 2.78
low-SES female AF3.AF4 2.10 0.165 52 1.77 2.44
Control male AF3.AF4 2.18 0.165 52 1.85 2.51
low-SES male AF3.AF4 2.08 0.197 52 1.68 2.47
Control female AF7.AF8 2.00 0.154 52 1.69 2.31
low-SES female AF7.AF8 2.14 0.118 52 1.90 2.37
Control male AF7.AF8 2.10 0.118 52 1.86 2.33
low-SES male AF7.AF8 1.94 0.140 52 1.66 2.22
Control female F1.F2 2.62 0.271 52 2.08 3.16
low-SES female F1.F2 2.74 0.208 52 2.32 3.16
Control male F1.F2 2.58 0.208 52 2.16 2.99
low-SES male F1.F2 2.39 0.247 52 1.89 2.89
Control female F3.F4 2.45 0.273 52 1.90 3.00
low-SES female F3.F4 2.54 0.210 52 2.12 2.97
Control male F3.F4 2.27 0.210 52 1.85 2.69
low-SES male F3.F4 2.25 0.250 52 1.75 2.76
Control female F5.F6 1.96 0.194 52 1.57 2.35
low-SES female F5.F6 2.40 0.148 52 2.11 2.70
Control male F5.F6 2.19 0.148 52 1.89 2.49
low-SES male F5.F6 2.04 0.177 52 1.69 2.40
Control female F7.F8 1.93 0.169 52 1.59 2.27
low-SES female F7.F8 2.13 0.129 52 1.87 2.39
Control male F7.F8 2.10 0.129 52 1.84 2.36
low-SES male F7.F8 2.07 0.154 52 1.76 2.38
Control female FC1.FC2 2.68 0.295 52 2.09 3.28
low-SES female FC1.FC2 2.94 0.226 52 2.48 3.39
Control male FC1.FC2 2.54 0.226 52 2.09 2.99
low-SES male FC1.FC2 2.46 0.269 52 1.92 3.00
Results are averaged over the levels of: hemisphere
Confidence level used: 0.95
hemisphere electrode_pair emmean SE df lower.CL upper.CL
Left Fp1.Fp2 2.02 0.0871 52 1.84 2.19
Right Fp1.Fp2 2.06 0.0826 52 1.89 2.23
Left AF3.AF4 2.19 0.1000 52 1.99 2.39
Right AF3.AF4 2.16 0.0937 52 1.98 2.35
Left AF7.AF8 2.02 0.0774 52 1.86 2.17
Right AF7.AF8 2.07 0.0748 52 1.92 2.22
Left F1.F2 2.59 0.1240 52 2.34 2.84
Right F1.F2 2.57 0.1240 52 2.32 2.82
Left F3.F4 2.36 0.1180 52 2.12 2.60
Right F3.F4 2.40 0.1250 52 2.15 2.65
Left F5.F6 2.07 0.0873 52 1.90 2.25
Right F5.F6 2.23 0.1020 52 2.02 2.43
Left F7.F8 2.01 0.0847 52 1.84 2.18
Right F7.F8 2.11 0.0827 52 1.94 2.27
Left FC1.FC2 2.64 0.1330 52 2.37 2.90
Right FC1.FC2 2.67 0.1300 52 2.41 2.93
Results are averaged over the levels of: Group, Sex
Confidence level used: 0.95
Group hemisphere electrode_pair emmean SE df lower.CL upper.CL
Control Left Fp1.Fp2 1.99 0.126 52 1.74 2.24
low-SES Left Fp1.Fp2 2.04 0.120 52 1.80 2.28
Control Right Fp1.Fp2 2.10 0.120 52 1.86 2.34
low-SES Right Fp1.Fp2 2.02 0.113 52 1.80 2.25
Control Left AF3.AF4 2.27 0.146 52 1.98 2.57
low-SES Left AF3.AF4 2.11 0.138 52 1.83 2.39
Control Right AF3.AF4 2.25 0.136 52 1.98 2.53
low-SES Right AF3.AF4 2.07 0.129 52 1.82 2.33
Control Left AF7.AF8 1.98 0.113 52 1.76 2.21
low-SES Left AF7.AF8 2.05 0.106 52 1.83 2.26
Control Right AF7.AF8 2.11 0.109 52 1.89 2.33
low-SES Right AF7.AF8 2.03 0.103 52 1.82 2.24
Control Left F1.F2 2.61 0.180 52 2.25 2.97
low-SES Left F1.F2 2.57 0.170 52 2.23 2.91
Control Right F1.F2 2.58 0.180 52 2.22 2.94
low-SES Right F1.F2 2.56 0.170 52 2.22 2.90
Control Left F3.F4 2.30 0.172 52 1.95 2.64
low-SES Left F3.F4 2.42 0.163 52 2.09 2.75
Control Right F3.F4 2.43 0.181 52 2.06 2.79
low-SES Right F3.F4 2.38 0.171 52 2.04 2.72
Control Left F5.F6 1.97 0.127 52 1.72 2.23
low-SES Left F5.F6 2.17 0.120 52 1.93 2.41
Control Right F5.F6 2.18 0.148 52 1.88 2.48
low-SES Right F5.F6 2.28 0.140 52 2.00 2.56
Control Left F7.F8 1.99 0.123 52 1.74 2.23
low-SES Left F7.F8 2.04 0.116 52 1.81 2.27
Control Right F7.F8 2.05 0.120 52 1.81 2.29
low-SES Right F7.F8 2.16 0.114 52 1.94 2.39
Control Left FC1.FC2 2.63 0.193 52 2.25 3.02
low-SES Left FC1.FC2 2.64 0.182 52 2.28 3.01
Control Right FC1.FC2 2.59 0.189 52 2.21 2.97
low-SES Right FC1.FC2 2.75 0.179 52 2.40 3.11
Results are averaged over the levels of: Sex
Confidence level used: 0.95
Sex hemisphere electrode_pair emmean SE df lower.CL upper.CL
female Left Fp1.Fp2 2.02 0.126 52 1.77 2.28
male Left Fp1.Fp2 2.01 0.120 52 1.77 2.25
female Right Fp1.Fp2 2.07 0.120 52 1.83 2.31
male Right Fp1.Fp2 2.05 0.113 52 1.82 2.28
female Left AF3.AF4 2.24 0.146 52 1.95 2.53
male Left AF3.AF4 2.14 0.138 52 1.87 2.42
female Right AF3.AF4 2.21 0.136 52 1.94 2.49
male Right AF3.AF4 2.11 0.129 52 1.86 2.37
female Left AF7.AF8 2.07 0.113 52 1.84 2.30
male Left AF7.AF8 1.96 0.106 52 1.75 2.17
female Right AF7.AF8 2.07 0.109 52 1.85 2.28
male Right AF7.AF8 2.07 0.103 52 1.87 2.28
female Left F1.F2 2.66 0.180 52 2.30 3.02
male Left F1.F2 2.52 0.170 52 2.18 2.87
female Right F1.F2 2.70 0.180 52 2.34 3.06
male Right F1.F2 2.44 0.170 52 2.10 2.78
female Left F3.F4 2.46 0.172 52 2.12 2.81
male Left F3.F4 2.25 0.163 52 1.93 2.58
female Right F3.F4 2.53 0.181 52 2.17 2.90
male Right F3.F4 2.27 0.171 52 1.93 2.62
female Left F5.F6 2.15 0.127 52 1.90 2.41
male Left F5.F6 1.99 0.120 52 1.75 2.23
female Right F5.F6 2.22 0.148 52 1.92 2.51
male Right F5.F6 2.24 0.140 52 1.96 2.52
female Left F7.F8 2.03 0.123 52 1.78 2.28
male Left F7.F8 2.00 0.116 52 1.76 2.23
female Right F7.F8 2.03 0.120 52 1.79 2.28
male Right F7.F8 2.18 0.114 52 1.95 2.40
female Left FC1.FC2 2.77 0.193 52 2.38 3.15
male Left FC1.FC2 2.51 0.182 52 2.14 2.88
female Right FC1.FC2 2.86 0.189 52 2.48 3.24
male Right FC1.FC2 2.49 0.179 52 2.13 2.85
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex hemisphere electrode_pair emmean SE df lower.CL upper.CL
Control female Left Fp1.Fp2 1.96 0.201 52 1.55 2.36
low-SES female Left Fp1.Fp2 2.09 0.154 52 1.78 2.40
Control male Left Fp1.Fp2 2.02 0.154 52 1.71 2.33
low-SES male Left Fp1.Fp2 2.00 0.183 52 1.63 2.37
Control female Right Fp1.Fp2 2.07 0.190 52 1.69 2.46
low-SES female Right Fp1.Fp2 2.07 0.146 52 1.77 2.36
Control male Right Fp1.Fp2 2.12 0.146 52 1.82 2.41
low-SES male Right Fp1.Fp2 1.98 0.174 52 1.63 2.33
Control female Left AF3.AF4 2.31 0.232 52 1.84 2.77
low-SES female Left AF3.AF4 2.17 0.178 52 1.81 2.53
Control male Left AF3.AF4 2.24 0.178 52 1.88 2.59
low-SES male Left AF3.AF4 2.05 0.212 52 1.63 2.47
Control female Right AF3.AF4 2.39 0.216 52 1.96 2.82
low-SES female Right AF3.AF4 2.04 0.166 52 1.71 2.37
Control male Right AF3.AF4 2.12 0.166 52 1.79 2.45
low-SES male Right AF3.AF4 2.11 0.197 52 1.71 2.50
Control female Left AF7.AF8 2.00 0.179 52 1.64 2.35
low-SES female Left AF7.AF8 2.14 0.137 52 1.87 2.42
Control male Left AF7.AF8 1.97 0.137 52 1.70 2.25
low-SES male Left AF7.AF8 1.95 0.163 52 1.62 2.28
Control female Right AF7.AF8 2.00 0.172 52 1.65 2.35
low-SES female Right AF7.AF8 2.13 0.132 52 1.87 2.40
Control male Right AF7.AF8 2.22 0.132 52 1.95 2.48
low-SES male Right AF7.AF8 1.93 0.157 52 1.61 2.24
Control female Left F1.F2 2.58 0.285 52 2.01 3.15
low-SES female Left F1.F2 2.73 0.219 52 2.29 3.17
Control male Left F1.F2 2.65 0.219 52 2.21 3.08
low-SES male Left F1.F2 2.40 0.260 52 1.88 2.93
Control female Right F1.F2 2.66 0.286 52 2.09 3.23
low-SES female Right F1.F2 2.75 0.219 52 2.31 3.18
Control male Right F1.F2 2.51 0.219 52 2.07 2.95
low-SES male Right F1.F2 2.37 0.261 52 1.85 2.90
Control female Left F3.F4 2.35 0.273 52 1.80 2.89
low-SES female Left F3.F4 2.58 0.209 52 2.16 3.00
Control male Left F3.F4 2.25 0.209 52 1.83 2.67
low-SES male Left F3.F4 2.26 0.249 52 1.76 2.76
Control female Right F3.F4 2.56 0.287 52 1.98 3.14
low-SES female Right F3.F4 2.51 0.220 52 2.07 2.95
Control male Right F3.F4 2.29 0.220 52 1.85 2.74
low-SES male Right F3.F4 2.25 0.262 52 1.72 2.78
Control female Left F5.F6 1.94 0.201 52 1.54 2.35
low-SES female Left F5.F6 2.36 0.154 52 2.05 2.67
Control male Left F5.F6 2.00 0.154 52 1.69 2.31
low-SES male Left F5.F6 1.98 0.184 52 1.61 2.35
Control female Right F5.F6 1.98 0.235 52 1.51 2.45
low-SES female Right F5.F6 2.45 0.180 52 2.09 2.81
Control male Right F5.F6 2.38 0.180 52 2.02 2.74
low-SES male Right F5.F6 2.11 0.214 52 1.68 2.54
Control female Left F7.F8 1.95 0.195 52 1.55 2.34
low-SES female Left F7.F8 2.11 0.150 52 1.81 2.41
Control male Left F7.F8 2.03 0.150 52 1.72 2.33
low-SES male Left F7.F8 1.97 0.178 52 1.61 2.32
Control female Right F7.F8 1.92 0.191 52 1.53 2.30
low-SES female Right F7.F8 2.15 0.146 52 1.86 2.45
Control male Right F7.F8 2.18 0.146 52 1.88 2.47
low-SES male Right F7.F8 2.18 0.174 52 1.83 2.52
Control female Left FC1.FC2 2.65 0.306 52 2.04 3.27
low-SES female Left FC1.FC2 2.88 0.234 52 2.41 3.35
Control male Left FC1.FC2 2.61 0.234 52 2.14 3.08
low-SES male Left FC1.FC2 2.41 0.279 52 1.85 2.97
Control female Right FC1.FC2 2.71 0.300 52 2.11 3.31
low-SES female Right FC1.FC2 3.00 0.230 52 2.54 3.46
Control male Right FC1.FC2 2.47 0.230 52 2.01 2.93
low-SES male Right FC1.FC2 2.51 0.274 52 1.96 3.06
Confidence level used: 0.95 5.4 Aperiodic Offset
options(width = 100)
aperiodic_offset_rep_anova = aov_ez("Subject", "Offset", alpha_power_data, between = c("Group", "Sex"), within = c("hemisphere", "electrode_pair"))G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = aperiodic_offset_rep_anova$data$long) / length(levels(aperiodic_offset_rep_anova$data$long$electrode_pair)) / length(levels(aperiodic_offset_rep_anova$data$long$hemisphere))
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 19 20 39
low-SES 24 14 38
Sum 43 34 77aperiodic_param_afex_plot <-
afex_plot(
aperiodic_offset_rep_anova,
x = "Group",
trace = "Sex",
panel = "hemisphere",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(aperiodic_param_afex_plot))
Anova Table (Type 3 tests)
Response: Offset
Effect df MSE F ges p.value
1 Group 1, 73 2.35 2.64 .025 .109
2 Sex 1, 73 2.35 5.40 * .050 .023
3 Group:Sex 1, 73 2.35 1.08 .010 .301
4 hemisphere 1, 73 0.09 0.23 <.001 .632
5 Group:hemisphere 1, 73 0.09 1.42 <.001 .238
6 Sex:hemisphere 1, 73 0.09 0.15 <.001 .703
7 Group:Sex:hemisphere 1, 73 0.09 0.00 <.001 .981
8 electrode_pair 2.37, 172.65 0.29 33.95 *** .087 <.001
9 Group:electrode_pair 2.37, 172.65 0.29 0.25 <.001 .816
10 Sex:electrode_pair 2.37, 172.65 0.29 2.71 + .007 .060
11 Group:Sex:electrode_pair 2.37, 172.65 0.29 0.43 .001 .688
12 hemisphere:electrode_pair 4.75, 346.40 0.04 0.77 <.001 .568
13 Group:hemisphere:electrode_pair 4.75, 346.40 0.04 1.02 <.001 .406
14 Sex:hemisphere:electrode_pair 4.75, 346.40 0.04 1.82 .001 .112
15 Group:Sex:hemisphere:electrode_pair 4.75, 346.40 0.04 1.40 .001 .225
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1
Sphericity correction method: GG G3;NOTE: Results may be misleading due to involvement in interactions
g Group emmean SE df lower.CL upper.CL
Control 1.091 0.0614 73 0.969 1.21
low-SES 0.946 0.0644 73 0.818 1.07
Results are averaged over the levels of: Sex, electrode_pair, hemisphere
Confidence level used: 0.95
________________________________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Sex emmean SE df lower.CL upper.CL
female 0.915 0.0588 73 0.798 1.03
male 1.122 0.0668 73 0.989 1.25
Results are averaged over the levels of: Group, electrode_pair, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
female - male -0.207 0.089 73 -2.323 0.0230
Results are averaged over the levels of: Group, electrode_pair, hemisphere
Group Sex emmean SE df lower.CL upper.CL
Control female 1.034 0.0879 73 0.859 1.209
low-SES female 0.797 0.0782 73 0.641 0.953
Control male 1.148 0.0857 73 0.977 1.319
low-SES male 1.096 0.1020 73 0.892 1.300
Results are averaged over the levels of: electrode_pair, hemisphere
Confidence level used: 0.95
hemisphere emmean SE df lower.CL upper.CL
Left 1.02 0.0465 73 0.930 1.12
Right 1.01 0.0442 73 0.926 1.10
Results are averaged over the levels of: Group, Sex, electrode_pair
Confidence level used: 0.95
Group hemisphere emmean SE df lower.CL upper.CL
Control Left 1.106 0.0642 73 0.978 1.23
low-SES Left 0.940 0.0673 73 0.806 1.07
Control Right 1.076 0.0610 73 0.955 1.20
low-SES Right 0.953 0.0640 73 0.825 1.08
Results are averaged over the levels of: Sex, electrode_pair
Confidence level used: 0.95
Sex hemisphere emmean SE df lower.CL upper.CL
female Left 0.916 0.0615 73 0.794 1.04
male Left 1.130 0.0698 73 0.991 1.27
female Right 0.914 0.0584 73 0.798 1.03
male Right 1.114 0.0663 73 0.982 1.25
Results are averaged over the levels of: Group, electrode_pair
Confidence level used: 0.95
Group Sex hemisphere emmean SE df lower.CL upper.CL
Control female Left 1.045 0.0919 73 0.862 1.229
low-SES female Left 0.787 0.0818 73 0.624 0.950
Control male Left 1.166 0.0896 73 0.987 1.344
low-SES male Left 1.093 0.1070 73 0.880 1.307
Control female Right 1.022 0.0873 73 0.848 1.196
low-SES female Right 0.807 0.0777 73 0.652 0.961
Control male Right 1.130 0.0851 73 0.960 1.300
low-SES male Right 1.099 0.1020 73 0.896 1.301
Results are averaged over the levels of: electrode_pair
Confidence level used: 0.95
________________________________________________________________________________
$emmeans
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 1.212 0.0625 73 1.087 1.337
AF3.AF4 1.127 0.0581 73 1.012 1.243
AF7.AF8 1.113 0.0526 73 1.008 1.217
F1.F2 1.065 0.0432 73 0.979 1.151
F3.F4 0.845 0.0469 73 0.751 0.938
F5.F6 0.775 0.0454 73 0.684 0.865
F7.F8 0.960 0.0456 73 0.870 1.051
FC1.FC2 1.052 0.0460 73 0.961 1.144
Results are averaged over the levels of: Group, Sex, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 0.0847 0.0176 73 4.814 0.0002
Fp1.Fp2 - AF7.AF8 0.0995 0.0315 73 3.156 0.0455
Fp1.Fp2 - F1.F2 0.1474 0.0471 73 3.128 0.0490
Fp1.Fp2 - F3.F4 0.3672 0.0499 73 7.355 <.0001
Fp1.Fp2 - F5.F6 0.4373 0.0463 73 9.450 <.0001
Fp1.Fp2 - F7.F8 0.2518 0.0419 73 6.004 <.0001
Fp1.Fp2 - FC1.FC2 0.1598 0.0548 73 2.915 0.0844
AF3.AF4 - AF7.AF8 0.0148 0.0303 73 0.488 0.9997
AF3.AF4 - F1.F2 0.0627 0.0395 73 1.588 0.7558
AF3.AF4 - F3.F4 0.2825 0.0413 73 6.836 <.0001
AF3.AF4 - F5.F6 0.3526 0.0390 73 9.031 <.0001
AF3.AF4 - F7.F8 0.1671 0.0377 73 4.435 0.0008
AF3.AF4 - FC1.FC2 0.0751 0.0473 73 1.590 0.7547
AF7.AF8 - F1.F2 0.0479 0.0405 73 1.183 0.9344
AF7.AF8 - F3.F4 0.2677 0.0402 73 6.661 <.0001
AF7.AF8 - F5.F6 0.3378 0.0321 73 10.534 <.0001
AF7.AF8 - F7.F8 0.1523 0.0258 73 5.911 <.0001
AF7.AF8 - FC1.FC2 0.0603 0.0467 73 1.291 0.8993
F1.F2 - F3.F4 0.2199 0.0182 73 12.083 <.0001
F1.F2 - F5.F6 0.2899 0.0253 73 11.477 <.0001
F1.F2 - F7.F8 0.1044 0.0316 73 3.301 0.0306
F1.F2 - FC1.FC2 0.0125 0.0146 73 0.851 0.9893
F3.F4 - F5.F6 0.0701 0.0203 73 3.459 0.0194
F3.F4 - F7.F8 -0.1155 0.0321 73 -3.595 0.0130
F3.F4 - FC1.FC2 -0.2074 0.0238 73 -8.700 <.0001
F5.F6 - F7.F8 -0.1855 0.0223 73 -8.317 <.0001
F5.F6 - FC1.FC2 -0.2775 0.0305 73 -9.093 <.0001
F7.F8 - FC1.FC2 -0.0919 0.0367 73 -2.507 0.2089
Results are averaged over the levels of: Group, Sex, hemisphere
P value adjustment: tukey method for comparing a family of 8 estimates
Group electrode_pair emmean SE df lower.CL upper.CL
Control Fp1.Fp2 1.296 0.0863 73 1.124 1.468
low-SES Fp1.Fp2 1.128 0.0906 73 0.947 1.308
Control AF3.AF4 1.205 0.0801 73 1.045 1.365
low-SES AF3.AF4 1.050 0.0841 73 0.882 1.218
Control AF7.AF8 1.193 0.0726 73 1.048 1.338
low-SES AF7.AF8 1.032 0.0762 73 0.880 1.184
Control F1.F2 1.135 0.0596 73 1.016 1.254
low-SES F1.F2 0.995 0.0626 73 0.870 1.119
Control F3.F4 0.916 0.0647 73 0.787 1.045
low-SES F3.F4 0.774 0.0680 73 0.638 0.909
Control F5.F6 0.826 0.0627 73 0.701 0.951
low-SES F5.F6 0.723 0.0658 73 0.592 0.854
Control F7.F8 1.047 0.0628 73 0.922 1.172
low-SES F7.F8 0.874 0.0660 73 0.742 1.005
Control FC1.FC2 1.109 0.0635 73 0.982 1.235
low-SES FC1.FC2 0.996 0.0666 73 0.863 1.129
Results are averaged over the levels of: Sex, hemisphere
Confidence level used: 0.95
Sex electrode_pair emmean SE df lower.CL upper.CL
female Fp1.Fp2 1.125 0.0827 73 0.960 1.290
male Fp1.Fp2 1.299 0.0938 73 1.112 1.486
female AF3.AF4 1.031 0.0768 73 0.878 1.184
male AF3.AF4 1.224 0.0871 73 1.050 1.397
female AF7.AF8 1.027 0.0696 73 0.888 1.166
male AF7.AF8 1.198 0.0790 73 1.041 1.355
female F1.F2 0.918 0.0571 73 0.804 1.032
male F1.F2 1.212 0.0648 73 1.082 1.341
female F3.F4 0.721 0.0621 73 0.597 0.845
male F3.F4 0.969 0.0704 73 0.828 1.109
female F5.F6 0.698 0.0601 73 0.578 0.818
male F5.F6 0.852 0.0682 73 0.716 0.988
female F7.F8 0.921 0.0602 73 0.801 1.041
male F7.F8 0.999 0.0684 73 0.863 1.136
female FC1.FC2 0.882 0.0608 73 0.760 1.003
male FC1.FC2 1.223 0.0690 73 1.085 1.360
Results are averaged over the levels of: Group, hemisphere
Confidence level used: 0.95
Group Sex electrode_pair emmean SE df lower.CL upper.CL
Control female Fp1.Fp2 1.221 0.1240 73 0.974 1.467
low-SES female Fp1.Fp2 1.029 0.1100 73 0.810 1.248
Control male Fp1.Fp2 1.372 0.1200 73 1.132 1.612
low-SES male Fp1.Fp2 1.227 0.1440 73 0.940 1.514
Control female AF3.AF4 1.156 0.1150 73 0.927 1.385
low-SES female AF3.AF4 0.906 0.1020 73 0.702 1.109
Control male AF3.AF4 1.254 0.1120 73 1.031 1.476
low-SES male AF3.AF4 1.194 0.1340 73 0.928 1.460
Control female AF7.AF8 1.164 0.1040 73 0.957 1.371
low-SES female AF7.AF8 0.890 0.0925 73 0.706 1.074
Control male AF7.AF8 1.222 0.1010 73 1.020 1.424
low-SES male AF7.AF8 1.174 0.1210 73 0.933 1.415
Control female F1.F2 1.027 0.0854 73 0.856 1.197
low-SES female F1.F2 0.809 0.0759 73 0.658 0.960
Control male F1.F2 1.243 0.0832 73 1.077 1.409
low-SES male F1.F2 1.180 0.0994 73 0.982 1.378
Control female F3.F4 0.844 0.0927 73 0.659 1.029
low-SES female F3.F4 0.598 0.0825 73 0.433 0.762
Control male F3.F4 0.988 0.0904 73 0.808 1.168
low-SES male F3.F4 0.950 0.1080 73 0.735 1.165
Control female F5.F6 0.789 0.0898 73 0.610 0.968
low-SES female F5.F6 0.606 0.0799 73 0.447 0.765
Control male F5.F6 0.863 0.0875 73 0.689 1.038
low-SES male F5.F6 0.840 0.1050 73 0.632 1.049
Control female F7.F8 1.068 0.0900 73 0.889 1.248
low-SES female F7.F8 0.774 0.0801 73 0.615 0.934
Control male F7.F8 1.026 0.0877 73 0.851 1.200
low-SES male F7.F8 0.973 0.1050 73 0.764 1.182
Control female FC1.FC2 1.001 0.0909 73 0.820 1.183
low-SES female FC1.FC2 0.762 0.0809 73 0.601 0.923
Control male FC1.FC2 1.216 0.0886 73 1.039 1.392
low-SES male FC1.FC2 1.230 0.1060 73 1.019 1.441
Results are averaged over the levels of: hemisphere
Confidence level used: 0.95
hemisphere electrode_pair emmean SE df lower.CL upper.CL
Left Fp1.Fp2 1.225 0.0626 73 1.100 1.349
Right Fp1.Fp2 1.199 0.0670 73 1.066 1.333
Left AF3.AF4 1.124 0.0607 73 1.003 1.245
Right AF3.AF4 1.130 0.0609 73 1.009 1.252
Left AF7.AF8 1.127 0.0555 73 1.016 1.238
Right AF7.AF8 1.098 0.0574 73 0.984 1.213
Left F1.F2 1.052 0.0461 73 0.960 1.144
Right F1.F2 1.077 0.0419 73 0.994 1.161
Left F3.F4 0.844 0.0512 73 0.742 0.946
Right F3.F4 0.846 0.0462 73 0.754 0.938
Left F5.F6 0.799 0.0495 73 0.701 0.898
Right F5.F6 0.750 0.0454 73 0.660 0.841
Left F7.F8 0.960 0.0492 73 0.862 1.058
Right F7.F8 0.960 0.0485 73 0.864 1.057
Left FC1.FC2 1.051 0.0483 73 0.955 1.148
Right FC1.FC2 1.053 0.0470 73 0.960 1.147
Results are averaged over the levels of: Group, Sex
Confidence level used: 0.95
Group hemisphere electrode_pair emmean SE df lower.CL upper.CL
Control Left Fp1.Fp2 1.310 0.0863 73 1.138 1.482
low-SES Left Fp1.Fp2 1.140 0.0906 73 0.959 1.320
Control Right Fp1.Fp2 1.283 0.0924 73 1.099 1.467
low-SES Right Fp1.Fp2 1.116 0.0970 73 0.923 1.310
Control Left AF3.AF4 1.213 0.0838 73 1.046 1.380
low-SES Left AF3.AF4 1.036 0.0879 73 0.861 1.211
Control Right AF3.AF4 1.197 0.0839 73 1.029 1.364
low-SES Right AF3.AF4 1.064 0.0881 73 0.889 1.240
Control Left AF7.AF8 1.248 0.0766 73 1.095 1.400
low-SES Left AF7.AF8 1.006 0.0804 73 0.846 1.166
Control Right AF7.AF8 1.139 0.0791 73 0.981 1.296
low-SES Right AF7.AF8 1.058 0.0831 73 0.892 1.223
Control Left F1.F2 1.131 0.0635 73 1.005 1.258
low-SES Left F1.F2 0.973 0.0667 73 0.840 1.106
Control Right F1.F2 1.139 0.0577 73 1.024 1.254
low-SES Right F1.F2 1.016 0.0606 73 0.895 1.137
Control Left F3.F4 0.933 0.0706 73 0.792 1.074
low-SES Left F3.F4 0.755 0.0741 73 0.607 0.902
Control Right F3.F4 0.899 0.0637 73 0.772 1.026
low-SES Right F3.F4 0.793 0.0669 73 0.660 0.926
Control Left F5.F6 0.853 0.0683 73 0.717 0.990
low-SES Left F5.F6 0.745 0.0717 73 0.602 0.888
Control Right F5.F6 0.799 0.0626 73 0.675 0.924
low-SES Right F5.F6 0.701 0.0658 73 0.570 0.832
Control Left F7.F8 1.054 0.0678 73 0.919 1.189
low-SES Left F7.F8 0.867 0.0712 73 0.725 1.009
Control Right F7.F8 1.040 0.0669 73 0.907 1.173
low-SES Right F7.F8 0.880 0.0702 73 0.741 1.020
Control Left FC1.FC2 1.104 0.0666 73 0.971 1.237
low-SES Left FC1.FC2 0.999 0.0699 73 0.860 1.138
Control Right FC1.FC2 1.113 0.0648 73 0.984 1.242
low-SES Right FC1.FC2 0.993 0.0680 73 0.857 1.128
Results are averaged over the levels of: Sex
Confidence level used: 0.95
Sex hemisphere electrode_pair emmean SE df lower.CL upper.CL
female Left Fp1.Fp2 1.119 0.0827 73 0.954 1.283
male Left Fp1.Fp2 1.331 0.0939 73 1.144 1.518
female Right Fp1.Fp2 1.131 0.0886 73 0.954 1.308
male Right Fp1.Fp2 1.268 0.1010 73 1.067 1.468
female Left AF3.AF4 1.010 0.0803 73 0.850 1.170
male Left AF3.AF4 1.239 0.0911 73 1.058 1.421
female Right AF3.AF4 1.052 0.0805 73 0.892 1.213
male Right AF3.AF4 1.208 0.0913 73 1.026 1.390
female Left AF7.AF8 1.065 0.0734 73 0.919 1.211
male Left AF7.AF8 1.189 0.0833 73 1.023 1.355
female Right AF7.AF8 0.989 0.0758 73 0.838 1.140
male Right AF7.AF8 1.207 0.0861 73 1.036 1.379
female Left F1.F2 0.901 0.0609 73 0.780 1.023
male Left F1.F2 1.203 0.0691 73 1.065 1.341
female Right F1.F2 0.934 0.0553 73 0.824 1.045
male Right F1.F2 1.220 0.0628 73 1.095 1.345
female Left F3.F4 0.716 0.0677 73 0.581 0.851
male Left F3.F4 0.971 0.0768 73 0.818 1.125
female Right F3.F4 0.726 0.0611 73 0.604 0.847
male Right F3.F4 0.966 0.0693 73 0.828 1.104
female Left F5.F6 0.719 0.0655 73 0.588 0.849
male Left F5.F6 0.880 0.0743 73 0.732 1.028
female Right F5.F6 0.677 0.0600 73 0.557 0.797
male Right F5.F6 0.823 0.0681 73 0.688 0.959
female Left F7.F8 0.945 0.0650 73 0.815 1.074
male Left F7.F8 0.976 0.0738 73 0.829 1.123
female Right F7.F8 0.898 0.0641 73 0.770 1.026
male Right F7.F8 1.022 0.0727 73 0.877 1.167
female Left FC1.FC2 0.856 0.0638 73 0.729 0.983
male Left FC1.FC2 1.247 0.0725 73 1.102 1.391
female Right FC1.FC2 0.907 0.0621 73 0.784 1.031
male Right FC1.FC2 1.199 0.0705 73 1.058 1.339
Results are averaged over the levels of: Group
Confidence level used: 0.95
Group Sex hemisphere electrode_pair emmean SE df lower.CL upper.CL
Control female Left Fp1.Fp2 1.210 0.1240 73 0.964 1.457
low-SES female Left Fp1.Fp2 1.027 0.1100 73 0.808 1.246
Control male Left Fp1.Fp2 1.409 0.1200 73 1.169 1.649
low-SES male Left Fp1.Fp2 1.253 0.1440 73 0.966 1.540
Control female Right Fp1.Fp2 1.231 0.1320 73 0.967 1.495
low-SES female Right Fp1.Fp2 1.031 0.1180 73 0.796 1.266
Control male Right Fp1.Fp2 1.335 0.1290 73 1.077 1.592
low-SES male Right Fp1.Fp2 1.201 0.1540 73 0.894 1.508
Control female Left AF3.AF4 1.158 0.1200 73 0.919 1.397
low-SES female Left AF3.AF4 0.861 0.1070 73 0.649 1.074
Control male Left AF3.AF4 1.268 0.1170 73 1.035 1.501
low-SES male Left AF3.AF4 1.210 0.1400 73 0.932 1.489
Control female Right AF3.AF4 1.155 0.1200 73 0.915 1.394
low-SES female Right AF3.AF4 0.950 0.1070 73 0.737 1.164
Control male Right AF3.AF4 1.239 0.1170 73 1.005 1.472
low-SES male Right AF3.AF4 1.178 0.1400 73 0.899 1.457
Control female Left AF7.AF8 1.247 0.1100 73 1.028 1.466
low-SES female Left AF7.AF8 0.883 0.0976 73 0.688 1.077
Control male Left AF7.AF8 1.248 0.1070 73 1.035 1.461
low-SES male Left AF7.AF8 1.130 0.1280 73 0.875 1.384
Control female Right AF7.AF8 1.081 0.1130 73 0.855 1.307
low-SES female Right AF7.AF8 0.897 0.1010 73 0.696 1.098
Control male Right AF7.AF8 1.196 0.1100 73 0.976 1.416
low-SES male Right AF7.AF8 1.218 0.1320 73 0.955 1.481
Control female Left F1.F2 1.026 0.0910 73 0.845 1.208
low-SES female Left F1.F2 0.776 0.0810 73 0.615 0.937
Control male Left F1.F2 1.236 0.0887 73 1.059 1.413
low-SES male Left F1.F2 1.170 0.1060 73 0.959 1.382
Control female Right F1.F2 1.027 0.0827 73 0.862 1.192
low-SES female Right F1.F2 0.842 0.0736 73 0.695 0.988
Control male Right F1.F2 1.250 0.0806 73 1.090 1.411
low-SES male Right F1.F2 1.190 0.0963 73 0.998 1.382
Control female Left F3.F4 0.856 0.1010 73 0.655 1.058
low-SES female Left F3.F4 0.576 0.0900 73 0.397 0.755
Control male Left F3.F4 1.010 0.0985 73 0.813 1.206
low-SES male Left F3.F4 0.933 0.1180 73 0.699 1.168
Control female Right F3.F4 0.832 0.0912 73 0.650 1.014
low-SES female Right F3.F4 0.619 0.0812 73 0.457 0.781
Control male Right F3.F4 0.966 0.0889 73 0.788 1.143
low-SES male Right F3.F4 0.967 0.1060 73 0.755 1.178
Control female Left F5.F6 0.818 0.0978 73 0.623 1.013
low-SES female Left F5.F6 0.619 0.0871 73 0.446 0.793
Control male Left F5.F6 0.889 0.0954 73 0.699 1.079
low-SES male Left F5.F6 0.871 0.1140 73 0.644 1.099
Control female Right F5.F6 0.761 0.0897 73 0.582 0.940
low-SES female Right F5.F6 0.593 0.0798 73 0.434 0.752
Control male Right F5.F6 0.838 0.0874 73 0.664 1.012
low-SES male Right F5.F6 0.809 0.1050 73 0.601 1.017
Control female Left F7.F8 1.062 0.0971 73 0.868 1.256
low-SES female Left F7.F8 0.827 0.0864 73 0.655 0.999
Control male Left F7.F8 1.046 0.0947 73 0.857 1.234
low-SES male Left F7.F8 0.907 0.1130 73 0.681 1.132
Control female Right F7.F8 1.075 0.0958 73 0.884 1.266
low-SES female Right F7.F8 0.721 0.0852 73 0.551 0.891
Control male Right F7.F8 1.005 0.0934 73 0.819 1.191
low-SES male Right F7.F8 1.040 0.1120 73 0.817 1.262
Control female Left FC1.FC2 0.986 0.0954 73 0.796 1.177
low-SES female Left FC1.FC2 0.725 0.0849 73 0.556 0.895
Control male Left FC1.FC2 1.221 0.0930 73 1.036 1.407
low-SES male Left FC1.FC2 1.272 0.1110 73 1.051 1.494
Control female Right FC1.FC2 1.016 0.0928 73 0.832 1.201
low-SES female Right FC1.FC2 0.798 0.0825 73 0.634 0.963
Control male Right FC1.FC2 1.210 0.0904 73 1.030 1.390
low-SES male Right FC1.FC2 1.187 0.1080 73 0.972 1.403
Confidence level used: 0.95 5.5 Aperiodic Slope
options(width = 100)
aperiodic_slope_rep_anova = aov_ez("Subject", "Slope", alpha_power_data, between = c("Group", "Sex"), within = c("hemisphere", "electrode_pair"))G3;Contrasts set to contr.sum for the following variables: Group, Sex
gmytable <- xtabs(~ Group + Sex, data = aperiodic_slope_rep_anova$data$long) / length(levels(aperiodic_slope_rep_anova$data$long$electrode_pair)) / length(levels(aperiodic_slope_rep_anova$data$long$hemisphere))
ftable(addmargins(mytable)) Sex female male Sum
Group
Control 19 20 39
low-SES 24 14 38
Sum 43 34 77aperiodic_param_afex_plot <-
afex_plot(
aperiodic_slope_rep_anova,
x = "Group",
trace = "Sex",
panel = "hemisphere",
error = "between",
error_arg = list(width = .1),
dodge = .3,
data_arg = list(
position =
position_jitterdodge(
jitter.width = .1,
jitter.height = .1,
dodge.width = .3 ## needs to be same as dodge
)),
mapping = c("color"),
point_arg = list(size = 4)
)
suppressWarnings(print(aperiodic_param_afex_plot))
Anova Table (Type 3 tests)
Response: Slope
Effect df MSE F ges p.value
1 Group 1, 73 1.70 5.43 * .054 .023
2 Sex 1, 73 1.70 9.19 ** .088 .003
3 Group:Sex 1, 73 1.70 0.32 .003 .574
4 hemisphere 1, 73 0.05 2.75 <.001 .102
5 Group:hemisphere 1, 73 0.05 0.00 <.001 .968
6 Sex:hemisphere 1, 73 0.05 0.37 <.001 .543
7 Group:Sex:hemisphere 1, 73 0.05 0.19 <.001 .664
8 electrode_pair 3.36, 245.05 0.11 56.12 *** .110 <.001
9 Group:electrode_pair 3.36, 245.05 0.11 0.23 <.001 .892
10 Sex:electrode_pair 3.36, 245.05 0.11 2.55 * .006 .049
11 Group:Sex:electrode_pair 3.36, 245.05 0.11 0.87 .002 .467
12 hemisphere:electrode_pair 4.92, 359.10 0.02 2.48 * .002 .033
13 Group:hemisphere:electrode_pair 4.92, 359.10 0.02 0.60 <.001 .698
14 Sex:hemisphere:electrode_pair 4.92, 359.10 0.02 1.00 <.001 .419
15 Group:Sex:hemisphere:electrode_pair 4.92, 359.10 0.02 0.81 <.001 .543
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1
Sphericity correction method: GG G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Group emmean SE df lower.CL upper.CL
Control 1.18 0.0522 73 1.077 1.28
low-SES 1.00 0.0547 73 0.896 1.11
Results are averaged over the levels of: Sex, electrode_pair, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Control - (low-SES) 0.176 0.0756 73 2.330 0.0226
Results are averaged over the levels of: Sex, electrode_pair, hemisphere ________________________________________________________________________________G3;NOTE: Results may be misleading due to involvement in interactions
g$emmeans
Sex emmean SE df lower.CL upper.CL
female 0.978 0.0500 73 0.879 1.08
male 1.207 0.0567 73 1.094 1.32
Results are averaged over the levels of: Group, electrode_pair, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
female - male -0.229 0.0756 73 -3.032 0.0034
Results are averaged over the levels of: Group, electrode_pair, hemisphere ________________________________________________________________________________$emmeans
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 1.112 0.0467 73 1.019 1.205
AF3.AF4 1.098 0.0475 73 1.004 1.193
AF7.AF8 1.020 0.0368 73 0.947 1.094
F1.F2 1.273 0.0395 73 1.194 1.352
F3.F4 1.024 0.0431 73 0.938 1.110
F5.F6 0.915 0.0404 73 0.834 0.996
F7.F8 0.987 0.0364 73 0.914 1.059
FC1.FC2 1.313 0.0405 73 1.232 1.393
Results are averaged over the levels of: Group, Sex, hemisphere
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 0.01347 0.0149 73 0.906 0.9846
Fp1.Fp2 - AF7.AF8 0.09137 0.0244 73 3.745 0.0082
Fp1.Fp2 - F1.F2 -0.16110 0.0301 73 -5.348 <.0001
Fp1.Fp2 - F3.F4 0.08735 0.0321 73 2.723 0.1324
Fp1.Fp2 - F5.F6 0.19673 0.0310 73 6.353 <.0001
Fp1.Fp2 - F7.F8 0.12508 0.0282 73 4.438 0.0008
Fp1.Fp2 - FC1.FC2 -0.20073 0.0359 73 -5.594 <.0001
AF3.AF4 - AF7.AF8 0.07790 0.0245 73 3.174 0.0434
AF3.AF4 - F1.F2 -0.17456 0.0292 73 -5.987 <.0001
AF3.AF4 - F3.F4 0.07389 0.0308 73 2.395 0.2593
AF3.AF4 - F5.F6 0.18326 0.0303 73 6.058 <.0001
AF3.AF4 - F7.F8 0.11162 0.0288 73 3.871 0.0055
AF3.AF4 - FC1.FC2 -0.21419 0.0343 73 -6.253 <.0001
AF7.AF8 - F1.F2 -0.25246 0.0274 73 -9.213 <.0001
AF7.AF8 - F3.F4 -0.00401 0.0269 73 -0.149 1.0000
AF7.AF8 - F5.F6 0.10536 0.0220 73 4.794 0.0002
AF7.AF8 - F7.F8 0.03372 0.0182 73 1.855 0.5857
AF7.AF8 - FC1.FC2 -0.29209 0.0305 73 -9.591 <.0001
F1.F2 - F3.F4 0.24845 0.0197 73 12.585 <.0001
F1.F2 - F5.F6 0.35782 0.0235 73 15.247 <.0001
F1.F2 - F7.F8 0.28618 0.0225 73 12.699 <.0001
F1.F2 - FC1.FC2 -0.03963 0.0108 73 -3.683 0.0099
F3.F4 - F5.F6 0.10937 0.0212 73 5.151 0.0001
F3.F4 - F7.F8 0.03773 0.0246 73 1.535 0.7862
F3.F4 - FC1.FC2 -0.28808 0.0236 73 -12.225 <.0001
F5.F6 - F7.F8 -0.07164 0.0168 73 -4.262 0.0015
F5.F6 - FC1.FC2 -0.39746 0.0255 73 -15.588 <.0001
F7.F8 - FC1.FC2 -0.32581 0.0254 73 -12.842 <.0001
Results are averaged over the levels of: Group, Sex, hemisphere
P value adjustment: tukey method for comparing a family of 8 estimates ________________________________________________________________________________$emmeans
Sex = female:
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 0.962 0.0617 73 0.839 1.085
AF3.AF4 0.936 0.0629 73 0.811 1.062
AF7.AF8 0.913 0.0486 73 0.816 1.010
F1.F2 1.154 0.0523 73 1.050 1.258
F3.F4 0.916 0.0570 73 0.803 1.030
F5.F6 0.834 0.0535 73 0.728 0.941
F7.F8 0.908 0.0482 73 0.812 1.004
FC1.FC2 1.201 0.0536 73 1.094 1.308
Sex = male:
electrode_pair emmean SE df lower.CL upper.CL
Fp1.Fp2 1.262 0.0700 73 1.122 1.402
AF3.AF4 1.260 0.0713 73 1.118 1.403
AF7.AF8 1.128 0.0552 73 1.018 1.238
F1.F2 1.392 0.0593 73 1.273 1.510
F3.F4 1.132 0.0647 73 1.003 1.261
F5.F6 0.996 0.0607 73 0.875 1.117
F7.F8 1.065 0.0547 73 0.956 1.174
FC1.FC2 1.424 0.0608 73 1.303 1.545
Results are averaged over the levels of: Group, hemisphere
Confidence level used: 0.95
$contrasts
Sex = female:
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 0.02540 0.0196 73 1.293 0.8987
Fp1.Fp2 - AF7.AF8 0.04855 0.0323 73 1.505 0.8022
Fp1.Fp2 - F1.F2 -0.19257 0.0398 73 -4.835 0.0002
Fp1.Fp2 - F3.F4 0.04520 0.0424 73 1.065 0.9619
Fp1.Fp2 - F5.F6 0.12741 0.0409 73 3.111 0.0512
Fp1.Fp2 - F7.F8 0.05335 0.0373 73 1.432 0.8396
Fp1.Fp2 - FC1.FC2 -0.23939 0.0474 73 -5.045 0.0001
AF3.AF4 - AF7.AF8 0.02315 0.0325 73 0.713 0.9963
AF3.AF4 - F1.F2 -0.21797 0.0386 73 -5.654 <.0001
AF3.AF4 - F3.F4 0.01980 0.0408 73 0.485 0.9997
AF3.AF4 - F5.F6 0.10200 0.0400 73 2.550 0.1917
AF3.AF4 - F7.F8 0.02795 0.0381 73 0.733 0.9957
AF3.AF4 - FC1.FC2 -0.26479 0.0453 73 -5.845 <.0001
AF7.AF8 - F1.F2 -0.24112 0.0362 73 -6.654 <.0001
AF7.AF8 - F3.F4 -0.00336 0.0355 73 -0.094 1.0000
AF7.AF8 - F5.F6 0.07885 0.0291 73 2.714 0.1351
AF7.AF8 - F7.F8 0.00480 0.0240 73 0.200 1.0000
AF7.AF8 - FC1.FC2 -0.28794 0.0403 73 -7.150 <.0001
F1.F2 - F3.F4 0.23777 0.0261 73 9.108 <.0001
F1.F2 - F5.F6 0.31997 0.0310 73 10.311 <.0001
F1.F2 - F7.F8 0.24592 0.0298 73 8.253 <.0001
F1.F2 - FC1.FC2 -0.04682 0.0142 73 -3.290 0.0315
F3.F4 - F5.F6 0.08221 0.0281 73 2.928 0.0817
F3.F4 - F7.F8 0.00816 0.0325 73 0.251 1.0000
F3.F4 - FC1.FC2 -0.28458 0.0312 73 -9.133 <.0001
F5.F6 - F7.F8 -0.07405 0.0222 73 -3.332 0.0281
F5.F6 - FC1.FC2 -0.36679 0.0337 73 -10.879 <.0001
F7.F8 - FC1.FC2 -0.29274 0.0335 73 -8.726 <.0001
Sex = male:
contrast estimate SE df t.ratio p.value
Fp1.Fp2 - AF3.AF4 0.00153 0.0223 73 0.069 1.0000
Fp1.Fp2 - AF7.AF8 0.13418 0.0366 73 3.665 0.0105
Fp1.Fp2 - F1.F2 -0.12962 0.0452 73 -2.868 0.0945
Fp1.Fp2 - F3.F4 0.12951 0.0481 73 2.690 0.1423
Fp1.Fp2 - F5.F6 0.26605 0.0465 73 5.726 <.0001
Fp1.Fp2 - F7.F8 0.19681 0.0423 73 4.654 0.0004
Fp1.Fp2 - FC1.FC2 -0.16207 0.0538 73 -3.010 0.0665
AF3.AF4 - AF7.AF8 0.13265 0.0368 73 3.602 0.0127
AF3.AF4 - F1.F2 -0.13115 0.0438 73 -2.998 0.0687
AF3.AF4 - F3.F4 0.12798 0.0463 73 2.765 0.1202
AF3.AF4 - F5.F6 0.26452 0.0454 73 5.827 <.0001
AF3.AF4 - F7.F8 0.19528 0.0433 73 4.514 0.0006
AF3.AF4 - FC1.FC2 -0.16360 0.0514 73 -3.183 0.0423
AF7.AF8 - F1.F2 -0.26380 0.0411 73 -6.416 <.0001
AF7.AF8 - F3.F4 -0.00467 0.0403 73 -0.116 1.0000
AF7.AF8 - F5.F6 0.13187 0.0330 73 3.999 0.0036
AF7.AF8 - F7.F8 0.06264 0.0273 73 2.296 0.3099
AF7.AF8 - FC1.FC2 -0.29625 0.0457 73 -6.482 <.0001
F1.F2 - F3.F4 0.25913 0.0296 73 8.748 <.0001
F1.F2 - F5.F6 0.39567 0.0352 73 11.236 <.0001
F1.F2 - F7.F8 0.32644 0.0338 73 9.654 <.0001
F1.F2 - FC1.FC2 -0.03245 0.0161 73 -2.010 0.4827
F3.F4 - F5.F6 0.13654 0.0319 73 4.286 0.0014
F3.F4 - F7.F8 0.06730 0.0369 73 1.824 0.6058
F3.F4 - FC1.FC2 -0.29158 0.0354 73 -8.246 <.0001
F5.F6 - F7.F8 -0.06924 0.0252 73 -2.745 0.1258
F5.F6 - FC1.FC2 -0.42812 0.0383 73 -11.190 <.0001
F7.F8 - FC1.FC2 -0.35889 0.0381 73 -9.427 <.0001
Results are averaged over the levels of: Group, hemisphere
P value adjustment: tukey method for comparing a family of 8 estimates ________________________________________________________________________________$emmeans
electrode_pair = Fp1.Fp2:
hemisphere emmean SE df lower.CL upper.CL
Left 1.139 0.0473 73 1.044 1.233
Right 1.085 0.0486 73 0.988 1.182
electrode_pair = AF3.AF4:
hemisphere emmean SE df lower.CL upper.CL
Left 1.106 0.0508 73 1.005 1.207
Right 1.091 0.0486 73 0.994 1.188
electrode_pair = AF7.AF8:
hemisphere emmean SE df lower.CL upper.CL
Left 1.038 0.0372 73 0.964 1.112
Right 1.003 0.0398 73 0.924 1.082
electrode_pair = F1.F2:
hemisphere emmean SE df lower.CL upper.CL
Left 1.258 0.0420 73 1.174 1.342
Right 1.288 0.0384 73 1.211 1.365
electrode_pair = F3.F4:
hemisphere emmean SE df lower.CL upper.CL
Left 1.022 0.0472 73 0.928 1.116
Right 1.027 0.0427 73 0.942 1.112
electrode_pair = F5.F6:
hemisphere emmean SE df lower.CL upper.CL
Left 0.947 0.0434 73 0.861 1.034
Right 0.883 0.0421 73 0.799 0.967
electrode_pair = F7.F8:
hemisphere emmean SE df lower.CL upper.CL
Left 1.006 0.0376 73 0.931 1.081
Right 0.967 0.0386 73 0.891 1.044
electrode_pair = FC1.FC2:
hemisphere emmean SE df lower.CL upper.CL
Left 1.309 0.0415 73 1.226 1.392
Right 1.316 0.0413 73 1.234 1.398
Results are averaged over the levels of: Group, Sex
Confidence level used: 0.95
$contrasts
electrode_pair = Fp1.Fp2:
contrast estimate SE df t.ratio p.value
Left - Right 0.05394 0.0222 73 2.426 0.0177
electrode_pair = AF3.AF4:
contrast estimate SE df t.ratio p.value
Left - Right 0.01468 0.0290 73 0.506 0.6145
electrode_pair = AF7.AF8:
contrast estimate SE df t.ratio p.value
Left - Right 0.03506 0.0229 73 1.532 0.1298
electrode_pair = F1.F2:
contrast estimate SE df t.ratio p.value
Left - Right -0.03011 0.0152 73 -1.979 0.0516
electrode_pair = F3.F4:
contrast estimate SE df t.ratio p.value
Left - Right -0.00484 0.0255 73 -0.190 0.8501
electrode_pair = F5.F6:
contrast estimate SE df t.ratio p.value
Left - Right 0.06459 0.0279 73 2.317 0.0233
electrode_pair = F7.F8:
contrast estimate SE df t.ratio p.value
Left - Right 0.03852 0.0223 73 1.727 0.0883
electrode_pair = FC1.FC2:
contrast estimate SE df t.ratio p.value
Left - Right -0.00697 0.0168 73 -0.415 0.6791
Results are averaged over the levels of: Group, Sex LS0tCnRpdGxlOiAiQWxwaGEgcG93ZXIgZXN0aW1hdGVkIGJ5IEZPT09GIgphdXRob3I6ICJBbHZhcm8gUml2ZXJhLVJlaSIKZGF0ZTogImByIGZvcm1hdChTeXMudGltZSgpLCAnJWQgJUIsICVZJylgIgpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vazoKICAgIGNvZGVfZm9sZGluZzogaGlkZQogICAgaGlnaGxpZ2h0OiB0YW5nbwogICAgbnVtYmVyX3NlY3Rpb25zOiB5ZXMKICAgIHRoZW1lOiBjZXJ1bGVhbgogICAgdG9jOiB5ZXMKICAgIHRvY19mbG9hdDoKICAgICAgY29sbGFwc2VkOiBubwogICAgICBzbW9vdGhfc2Nyb2xsOiBubwogIHBkZl9kb2N1bWVudDoKICAgIHRvYzogeWVzCnN1YnRpdGxlOiBGb25kZWN5dCAyMDE3Ci0tLQoKPCEtLSBgYGB7ciBzZXR1cCwgaW5jbHVkZSA9IEZBTFNFfSAgLS0+CjwhLS0ga25pdHI6Om9wdHNfY2h1bmskc2V0KHdhcm5pbmcgPSBGQUxTRSwgbWVzc2FnZSA9IEZBTFNFKSAgLS0+CjwhLS0gYGBgIC0tPgoKYGBge3IgQ2xlYW4gYW5kIExvYWQgTGlicmFyaWVzLCBtZXNzYWdlPUZBTFNFfQpjYXQoIlwwMTQiKSAgICAgIyBjbGVhbiB0ZXJtaW5hbApybShsaXN0ID0gbHMoKSkgIyBjbGVhbiB3b3Jrc3BhY2UKdHJ5KGRldi5vZmYoKSwgc2lsZW50ID0gVFJVRSkgIyBjbG9zZSBhbGwgcGxvdHMKbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkoYWZleCkKbGlicmFyeShsbWVyVGVzdCkKbGlicmFyeShlbW1lYW5zKQpsaWJyYXJ5KEdHYWxseSkKbGlicmFyeShlYXN5c3RhdHMpCmBgYAoKYGBge3IgU2V0IERlZmF1bHRzfQp0aGVtZV9zZXQoCiAgdGhlbWVfbWluaW1hbCgpCikKCmFfcG9zdGVyaW9yaSA8LSBmdW5jdGlvbihhZmV4X2Fvdiwgc2lnX2xldmVsID0gLjA1KSB7CiAgZmFjdG9ycyAgPC0gYXMubGlzdChyb3duYW1lcyhhZmV4X2FvdiRhbm92YV90YWJsZSkpCiAgZm9yIChqIGluIDE6bGVuZ3RoKGZhY3RvcnMpKSB7CiAgICBpZiAoZ3JlcGwoIjoiLCBmYWN0b3JzW1tqXV0pKSB7CiAgICAgIGZhY3RvcnNbW2pdXSA8LSB1bmxpc3Qoc3Ryc3BsaXQoZmFjdG9yc1tbal1dLCAiOiIpKQogICAgfQogIH0KICBwX3ZhbHVlcyA8LSBhZmV4X2FvdiRhbm92YV90YWJsZSRgUHIoPkYpYAogIGZvciAoaSBpbiAxOmxlbmd0aChwX3ZhbHVlcykpIHsKICAgIGlmIChwX3ZhbHVlc1tpXSA8PSBzaWdfbGV2ZWwpIHsKICAgICAgY2F0KHJlcCgiXyIsIDgwKSwgJ1xuJywgc2VwID0gIiIpCiAgICAgIHByaW50KGVtbWVhbnMoYWZleF9hb3YsIGZhY3RvcnNbW2ldXSwgY29udHIgPSAicGFpcndpc2UiKSkKICAgIH0gZWxzZSB7CiAgICAgIHByaW50KGVtbWVhbnMoYWZleF9hb3YsIGZhY3RvcnNbW2ldXSkpCiAgICB9CiAgfQp9CgphX3Bvc3RlcmlvcmlfbG1lciA8LSBmdW5jdGlvbihsbWVyX29iaiwgc2lnX2xldmVsID0gLjA1KSB7CiAgYW5vdmFfbG1lciA8LSBhbm92YShsbWVyX29iaikKICBmYWN0b3JzICA8LSBhcy5saXN0KHJvdy5uYW1lcyhhbm92YV9sbWVyKSkKICBmb3IgKGogaW4gMTpsZW5ndGgoZmFjdG9ycykpIHsKICAgIGlmIChncmVwbCgnOicsIGZhY3RvcnNbW2pdXSkpIHsKICAgICAgZmFjdG9yc1tbal1dIDwtIHVubGlzdChzdHJzcGxpdChmYWN0b3JzW1tqXV0sICc6JykpCiAgICB9CiAgfQogIHBfdmFsdWVzIDwtIGFub3ZhX2xtZXIkYFByKD5GKWAKICBmb3IgKGkgaW4gMTpsZW5ndGgocF92YWx1ZXMpKSB7CiAgICBpZiAocF92YWx1ZXNbaV0gPD0gc2lnX2xldmVsKSB7CiAgICAgIGNhdChyZXAoJ18nLCA2MCksICdcbicsIHNlcCA9ICcnKQogICAgICBwcmludChlbW1lYW5zKGxtZXJfb2JqLCBmYWN0b3JzW1tpXV0sIGNvbnRyID0gJ3BhaXJ3aXNlJykpCiAgICB9CiAgfQp9CmBgYAoKYGBge3IgTG9hZCBEYXRhfQojIG1hc3Rlcl9kaXIgICAgICAgICAgICAgICAgIDwtICcvaG9tZS9hbHZhcm8vSW5zeW5jL0RyaXZlLzAwTEFCV09SS1NfZHJpdmUvSHVlcGUvRmRjeXRfMjAxNy9yZXN0aW5nX2h1ZXBlJwojIGRhdGFfZGlyICAgICAgICAgICAgICAgICAgIDwtIHBhc3RlKG1hc3Rlcl9kaXIsICdGQUFfcmVzdWx0cycsICBzZXAgPSAnLycpCiMgYWxwaGFfcG93ZXJfZGF0YV9uYW1lICAgICAgPC0gcGFzdGUoZGF0YV9kaXIsICdmb29mX2RhdGFfMl90b180OF9Iei5jc3YnLCBzZXA9Jy8nKQphbHBoYV9wb3dlcl9kYXRhICAgICAgICAgICA8LSByZWFkLnRhYmxlKCdmb29mX2RhdGFfMl90b180OF9Iei5jc3YnLCBoZWFkZXIgPSBUUlVFLCBzdHJpcC53aGl0ZSA9IFRSVUUsIHNlcCA9ICIsIikKYWxwaGFfcG93ZXJfZGF0YSREYXRhc2V0ICAgPC0gZmFjdG9yKGFscGhhX3Bvd2VyX2RhdGEkRGF0YXNldCkKYWxwaGFfcG93ZXJfZGF0YSRFbGVjdHJvZGUgPC0gZmFjdG9yKGFscGhhX3Bvd2VyX2RhdGEkRWxlY3Ryb2RlKQphbHBoYV9wb3dlcl9kYXRhJFN1YmplY3QgICA8LSBmYWN0b3IoYWxwaGFfcG93ZXJfZGF0YSRTdWJqZWN0KQphbHBoYV9wb3dlcl9kYXRhJGhlbWlzcGhlcmVbYWxwaGFfcG93ZXJfZGF0YSRFbGVjdHJvZGUgJWluJSBjKCdFMDkzLUZwMScsICAnRTA5Mi1BRjNhJywgICdFMDk0LUFGNycsICAnRTA4OS1GMWEnLCAgJ0UxMDAtRjNhJywgICdFMTAxLUY1YScsICAnRTEwMy1GNycsICAnRTA4OC1GQzFhJyldIDwtICdMZWZ0JwphbHBoYV9wb3dlcl9kYXRhJGhlbWlzcGhlcmVbYWxwaGFfcG93ZXJfZGF0YSRFbGVjdHJvZGUgJWluJSBjKCdFMDgwLUZwMicsICAnRTA3OS1BRjRhJywgICdFMDcyLUFGOCcsICAnRTA3Ni1GMmEnLCAgJ0UwNjgtRjRhJywgICdFMDY5LUY2YScsICAnRTA3MS1GOCcsICAnRTA3NS1GQzJhJyldIDwtICdSaWdodCcKYWxwaGFfcG93ZXJfZGF0YSRoZW1pc3BoZXJlIDwtIGZhY3RvcihhbHBoYV9wb3dlcl9kYXRhJGhlbWlzcGhlcmUpCmFscGhhX3Bvd2VyX2RhdGEkZWxlY3Ryb2RlX3BhaXJbYWxwaGFfcG93ZXJfZGF0YSRFbGVjdHJvZGUgJWluJSBjKCdFMDkzLUZwMScgLCAnRTA4MC1GcDInKV0gIDwtICdGcDEtRnAyJwphbHBoYV9wb3dlcl9kYXRhJGVsZWN0cm9kZV9wYWlyW2FscGhhX3Bvd2VyX2RhdGEkRWxlY3Ryb2RlICVpbiUgYygnRTA5Mi1BRjNhJywgJ0UwNzktQUY0YScpXSA8LSAnQUYzLUFGNCcKYWxwaGFfcG93ZXJfZGF0YSRlbGVjdHJvZGVfcGFpclthbHBoYV9wb3dlcl9kYXRhJEVsZWN0cm9kZSAlaW4lIGMoJ0UwOTQtQUY3JyAsICdFMDcyLUFGOCcpXSAgPC0gJ0FGNy1BRjgnCmFscGhhX3Bvd2VyX2RhdGEkZWxlY3Ryb2RlX3BhaXJbYWxwaGFfcG93ZXJfZGF0YSRFbGVjdHJvZGUgJWluJSBjKCdFMDg5LUYxYScgLCAnRTA3Ni1GMmEnKV0gIDwtICdGMS1GMicKYWxwaGFfcG93ZXJfZGF0YSRlbGVjdHJvZGVfcGFpclthbHBoYV9wb3dlcl9kYXRhJEVsZWN0cm9kZSAlaW4lIGMoJ0UxMDAtRjNhJyAsICdFMDY4LUY0YScpXSAgPC0gJ0YzLUY0JwphbHBoYV9wb3dlcl9kYXRhJGVsZWN0cm9kZV9wYWlyW2FscGhhX3Bvd2VyX2RhdGEkRWxlY3Ryb2RlICVpbiUgYygnRTEwMS1GNWEnICwgJ0UwNjktRjZhJyldICA8LSAnRjUtRjYnCmFscGhhX3Bvd2VyX2RhdGEkZWxlY3Ryb2RlX3BhaXJbYWxwaGFfcG93ZXJfZGF0YSRFbGVjdHJvZGUgJWluJSBjKCdFMTAzLUY3JyAgLCAnRTA3MS1GOCcpXSAgIDwtICdGNy1GOCcKYWxwaGFfcG93ZXJfZGF0YSRlbGVjdHJvZGVfcGFpclthbHBoYV9wb3dlcl9kYXRhJEVsZWN0cm9kZSAlaW4lIGMoJ0UwODgtRkMxYScsICdFMDc1LUZDMmEnKV0gPC0gJ0ZDMS1GQzInCmFscGhhX3Bvd2VyX2RhdGEkZWxlY3Ryb2RlX3BhaXIgPC0gZmFjdG9yKGFscGhhX3Bvd2VyX2RhdGEkZWxlY3Ryb2RlX3BhaXIsIGxldmVscyA9IGMoJ0ZwMS1GcDInLCAnQUYzLUFGNCcsICdBRjctQUY4JywgJ0YxLUYyJywgJ0YzLUY0JywnRjUtRjYnLCAnRjctRjgnLCAnRkMxLUZDMicpKQpncm91cF9pZCAgIDwtIHJlYWQudGFibGUoIi9ob21lL2FsdmFyby9JbnN5bmMvRHJpdmUvMDBFRUcvUHJveWVjdG9zL0h1ZXBlL2ZkY3l0XzIwMTcvUmVnaXN0cm8tRXZhbHVhY2lvbmVzLUZEQ1lULURILTIwMTcgLSBHZW5lcmFsIEFOT05JTU8uY3N2IiwKICAgICAgICAgICAgICAgICAgICAgICAgIHNlcCA9ICIsIiwgaGVhZGVyID0gVFJVRSwgY29sLm5hbWVzID0gYygiZnVsbC5pZCIsICJTdWJqZWN0IiwgIlNleCIsICJHcm91cCIsICJTdHJlc3MiKSkKZ3JvdXBfaWQkU2V4ICAgICAgICAgICA8LSBmYWN0b3IoZ3JvdXBfaWQkU2V4KQpncm91cF9pZCRHcm91cCAgICAgICAgIDwtIGZhY3Rvcihncm91cF9pZCRHcm91cCkKbGV2ZWxzKGdyb3VwX2lkJFNleCkgICA8LSBsaXN0KGZlbWFsZSA9ICJGIiwgbWFsZSA9ICJNIikKbGV2ZWxzKGdyb3VwX2lkJEdyb3VwKSA8LSBsaXN0KENvbnRyb2wgPSAiQ04iLCAibG93LVNFUyIgPSAiRVgiKQpncm91cF9pZCAgICAgICAgICAgICAgIDwtIGdyb3VwX2lkW2MoJ1N1YmplY3QnLCAnR3JvdXAnLCAnU2V4JyldCmFscGhhX3Bvd2VyX2RhdGEgICAgICAgPC0gbWVyZ2UoYWxwaGFfcG93ZXJfZGF0YSwgZ3JvdXBfaWQsIGJ5ID0gJ1N1YmplY3QnKQp3cml0ZS5jc3YoYWxwaGFfcG93ZXJfZGF0YSwgICdhbHBoYV9wb3dlcl9kYXRhX2NsZWFuLmNzdicsICByb3cubmFtZXMgPSBGQUxTRSkKYXN5bW1ldHJ5X0ZwMl9GcDEgPC0gYygpCmFzeW1tZXRyeV9BRjRfQUYzIDwtIGMoKQphc3ltbWV0cnlfQUY4X0FGNyA8LSBjKCkKYXN5bW1ldHJ5X0YyX0YxICAgPC0gYygpCmFzeW1tZXRyeV9GNF9GMyAgIDwtIGMoKQphc3ltbWV0cnlfRjZfRjUgICA8LSBjKCkKYXN5bW1ldHJ5X0Y4X0Y3ICAgPC0gYygpCmFzeW1tZXRyeV9GQzJfRkMxIDwtIGMoKQpTdWJqZWN0IDwtIHVuaXF1ZShhbHBoYV9wb3dlcl9kYXRhJFN1YmplY3QpCmZvciAoc3ViaiBpbiBTdWJqZWN0KSB7CiAgc3ViamVjdF9kYXRhIDwtIHN1YnNldChhbHBoYV9wb3dlcl9kYXRhLCBTdWJqZWN0ID09IHN1YmopCiAgYXN5bW1ldHJ5X0ZwMl9GcDEgPC0gYyhhc3ltbWV0cnlfRnAyX0ZwMSwgc3ViamVjdF9kYXRhW3doaWNoKHN1YmplY3RfZGF0YSRFbGVjdHJvZGUgPT0gJ0UwODAtRnAyJykgLCA1XSAtIHN1YmplY3RfZGF0YVt3aGljaChzdWJqZWN0X2RhdGEkRWxlY3Ryb2RlPT0nRTA5My1GcDEnKSAsIDVdKQogIGFzeW1tZXRyeV9BRjRfQUYzIDwtIGMoYXN5bW1ldHJ5X0FGNF9BRjMsIHN1YmplY3RfZGF0YVt3aGljaChzdWJqZWN0X2RhdGEkRWxlY3Ryb2RlID09ICdFMDc5LUFGNGEnKSwgNV0gLSBzdWJqZWN0X2RhdGFbd2hpY2goc3ViamVjdF9kYXRhJEVsZWN0cm9kZT09J0UwOTItQUYzYScpLCA1XSkKICBhc3ltbWV0cnlfQUY4X0FGNyA8LSBjKGFzeW1tZXRyeV9BRjhfQUY3LCBzdWJqZWN0X2RhdGFbd2hpY2goc3ViamVjdF9kYXRhJEVsZWN0cm9kZSA9PSAnRTA3Mi1BRjgnKSAsIDVdIC0gc3ViamVjdF9kYXRhW3doaWNoKHN1YmplY3RfZGF0YSRFbGVjdHJvZGU9PSdFMDk0LUFGNycpICwgNV0pCiAgYXN5bW1ldHJ5X0YyX0YxICAgPC0gYyhhc3ltbWV0cnlfRjJfRjEgICwgc3ViamVjdF9kYXRhW3doaWNoKHN1YmplY3RfZGF0YSRFbGVjdHJvZGUgPT0gJ0UwNzYtRjJhJykgLCA1XSAtIHN1YmplY3RfZGF0YVt3aGljaChzdWJqZWN0X2RhdGEkRWxlY3Ryb2RlPT0nRTA4OS1GMWEnKSAsIDVdKQogIGFzeW1tZXRyeV9GNF9GMyAgIDwtIGMoYXN5bW1ldHJ5X0Y0X0YzICAsIHN1YmplY3RfZGF0YVt3aGljaChzdWJqZWN0X2RhdGEkRWxlY3Ryb2RlID09ICdFMDY4LUY0YScpICwgNV0gLSBzdWJqZWN0X2RhdGFbd2hpY2goc3ViamVjdF9kYXRhJEVsZWN0cm9kZT09J0UxMDAtRjNhJykgLCA1XSkKICBhc3ltbWV0cnlfRjZfRjUgICA8LSBjKGFzeW1tZXRyeV9GNl9GNSAgLCBzdWJqZWN0X2RhdGFbd2hpY2goc3ViamVjdF9kYXRhJEVsZWN0cm9kZSA9PSAnRTA2OS1GNmEnKSAsIDVdIC0gc3ViamVjdF9kYXRhW3doaWNoKHN1YmplY3RfZGF0YSRFbGVjdHJvZGU9PSdFMTAxLUY1YScpICwgNV0pCiAgYXN5bW1ldHJ5X0Y4X0Y3ICAgPC0gYyhhc3ltbWV0cnlfRjhfRjcgICwgc3ViamVjdF9kYXRhW3doaWNoKHN1YmplY3RfZGF0YSRFbGVjdHJvZGUgPT0gJ0UwNzEtRjgnKSAgLCA1XSAtIHN1YmplY3RfZGF0YVt3aGljaChzdWJqZWN0X2RhdGEkRWxlY3Ryb2RlPT0nRTEwMy1GNycpICAsIDVdKQogIGFzeW1tZXRyeV9GQzJfRkMxIDwtIGMoYXN5bW1ldHJ5X0ZDMl9GQzEsIHN1YmplY3RfZGF0YVt3aGljaChzdWJqZWN0X2RhdGEkRWxlY3Ryb2RlID09ICdFMDc1LUZDMmEnKSwgNV0gLSBzdWJqZWN0X2RhdGFbd2hpY2goc3ViamVjdF9kYXRhJEVsZWN0cm9kZT09J0UwODgtRkMxYScpLCA1XSkKfQphbHBoYV9hc3ltbWV0cnlfZGF0YSA8LSBkYXRhLmZyYW1lKFN1YmplY3QsIGFzeW1tZXRyeV9GcDJfRnAxLCBhc3ltbWV0cnlfQUY0X0FGMywgYXN5bW1ldHJ5X0FGOF9BRjcsIGFzeW1tZXRyeV9GMl9GMSwgYXN5bW1ldHJ5X0Y0X0YzLCBhc3ltbWV0cnlfRjZfRjUsIGFzeW1tZXRyeV9GOF9GNywgYXN5bW1ldHJ5X0ZDMl9GQzEpCmFscGhhX2FzeW1tZXRyeV9kYXRhIDwtIG1lcmdlKGFscGhhX2FzeW1tZXRyeV9kYXRhLCBncm91cF9pZCwgYnkgPSAnU3ViamVjdCcpCmFscGhhX2FzeW1tZXRyeV9kYXRhX2xvbmcgPC0gcGl2b3RfbG9uZ2VyKGFscGhhX2FzeW1tZXRyeV9kYXRhLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2xzICAgICAgICAgPSBhc3ltbWV0cnlfRnAyX0ZwMTphc3ltbWV0cnlfRkMyX0ZDMSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5hbWVzX3RvICAgICA9ICdFbGVjdHJvZGVzJyAsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBuYW1lc19wcmVmaXggPSAnYXN5bW1ldHJ5XycsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBuYW1lc19wdHlwZXMgPSBmYWN0b3IoKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgdmFsdWVzX3RvICAgID0gJ0FzaW1tZXRyeScpCndyaXRlLmNzdihhbHBoYV9hc3ltbWV0cnlfZGF0YSAgICAgLCAnYWxwaGFfYXN5bW1ldHJ5X2RhdGFfY2xlYW4uY3N2JyAgICAgLCByb3cubmFtZXMgPSBGQUxTRSkKd3JpdGUuY3N2KGFscGhhX2FzeW1tZXRyeV9kYXRhX2xvbmcsICdhbHBoYV9hc3ltbWV0cnlfZGF0YV9jbGVhbl9sb25nLmNzdicsIHJvdy5uYW1lcyA9IEZBTFNFKQpgYGAKCiMgU3BlY3RyYWwgZGVjb21wb3NpdGlvbgoKYGBge3IgcGFydGljaXBhbnRzLCBmaWcud2lkdGggPSAxMn0Kb3B0aW9ucyh3aWR0aCA9IDEwMCkKbXl0YWJsZSA8LSB4dGFicyh+IEdyb3VwICsgU2V4LCBkYXRhID0gYWxwaGFfYXN5bW1ldHJ5X2RhdGEpCmZ0YWJsZShhZGRtYXJnaW5zKG15dGFibGUpKQpgYGAKCi0gKipJbmZpbml0eSBSZWZlcmVuY2Ugb3IgUmVmZXJlbmNlIEVsZWN0cm9kZSBTdGFuZGFyZGl6YXRpb24gVGVjaG5pcXVlIChSRVNUKSoqLlwKPGRpdiBjbGFzcz0iY3NsLWVudHJ5Ij5Eb25nLCBMLiwgTGksIEYuLCBMaXUsIFEuLCBXZW4sIFguLCBMYWksIFkuLCBYdSwgUC4sICYjMzg7IFlhbywgRC4gKDIwMTcpLiBNQVRMQUIgVG9vbGJveGVzIGZvciBSZWZlcmVuY2UgRWxlY3Ryb2RlIFN0YW5kYXJkaXphdGlvbiBUZWNobmlxdWUgKFJFU1QpIG9mIFNjYWxwIEVFRy4gPGk+RnJvbnRpZXJzIGluIE5ldXJvc2NpZW5jZTwvaT4sIDxpPjExPC9pPihPQ1QpLCA2MDEuIGh0dHBzOi8vZG9pLm9yZy8xMC4zMzg5L2ZuaW5zLjIwMTcuMDA2MDE8L2Rpdj4KLSAqKjEyMCBjb25zZWN1dGl2ZSBzZWdtZW50cywgNSBzZWNvbmRzIGVhY2gqKi5cCi0gKipQU0QgY29tcHV0ZWQgd2l0aCBXZWxjaCdzIG1ldGhvZCoqLgoKIyMgU2NhbHAgTWFwLCBtZWFuIHBvd2VyIDktMTEgSHoKIVtdKFJFU1RfUFNEX3NjYWxwLnBuZykKCiMjIFBTRCB0b3BvZ3JhcGh5LCAxIHRvIDU1IEh6LCBncmFuZCBhdmVyYWdlCiFbXShSRVNUX1BTRF90b3BvZ3JhcGh5LnBuZykKCiMjIFBTRCB0b3BvZ3JhcGh5LCA0IHRvIDMwIEh6LCBncmFuZCBhdmVyYWdlCiFbXShSRVNUX1BTRF90b3BvZ3JhcGh5XzRfdG9fMzBfSHoucG5nKQoKIyMgRnJvbnRhbCBFbGVjdHJvZGVzLCA0IHRvIDMwIEh6LCBncmFuZCBhdmVyYWdlCi0gMSBzdGFuZGFyZCBlcnJvciAgICAgYmFuZHMKIVtGcm9udGFsIGVsZWN0cm9kZXNdKFJFU1RfUFNEXzRfdG9fMzBfSHoucG5nKQoKIyBQYXJhbWV0ZXJpemF0aW9uIG9mIG5ldXJhbCBwb3dlciBzcGVjdHJhCiogRWFjaCBpbmRpdmlkdWFsIFBTRCBpcyByZWdhcmRlZCBhcyBhIGNvbWJpbmF0aW9uIG9mIGFuIGFwZXJpb2RpYyBjb21wb25lbnQsIGNoYXJhY3Rlcml6ZWQgYnkgYW4gKipvZmZzZXQqKiBhbmQgYSAqKnNsb3BlKiosIGFuZCBwdXRhdGl2ZSBwZXJpb2RpYyBvc2NpbGxhdG9yeSBwZWFrczoKIVtdKHBlcmlvZGljX2FuZF9hcGVyaW9kaWMucG5nKQoKCiogVGhlIHByb2Nlc3MgZ29lcyBsaWtlIHRoaXMuLi4KICArIFBvd2VyIFNwZWN0cnVtIFAgRml0OgokJFAgPSBMICsgXHN1bV97bj0wfV5OIEdfbiQkCiAgKyBXaGVyZToKJCRMID0gYXBlcmlvZGljXCBzaWduYWwkJAokJEdfbiA9IEdhdXNzaWFuXCBmaXRcIHRvXCBhXCBwZWFrJCQKJCROID0gdG90YWxcIHBlYWtzXCBleHRyYWN0ZWQkJAoKICArIEFwZXJpb2RpYyAoYmFja2dyb3VuZCkgRml0IEw6CiQkTCA9IGIgLSBcbG9nX3sxMH0oayArIEZeXGNoaSkkJAogICsgV2hlcmU6CiQkYiA9IGJyb2FkYmFuZFwgb2Zmc2V0JCQKJCRrID0gYGBrbmVlIlwgcGFyYW1ldGVyJCQKJCRGID0gdmVjdG9yXCBvZlwgZnJlcXVlbmNpZXMkJAokJFxjaGkgPSBhcGVyaW9kaWNcIHNsb3BlJCQKIVtdKGFsZ29yaXRobS5wbmcpCgo8ZGl2IGNsYXNzPSJjc2wtZW50cnkiPkRvbm9naHVlLCBULiwgSGFsbGVyLCBNLiwgUGV0ZXJzb24sIEUuIEouLCBWYXJtYSwgUC4sIFNlYmFzdGlhbiwgUC4sIEdhbywgUi4sIE5vdG8sIFQuLCBMYXJhLCBBLiBILiwgV2FsbGlzLCBKLiBELiwgS25pZ2h0LCBSLiBULiwgU2hlc3R5dWssIEEuLCAmIzM4OyBWb3l0ZWssIEIuICgyMDIwKS4gUGFyYW1ldGVyaXppbmcgbmV1cmFsIHBvd2VyIHNwZWN0cmEgaW50byBwZXJpb2RpYyBhbmQgYXBlcmlvZGljIGNvbXBvbmVudHMuIDxpPk5hdHVyZSBOZXVyb3NjaWVuY2U8L2k+LCA8aT4yMzwvaT4oMTIpLCAxNjU14oCTMTY2NS4gaHR0cHM6Ly9kb2kub3JnLzEwLjEwMzgvczQxNTkzLTAyMC0wMDc0NC14PC9kaXY+CgojIyBJbmRpdmlkdWFsIEVsZWN0cm9kZXMKIVtdKC4uL0Zyb250YWxfZm9vb2ZfcGxvdHNfYnlfZWxlY3Ryb2RlL3Jlc3RfMDU2X3Jlc3RpbmdfRnAxLnBuZykKIVtdKC4uL0Zyb250YWxfZm9vb2ZfcGxvdHNfYnlfZWxlY3Ryb2RlL3Jlc3RfMDU2X3Jlc3RpbmdfRnAyLnBuZykKCiMjIEVudGlyZSBTdWJqZWN0CiFbXSguLi9Gcm9udGFsX2Zvb29mX3Bsb3RzX2J5X3N1YmplY3QvMDU2X3Jlc3RpbmcucG5nKQoKIyBHZW5lcmFsIEFzeW1tZXRyeSBEZXNjcmlwdGlvbgpgYGB7ciBnZW5lcmFsLCBmaWcud2lkdGggPSAxMn0Kb3B0aW9ucyh3aWR0aCA9IDEwMCkKc3VtbWFyeShhbHBoYV9hc3ltbWV0cnlfZGF0YSkKYXN5bW1ldHJ5X3BhaXJzIDwtIGMoJ2FzeW1tZXRyeV9GcDJfRnAxJywgJ2FzeW1tZXRyeV9BRjRfQUYzJywgJ2FzeW1tZXRyeV9BRjhfQUY3JywgJ2FzeW1tZXRyeV9GMl9GMScsICdhc3ltbWV0cnlfRjRfRjMnLCAnYXN5bW1ldHJ5X0Y2X0Y1JywgJ2FzeW1tZXRyeV9GOF9GNycsICdhc3ltbWV0cnlfRkMyX0ZDMScpCmFzeW1tZXRyeV9wYWlyc19wYWlycyA8LSBnZ3BhaXJzKGFscGhhX2FzeW1tZXRyeV9kYXRhLAogICAgICAgICAgICAgICAgICAgICAgIGNvbHVtbnMgPSBhc3ltbWV0cnlfcGFpcnMsCiAgICAgICAgICAgICAgICAgICAgICAgYWVzKGNvbG91ciA9IEdyb3VwLCBhbHBoYSA9IC4yNSksCiAgICAgICAgICAgICAgICAgICAgICAgcHJvZ3Jlc3MgPSBGQUxTRSwKICAgICAgICAgICAgICAgICAgICAgICBsb3dlciA9IGxpc3QoY29udGludW91cyA9IHdyYXAoInBvaW50cyIpKSkKc3VwcHJlc3NXYXJuaW5ncyhwcmludChhc3ltbWV0cnlfcGFpcnNfcGFpcnMpKQpgYGAKCjwhLS0gIyBBbHBoYSBBc3ltbWV0cnkgPSAgJFxsb2dfezEwfShcYWxwaGFfe3JpZ2h0fSkgLSBcbG9nX3sxMH0oXGFscGhhX3tsZWZ0fSkkIC0tPgojIEFscGhhIEFzeW1tZXRyeSA9ICBwZWFrX3Bvd2VyKHJpZ2h0KSAtIHBlYWtfcG93ZXIobGVmdCkKCiMjIE1peGVkIE1vZGVsCgpgYGB7ciBtaXhlZCwgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCnN1bW1hcnkoYWxwaGFfYXN5bW1ldHJ5X2RhdGFfbG9uZykKYWxwaGFfYXN5bW1ldHJ5X2xtZXIgPSBsbWVyKEFzaW1tZXRyeSB+IEdyb3VwKlNleCArICgxfFN1YmplY3QpICArICgxfEVsZWN0cm9kZXMpLCBhbHBoYV9hc3ltbWV0cnlfZGF0YV9sb25nKQphZmV4X3Bsb3QoCiAgICBhbHBoYV9hc3ltbWV0cnlfbG1lciwKICAgIHggICAgID0gIkdyb3VwIiwKICAgIHRyYWNlID0gIlNleCIsCiAgICBpZCAgICA9ICJTdWJqZWN0IiwKICAgIGVycm9yX2FyZyA9IGxpc3Qod2lkdGggPSAuMSksCiAgICBkb2RnZSA9IC4zLAogICAgZGF0YV9hcmcgID0gbGlzdCgKICAgICAgcG9zaXRpb24gPSAKICAgICAgICBwb3NpdGlvbl9qaXR0ZXJkb2RnZSgKICAgICAgICAgIGppdHRlci53aWR0aCA9IC4xLCAKICAgICAgICAgIGRvZGdlLndpZHRoICA9IC4zICAjIyBuZWVkcyB0byBiZSBzYW1lIGFzIGRvZGdlCiAgICAgICAgKSksCiAgICBtYXBwaW5nID0gYygiY29sb3IiKSwKICAgIHBvaW50X2FyZyA9IGxpc3Qoc2l6ZSA9IDQpCiAgKQpzdW1tYXJ5KGFscGhhX2FzeW1tZXRyeV9sbWVyKQptb2RlbF9wZXJmb3JtYW5jZShhbHBoYV9hc3ltbWV0cnlfbG1lcikKcmFub3ZhKGFscGhhX2FzeW1tZXRyeV9sbWVyKQphbm92YShhbHBoYV9hc3ltbWV0cnlfbG1lcikKaW50ZXJwcmV0KG9tZWdhX3NxdWFyZWQoYWxwaGFfYXN5bW1ldHJ5X2xtZXIsIGFsdGVybmF0aXZlID0gJ3R3by5zaWRlZCcpLCBydWxlcyA9ICdmaWVsZDIwMTMnKQphX3Bvc3RlcmlvcmlfbG1lcihhbHBoYV9hc3ltbWV0cnlfbG1lcikKYGBgCgojIyBGcDItRnAxIHBhaXIKCmBgYHtyIGFzeW1tZXRyeV9GcDJfRnAxLCBmaWcud2lkdGggPSAxMn0Kb3B0aW9ucyh3aWR0aCA9IDEwMCkKYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSA9IGFvdl9leigiU3ViamVjdCIsICJhc3ltbWV0cnlfRnAyX0ZwMSIsIGFscGhhX2FzeW1tZXRyeV9kYXRhLCBiZXR3ZWVuID0gYygiR3JvdXAiLCAiU2V4IikpCm15dGFibGUgPC0geHRhYnMofiBHcm91cCArIFNleCwgZGF0YSA9IGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEkZGF0YSRsb25nKQpmdGFibGUoYWRkbWFyZ2lucyhteXRhYmxlKSkKIyBhbHBoYV9hc3ltbWV0cnlfcmFpbiA8LSBnZ3Bsb3QoYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSRkYXRhJGxvbmcsIGFlcyh5ID0gYXN5bW1ldHJ5X0ZwMl9GcDEsIHggPSBHcm91cCwgY29sb3IgPSBTZXgsIGZpbGwgPSBTZXgpKSArCiMgICAjIGdndGl0bGUoImFscGhhX2FzeW1tZXRyeSIpICsKIyAgIHlsYWIoInBvd2VyIikgKwojICAgc3RhdF9oYWxmZXllKAojICAgICB0cmltICAgPSBGQUxTRSwgCiMgICAgIGFkanVzdCA9IDEsIAojICAgICAud2lkdGggPSAwLCAKIyAgICAganVzdGlmaWNhdGlvbiA9IC0uMTUsIAojICAgICBhbHBoYSAgPSAuNSwKIyAgICAgcG9pbnRfY29sb3VyID0gTkEpICsgCiMgICAjIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0nbm9uZScpCiMgICAjIGdlb21fYm94cGxvdCh3aWR0aCA9IC4xNSwgYWxwaGEgPSAuMiwgb3V0bGllci5zaGFwZSA9IE5BKSArCiMgICBnZW9tX3BvaW50KHNpemUgPSAyLCBhbHBoYSA9IC40LCBwb3NpdGlvbiA9IHBvc2l0aW9uX2ppdHRlcih3aWR0aCA9IC4wNSwgaGVpZ2h0ID0gMCkpIAojIHN1cHByZXNzV2FybmluZ3MocHJpbnQoYWxwaGFfYXN5bW1ldHJ5X3JhaW4pKQphbHBoYV9hc3ltbWV0cnlfYWZleF9wbG90IDwtCiAgYWZleF9wbG90KAogICAgYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSwKICAgIHggPSAiR3JvdXAiLAogICAgdHJhY2UgPSAiU2V4IiwKICAgIGVycm9yID0gImJldHdlZW4iLAogICAgZXJyb3JfYXJnID0gbGlzdCh3aWR0aCA9IC4xKSwKICAgIGRvZGdlID0gLjMsCiAgICBkYXRhX2FyZyAgPSBsaXN0KAogICAgICBwb3NpdGlvbiA9IAogICAgICAgIHBvc2l0aW9uX2ppdHRlcmRvZGdlKAogICAgICAgICAgaml0dGVyLndpZHRoICA9IC4xLCAKICAgICAgICAgIGppdHRlci5oZWlnaHQgPSAuMSwgCiAgICAgICAgICBkb2RnZS53aWR0aCAgID0gLjMgICMjIG5lZWRzIHRvIGJlIHNhbWUgYXMgZG9kZ2UKICAgICAgICApKSwKICAgIG1hcHBpbmcgPSBjKCJjb2xvciIpLAogICAgcG9pbnRfYXJnID0gbGlzdChzaXplID0gNCkKICApCnN1cHByZXNzV2FybmluZ3MocHJpbnQoYWxwaGFfYXN5bW1ldHJ5X2FmZXhfcGxvdCkpCm5pY2UoYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSkKYV9wb3N0ZXJpb3JpKGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEpCmBgYAoKIyMgQUY0LUFGMyBwYWlyCgpgYGB7ciBhc3ltbWV0cnlfQUY0X0FGMywgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCmFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEgPSBhb3ZfZXooIlN1YmplY3QiLCAiYXN5bW1ldHJ5X0FGNF9BRjMiLCBhbHBoYV9hc3ltbWV0cnlfZGF0YSwgYmV0d2VlbiA9IGMoIkdyb3VwIiwgIlNleCIpKQpteXRhYmxlIDwtIHh0YWJzKH4gR3JvdXAgKyBTZXgsIGRhdGEgPSBhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhJGRhdGEkbG9uZykKZnRhYmxlKGFkZG1hcmdpbnMobXl0YWJsZSkpCmFscGhhX2FzeW1tZXRyeV9hZmV4X3Bsb3QgPC0KICBhZmV4X3Bsb3QoCiAgICBhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhLAogICAgeCA9ICJHcm91cCIsCiAgICB0cmFjZSA9ICJTZXgiLAogICAgZXJyb3IgPSAiYmV0d2VlbiIsCiAgICBlcnJvcl9hcmcgPSBsaXN0KHdpZHRoID0gLjEpLAogICAgZG9kZ2UgPSAuMywKICAgIGRhdGFfYXJnICA9IGxpc3QoCiAgICAgIHBvc2l0aW9uID0gCiAgICAgICAgcG9zaXRpb25faml0dGVyZG9kZ2UoCiAgICAgICAgICBqaXR0ZXIud2lkdGggID0gLjEsIAogICAgICAgICAgaml0dGVyLmhlaWdodCA9IC4xLCAKICAgICAgICAgIGRvZGdlLndpZHRoICAgPSAuMyAgIyMgbmVlZHMgdG8gYmUgc2FtZSBhcyBkb2RnZQogICAgICAgICkpLAogICAgbWFwcGluZyA9IGMoImNvbG9yIiksCiAgICBwb2ludF9hcmcgPSBsaXN0KHNpemUgPSA0KQogICkKc3VwcHJlc3NXYXJuaW5ncyhwcmludChhbHBoYV9hc3ltbWV0cnlfYWZleF9wbG90KSkKbmljZShhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhKQphX3Bvc3RlcmlvcmkoYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSkKYGBgCgojIyBBRjgtQUY3IHBhaXIKCmBgYHtyIGFzeW1tZXRyeV9BRjhfQUY3LCBmaWcud2lkdGggPSAxMn0Kb3B0aW9ucyh3aWR0aCA9IDEwMCkKYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSA9IGFvdl9leigiU3ViamVjdCIsICJhc3ltbWV0cnlfQUY4X0FGNyIsIGFscGhhX2FzeW1tZXRyeV9kYXRhLCBiZXR3ZWVuID0gYygiR3JvdXAiLCAiU2V4IikpCm15dGFibGUgPC0geHRhYnMofiBHcm91cCArIFNleCwgZGF0YSA9IGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEkZGF0YSRsb25nKQpmdGFibGUoYWRkbWFyZ2lucyhteXRhYmxlKSkKYWxwaGFfYXN5bW1ldHJ5X2FmZXhfcGxvdCA8LQogIGFmZXhfcGxvdCgKICAgIGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEsCiAgICB4ID0gIkdyb3VwIiwKICAgIHRyYWNlID0gIlNleCIsCiAgICBlcnJvciA9ICJiZXR3ZWVuIiwKICAgIGVycm9yX2FyZyA9IGxpc3Qod2lkdGggPSAuMSksCiAgICBkb2RnZSA9IC4zLAogICAgZGF0YV9hcmcgID0gbGlzdCgKICAgICAgcG9zaXRpb24gPSAKICAgICAgICBwb3NpdGlvbl9qaXR0ZXJkb2RnZSgKICAgICAgICAgIGppdHRlci53aWR0aCAgPSAuMSwgCiAgICAgICAgICBqaXR0ZXIuaGVpZ2h0ID0gLjEsIAogICAgICAgICAgZG9kZ2Uud2lkdGggICA9IC4zICAjIyBuZWVkcyB0byBiZSBzYW1lIGFzIGRvZGdlCiAgICAgICAgKSksCiAgICBtYXBwaW5nID0gYygiY29sb3IiKSwKICAgIHBvaW50X2FyZyA9IGxpc3Qoc2l6ZSA9IDQpCiAgKQpzdXBwcmVzc1dhcm5pbmdzKHByaW50KGFscGhhX2FzeW1tZXRyeV9hZmV4X3Bsb3QpKQpuaWNlKGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEpCmFfcG9zdGVyaW9yaShhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhKQpgYGAKIyMgRjItRjEgcGFpcgoKYGBge3IgYXN5bW1ldHJ5X0YyX0YxLCBmaWcud2lkdGggPSAxMn0Kb3B0aW9ucyh3aWR0aCA9IDEwMCkKYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSA9IGFvdl9leigiU3ViamVjdCIsICJhc3ltbWV0cnlfRjJfRjEiLCBhbHBoYV9hc3ltbWV0cnlfZGF0YSwgYmV0d2VlbiA9IGMoIkdyb3VwIiwgIlNleCIpKQpteXRhYmxlIDwtIHh0YWJzKH4gR3JvdXAgKyBTZXgsIGRhdGEgPSBhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhJGRhdGEkbG9uZykKZnRhYmxlKGFkZG1hcmdpbnMobXl0YWJsZSkpCmFscGhhX2FzeW1tZXRyeV9hZmV4X3Bsb3QgPC0KICBhZmV4X3Bsb3QoCiAgICBhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhLAogICAgeCA9ICJHcm91cCIsCiAgICB0cmFjZSA9ICJTZXgiLAogICAgZXJyb3IgPSAiYmV0d2VlbiIsCiAgICBlcnJvcl9hcmcgPSBsaXN0KHdpZHRoID0gLjEpLAogICAgZG9kZ2UgPSAuMywKICAgIGRhdGFfYXJnICA9IGxpc3QoCiAgICAgIHBvc2l0aW9uID0gCiAgICAgICAgcG9zaXRpb25faml0dGVyZG9kZ2UoCiAgICAgICAgICBqaXR0ZXIud2lkdGggID0gLjEsIAogICAgICAgICAgaml0dGVyLmhlaWdodCA9IC4xLCAKICAgICAgICAgIGRvZGdlLndpZHRoICAgPSAuMyAgIyMgbmVlZHMgdG8gYmUgc2FtZSBhcyBkb2RnZQogICAgICAgICkpLAogICAgbWFwcGluZyA9IGMoImNvbG9yIiksCiAgICBwb2ludF9hcmcgPSBsaXN0KHNpemUgPSA0KQogICkKc3VwcHJlc3NXYXJuaW5ncyhwcmludChhbHBoYV9hc3ltbWV0cnlfYWZleF9wbG90KSkKbmljZShhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhKQphX3Bvc3RlcmlvcmkoYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSkKYGBgCgojIyBGNC1GMyBwYWlyCgpgYGB7ciBhc3ltbWV0cnlfRjRfRjMsIGZpZy53aWR0aCA9IDEyfQpvcHRpb25zKHdpZHRoID0gMTAwKQphbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhID0gYW92X2V6KCJTdWJqZWN0IiwgImFzeW1tZXRyeV9GNF9GMyIsIGFscGhhX2FzeW1tZXRyeV9kYXRhLCBiZXR3ZWVuID0gYygiR3JvdXAiLCAiU2V4IikpCm15dGFibGUgPC0geHRhYnMofiBHcm91cCArIFNleCwgZGF0YSA9IGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEkZGF0YSRsb25nKQpmdGFibGUoYWRkbWFyZ2lucyhteXRhYmxlKSkKYWxwaGFfYXN5bW1ldHJ5X2FmZXhfcGxvdCA8LQogIGFmZXhfcGxvdCgKICAgIGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEsCiAgICB4ID0gIkdyb3VwIiwKICAgIHRyYWNlID0gIlNleCIsCiAgICBlcnJvciA9ICJiZXR3ZWVuIiwKICAgIGVycm9yX2FyZyA9IGxpc3Qod2lkdGggPSAuMSksCiAgICBkb2RnZSA9IC4zLAogICAgZGF0YV9hcmcgID0gbGlzdCgKICAgICAgcG9zaXRpb24gPSAKICAgICAgICBwb3NpdGlvbl9qaXR0ZXJkb2RnZSgKICAgICAgICAgIGppdHRlci53aWR0aCAgPSAuMSwgCiAgICAgICAgICBqaXR0ZXIuaGVpZ2h0ID0gLjEsIAogICAgICAgICAgZG9kZ2Uud2lkdGggICA9IC4zICAjIyBuZWVkcyB0byBiZSBzYW1lIGFzIGRvZGdlCiAgICAgICAgKSksCiAgICBtYXBwaW5nID0gYygiY29sb3IiKSwKICAgIHBvaW50X2FyZyA9IGxpc3Qoc2l6ZSA9IDQpCiAgKQpzdXBwcmVzc1dhcm5pbmdzKHByaW50KGFscGhhX2FzeW1tZXRyeV9hZmV4X3Bsb3QpKQpuaWNlKGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEpCmFfcG9zdGVyaW9yaShhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhKQpgYGAKCiMjIEY2LUY1IHBhaXIKCmBgYHtyIGFzeW1tZXRyeV9GNl9GNSwgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCmFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEgPSBhb3ZfZXooIlN1YmplY3QiLCAiYXN5bW1ldHJ5X0Y2X0Y1IiwgYWxwaGFfYXN5bW1ldHJ5X2RhdGEsIGJldHdlZW4gPSBjKCJHcm91cCIsICJTZXgiKSkKbXl0YWJsZSA8LSB4dGFicyh+IEdyb3VwICsgU2V4LCBkYXRhID0gYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSRkYXRhJGxvbmcpCmZ0YWJsZShhZGRtYXJnaW5zKG15dGFibGUpKQphbHBoYV9hc3ltbWV0cnlfYWZleF9wbG90IDwtCiAgYWZleF9wbG90KAogICAgYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSwKICAgIHggPSAiR3JvdXAiLAogICAgdHJhY2UgPSAiU2V4IiwKICAgIGVycm9yID0gImJldHdlZW4iLAogICAgZXJyb3JfYXJnID0gbGlzdCh3aWR0aCA9IC4xKSwKICAgIGRvZGdlID0gLjMsCiAgICBkYXRhX2FyZyAgPSBsaXN0KAogICAgICBwb3NpdGlvbiA9IAogICAgICAgIHBvc2l0aW9uX2ppdHRlcmRvZGdlKAogICAgICAgICAgaml0dGVyLndpZHRoICA9IC4xLCAKICAgICAgICAgIGppdHRlci5oZWlnaHQgPSAuMSwgCiAgICAgICAgICBkb2RnZS53aWR0aCAgID0gLjMgICMjIG5lZWRzIHRvIGJlIHNhbWUgYXMgZG9kZ2UKICAgICAgICApKSwKICAgIG1hcHBpbmcgPSBjKCJjb2xvciIpLAogICAgcG9pbnRfYXJnID0gbGlzdChzaXplID0gNCkKICApCnN1cHByZXNzV2FybmluZ3MocHJpbnQoYWxwaGFfYXN5bW1ldHJ5X2FmZXhfcGxvdCkpCm5pY2UoYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSkKYV9wb3N0ZXJpb3JpKGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEpCmBgYAoKIyMgRjgtRjcgcGFpcgoKYGBge3IgYWFzeW1tZXRyeV9GOF9GNywgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCmFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEgPSBhb3ZfZXooIlN1YmplY3QiLCAiYXN5bW1ldHJ5X0Y4X0Y3IiwgYWxwaGFfYXN5bW1ldHJ5X2RhdGEsIGJldHdlZW4gPSBjKCJHcm91cCIsICJTZXgiKSkKbXl0YWJsZSA8LSB4dGFicyh+IEdyb3VwICsgU2V4LCBkYXRhID0gYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSRkYXRhJGxvbmcpCmZ0YWJsZShhZGRtYXJnaW5zKG15dGFibGUpKQphbHBoYV9hc3ltbWV0cnlfYWZleF9wbG90IDwtCiAgYWZleF9wbG90KAogICAgYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSwKICAgIHggPSAiR3JvdXAiLAogICAgdHJhY2UgPSAiU2V4IiwKICAgIGVycm9yID0gImJldHdlZW4iLAogICAgZXJyb3JfYXJnID0gbGlzdCh3aWR0aCA9IC4xKSwKICAgIGRvZGdlID0gLjMsCiAgICBkYXRhX2FyZyAgPSBsaXN0KAogICAgICBwb3NpdGlvbiA9IAogICAgICAgIHBvc2l0aW9uX2ppdHRlcmRvZGdlKAogICAgICAgICAgaml0dGVyLndpZHRoICA9IC4xLCAKICAgICAgICAgIGppdHRlci5oZWlnaHQgPSAuMSwgCiAgICAgICAgICBkb2RnZS53aWR0aCAgID0gLjMgICMjIG5lZWRzIHRvIGJlIHNhbWUgYXMgZG9kZ2UKICAgICAgICApKSwKICAgIG1hcHBpbmcgPSBjKCJjb2xvciIpLAogICAgcG9pbnRfYXJnID0gbGlzdChzaXplID0gNCkKICApCnN1cHByZXNzV2FybmluZ3MocHJpbnQoYWxwaGFfYXN5bW1ldHJ5X2FmZXhfcGxvdCkpCm5pY2UoYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSkKYV9wb3N0ZXJpb3JpKGFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEpCmBgYAoKIyMgRkMyLUZDMSBwYWlyCgpgYGB7ciBhc3ltbWV0cnlfRkMyX0ZDMSwgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCmFscGhhX2FzeW1tZXRyeV9yZXBfYW5vdmEgPSBhb3ZfZXooIlN1YmplY3QiLCAiYXN5bW1ldHJ5X0ZDMl9GQzEiLCBhbHBoYV9hc3ltbWV0cnlfZGF0YSwgYmV0d2VlbiA9IGMoIkdyb3VwIiwgIlNleCIpKQpteXRhYmxlIDwtIHh0YWJzKH4gR3JvdXAgKyBTZXgsIGRhdGEgPSBhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhJGRhdGEkbG9uZykKZnRhYmxlKGFkZG1hcmdpbnMobXl0YWJsZSkpCmFscGhhX2FzeW1tZXRyeV9hZmV4X3Bsb3QgPC0KICBhZmV4X3Bsb3QoCiAgICBhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhLAogICAgeCA9ICJHcm91cCIsCiAgICB0cmFjZSA9ICJTZXgiLAogICAgZXJyb3IgPSAiYmV0d2VlbiIsCiAgICBlcnJvcl9hcmcgPSBsaXN0KHdpZHRoID0gLjEpLAogICAgZG9kZ2UgPSAuMywKICAgIGRhdGFfYXJnICA9IGxpc3QoCiAgICAgIHBvc2l0aW9uID0gCiAgICAgICAgcG9zaXRpb25faml0dGVyZG9kZ2UoCiAgICAgICAgICBqaXR0ZXIud2lkdGggID0gLjEsIAogICAgICAgICAgaml0dGVyLmhlaWdodCA9IC4xLCAKICAgICAgICAgIGRvZGdlLndpZHRoICAgPSAuMyAgIyMgbmVlZHMgdG8gYmUgc2FtZSBhcyBkb2RnZQogICAgICAgICkpLAogICAgbWFwcGluZyA9IGMoImNvbG9yIiksCiAgICBwb2ludF9hcmcgPSBsaXN0KHNpemUgPSA0KQogICkKc3VwcHJlc3NXYXJuaW5ncyhwcmludChhbHBoYV9hc3ltbWV0cnlfYWZleF9wbG90KSkKbmljZShhbHBoYV9hc3ltbWV0cnlfcmVwX2Fub3ZhKQphX3Bvc3RlcmlvcmkoYWxwaGFfYXN5bW1ldHJ5X3JlcF9hbm92YSkKYGBgCgojIEFscGhhIHBlYWsgcGFyYW1ldGVycyBhbmQgYXBlcmlvZGljIGFjdGl2aXR5CgpgYGB7ciBwZWFrX3BhcmFtcywgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCnN1bW1hcnkoYWxwaGFfcG93ZXJfZGF0YSkKc3BlY19wYXJhbXMgPC0gYygnQW1wbGl0dWRlJywgJ0ZyZXF1ZW5jeScsICdXaWR0aCcsICdPZmZzZXQnLCAnU2xvcGUnKQpzcGVjX3BhcmFtc19wYWlycyA8LSBnZ3BhaXJzKGFscGhhX3Bvd2VyX2RhdGEsCiAgICAgICAgICAgICAgICAgICAgICAgY29sdW1ucyA9IHNwZWNfcGFyYW1zLAogICAgICAgICAgICAgICAgICAgICAgIGFlcyhjb2xvdXIgPSBHcm91cCwgYWxwaGEgPSAuMjUpLAogICAgICAgICAgICAgICAgICAgICAgIHByb2dyZXNzID0gRkFMU0UsCiAgICAgICAgICAgICAgICAgICAgICAgbG93ZXIgPSBsaXN0KGNvbnRpbnVvdXMgPSB3cmFwKCJwb2ludHMiKSkpCnN1cHByZXNzV2FybmluZ3MocHJpbnQoc3BlY19wYXJhbXNfcGFpcnMpKQpgYGAKCiMjIEFscGhhIFBvd2VyCgpgYGB7ciBhbHBoYV9wb3dlciwgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCmFscGhhX3Bvd19yZXBfYW5vdmEgPSBhb3ZfZXooIlN1YmplY3QiLCAiQW1wbGl0dWRlIiwgYWxwaGFfcG93ZXJfZGF0YSwgYmV0d2VlbiA9IGMoIkdyb3VwIiwgIlNleCIpLCB3aXRoaW4gPSBjKCJoZW1pc3BoZXJlIiwgImVsZWN0cm9kZV9wYWlyIikpCm15dGFibGUgPC0geHRhYnMofiBHcm91cCArIFNleCwgZGF0YSA9IGFscGhhX3Bvd19yZXBfYW5vdmEkZGF0YSRsb25nKSAvIGxlbmd0aChsZXZlbHMoYWxwaGFfcG93X3JlcF9hbm92YSRkYXRhJGxvbmckZWxlY3Ryb2RlX3BhaXIpKSAvIGxlbmd0aChsZXZlbHMoYWxwaGFfcG93X3JlcF9hbm92YSRkYXRhJGxvbmckaGVtaXNwaGVyZSkpCmZ0YWJsZShhZGRtYXJnaW5zKG15dGFibGUpKQphbHBoYV9wb3dfYWZleF9wbG90IDwtCiAgYWZleF9wbG90KAogICAgYWxwaGFfcG93X3JlcF9hbm92YSwKICAgIHggPSAiR3JvdXAiLAogICAgdHJhY2UgPSAiU2V4IiwKICAgIHBhbmVsID0gImhlbWlzcGhlcmUiLAogICAgZXJyb3IgPSAiYmV0d2VlbiIsCiAgICBlcnJvcl9hcmcgPSBsaXN0KHdpZHRoID0gLjEpLAogICAgZG9kZ2UgPSAuMywKICAgIGRhdGFfYXJnICA9IGxpc3QoCiAgICAgIHBvc2l0aW9uID0gCiAgICAgICAgcG9zaXRpb25faml0dGVyZG9kZ2UoCiAgICAgICAgICBqaXR0ZXIud2lkdGggID0gLjEsIAogICAgICAgICAgaml0dGVyLmhlaWdodCA9IC4xLCAKICAgICAgICAgIGRvZGdlLndpZHRoICAgPSAuMyAgIyMgbmVlZHMgdG8gYmUgc2FtZSBhcyBkb2RnZQogICAgICAgICkpLAogICAgbWFwcGluZyA9IGMoImNvbG9yIiksCiAgICBwb2ludF9hcmcgPSBsaXN0KHNpemUgPSA0KQogICkKc3VwcHJlc3NXYXJuaW5ncyhwcmludChhbHBoYV9wb3dfYWZleF9wbG90KSkKbmljZShhbHBoYV9wb3dfcmVwX2Fub3ZhKQojIGFfcG9zdGVyaW9yaShhbHBoYV9wb3dfcmVwX2Fub3ZhKQplbW1lYW5zKGFscGhhX3Bvd19yZXBfYW5vdmEsIHBhaXJ3aXNlIH4gaGVtaXNwaGVyZSB8IEdyb3VwKQpjYXQocmVwKCJfIiwgODApLCAnXG4nLCBzZXAgPSAiIikKZW1tZWFucyhhbHBoYV9wb3dfcmVwX2Fub3ZhLCBwYWlyd2lzZSB+IGVsZWN0cm9kZV9wYWlyKQpjYXQocmVwKCJfIiwgODApLCAnXG4nLCBzZXAgPSAiIikKZW1tZWFucyhhbHBoYV9wb3dfcmVwX2Fub3ZhLCBwYWlyd2lzZSB+IGVsZWN0cm9kZV9wYWlyIHwgaGVtaXNwaGVyZSkKYGBgCgojIyBBbHBoYSBGcmVxdWVuY3kKCmBgYHtyIGFscGhhX2ZyZXF1ZW5jeSwgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCmFscGhhX2ZyZXFfcmVwX2Fub3ZhID0gYW92X2V6KCJTdWJqZWN0IiwgIkZyZXF1ZW5jeSIsIGFscGhhX3Bvd2VyX2RhdGEsIGJldHdlZW4gPSBjKCJHcm91cCIsICJTZXgiKSwgd2l0aGluID0gYygiaGVtaXNwaGVyZSIsICJlbGVjdHJvZGVfcGFpciIpKQpteXRhYmxlIDwtIHh0YWJzKH4gR3JvdXAgKyBTZXgsIGRhdGEgPSBhbHBoYV9mcmVxX3JlcF9hbm92YSRkYXRhJGxvbmcpIC8gbGVuZ3RoKGxldmVscyhhbHBoYV9mcmVxX3JlcF9hbm92YSRkYXRhJGxvbmckZWxlY3Ryb2RlX3BhaXIpKSAvIGxlbmd0aChsZXZlbHMoYWxwaGFfZnJlcV9yZXBfYW5vdmEkZGF0YSRsb25nJGhlbWlzcGhlcmUpKQpmdGFibGUoYWRkbWFyZ2lucyhteXRhYmxlKSkKYWxwaGFfZnJlcV9hZmV4X3Bsb3QgPC0KICBhZmV4X3Bsb3QoCiAgICBhbHBoYV9mcmVxX3JlcF9hbm92YSwKICAgIHggPSAiR3JvdXAiLAogICAgdHJhY2UgPSAiU2V4IiwKICAgIHBhbmVsID0gImhlbWlzcGhlcmUiLAogICAgZXJyb3IgPSAiYmV0d2VlbiIsCiAgICBlcnJvcl9hcmcgPSBsaXN0KHdpZHRoID0gLjEpLAogICAgZG9kZ2UgPSAuMywKICAgIGRhdGFfYXJnICA9IGxpc3QoCiAgICAgIHBvc2l0aW9uID0gCiAgICAgICAgcG9zaXRpb25faml0dGVyZG9kZ2UoCiAgICAgICAgICBqaXR0ZXIud2lkdGggID0gLjEsIAogICAgICAgICAgaml0dGVyLmhlaWdodCA9IC4xLCAKICAgICAgICAgIGRvZGdlLndpZHRoICAgPSAuMyAgIyMgbmVlZHMgdG8gYmUgc2FtZSBhcyBkb2RnZQogICAgICAgICkpLAogICAgbWFwcGluZyA9IGMoImNvbG9yIiksCiAgICBwb2ludF9hcmcgPSBsaXN0KHNpemUgPSA0KQogICkKc3VwcHJlc3NXYXJuaW5ncyhwcmludChhbHBoYV9mcmVxX2FmZXhfcGxvdCkpCm5pY2UoYWxwaGFfZnJlcV9yZXBfYW5vdmEpCmFfcG9zdGVyaW9yaShhbHBoYV9mcmVxX3JlcF9hbm92YSkKYGBgCiMjIEFscGhhIEJhbmR3aWR0aAoKYGBge3IgYWxwaGFfYmFuZHdpZHRoLCBmaWcud2lkdGggPSAxMn0Kb3B0aW9ucyh3aWR0aCA9IDEwMCkKYWxwaGFfYndkX3JlcF9hbm92YSA9IGFvdl9leigiU3ViamVjdCIsICJXaWR0aCIsIGFscGhhX3Bvd2VyX2RhdGEsIGJldHdlZW4gPSBjKCJHcm91cCIsICJTZXgiKSwgd2l0aGluID0gYygiaGVtaXNwaGVyZSIsICJlbGVjdHJvZGVfcGFpciIpKQpteXRhYmxlIDwtIHh0YWJzKH4gR3JvdXAgKyBTZXgsIGRhdGEgPSBhbHBoYV9id2RfcmVwX2Fub3ZhJGRhdGEkbG9uZykgLyBsZW5ndGgobGV2ZWxzKGFscGhhX2J3ZF9yZXBfYW5vdmEkZGF0YSRsb25nJGVsZWN0cm9kZV9wYWlyKSkgLyBsZW5ndGgobGV2ZWxzKGFscGhhX2J3ZF9yZXBfYW5vdmEkZGF0YSRsb25nJGhlbWlzcGhlcmUpKQpmdGFibGUoYWRkbWFyZ2lucyhteXRhYmxlKSkKYWxwaGFfYndkX2FmZXhfcGxvdCA8LQogIGFmZXhfcGxvdCgKICAgIGFscGhhX2J3ZF9yZXBfYW5vdmEsCiAgICB4ID0gIkdyb3VwIiwKICAgIHRyYWNlID0gIlNleCIsCiAgICBwYW5lbCA9ICJoZW1pc3BoZXJlIiwKICAgIGVycm9yID0gImJldHdlZW4iLAogICAgZXJyb3JfYXJnID0gbGlzdCh3aWR0aCA9IC4xKSwKICAgIGRvZGdlID0gLjMsCiAgICBkYXRhX2FyZyAgPSBsaXN0KAogICAgICBwb3NpdGlvbiA9IAogICAgICAgIHBvc2l0aW9uX2ppdHRlcmRvZGdlKAogICAgICAgICAgaml0dGVyLndpZHRoICA9IC4xLCAKICAgICAgICAgIGppdHRlci5oZWlnaHQgPSAuMSwgCiAgICAgICAgICBkb2RnZS53aWR0aCAgID0gLjMgICMjIG5lZWRzIHRvIGJlIHNhbWUgYXMgZG9kZ2UKICAgICAgICApKSwKICAgIG1hcHBpbmcgPSBjKCJjb2xvciIpLAogICAgcG9pbnRfYXJnID0gbGlzdChzaXplID0gNCkKICApCnN1cHByZXNzV2FybmluZ3MocHJpbnQoYWxwaGFfYndkX2FmZXhfcGxvdCkpCm5pY2UoYWxwaGFfYndkX3JlcF9hbm92YSkKYV9wb3N0ZXJpb3JpKGFscGhhX2J3ZF9yZXBfYW5vdmEpCmBgYAoKIyMgQXBlcmlvZGljIE9mZnNldAoKYGBge3IgYXBlcmlvZGljX29mZnNldCwgZmlnLndpZHRoID0gMTJ9Cm9wdGlvbnMod2lkdGggPSAxMDApCmFwZXJpb2RpY19vZmZzZXRfcmVwX2Fub3ZhID0gYW92X2V6KCJTdWJqZWN0IiwgIk9mZnNldCIsIGFscGhhX3Bvd2VyX2RhdGEsIGJldHdlZW4gPSBjKCJHcm91cCIsICJTZXgiKSwgd2l0aGluID0gYygiaGVtaXNwaGVyZSIsICJlbGVjdHJvZGVfcGFpciIpKQpteXRhYmxlIDwtIHh0YWJzKH4gR3JvdXAgKyBTZXgsIGRhdGEgPSBhcGVyaW9kaWNfb2Zmc2V0X3JlcF9hbm92YSRkYXRhJGxvbmcpIC8gbGVuZ3RoKGxldmVscyhhcGVyaW9kaWNfb2Zmc2V0X3JlcF9hbm92YSRkYXRhJGxvbmckZWxlY3Ryb2RlX3BhaXIpKSAvIGxlbmd0aChsZXZlbHMoYXBlcmlvZGljX29mZnNldF9yZXBfYW5vdmEkZGF0YSRsb25nJGhlbWlzcGhlcmUpKQpmdGFibGUoYWRkbWFyZ2lucyhteXRhYmxlKSkKYXBlcmlvZGljX3BhcmFtX2FmZXhfcGxvdCA8LQogIGFmZXhfcGxvdCgKICAgIGFwZXJpb2RpY19vZmZzZXRfcmVwX2Fub3ZhLAogICAgeCA9ICJHcm91cCIsCiAgICB0cmFjZSA9ICJTZXgiLAogICAgcGFuZWwgPSAiaGVtaXNwaGVyZSIsCiAgICBlcnJvciA9ICJiZXR3ZWVuIiwKICAgIGVycm9yX2FyZyA9IGxpc3Qod2lkdGggPSAuMSksCiAgICBkb2RnZSA9IC4zLAogICAgZGF0YV9hcmcgID0gbGlzdCgKICAgICAgcG9zaXRpb24gPSAKICAgICAgICBwb3NpdGlvbl9qaXR0ZXJkb2RnZSgKICAgICAgICAgIGppdHRlci53aWR0aCAgPSAuMSwgCiAgICAgICAgICBqaXR0ZXIuaGVpZ2h0ID0gLjEsIAogICAgICAgICAgZG9kZ2Uud2lkdGggICA9IC4zICAjIyBuZWVkcyB0byBiZSBzYW1lIGFzIGRvZGdlCiAgICAgICAgKSksCiAgICBtYXBwaW5nID0gYygiY29sb3IiKSwKICAgIHBvaW50X2FyZyA9IGxpc3Qoc2l6ZSA9IDQpCiAgKQpzdXBwcmVzc1dhcm5pbmdzKHByaW50KGFwZXJpb2RpY19wYXJhbV9hZmV4X3Bsb3QpKQpuaWNlKGFwZXJpb2RpY19vZmZzZXRfcmVwX2Fub3ZhKQphX3Bvc3RlcmlvcmkoYXBlcmlvZGljX29mZnNldF9yZXBfYW5vdmEpCmBgYAoKIyMgQXBlcmlvZGljIFNsb3BlCgpgYGB7ciBhcGVyaW9kaWNfc2xvcGUsIGZpZy53aWR0aCA9IDEyfQpvcHRpb25zKHdpZHRoID0gMTAwKQphcGVyaW9kaWNfc2xvcGVfcmVwX2Fub3ZhID0gYW92X2V6KCJTdWJqZWN0IiwgIlNsb3BlIiwgYWxwaGFfcG93ZXJfZGF0YSwgYmV0d2VlbiA9IGMoIkdyb3VwIiwgIlNleCIpLCB3aXRoaW4gPSBjKCJoZW1pc3BoZXJlIiwgImVsZWN0cm9kZV9wYWlyIikpCm15dGFibGUgPC0geHRhYnMofiBHcm91cCArIFNleCwgZGF0YSA9IGFwZXJpb2RpY19zbG9wZV9yZXBfYW5vdmEkZGF0YSRsb25nKSAvIGxlbmd0aChsZXZlbHMoYXBlcmlvZGljX3Nsb3BlX3JlcF9hbm92YSRkYXRhJGxvbmckZWxlY3Ryb2RlX3BhaXIpKSAvIGxlbmd0aChsZXZlbHMoYXBlcmlvZGljX3Nsb3BlX3JlcF9hbm92YSRkYXRhJGxvbmckaGVtaXNwaGVyZSkpCmZ0YWJsZShhZGRtYXJnaW5zKG15dGFibGUpKQphcGVyaW9kaWNfcGFyYW1fYWZleF9wbG90IDwtCiAgYWZleF9wbG90KAogICAgYXBlcmlvZGljX3Nsb3BlX3JlcF9hbm92YSwKICAgIHggPSAiR3JvdXAiLAogICAgdHJhY2UgPSAiU2V4IiwKICAgIHBhbmVsID0gImhlbWlzcGhlcmUiLAogICAgZXJyb3IgPSAiYmV0d2VlbiIsCiAgICBlcnJvcl9hcmcgPSBsaXN0KHdpZHRoID0gLjEpLAogICAgZG9kZ2UgPSAuMywKICAgIGRhdGFfYXJnICA9IGxpc3QoCiAgICAgIHBvc2l0aW9uID0gCiAgICAgICAgcG9zaXRpb25faml0dGVyZG9kZ2UoCiAgICAgICAgICBqaXR0ZXIud2lkdGggID0gLjEsIAogICAgICAgICAgaml0dGVyLmhlaWdodCA9IC4xLCAKICAgICAgICAgIGRvZGdlLndpZHRoICAgPSAuMyAgIyMgbmVlZHMgdG8gYmUgc2FtZSBhcyBkb2RnZQogICAgICAgICkpLAogICAgbWFwcGluZyA9IGMoImNvbG9yIiksCiAgICBwb2ludF9hcmcgPSBsaXN0KHNpemUgPSA0KQogICkKc3VwcHJlc3NXYXJuaW5ncyhwcmludChhcGVyaW9kaWNfcGFyYW1fYWZleF9wbG90KSkKbmljZShhcGVyaW9kaWNfc2xvcGVfcmVwX2Fub3ZhKQojIGFfcG9zdGVyaW9yaShhcGVyaW9kaWNfc2xvcGVfcmVwX2Fub3ZhKQplbW1lYW5zKGFwZXJpb2RpY19zbG9wZV9yZXBfYW5vdmEsIHBhaXJ3aXNlIH4gR3JvdXApCmNhdChyZXAoIl8iLCA4MCksICdcbicsIHNlcCA9ICIiKQplbW1lYW5zKGFwZXJpb2RpY19zbG9wZV9yZXBfYW5vdmEsIHBhaXJ3aXNlIH4gU2V4KQpjYXQocmVwKCJfIiwgODApLCAnXG4nLCBzZXAgPSAiIikKZW1tZWFucyhhcGVyaW9kaWNfc2xvcGVfcmVwX2Fub3ZhLCBwYWlyd2lzZSB+IGVsZWN0cm9kZV9wYWlyKQpjYXQocmVwKCJfIiwgODApLCAnXG4nLCBzZXAgPSAiIikKZW1tZWFucyhhcGVyaW9kaWNfc2xvcGVfcmVwX2Fub3ZhLCBwYWlyd2lzZSB+IGVsZWN0cm9kZV9wYWlyIHwgU2V4KQpjYXQocmVwKCJfIiwgODApLCAnXG4nLCBzZXAgPSAiIikKZW1tZWFucyhhcGVyaW9kaWNfc2xvcGVfcmVwX2Fub3ZhLCBwYWlyd2lzZSB+IGhlbWlzcGhlcmUgfCBlbGVjdHJvZGVfcGFpcikKYGBgCg==