10  What other gender disparities do we observe in oral sessions?

Author

Rebecca S. Chen & Petroula Botsidou

Note: this analysis is only part of the supplementary information of the manuscript

What other gender disparities can we observe in oral sessions?

More specifically, we also collected data on whether positive appraisal was given by a questioner, whether a person asked a question without raising their hand or being chosen by the session host to do so (a jumper), whether a speaker talked longer than their allocated time, and if a question was criticizing the speaker. Depending on the question, we investigated the effect of the question asker’s gender, speaker’s gender, host’s gender while correcting for appropriate confounders such as age or talk number.

For all analysis, we only excluded questions that were follow-up questions, and those asked by the host. This means that these data include manipulated and unmanipulated sessions, as we do not expect our manipulation to interfere with what a questioner says exactly.

10.1 Positive appraisal

Positive appraisal was defined as any positive words towards the speaker, but excluding when simply “Thank you for your talk” was said. However, if this comment included a compliment such as “Thank you for your nice talk”, this was counted as positive appraisal.

We used binomial GLMMs to address the effect of the gender and age of the questioner and gender and career stage of the speaker, and corrected for the question number within that Q&A session. We assessed the fit of each variable using a likelihood ratio test and included only the variables that explained significant variation in the final model.

10.1.1 Receiving

# explore the data 
table(data_analysis$compliment, data_analysis$questioner_gender) %>% kbl() %>%
  kable_classic_2()
F M
No 275 325
Yes 216 214
table(data_analysis$compliment, data_analysis$speaker_gender) %>% kbl() %>%
  kable_classic_2()
F M
No 360 235
Yes 274 147
# build initial models

## speaker gender

model_compliment_s_gender <- glmer(compliment ~ speaker_gender + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_gender)), family = "binomial") 

model_compliment_s_gender_null <- glmer(compliment ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_gender)), family = "binomial") 

anova(model_compliment_s_gender, model_compliment_s_gender_null) # not significant
Data: subset(data_analysis, !is.na(speaker_gender))
Models:
model_compliment_s_gender_null: compliment ~ 1 + (1 | session_id/talk_id)
model_compliment_s_gender: compliment ~ speaker_gender + (1 | session_id/talk_id)
                               npar    AIC    BIC  logLik deviance  Chisq Df
model_compliment_s_gender_null    3 1351.1 1365.8 -672.53   1345.1          
model_compliment_s_gender         4 1351.0 1370.6 -671.47   1343.0 2.1136  1
                               Pr(>Chisq)
model_compliment_s_gender_null           
model_compliment_s_gender           0.146
## speaker career stage

model_compliment_s_career <- glmer(compliment ~ speaker_career_short + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_career_short)), family = "binomial") 

model_compliment_s_career_null <- glmer(compliment ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_career_short)), family = "binomial") 

anova(model_compliment_s_career, model_compliment_s_career_null) # almost significant
Data: subset(data_analysis, !is.na(speaker_career_short))
Models:
model_compliment_s_career_null: compliment ~ 1 + (1 | session_id/talk_id)
model_compliment_s_career: compliment ~ speaker_career_short + (1 | session_id/talk_id)
                               npar    AIC    BIC  logLik deviance  Chisq Df
model_compliment_s_career_null    3 1254.6 1269.2 -624.31   1248.6          
model_compliment_s_career         5 1254.0 1278.2 -621.98   1244.0 4.6456  2
                               Pr(>Chisq)  
model_compliment_s_career_null             
model_compliment_s_career           0.098 .
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## final model: nothing significant

10.1.2 Giving

## question nr 

model_compliment_q_nr <- glmer(compliment ~ question_nr + (1|session_id/talk_id), data = subset(data_analysis, !is.na(question_nr)), family = "binomial") 

model_compliment_q_nr_null <- glmer(compliment ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(question_nr)), family = "binomial") 

anova(model_compliment_q_nr, model_compliment_q_nr_null) # significant
Data: subset(data_analysis, !is.na(question_nr))
Models:
model_compliment_q_nr_null: compliment ~ 1 + (1 | session_id/talk_id)
model_compliment_q_nr: compliment ~ question_nr + (1 | session_id/talk_id)
                           npar    AIC    BIC  logLik deviance Chisq Df
model_compliment_q_nr_null    3 1398.7 1413.6 -696.36   1392.7         
model_compliment_q_nr         4 1333.9 1353.8 -662.97   1325.9 66.79  1
                           Pr(>Chisq)    
model_compliment_q_nr_null               
model_compliment_q_nr       3.021e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## questioner gender

model_compliment_q_gender <- glmer(compliment ~ questioner_gender + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_gender)), family = "binomial") 

model_compliment_q_gender_null <- glmer(compliment ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_gender)), family = "binomial") 

anova(model_compliment_q_gender, model_compliment_q_gender_null) # not significant
Data: subset(data_analysis, !is.na(questioner_gender))
Models:
model_compliment_q_gender_null: compliment ~ 1 + (1 | session_id/talk_id)
model_compliment_q_gender: compliment ~ questioner_gender + (1 | session_id/talk_id)
                               npar    AIC    BIC  logLik deviance  Chisq Df
model_compliment_q_gender_null    3 1364.9 1379.7 -679.45   1358.9          
model_compliment_q_gender         4 1365.5 1385.2 -678.73   1357.5 1.4431  1
                               Pr(>Chisq)
model_compliment_q_gender_null           
model_compliment_q_gender          0.2296
## questioner age

model_compliment_q_age <- glmer(compliment ~ questioner_age + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_age)), family = "binomial") 

model_compliment_q_age_null <- glmer(compliment ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_age)), family = "binomial") 

anova(model_compliment_q_age, model_compliment_q_age_null) # significant
Data: subset(data_analysis, !is.na(questioner_age))
Models:
model_compliment_q_age_null: compliment ~ 1 + (1 | session_id/talk_id)
model_compliment_q_age: compliment ~ questioner_age + (1 | session_id/talk_id)
                            npar    AIC    BIC  logLik deviance  Chisq Df
model_compliment_q_age_null    3 1085.8 1099.9 -539.87   1079.8          
model_compliment_q_age         5 1081.9 1105.5 -535.96   1071.9 7.8264  2
                            Pr(>Chisq)  
model_compliment_q_age_null             
model_compliment_q_age         0.01998 *
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
# only question number and age questioner significant (LRT p < 0.05)

# build final model

model_compliment_null <- glmer(compliment ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(question_nr) & !is.na(questioner_age)), family = "binomial") 

model_compliment <- glmer(compliment ~ question_nr + questioner_age + (1|session_id/talk_id), data = subset(data_analysis, !is.na(question_nr) & !is.na(questioner_age)), family = "binomial") 

# use helper function to collect output

model_compliment_out <- collect_out(model = model_compliment, null = model_compliment_null, n_factors = 3, name = "m_compliment", type = "qa", save = "yes", dir = "../results/question-asking/")

model_compliment_out %>% t() %>% kbl() %>%
  kable_classic_2()
model_name m_compliment
AIC 1032.187
n_obs 822
lrt_pval 0
lrt_chisq 59.56
intercept_estimate 0.921
intercept_estimate_prop 0.715
intercept_pval 0
intercept_ci_lower 0.497
intercept_ci_higher 1.344
n_factors 3
est_question_nr -0.365
lowerCI_question_nr -0.473
higherCI_question_nr -0.257
pval_question_nr 0
zval_question_nr -6.607
est_questioner_age2 -0.607
lowerCI_questioner_age2 -0.976
higherCI_questioner_age2 -0.238
pval_questioner_age2 0.001
zval_questioner_age2 -3.222
est_questioner_age3 -0.493
lowerCI_questioner_age3 -0.99
higherCI_questioner_age3 0.004
pval_questioner_age3 0.052
zval_questioner_age3 -1.945

These results show is that the likelihood of a questioner giving a words of positive appraisal is higher for questions asked earlier in the Q&A, and that especially younger audience members (age < 35 years) are more likely to give words of positive appraisal compared to older age categories. Gender of either the question asker or speaker were not significant.

10.2 Jumpers

Next, we asked whether men or women are more likely to ask a question without being allocated to do so (i.e. chosen by the session host to ask your question). Since this might have to do with the perceived ‘authority’ of the session host, we controlled for the age and gender of the session host as well.

We used binomial GLMMs to address the effect of the gender and age of the questioner, and corrected for question number and age and gender of the session host. We assessed the fit of each variable using a likelihood ratio test and included only the variables that explained significant variation in the final model.

# recode NA to 'no jumper' = 0
jumperdata <- data_analysis %>% mutate(jumper = as.factor(case_when(
  is.na(jumper) ~ "0",
  jumper == "1" ~ "1" )))

# explore the data 
table(jumperdata$jumper, jumperdata$questioner_gender) %>% kbl() %>%
  kable_classic_2()
F M
0 486 526
1 5 13
table(jumperdata$jumper, jumperdata$questioner_age) %>% kbl() %>%
  kable_classic_2()
1 2 3
0 391 302 119
1 4 5 1
# build initial models

## question nr

model_jumper_q_nr <- glmer(jumper ~ question_nr + (1|session_id/talk_id), data = subset(jumperdata, !is.na(question_nr)), family = "binomial") 

model_jumper_q_nr_null <- glmer(jumper ~ 1 + (1|session_id/talk_id), data = subset(jumperdata, !is.na(question_nr)), family = "binomial") 

anova(model_jumper_q_nr, model_jumper_q_nr_null) # not significant
Data: subset(jumperdata, !is.na(question_nr))
Models:
model_jumper_q_nr_null: jumper ~ 1 + (1 | session_id/talk_id)
model_jumper_q_nr: jumper ~ question_nr + (1 | session_id/talk_id)
                       npar    AIC    BIC  logLik deviance  Chisq Df Pr(>Chisq)
model_jumper_q_nr_null    3 188.24 203.13 -91.122   182.24                     
model_jumper_q_nr         4 188.12 207.97 -90.062   180.12 2.1197  1     0.1454
## questioner age

model_jumper_q_age <- glmer(jumper ~ questioner_age + (1|session_id/talk_id), data = subset(jumperdata, !is.na(questioner_age)), family = "binomial") 

model_jumper_q_age_null <- glmer(jumper ~ 1 + (1|session_id/talk_id), data = subset(jumperdata, !is.na(questioner_age)), family = "binomial") 

anova(model_jumper_q_age, model_jumper_q_age_null) # not significant
Data: subset(jumperdata, !is.na(questioner_age))
Models:
model_jumper_q_age_null: jumper ~ 1 + (1 | session_id/talk_id)
model_jumper_q_age: jumper ~ questioner_age + (1 | session_id/talk_id)
                        npar    AIC    BIC  logLik deviance Chisq Df Pr(>Chisq)
model_jumper_q_age_null    3 114.06 128.20 -54.030   108.06                    
model_jumper_q_age         5 117.36 140.92 -53.679   107.36 0.702  2      0.704
## questioner gender

model_jumper_q_gender <- glmer(jumper ~ questioner_gender + (1|session_id/talk_id), data = subset(jumperdata, !is.na(questioner_gender)), family = "binomial") 

model_jumper_q_gender_null <- glmer(jumper ~ 1 + (1|session_id/talk_id), data = subset(jumperdata, !is.na(questioner_gender)), family = "binomial") 

anova(model_jumper_q_gender, model_jumper_q_gender_null) # almost significant
Data: subset(jumperdata, !is.na(questioner_gender))
Models:
model_jumper_q_gender_null: jumper ~ 1 + (1 | session_id/talk_id)
model_jumper_q_gender: jumper ~ questioner_gender + (1 | session_id/talk_id)
                           npar    AIC    BIC  logLik deviance  Chisq Df
model_jumper_q_gender_null    3 187.37 202.19 -90.687   181.37          
model_jumper_q_gender         4 186.35 206.10 -89.173   178.35 3.0267  1
                           Pr(>Chisq)  
model_jumper_q_gender_null             
model_jumper_q_gender          0.0819 .
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(model_jumper_q_gender)
Generalized linear mixed model fit by maximum likelihood (Laplace
  Approximation) [glmerMod]
 Family: binomial  ( logit )
Formula: jumper ~ questioner_gender + (1 | session_id/talk_id)
   Data: subset(jumperdata, !is.na(questioner_gender))

     AIC      BIC   logLik deviance df.resid 
   186.3    206.1    -89.2    178.3     1026 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-0.1572 -0.1572 -0.1572 -0.1014  9.8590 

Random effects:
 Groups             Name        Variance Std.Dev.
 talk_id:session_id (Intercept) 0        0       
 session_id         (Intercept) 0        0       
Number of obs: 1030, groups:  talk_id:session_id, 339; session_id, 64

Fixed effects:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)         -4.5768     0.4495 -10.182   <2e-16 ***
questioner_genderM   0.8764     0.5300   1.654   0.0982 .  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Correlation of Fixed Effects:
            (Intr)
qstnr_gndrM -0.848
optimizer (Nelder_Mead) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
## host age

model_jumper_h_age <- glmer(jumper ~ host_1_age + (1|session_id/talk_id), data = subset(jumperdata, !is.na(host_1_age)), family = "binomial") 

model_jumper_h_age_null <- glmer(jumper ~ 1 + (1|session_id/talk_id), data = subset(jumperdata, !is.na(host_1_age)), family = "binomial") 

anova(model_jumper_h_age, model_jumper_h_age_null) # not significant
Data: subset(jumperdata, !is.na(host_1_age))
Models:
model_jumper_h_age_null: jumper ~ 1 + (1 | session_id/talk_id)
model_jumper_h_age: jumper ~ host_1_age + (1 | session_id/talk_id)
                        npar    AIC    BIC  logLik deviance  Chisq Df
model_jumper_h_age_null    3 129.37 143.27 -61.686   123.37          
model_jumper_h_age         5 130.14 153.30 -60.068   120.14 3.2346  2
                        Pr(>Chisq)
model_jumper_h_age_null           
model_jumper_h_age          0.1984
## host gender
 
model_jumper_h_gender <- glmer(jumper ~ host_1_gender + (1|session_id/talk_id), data = subset(jumperdata, !is.na(host_1_gender)), family = "binomial") 

model_jumper_h_gender_null <- glmer(jumper ~ 1 + (1|session_id/talk_id), data = subset(jumperdata, !is.na(host_1_gender)), family = "binomial") 

anova(model_jumper_h_gender, model_jumper_h_gender_null) # not significant
Data: subset(jumperdata, !is.na(host_1_gender))
Models:
model_jumper_h_gender_null: jumper ~ 1 + (1 | session_id/talk_id)
model_jumper_h_gender: jumper ~ host_1_gender + (1 | session_id/talk_id)
                           npar    AIC    BIC  logLik deviance  Chisq Df
model_jumper_h_gender_null    3 169.48 184.15 -81.741   163.48          
model_jumper_h_gender         4 171.16 190.72 -81.579   163.16 0.3243  1
                           Pr(>Chisq)
model_jumper_h_gender_null           
model_jumper_h_gender           0.569
# nothing is significant, trend for gender questioner

These results show is that the likelihood of a person jumping a question not significantly affected by any of the variables. There was however a tendency for men to be more likely to jump a question compared to women. Note that the inference of these models is however limited, since jumpers were rare (N = 18).

10.3 Speaker over time

Next, we investigated the probability that a speaker talks for longer than their allocated speaking time is affected by speaker gender or career stage. We did not expect that any other confounding variables would explain variation in speaking overtime, since this is something that was prepared by only the speaker.

# have to summarize by talk not per question
data_overtime <- data_analysis %>% select(talk_id, session_id, overtime, speaker_gender, speaker_career_short) %>% unique()

# N = 
nrow(data_overtime)
[1] 340
# explore data
table(data_overtime$overtime, data_overtime$speaker_gender) %>% kbl() %>%
  kable_classic_2()
F M
N 181 97
Y 29 21
table(data_overtime$overtime, data_overtime$speaker_career_short) %>% kbl() %>%
  kable_classic_2()
Early career Late career Mid career
N 113 18 123
Y 22 9 14
# build initial models

## speaker gender
model_overtime_s_gender <- glmer(overtime ~ speaker_gender + (1|session_id/talk_id), data = subset(data_overtime, !is.na(speaker_gender)), family = "binomial") 

model_overtime_s_gender_null <- glmer(overtime ~ 1 + (1|session_id/talk_id), data = subset(data_overtime, !is.na(speaker_gender)), family = "binomial") 

anova(model_overtime_s_gender, model_overtime_s_gender_null) # not significant
Data: subset(data_overtime, !is.na(speaker_gender))
Models:
model_overtime_s_gender_null: overtime ~ 1 + (1 | session_id/talk_id)
model_overtime_s_gender: overtime ~ speaker_gender + (1 | session_id/talk_id)
                             npar    AIC    BIC  logLik deviance  Chisq Df
model_overtime_s_gender_null    3 160.29 171.67 -77.145   154.29          
model_overtime_s_gender         4 162.21 177.38 -77.105   154.21 0.0795  1
                             Pr(>Chisq)
model_overtime_s_gender_null           
model_overtime_s_gender           0.778
## speaker career 
model_overtime_s_career <- glmer(overtime ~ speaker_career_short + (1|session_id/talk_id), data = subset(data_overtime, !is.na(speaker_career_short)), family = "binomial") 

model_overtime_s_career_null <- glmer(overtime ~ 1 + (1|session_id/talk_id), data = subset(data_overtime, !is.na(speaker_career_short)), family = "binomial") 

anova(model_overtime_s_career, model_overtime_s_career_null) # significant
Data: subset(data_overtime, !is.na(speaker_career_short))
Models:
model_overtime_s_career_null: overtime ~ 1 + (1 | session_id/talk_id)
model_overtime_s_career: overtime ~ speaker_career_short + (1 | session_id/talk_id)
                             npar    AIC    BIC   logLik deviance  Chisq Df
model_overtime_s_career_null    3 253.24 264.34 -123.621   247.24          
model_overtime_s_career         5 148.76 167.26  -69.381   138.76 108.48  2
                             Pr(>Chisq)    
model_overtime_s_career_null               
model_overtime_s_career       < 2.2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(model_overtime_s_career)
Generalized linear mixed model fit by maximum likelihood (Laplace
  Approximation) [glmerMod]
 Family: binomial  ( logit )
Formula: overtime ~ speaker_career_short + (1 | session_id/talk_id)
   Data: subset(data_overtime, !is.na(speaker_career_short))

     AIC      BIC   logLik deviance df.resid 
   148.8    167.3    -69.4    138.8      294 

Scaled residuals: 
      Min        1Q    Median        3Q       Max 
-0.005965 -0.002171 -0.001696 -0.001311  0.078898 

Random effects:
 Groups             Name        Variance Std.Dev.
 talk_id:session_id (Intercept) 2882.12  53.685  
 session_id         (Intercept)   83.74   9.151  
Number of obs: 299, groups:  talk_id:session_id, 299; session_id, 63

Fixed effects:
                                  Estimate Std. Error z value Pr(>|z|)    
(Intercept)                     -12.700754   0.001747 -7271.0   <2e-16 ***
speaker_career_shortLate career   1.103260   0.001747   631.6   <2e-16 ***
speaker_career_shortMid career   -0.566997   0.001747  -324.6   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Correlation of Fixed Effects:
            (Intr) sp__Lc
spkr_crr_Lc 0.000        
spkr_crr_Mc 0.000  0.000 
optimizer (Nelder_Mead) convergence code: 0 (OK)
Model failed to converge with max|grad| = 0.0404888 (tol = 0.002, component 1)
# use helper function to collect output

model_overtime_out <- collect_out(model = model_overtime_s_career, null = model_overtime_s_career_null, n_factors = 2, name = "m_overtime", type = "qa", save = "yes", dir = "../results/question-asking/")

model_overtime_out %>% t() %>% kbl() %>%
  kable_classic_2()
model_name m_overtime
AIC 148.762
n_obs 299
lrt_pval 0
lrt_chisq 108.48
intercept_estimate -12.701
intercept_estimate_prop 0
intercept_pval 0
intercept_ci_lower -12.704
intercept_ci_higher -12.697
n_factors 2
est_speaker_career_shortLate career 1.103
lowerCI_speaker_career_shortLate career 1.1
higherCI_speaker_career_shortLate career 1.107
pval_speaker_career_shortLate career 0
zval_speaker_career_shortLate career 631.607
est_speaker_career_shortMid career -0.567
lowerCI_speaker_career_shortMid career -0.57
higherCI_speaker_career_shortMid career -0.564
pval_speaker_career_shortMid career 0
zval_speaker_career_shortMid career -324.618

The model output shows that late career speakers were more likely to speak overtime and mid-career speakers were the least likely to speak overtime.

10.4 Critical question

Next, we assessed whether the likelihood of receiving a ‘critical’ question is affected by the gender and age of the questioner or the gender and career stage of the speaker.

10.4.1 Receiving

# explore data
table(data_analysis$question_type_e, data_analysis$questioner_gender) %>% kbl() %>%  kable_classic_2()
F M
0 458 499
1 19 23
table(data_analysis$question_type_e, data_analysis$questioner_age) %>% kbl() %>%  kable_classic_2()
1 2 3
0 381 294 102
1 5 10 13
table(data_analysis$question_type_e, data_analysis$speaker_gender) %>% kbl() %>%
  kable_classic_2()
F M
0 593 349
1 19 23
table(data_analysis$question_type_e, data_analysis$speaker_career_short) %>% kbl() %>%  kable_classic_2()
Early career Late career Mid career
0 434 82 361
1 13 3 18
# build initial models

## speaker gender
model_critical_s_gender <- glmer(question_type_e ~ speaker_gender + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_gender)), family = "binomial") 

model_critical_s_gender_null <- glmer(question_type_e ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_gender)), family = "binomial") 

anova(model_critical_s_gender, model_critical_s_gender_null) # not significant, sort of trend
Data: subset(data_analysis, !is.na(speaker_gender))
Models:
model_critical_s_gender_null: question_type_e ~ 1 + (1 | session_id/talk_id)
model_critical_s_gender: question_type_e ~ speaker_gender + (1 | session_id/talk_id)
                             npar    AIC    BIC  logLik deviance  Chisq Df
model_critical_s_gender_null    3 317.43 332.11 -155.72   311.43          
model_critical_s_gender         4 318.13 337.69 -155.06   310.12 1.3065  1
                             Pr(>Chisq)
model_critical_s_gender_null           
model_critical_s_gender           0.253
## speaker career stage
model_critical_s_career <- glmer(question_type_e ~ speaker_career_short + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_career_short)), family = "binomial") 

model_critical_s_career_null <- glmer(question_type_e ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(speaker_career_short)), family = "binomial") 

anova(model_critical_s_career, model_critical_s_career_null) # not significant
Data: subset(data_analysis, !is.na(speaker_career_short))
Models:
model_critical_s_career_null: question_type_e ~ 1 + (1 | session_id/talk_id)
model_critical_s_career: question_type_e ~ speaker_career_short + (1 | session_id/talk_id)
                             npar    AIC    BIC  logLik deviance  Chisq Df
model_critical_s_career_null    3 263.75 278.20 -128.88   257.75          
model_critical_s_career         5 267.34 291.42 -128.67   257.34 0.4086  2
                             Pr(>Chisq)
model_critical_s_career_null           
model_critical_s_career          0.8152
## nothing significant

10.4.2 Giving

## questioner gender
model_critical_q_gender <- glmer(question_type_e ~ questioner_gender + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_gender)), family = "binomial") 

model_critical_q_gender_null <- glmer(question_type_e ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_gender)), family = "binomial") 

anova(model_critical_q_gender, model_critical_q_gender_null) # not significant
Data: subset(data_analysis, !is.na(questioner_gender))
Models:
model_critical_q_gender_null: question_type_e ~ 1 + (1 | session_id/talk_id)
model_critical_q_gender: question_type_e ~ questioner_gender + (1 | session_id/talk_id)
                             npar    AIC    BIC  logLik deviance Chisq Df
model_critical_q_gender_null    3 317.28 332.00 -155.64   311.28         
model_critical_q_gender         4 318.42 338.04 -155.21   310.42 0.863  1
                             Pr(>Chisq)
model_critical_q_gender_null           
model_critical_q_gender          0.3529
## questioner age
model_critical_q_age <- glmer(question_type_e ~ questioner_age + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_age)), family = "binomial") 

model_critical_q_age_null <- glmer(question_type_e ~ 1 + (1|session_id/talk_id), data = subset(data_analysis, !is.na(questioner_age)), family = "binomial") 

anova(model_critical_q_age, model_critical_q_age_null) # significant, older higher probability
Data: subset(data_analysis, !is.na(questioner_age))
Models:
model_critical_q_age_null: question_type_e ~ 1 + (1 | session_id/talk_id)
model_critical_q_age: question_type_e ~ questioner_age + (1 | session_id/talk_id)
                          npar    AIC    BIC  logLik deviance  Chisq Df
model_critical_q_age_null    3 203.28 217.36 -98.642   197.28          
model_critical_q_age         5 197.07 220.53 -93.536   187.07 10.212  2
                          Pr(>Chisq)   
model_critical_q_age_null              
model_critical_q_age         0.00606 **
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
model_critical_q_age_out <- collect_out(model = model_critical_q_age, null = model_critical_q_age_null, n_factors = 2, name = "m_critical_q_age", type = "qa", save = "yes", dir = "../results/question-asking/")

model_critical_q_age_out %>% t() %>% kbl() %>%  kable_classic_2()
model_name m_critical_q_age
AIC 197.072
n_obs 805
lrt_pval 0.006
lrt_chisq 10.212
intercept_estimate -10.888
intercept_estimate_prop 0
intercept_pval 0
intercept_ci_lower -10.893
intercept_ci_higher -10.883
n_factors 2
est_questioner_age2 2.188
lowerCI_questioner_age2 2.183
higherCI_questioner_age2 2.193
pval_questioner_age2 0
zval_questioner_age2 826.121
est_questioner_age3 3.068
lowerCI_questioner_age3 3.063
higherCI_questioner_age3 3.073
pval_questioner_age3 0
zval_questioner_age3 1158.491
# only age questioner significant

These results show is that the likelihood of a person asking or receiving a critical question was affected by the age category of the question asker. Note that the inference of these models is however limited, since critical questions were rare and subject to observer bias (N = 43).