Prep

Data & libraries

library(tidyverse)
library(psych)
library(ggplot2)
library(lmerTest)
library(lme4)
library(knitr)
library(interactions)
library(patchwork)
library(sjPlot)
library(cowplot)
library(hrbrthemes)
library(ordinal)
library(marginaleffects)

d1 <- read.csv("Sexism_PostElection_Study1comb.csv", header = T)
d2 <- read.csv("Sexism_PostElection_Study2comb.csv", header = T)
d3 <- read.csv("Sexism_PostElection_Study3comb.csv", header = T)

Study 1

VotedFor: If you could have voted in the 2024 U.S. Presidential Election/were an American citizen this week, who would you have voted for?

Var Construction

d1$ASI_3.r <- 5 - d1$ASI_3
d1$ASI_6.r <- 5 - d1$ASI_6
d1$ASI_7.r <- 5 - d1$ASI_7
d1$ASI_13.r <- 5 - d1$ASI_13
d1$ASI_18.r <- 5 - d1$ASI_18
d1$ASI_21.r <- 5 - d1$ASI_21

d1$ASI <- rowMeans(d1[,c("ASI_1", "ASI_2", "ASI_3.r", "ASI_4", "ASI_5", "ASI_6.r", "ASI_7.r", "ASI_8", "ASI_9", "ASI_10", "ASI_11", "ASI_12", "ASI_13.r", "ASI_14", "ASI_15", "ASI_16", "ASI_17", "ASI_18.r", "ASI_19", "ASI_20", "ASI_21.r", "ASI_22")], na.rm = T) # ambivalent sexism inventory
d1$BS <- rowMeans(d1[,c("ASI_1",  "ASI_3.r", "ASI_6.r", "ASI_8", "ASI_9", "ASI_12", "ASI_13.r", "ASI_17", "ASI_19", "ASI_20", "ASI_22")], na.rm = T) # benevolent sexism
d1$HS <- rowMeans(d1[,c("ASI_2",  "ASI_4", "ASI_5", "ASI_7.r", "ASI_10", "ASI_11",   "ASI_14", "ASI_15", "ASI_16", "ASI_18.r", "ASI_21.r")], na.rm = T) # hostile sexism

d1$SRQ_1.r <- 100 - d1$SRQ_1
d1$SRQ_2.r <- 100 - d1$SRQ_2
d1$SRQ_3.r <- 100 - d1$SRQ_3
d1$SRQ_4.r <- 100 - d1$SRQ_4
d1$SRQ_5.r <- 100 - d1$SRQ_5
d1$SRQ <- rowMeans(d1[,c("SRQ_1.r", "SRQ_2.r", "SRQ_3.r", "SRQ_4.r", "SRQ_5.r", "SRQ_6", "SRQ_7", "SRQ_8", "SRQ_9", "SRQ_10", "SRQ_11", "SRQ_12", "SRQ_13")], na.rm = T) # Social roles questionnaire

d1$SDO7_3.r <- 8 - d1$SDO7_3
d1$SDO7_4.r <- 8 - d1$SDO7_4
d1$SDO7_7.r <- 8 - d1$SDO7_7
d1$SDO7_8.r <- 8 - d1$SDO7_8

d1$SDO <- rowMeans(d1[,c("SDO7_1", "SDO7_2", "SDO7_3.r", "SDO7_4.r", "SDO7_5", "SDO7_6", "SDO7_7.r", "SDO7_8.r")], na.rm = T) # social dominance orientation


d1$VotedHarris_num <- ifelse(d1$VotedFor == 1, 1, 0)
d1$VotedHarris <- as.factor(d1$VotedHarris_num)


psych::alpha(d1[,c("Mood", "MentalHealth", "StressAnxiety", "Safety")])
## 
## Reliability analysis   
## Call: psych::alpha(x = d1[, c("Mood", "MentalHealth", "StressAnxiety", 
##     "Safety")])
## 
##   raw_alpha std.alpha G6(smc) average_r S/N    ase mean  sd median_r
##       0.92      0.92     0.9      0.74  11 0.0084  4.4 1.8     0.74
## 
##     95% confidence boundaries 
##          lower alpha upper
## Feldt      0.9  0.92  0.93
## Duhachek   0.9  0.92  0.94
## 
##  Reliability if an item is dropped:
##               raw_alpha std.alpha G6(smc) average_r S/N alpha se   var.r med.r
## Mood               0.91      0.91    0.87      0.76 9.7    0.011 0.00251  0.78
## MentalHealth       0.87      0.87    0.82      0.69 6.8    0.014 0.00022  0.70
## StressAnxiety      0.90      0.90    0.86      0.75 9.1    0.011 0.00231  0.78
## Safety             0.90      0.90    0.87      0.75 9.2    0.011 0.00447  0.78
## 
##  Item statistics 
##                 n raw.r std.r r.cor r.drop mean  sd
## Mood          244  0.87  0.88  0.82   0.78  4.9 1.8
## MentalHealth  244  0.94  0.94  0.92   0.88  4.1 2.0
## StressAnxiety 244  0.90  0.89  0.84   0.80  4.1 2.1
## Safety        244  0.89  0.89  0.83   0.80  4.3 2.0
## 
## Non missing response frequency for each item
##                  1    2    3    4    5    6    7 miss
## Mood          0.07 0.06 0.09 0.16 0.23 0.15 0.25    0
## MentalHealth  0.16 0.11 0.09 0.16 0.21 0.12 0.16    0
## StressAnxiety 0.18 0.11 0.11 0.12 0.16 0.13 0.19    0
## Safety        0.14 0.08 0.12 0.16 0.17 0.11 0.21    0
d1$Impact_index <- rowMeans(d1[,c("Mood", "MentalHealth", "StressAnxiety", "Safety")], na.rm = T)
  • Here & Study 2: Combining “Mood”, “MentalHealth”, “StressAnxiety”, “Safety” into a single impact index b/c of high internal reliability

Setting Classes & centering

# classes
d1$Demo_Age <- as.numeric(d1$Demo_Age) #Participant age
d1$Demo_Gender_num <- as.factor(d1$Demo_Gender) #0 = male, 1 = female
d1$Demo_Gender <- recode_factor(d1$Demo_Gender_num,
                                `0` = "Man",
                                `1` = "Woman") #0 = male, 1 = female
d1$Demo_Ethnicity <- as.factor(d1$Demo_Ethnicity) #Participant ethnicity (0 = non-Hispanic, 1 = Hispanic)
d1$Demo_Race <- as.factor(d1$Demo_Race) #Participant race
d1$Demo_Education <- as.numeric(d1$Demo_Education, na.rm = TRUE)
d1$BS <- as.numeric(d1$BS) #Benevolent sexism
d1$HS <- as.numeric(d1$HS) #Hostile sexism
d1$SRQ <- as.numeric(d1$SRQ) #Traditional gender role beliefs (higher = more traditional)
d1$SDO <- as.numeric(d1$SDO)
d1$CanadianCitizen <- as.factor(d1$CanadianCitizen)
d1$OftenThought <- as.numeric(d1$OftenThought)
d1$InfluenceCanada <- as.numeric(d1$InfluenceCanada)
d1$VotedFor <- as.factor(d1$VotedFor)

# centering
d1$age.c <- d1$Demo_Age - mean(d1$Demo_Age, na.rm = T)
d1$edu.c <- d1$Demo_Education - mean(d1$Demo_Education, na.rm = T)
d1$BS.c <- d1$BS - mean(d1$BS, na.rm = T)
d1$HS.c <- d1$HS - mean(d1$HS, na.rm = T)
d1$SRQ.c <- d1$SRQ - mean(d1$SRQ, na.rm = T)
d1$SDO.c <- d1$SDO - mean(d1$SDO, na.rm = T)
d1$Impact.c <- d1$Impact_index - mean(d1$Impact_index, na.rm = T)
d1$InfluenceCanada.c <- d1$InfluenceCanada - mean(d1$InfluenceCanada, na.rm = T)

Assessing between country metrics, normality, etc.

ASI & factors

describeBy(d1$ASI, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 121 1.62 0.86    1.5    1.59 1.01   0 3.68  3.68 0.26     -0.9 0.08
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad  min  max range  skew kurtosis   se
## X1    1 123 1.83 0.88   1.91    1.84 1.08 0.05 3.68  3.64 -0.12    -0.94 0.08
d1 %>%
  ggplot( aes(x=ASI, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Ambivalent Sexism Inventory") +
  facet_grid(Country~.)

describeBy(d1$BS, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 121 1.94 0.93   1.91     1.9 0.94   0 4.45  4.45 0.31    -0.32 0.08
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min  max range  skew kurtosis   se
## X1    1 123 2.18 1.03   2.18     2.2 1.21   0 4.09  4.09 -0.17    -0.84 0.09
d1 %>%
  ggplot( aes(x=BS, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Benevolent Sexism")+
  facet_grid(Country~.)

describeBy(d1$HS, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis  se
## X1    1 121  1.3 1.09   1.09    1.17 1.08   0 4.36  4.36 0.85    -0.19 0.1
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis  se
## X1    1 123 1.48 1.13   1.27    1.42 1.35   0 4.73  4.73 0.46    -0.77 0.1
d1 %>%
  ggplot( aes(x=HS, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Hostile Sexism")+
  facet_grid(Country~.)

SRQ

describeBy(d1$SRQ, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n  mean    sd median trimmed   mad min   max range skew kurtosis   se
## X1    1 121 26.15 15.81  23.08   25.28 17.11   0 64.62 64.62 0.43    -0.67 1.44
## ------------------------------------------------------------ 
## group: US
##    vars   n  mean    sd median trimmed   mad min   max range skew kurtosis   se
## X1    1 123 31.34 17.47     30   31.22 21.67   0 72.31 72.31 0.11    -1.01 1.58
d1 %>%
  ggplot( aes(x=SRQ, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 5) +
#  coord_cartesian(xlim = c(0,100)) +
  scale_x_continuous(breaks = seq(0,70,5)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Social Roles Questionnaire")+
  facet_grid(Country~.)

SDO

describeBy(d1$SDO, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 121 2.34 1.17      2     2.2 1.11   1 5.75  4.75    1     0.43 0.11
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 123 2.43 1.26   2.12     2.3 1.48   1 6.62  5.62 0.79    -0.02 0.11
d1 %>%
  ggplot( aes(x=SDO, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .5) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Social Dominance Orientation")+
  facet_grid(Country~.)

Personal Impact Measures

Mood = “How much has the U.S. election result impacted your overall mood this week?”

OftenThought = “How often have you thought about the U.S. election results in the past few days?”

StressAnxiety = “How much stress or anxiety do you feel regarding the outcome of the U.S. election?”

Safety = “Do you feel the U.S. election results have affected your sense of safety or well-being?”

MentalHealth = “How have the election results affected your overall mental health this week?”

describeBy(d1$Mood, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 121 4.64 1.71      5    4.76 1.48   1   7     6 -0.53    -0.46 0.16
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 123 5.11 1.87      5    5.33 2.97   1   7     6 -0.68    -0.64 0.17
d1 %>%
  ggplot( aes(x=Mood, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Mood")+
  facet_grid(Country~.)

describeBy(d1$StressAnxiety, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 121 3.99 2.06      4    3.99 2.97   1   7     6 -0.11    -1.31 0.19
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis  se
## X1    1 123 4.26 2.23      4    4.32 2.97   1   7     6 -0.17    -1.44 0.2
d1 %>%
  ggplot( aes(x=StressAnxiety, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Stress & Anxiety")+
  facet_grid(Country~.)

describeBy(d1$Safety, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 121 4.16 1.89      4     4.2 1.48   1   7     6 -0.23    -0.98 0.17
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 123  4.5 2.16      5    4.63 2.97   1   7     6 -0.27    -1.32 0.19
d1 %>%
  ggplot( aes(x=Safety, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Stress & Anxiety")+
  facet_grid(Country~.)

describeBy(d1$MentalHealth, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 121  3.9 1.93      4    3.88 1.48   1   7     6 -0.15     -1.1 0.18
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 123 4.39 2.07      5    4.48 2.97   1   7     6 -0.3    -1.25 0.19
d1 %>%
  ggplot( aes(x=MentalHealth, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Mental Health")+
  facet_grid(Country~.)

describeBy(d1$OftenThought, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean  sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 121 5.48 1.4      6    5.66 1.48   2   7     5 -0.83    -0.04 0.13
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 123  5.6 1.61      6    5.83 1.48   1   7     6 -0.93    -0.33 0.15
d1 %>%
  ggplot( aes(x=OftenThought, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "How often have you thought about the U.S. election this week?")+
  facet_grid(Country~.)

describeBy(d1$MentalHealth, d1$Country)
## 
##  Descriptive statistics by group 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 121  3.9 1.93      4    3.88 1.48   1   7     6 -0.15     -1.1 0.18
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 123 4.39 2.07      5    4.48 2.97   1   7     6 -0.3    -1.25 0.19
d1 %>%
  ggplot( aes(x=MentalHealth, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Mental Health")+
  facet_grid(Country~.)

Study 2

VotedFor: If you could have voted in the 2024 U.S. Presidential Election/were an American citizen this week, who would you have voted for?

Var Construction

d2$ASI_3.r <- 5 - d2$ASI_3
d2$ASI_6.r <- 5 - d2$ASI_6
d2$ASI_7.r <- 5 - d2$ASI_7
d2$ASI_13.r <- 5 - d2$ASI_13
d2$ASI_18.r <- 5 - d2$ASI_18
d2$ASI_21.r <- 5 - d2$ASI_21

d2$ASI <- rowMeans(d2[,c("ASI_1", "ASI_2", "ASI_3.r", "ASI_4", "ASI_5", "ASI_6.r", "ASI_7.r", "ASI_8", "ASI_9", "ASI_10", "ASI_11", "ASI_12", "ASI_13.r", "ASI_14", "ASI_15", "ASI_16", "ASI_17", "ASI_18.r", "ASI_19", "ASI_20", "ASI_21.r", "ASI_22")], na.rm = T) # ambivalent sexism inventory
d2$BS <- rowMeans(d2[,c("ASI_1",  "ASI_3.r", "ASI_6.r", "ASI_8", "ASI_9", "ASI_12", "ASI_13.r", "ASI_17", "ASI_19", "ASI_20", "ASI_22")], na.rm = T) # benevolent sexism
d2$HS <- rowMeans(d2[,c("ASI_2",  "ASI_4", "ASI_5", "ASI_7.r", "ASI_10", "ASI_11",   "ASI_14", "ASI_15", "ASI_16", "ASI_18.r", "ASI_21.r")], na.rm = T) # hostile sexism

d2$SRQ_1.r <- 100 - d2$SRQ_1
d2$SRQ_2.r <- 100 - d2$SRQ_2
d2$SRQ_3.r <- 100 - d2$SRQ_3
d2$SRQ_4.r <- 100 - d2$SRQ_4
d2$SRQ_5.r <- 100 - d2$SRQ_5
d2$SRQ <- rowMeans(d2[,c("SRQ_1.r", "SRQ_2.r", "SRQ_3.r", "SRQ_4.r", "SRQ_5.r", "SRQ_6", "SRQ_7", "SRQ_8", "SRQ_9", "SRQ_10", "SRQ_11", "SRQ_12", "SRQ_13")], na.rm = T) # Social roles questionnaire

d2$SDO7_3.r <- 8 - d2$SDO7_3
d2$SDO7_4.r <- 8 - d2$SDO7_4
d2$SDO7_7.r <- 8 - d2$SDO7_7
d2$SDO7_8.r <- 8 - d2$SDO7_8

d2$SDO <- rowMeans(d2[,c("SDO7_1", "SDO7_2", "SDO7_3.r", "SDO7_4.r", "SDO7_5", "SDO7_6", "SDO7_7.r", "SDO7_8.r")], na.rm = T) # social dominance orientation

d2$VotedHarris_num <- ifelse(d2$VotedFor == 1, 1, 0)
d2$VotedHarris <- as.factor(d2$VotedHarris_num)

psych::alpha(d2[,c("Mood", "StressAnxiety", "Safety")])
## 
## Reliability analysis   
## Call: psych::alpha(x = d2[, c("Mood", "StressAnxiety", "Safety")])
## 
##   raw_alpha std.alpha G6(smc) average_r S/N   ase mean  sd median_r
##       0.84      0.84    0.79      0.63 5.1 0.012  3.7 1.7     0.57
## 
##     95% confidence boundaries 
##          lower alpha upper
## Feldt     0.81  0.84  0.86
## Duhachek  0.81  0.84  0.86
## 
##  Reliability if an item is dropped:
##               raw_alpha std.alpha G6(smc) average_r S/N alpha se var.r med.r
## Mood               0.86      0.86    0.76      0.76 6.3    0.011    NA  0.76
## StressAnxiety      0.72      0.72    0.56      0.56 2.5    0.023    NA  0.56
## Safety             0.73      0.73    0.57      0.57 2.7    0.022    NA  0.57
## 
##  Item statistics 
##                 n raw.r std.r r.cor r.drop mean  sd
## Mood          593  0.81  0.82  0.65   0.60  3.9 1.9
## StressAnxiety 593  0.90  0.89  0.84   0.76  3.6 2.0
## Safety        593  0.89  0.89  0.83   0.75  3.7 2.0
## 
## Non missing response frequency for each item
##                  1    2    3    4    5    6    7 miss
## Mood          0.16 0.13 0.11 0.18 0.23 0.09 0.10    0
## StressAnxiety 0.22 0.14 0.13 0.16 0.16 0.11 0.09    0
## Safety        0.20 0.15 0.08 0.15 0.20 0.12 0.09    0
d2$Impact_index <- rowMeans(d2[,c("Mood", "StressAnxiety", "Safety")], na.rm = T)


d2$ideology <- rowMeans(d2[,grep("Ideology", colnames(d2))], na.rm = T)
d2$Ideology <- 8 - d2$ideology

describe(d2[,c("Country","Ideology","Ideology_1","Ideology_2","Ideology_3","Ideology_4")])
##            vars   n mean   sd median trimmed  mad min max range  skew kurtosis
## Country*      1 593 2.50 1.12   3.00    2.50 1.48   1   4     3  0.00    -1.36
## Ideology      2 593 3.36 1.46   3.25    3.28 1.11   1   7     6  0.44    -0.16
## Ideology_1    3 593 4.51 1.64   5.00    4.57 1.48   1   7     6 -0.33    -0.62
## Ideology_2    4 593 4.89 1.60   5.00    5.02 1.48   1   7     6 -0.58    -0.40
## Ideology_3    5 593 4.45 1.58   4.00    4.49 1.48   1   7     6 -0.31    -0.45
## Ideology_4    6 593 4.74 1.50   5.00    4.83 1.48   1   7     6 -0.48    -0.16
##              se
## Country*   0.05
## Ideology   0.06
## Ideology_1 0.07
## Ideology_2 0.07
## Ideology_3 0.06
## Ideology_4 0.06

Setting Classes

# classes
d2$Demo_Age <- as.numeric(d2$Demo_Age) #Participant age
d2$Demo_Gender_num <- as.factor(d2$Demo_Gender) #0 = male, 1 = female
d2$Demo_Gender <- recode_factor(d2$Demo_Gender_num,
                                `0` = "Man",
                                `1` = "Woman") #0 = male, 1 = female
d2$Demo_Ethnicity <- as.factor(d2$Demo_Ethnicity) #Participant ethnicity (0 = non-Hispanic, 1 = Hispanic)
d2$Demo_Race <- as.factor(d2$Demo_Race) #Participant race
d2$Demo_Education <- as.numeric(d2$Demo_Education, na.rm = TRUE)
d2$BS <- as.numeric(d2$BS) #Benevolent sexism
d2$HS <- as.numeric(d2$HS) #Hostile sexism
d2$SRQ <- as.numeric(d2$SRQ) #Traditional gender role beliefs (higher = more traditional)
d2$SDO <- as.numeric(d2$SDO)
d2$CanadianCitizen <- as.factor(d2$CanadianCitizen)
d2$OftenThought <- as.numeric(d2$OftenThought)
d2$InfluenceCanada <- as.numeric(d2$InfluenceCanada)
d2$VotedFor <- as.factor(d2$VotedFor)
d2$Country <- as.factor(d2$Country)

# centering
d2$age.c <- d2$Demo_Age - mean(d2$Demo_Age, na.rm = T)
d2$edu.c <- d2$Demo_Education - mean(d2$Demo_Education, na.rm = T)
d2$BS.c <- d2$BS - mean(d2$BS, na.rm = T)
d2$HS.c <- d2$HS - mean(d2$HS, na.rm = T)
d2$SRQ.c <- d2$SRQ - mean(d2$SRQ, na.rm = T)
d2$SDO.c <- d2$SDO - mean(d2$SDO, na.rm = T)
d2$Impact.c <- d2$Impact_index - mean(d2$Impact_index, na.rm = T)
d2$InfluenceCanada.c <- d2$InfluenceCanada - mean(d2$InfluenceCanada, na.rm = T)
d2$Ideology.c <- d2$Ideology - 4 # meaningful midpoint--midpoint centering

Assessing between country metrics, normality, etc.

ASI & factors

describeBy(d2$ASI, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad  min  max range skew kurtosis   se
## X1    1 147 1.77 0.94   1.82    1.75 1.01 0.05 4.05     4 0.19    -0.72 0.08
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range  skew kurtosis   se
## X1    1 149 1.73 0.98   1.86    1.73 1.21   0 3.73  3.73 -0.07    -1.11 0.08
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min  max range  skew kurtosis   se
## X1    1 150 2.01 0.95   2.07    2.03 1.11   0 4.09  4.09 -0.15    -0.86 0.08
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 147 1.89 0.93   1.86    1.89 1.01   0 4.45  4.45 0.01    -0.54 0.08
d2 %>%
  ggplot( aes(x=ASI, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Ambivalent Sexism Inventory") +
  facet_wrap(~Country)

describeBy(d2$BS, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 147 2.03 0.99   2.18    2.02 0.94   0 4.73  4.73 0.04    -0.55 0.08
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 149 1.92 1.07      2    1.91 1.21   0 4.73  4.73 0.05    -0.73 0.09
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min  max range  skew kurtosis   se
## X1    1 150 2.25 0.97   2.32    2.27 0.88   0 4.45  4.45 -0.14    -0.24 0.08
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 147 2.19 1.09   2.27    2.19 1.08   0   5     5 0.05    -0.26 0.09
d2 %>%
  ggplot( aes(x=BS, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Benevolent Sexism")+
  facet_wrap(~Country)

describeBy(d2$HS, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 147 1.51 1.13   1.45    1.45 1.35   0 4.36  4.36 0.39    -0.95 0.09
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 149 1.53 1.16   1.45    1.47 1.48   0 4.36  4.36 0.29    -1.03 0.09
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis  se
## X1    1 150 1.76 1.17   1.64    1.71 1.35   0 4.64  4.64 0.31    -0.92 0.1
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis  se
## X1    1 147 1.58 1.17   1.45     1.5 1.48   0 4.18  4.18 0.42    -0.87 0.1
d2 %>%
  ggplot( aes(x=HS, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Hostile Sexism") +
  facet_wrap(~Country)

SRQ

describeBy(d2$SRQ, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n  mean    sd median trimmed   mad min max range skew kurtosis   se
## X1    1 147 27.26 16.33  26.15   26.48 17.11   0  70    70 0.39    -0.65 1.35
## ------------------------------------------------------------ 
## group: Canada
##    vars   n  mean    sd median trimmed   mad min   max range skew kurtosis   se
## X1    1 149 25.74 17.82  25.38    25.2 23.95   0 66.92 66.92 0.22    -1.14 1.46
## ------------------------------------------------------------ 
## group: UK
##    vars   n  mean    sd median trimmed   mad min   max range skew kurtosis   se
## X1    1 150 31.13 16.05     30   31.04 19.39   0 66.92 66.92 0.09    -1.03 1.31
## ------------------------------------------------------------ 
## group: US
##    vars   n  mean    sd median trimmed   mad min   max range skew kurtosis   se
## X1    1 147 31.21 18.62  30.77   31.03 23.95   0 76.15 76.15 0.07    -1.05 1.54
d2 %>%
  ggplot( aes(x=SRQ, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 5) +
#  coord_cartesian(xlim = c(0,100)) +
  scale_x_continuous(breaks = seq(0,70,5)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Social Roles Questionnaire")+
  facet_wrap(~Country)

SDO

describeBy(d2$SDO, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 147  2.6 1.12    2.5    2.53 1.11   1 6.25  5.25 0.59     0.06 0.09
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 149 2.51 1.25   2.38    2.41 1.48   1   7     6 0.81     0.59 0.1
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed mad min  max range skew kurtosis  se
## X1    1 150 2.77 1.17   2.88    2.74 1.3   1 5.75  4.75 0.16    -0.79 0.1
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 147 2.44 1.39   2.25     2.3 1.85   1 5.88  4.88 0.57    -0.85 0.11
d2 %>%
  ggplot( aes(x=SDO, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .5) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Social Dominance Orientation")+
  facet_wrap(~Country)

Personal Impact Measures

Mood = “How much has the U.S. election result impacted your overall mood this week?”

StressAnxiety = “How much stress or anxiety do you feel regarding the outcome of the U.S. election?”

Safety = “Do you feel the U.S. election results have affected your sense of safety or well-being?”

describeBy(d2$Mood, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 147 3.91 1.78      4    3.89 1.48   1   7     6 -0.09    -0.85 0.15
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 149 3.79 1.82      4    3.77 1.48   1   7     6 -0.06    -1.12 0.15
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 150 3.53 1.82      4    3.46 1.48   1   7     6  0.1    -1.05 0.15
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 147 4.28 2.07      5    4.34 2.97   1   7     6 -0.29    -1.21 0.17
d2 %>%
  ggplot( aes(x=Mood, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Mood") +
  facet_wrap(~Country)

describeBy(d2$StressAnxiety, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 147 3.52 1.79      4    3.46 1.48   1   7     6 0.13    -1.02 0.15
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 149 3.66 1.92      4    3.61 2.97   1   7     6 0.06    -1.25 0.16
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 150 3.43 1.85      3    3.34 2.97   1   7     6 0.21     -1.1 0.15
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 147 3.73 2.28      4    3.66 2.97   1   7     6 0.14    -1.51 0.19
d2 %>%
  ggplot( aes(x=StressAnxiety, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Stress & Anxiety")+
  facet_wrap(~Country)

describeBy(d2$Safety, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 147 3.73 1.86      4    3.71 2.97   1   7     6 -0.04    -1.19 0.15
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 149 3.62 1.95      4    3.57 2.97   1   7     6 0.05    -1.31 0.16
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 150 3.53 1.86      4    3.46 2.97   1   7     6 0.15    -1.18 0.15
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 147 4.06 2.25      4    4.08 2.97   1   7     6 -0.15    -1.47 0.19
d2 %>%
  ggplot( aes(x=Safety, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Stress & Anxiety")+
  facet_wrap(~Country)

describeBy(d2$OftenThought, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 147  4.5 1.69      5    4.56 1.48   1   7     6 -0.33     -0.8 0.14
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 149 4.34 1.59      4    4.32 1.48   1   7     6 -0.06    -0.92 0.13
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean  sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 150 4.03 1.8      4    4.03 1.48   1   7     6 -0.09    -0.96 0.15
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 147 4.59 1.86      5    4.72 1.48   1   7     6 -0.45    -0.81 0.15
d2 %>%
  ggplot( aes(x=OftenThought, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 1) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "How often have you thought about the U.S. election this week?")+
  facet_wrap(~Country)

Ideology

describeBy(d2$Ideology, d2$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 147 3.12 1.17      3    3.08 1.48   1   7     6 0.41     0.48 0.1
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 149 3.36 1.49      3    3.29 1.48   1   7     6 0.48    -0.29 0.12
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 150 3.41 1.26    3.5    3.35 0.93   1   7     6 0.39     0.27 0.1
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 147 3.54 1.81   3.75    3.46 2.22   1   7     6 0.23       -1 0.15
d2 %>%
  ggplot( aes(x=Ideology, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .5) +
  scale_x_continuous(breaks = seq(1,7,1)) +
  coord_cartesian(xlim = c(1,7)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Ideology (Liberal to Conservative)")+
  facet_wrap(~Country)

Study 3

VotedFor: If you could have voted in the 2024 U.S. Presidential Election/were an American citizen this week, who would you have voted for?

Var Construction

d3$ASI_3.r <- 5 - d3$ASI_3
d3$ASI_6.r <- 5 - d3$ASI_6
d3$ASI_7.r <- 5 - d3$ASI_7
d3$ASI_13.r <- 5 - d3$ASI_13
d3$ASI_18.r <- 5 - d3$ASI_18
d3$ASI_21.r <- 5 - d3$ASI_21

d3$ASI <- rowMeans(d3[,c("ASI_1", "ASI_2", "ASI_3.r", "ASI_4", "ASI_5", "ASI_6.r", "ASI_7.r", "ASI_8", "ASI_9", "ASI_10", "ASI_11", "ASI_12", "ASI_13.r", "ASI_14", "ASI_15", "ASI_16", "ASI_17", "ASI_18.r", "ASI_19", "ASI_20", "ASI_21.r", "ASI_22")], na.rm = T) # ambivalent sexism inventory
d3$BS <- rowMeans(d3[,c("ASI_1",  "ASI_3.r", "ASI_6.r", "ASI_8", "ASI_9", "ASI_12", "ASI_13.r", "ASI_17", "ASI_19", "ASI_20", "ASI_22")], na.rm = T) # benevolent sexism
d3$HS <- rowMeans(d3[,c("ASI_2",  "ASI_4", "ASI_5", "ASI_7.r", "ASI_10", "ASI_11",   "ASI_14", "ASI_15", "ASI_16", "ASI_18.r", "ASI_21.r")], na.rm = T) # hostile sexism

d3$SRQ_1.r <- 100 - d3$SRQ_1
d3$SRQ_2.r <- 100 - d3$SRQ_2
d3$SRQ_3.r <- 100 - d3$SRQ_3
d3$SRQ_4.r <- 100 - d3$SRQ_4
d3$SRQ_5.r <- 100 - d3$SRQ_5
d3$SRQ <- rowMeans(d3[,c("SRQ_1.r", "SRQ_2.r", "SRQ_3.r", "SRQ_4.r", "SRQ_5.r", "SRQ_6", "SRQ_7", "SRQ_8", "SRQ_9", "SRQ_10", "SRQ_11", "SRQ_12", "SRQ_13")], na.rm = T) # Social roles questionnaire

d3$SDO7_3.r <- 8 - d3$SDO7_3
d3$SDO7_4.r <- 8 - d3$SDO7_4
d3$SDO7_7.r <- 8 - d3$SDO7_7
d3$SDO7_8.r <- 8 - d3$SDO7_8

d3$SDO <- rowMeans(d3[,c("SDO7_1", "SDO7_2", "SDO7_3.r", "SDO7_4.r", "SDO7_5", "SDO7_6", "SDO7_7.r", "SDO7_8.r")], na.rm = T) # social dominance orientation


d3$VotedHarris_num <- ifelse(d3$VotedFor == 1, 1, 0)
d3$VotedHarris <- as.factor(d3$VotedHarris_num)


d3$ideology <- rowMeans(d3[,grep("Ideology", colnames(d3))], na.rm = T)
d3$Ideology <- 8 - d3$ideology

describe(d3[,c("Country","Ideology","Ideology_1","Ideology_2","Ideology_3","Ideology_4")])
##            vars   n mean   sd median trimmed  mad min max range  skew kurtosis
## Country*      1 995 2.50 1.12   3.00    2.50 1.48   1   4     3  0.00    -1.36
## Ideology      2 995 3.49 1.55   3.25    3.43 1.85   1   7     6  0.32    -0.59
## Ideology_1    3 995 4.36 1.71   4.00    4.39 1.48   1   7     6 -0.23    -0.88
## Ideology_2    4 995 4.79 1.72   5.00    4.92 1.48   1   7     6 -0.54    -0.63
## Ideology_3    5 995 4.32 1.66   4.00    4.35 1.48   1   7     6 -0.19    -0.77
## Ideology_4    6 995 4.59 1.59   5.00    4.66 1.48   1   7     6 -0.38    -0.63
##              se
## Country*   0.04
## Ideology   0.05
## Ideology_1 0.05
## Ideology_2 0.05
## Ideology_3 0.05
## Ideology_4 0.05
psych::alpha(d3[,c("WomanTraits_1", "WomanTraits_2", "WomanTraits_3", "WomanTraits_4", "WomanTraits_5")])
## 
## Reliability analysis   
## Call: psych::alpha(x = d3[, c("WomanTraits_1", "WomanTraits_2", "WomanTraits_3", 
##     "WomanTraits_4", "WomanTraits_5")])
## 
##   raw_alpha std.alpha G6(smc) average_r S/N   ase mean   sd median_r
##       0.77      0.78    0.75      0.41 3.5 0.011  6.2 0.67      0.4
## 
##     95% confidence boundaries 
##          lower alpha upper
## Feldt     0.75  0.77   0.8
## Duhachek  0.75  0.77   0.8
## 
##  Reliability if an item is dropped:
##               raw_alpha std.alpha G6(smc) average_r S/N alpha se  var.r med.r
## WomanTraits_1      0.78      0.78    0.74      0.48 3.6    0.012 0.0074  0.46
## WomanTraits_2      0.75      0.75    0.71      0.43 3.0    0.012 0.0195  0.42
## WomanTraits_3      0.70      0.71    0.65      0.38 2.4    0.015 0.0056  0.36
## WomanTraits_4      0.71      0.72    0.68      0.39 2.6    0.014 0.0122  0.36
## WomanTraits_5      0.71      0.71    0.66      0.38 2.5    0.014 0.0099  0.37
## 
##  Item statistics 
##                 n raw.r std.r r.cor r.drop mean   sd
## WomanTraits_1 994  0.56  0.62  0.45   0.40  6.7 0.67
## WomanTraits_2 994  0.73  0.70  0.57   0.51  6.0 1.08
## WomanTraits_3 994  0.79  0.78  0.73   0.63  6.1 0.98
## WomanTraits_4 994  0.78  0.76  0.68   0.60  5.9 1.03
## WomanTraits_5 994  0.76  0.78  0.71   0.63  6.4 0.80
## 
## Non missing response frequency for each item
##               1    2    3    4    5    6    7 miss
## WomanTraits_1 0 0.00 0.00 0.01 0.04 0.19 0.75    0
## WomanTraits_2 0 0.00 0.02 0.07 0.16 0.35 0.40    0
## WomanTraits_3 0 0.01 0.01 0.04 0.18 0.37 0.39    0
## WomanTraits_4 0 0.00 0.01 0.08 0.21 0.34 0.36    0
## WomanTraits_5 0 0.00 0.00 0.02 0.11 0.34 0.53    0
psych::alpha(d3[,c("WomanTraits_6", "WomanTraits_7", "WomanTraits_8", "WomanTraits_9", "WomanTraits_10")])
## 
## Reliability analysis   
## Call: psych::alpha(x = d3[, c("WomanTraits_6", "WomanTraits_7", "WomanTraits_8", 
##     "WomanTraits_9", "WomanTraits_10")])
## 
##   raw_alpha std.alpha G6(smc) average_r S/N    ase mean  sd median_r
##       0.91      0.91     0.9      0.67  10 0.0045  5.4 1.1     0.67
## 
##     95% confidence boundaries 
##          lower alpha upper
## Feldt      0.9  0.91  0.92
## Duhachek   0.9  0.91  0.92
## 
##  Reliability if an item is dropped:
##                raw_alpha std.alpha G6(smc) average_r S/N alpha se  var.r med.r
## WomanTraits_6       0.89      0.89    0.87      0.67 8.1   0.0058 0.0086  0.68
## WomanTraits_7       0.89      0.89    0.87      0.67 8.2   0.0057 0.0046  0.67
## WomanTraits_8       0.91      0.91    0.88      0.71 9.7   0.0049 0.0026  0.71
## WomanTraits_9       0.88      0.88    0.87      0.65 7.6   0.0061 0.0056  0.65
## WomanTraits_10      0.89      0.89    0.88      0.67 8.2   0.0057 0.0062  0.67
## 
##  Item statistics 
##                  n raw.r std.r r.cor r.drop mean  sd
## WomanTraits_6  994  0.87  0.87  0.83   0.79  5.2 1.3
## WomanTraits_7  994  0.86  0.86  0.83   0.78  5.6 1.3
## WomanTraits_8  994  0.82  0.81  0.75   0.71  5.4 1.3
## WomanTraits_9  994  0.89  0.89  0.86   0.82  5.5 1.3
## WomanTraits_10 994  0.86  0.86  0.82   0.78  5.4 1.3
## 
## Non missing response frequency for each item
##                   1    2    3    4    5    6    7 miss
## WomanTraits_6  0.01 0.03 0.06 0.18 0.31 0.23 0.19    0
## WomanTraits_7  0.01 0.02 0.03 0.12 0.25 0.28 0.29    0
## WomanTraits_8  0.01 0.02 0.04 0.17 0.27 0.28 0.22    0
## WomanTraits_9  0.01 0.02 0.04 0.13 0.25 0.30 0.26    0
## WomanTraits_10 0.01 0.03 0.03 0.15 0.27 0.28 0.25    0
d3$WomanTraits_competence <- rowMeans(d3[,c("WomanTraits_1", "WomanTraits_2", "WomanTraits_3", "WomanTraits_4", "WomanTraits_5")], na.rm = T)

d3$WomanTraits_warmth <- rowMeans(d3[,c("WomanTraits_6", "WomanTraits_7", "WomanTraits_8", "WomanTraits_9", "WomanTraits_10")], na.rm = T)

Setting Classes & Centering

# classes
d3$Demo_Age <- as.numeric(d3$Demo_Age) #Participant age
d3$Demo_Gender_num <- as.factor(d3$Demo_Gender) #0 = male, 1 = female
d3$Demo_Gender <- recode_factor(d3$Demo_Gender_num,
                                `0` = "Man",
                                `1` = "Woman") #0 = male, 1 = female
d3$Demo_Ethnicity <- as.factor(d3$Demo_Ethnicity) #Participant ethnicity (0 = non-Hispanic, 1 = Hispanic)
d3$Demo_Race <- as.factor(d3$Demo_Race) #Participant race
d3$Demo_Education <- as.numeric(d3$Demo_Education, na.rm = TRUE)
d3$BS <- as.numeric(d3$BS) #Benevolent sexism
d3$HS <- as.numeric(d3$HS) #Hostile sexism
d3$SRQ <- as.numeric(d3$SRQ) #Traditional gender role beliefs (higher = more traditional)
d3$SDO <- as.numeric(d3$SDO)
d3$Citizen <- as.factor(d3$Citizen)
d3$VotedFor <- as.factor(d3$VotedFor)
d3$Country <- as.factor(d3$Country)

# centering
d3$age.c <- d3$Demo_Age - mean(d3$Demo_Age, na.rm = T)
d3$edu.c <- d3$Demo_Education - mean(d3$Demo_Education, na.rm = T)
d3$BS.c <- d3$BS - mean(d3$BS, na.rm = T)
d3$HS.c <- d3$HS - mean(d3$HS, na.rm = T)
d3$SRQ.c <- d3$SRQ - mean(d3$SRQ, na.rm = T)
d3$SDO.c <- d3$SDO - mean(d3$SDO, na.rm = T)
d3$Ideology.c <- d3$Ideology - 4 # meaningful midpoint--midpoint centering
d3$warmth.c <- d3$WomanTraits_warmth - mean(d3$WomanTraits_warmth, na.rm = T)
d3$competence.c <- d3$WomanTraits_competence - mean(d3$WomanTraits_competence, na.rm = T)

Assessing between country metrics, normality, etc.

ASI & factors

describeBy(d3$ASI, d3$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean  sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 247  1.7 0.9   1.68    1.67 1.01   0 3.95  3.95 0.27    -0.56 0.06
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range  skew kurtosis   se
## X1    1 250 1.85 0.86   1.86    1.86 0.98   0 4.05  4.05 -0.01    -0.71 0.05
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 249 1.87 0.91      2    1.88 1.01   0   4     4 -0.15    -0.73 0.06
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 249 1.71 1.05   1.64    1.68 1.21   0 4.45  4.45 0.25    -0.84 0.07
d3 %>%
  ggplot( aes(x=ASI, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Ambivalent Sexism Inventory") +
  facet_wrap(~Country)

describeBy(d3$BS, d3$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 247    2 0.94      2    1.99 0.94   0 4.45  4.45 0.09    -0.45 0.06
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min  max range  skew kurtosis   se
## X1    1 250 2.09 0.92   2.09     2.1 0.94   0 4.64  4.64 -0.08     -0.3 0.06
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean sd median trimmed  mad min  max range  skew kurtosis   se
## X1    1 249 2.06  1   2.09    2.08 1.08   0 4.55  4.55 -0.11    -0.57 0.06
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 249 1.94 1.13      2    1.93 1.35   0   5     5 0.08    -0.76 0.07
d3 %>%
  ggplot( aes(x=BS, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Benevolent Sexism")+
  facet_wrap(~Country)

describeBy(d3$HS, d3$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 247  1.4 1.12   1.27     1.3 1.21   0   5     5 0.71    -0.05 0.07
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 250 1.62 1.14   1.55    1.54 1.21   0   5     5  0.5    -0.29 0.07
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 249 1.67 1.12   1.64    1.63 1.35   0 4.64  4.64 0.27    -0.72 0.07
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 249 1.48 1.28   1.18    1.35 1.35   0   5     5 0.69    -0.47 0.08
d3 %>%
  ggplot( aes(x=HS, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .25) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
    theme_bw() +
    labs(fill="Gender",
         x = "Hostile Sexism") +
  facet_wrap(~Country)

SRQ

describeBy(d3$SRQ, d3$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n  mean    sd median trimmed   mad min   max range skew kurtosis   se
## X1    1 247 27.93 16.32  26.15   27.53 20.53   0 66.15 66.15 0.21    -0.95 1.04
## ------------------------------------------------------------ 
## group: Canada
##    vars   n  mean    sd median trimmed   mad min   max range skew kurtosis   se
## X1    1 250 29.94 16.38  29.23   29.51 18.25   0 76.15 76.15 0.23     -0.6 1.04
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean    sd median trimmed   mad min max range skew kurtosis   se
## X1    1 249 29.9 17.51  29.23   29.23 19.86   0  80    80 0.34    -0.52 1.11
## ------------------------------------------------------------ 
## group: US
##    vars   n  mean    sd median trimmed   mad min max range skew kurtosis   se
## X1    1 249 29.38 20.12  27.69   28.56 26.23   0  90    90 0.34    -0.77 1.28
d3 %>%
  ggplot( aes(x=SRQ, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = 5) +
  coord_cartesian(xlim = c(0,75)) +
  scale_x_continuous(breaks = seq(0,75,5)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Social Roles Questionnaire")+
  facet_wrap(~Country)

SDO

describeBy(d3$SDO, d3$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean   sd median trimmed mad min max range skew kurtosis   se
## X1    1 247 2.43 1.17   2.25    2.33 1.3   1 6.5   5.5 0.66    -0.05 0.07
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 250 2.55 1.26    2.5    2.45 1.48   1   7     6 0.62    -0.19 0.08
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 249 2.74 1.37   2.62    2.66 1.85   1   7     6 0.44    -0.57 0.09
## ------------------------------------------------------------ 
## group: US
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 249  2.4 1.58   1.75    2.17 1.11   1   7     6 0.98    -0.05 0.1
d3 %>%
  ggplot( aes(x=SDO, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .5) +
  scale_x_continuous(breaks = seq(1,7,1)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Social Dominance Orientation")+
  facet_wrap(~Country)

Ideology

describeBy(d3$Ideology, d3$Country)
## 
##  Descriptive statistics by group 
## group: Aus
##    vars   n mean  sd median trimmed  mad min max range skew kurtosis   se
## X1    1 247 3.22 1.3      3    3.15 1.48   1   7     6 0.45    -0.11 0.08
## ------------------------------------------------------------ 
## group: Canada
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 250 3.59 1.52   3.75    3.57 1.85   1   7     6 0.16    -0.64 0.1
## ------------------------------------------------------------ 
## group: UK
##    vars   n mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 249 3.68 1.51    3.5    3.64 1.85   1   7     6 0.23    -0.66 0.1
## ------------------------------------------------------------ 
## group: US
##    vars   n mean  sd median trimmed  mad min max range skew kurtosis   se
## X1    1 249 3.45 1.8   3.25    3.35 1.85   1   7     6 0.36    -0.92 0.11
d3 %>%
  ggplot( aes(x=Ideology, fill=Demo_Gender)) +
    geom_histogram( color="#e9ecef", alpha=0.6, position = 'identity',
                    binwidth = .5) +
  scale_x_continuous(breaks = seq(1,7,1)) +
  coord_cartesian(xlim = c(1,7)) +
    scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
    labs(fill="Gender",
         x = "Ideology (Liberal to Conservative)")+
  facet_wrap(~Country)

Analyses: Voting for Harris

  • Note on analyses: These models are binomial logistic regressions, evaluating how predictors affect likelihood of voting for Harris (vs. not)
  • Model results will be presented as odds ratios
    • ORs greater than 1 indicate positive association with likelihood to vote for Harris
    • ORs less than 1 indicate negative association with likelihood to vote for Harris
    • (the “Statistic” column in each model table is the Z-score assoc. with each OR)
  • Models control for demographic info
  • Continuous predictors are mean-centered

Study 1

Codes

d1$US_Can <- NA
d1$US_Can[d1$Country == "US"] <- -1/2
d1$US_Can[d1$Country == "Canada"] <- 1/2

d1$Country <- factor(d1$Country, levels = c("US","Canada"))

d1$man_woman <- NA
d1$man_woman[d1$Demo_Gender == "Woman"] <- 1/2
d1$man_woman[d1$Demo_Gender == "Man"] <- -1/2


contrasts(d1$Demo_Race) <- contr.poly(nlevels(d1$Demo_Race))

Models

(Dropping interaction terms with country because of model instability)

Sexist Attitudes x Gender

s1_glm1 <- glm(VotedHarris ~ (HS.c + BS.c + SRQ.c + SDO.c) * man_woman + Country
            + (edu.c + age.c + Demo_Race + Demo_Ethnicity + InfluenceCanada.c),
              data = d1,
              family = binomial)

# summary(s1_glm1)
tab_model(s1_glm1,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 0.79 0.43 0.27 – 2.29 -0.44 0.657
HS c 0.33 0.10 0.17 – 0.58 -3.58 <0.001
BS c 0.80 0.21 0.47 – 1.35 -0.83 0.404
SRQ c 0.99 0.02 0.96 – 1.03 -0.32 0.748
SDO c 0.68 0.15 0.43 – 1.05 -1.72 0.085
man woman 0.80 0.36 0.32 – 1.89 -0.50 0.614
Country [Canada] 1.54 0.66 0.66 – 3.61 1.01 0.314
edu c 0.84 0.07 0.71 – 1.00 -1.96 0.050
age c 1.02 0.02 0.98 – 1.06 0.98 0.330
Demo Race [.L] 0.85 1.15 0.06 – 14.43 -0.12 0.901
Demo Race [.Q] 0.38 0.48 0.03 – 4.74 -0.76 0.447
Demo Race [.C] 3.39 3.35 0.51 – 26.05 1.24 0.215
Demo Race [^4] 0.73 0.54 0.16 – 3.08 -0.43 0.668
Demo Race [^5] 2.34 1.26 0.82 – 6.85 1.58 0.114
Demo Ethnicity [1] 5.46 4.74 1.01 – 31.22 1.95 0.051
InfluenceCanada c 1.57 0.18 1.26 – 1.99 3.89 <0.001
HS c × man woman 0.12 0.08 0.03 – 0.39 -3.31 0.001
BS c × man woman 2.16 1.14 0.78 – 6.26 1.45 0.147
SRQ c × man woman 0.97 0.04 0.90 – 1.04 -0.84 0.403
SDO c × man woman 3.91 1.94 1.52 – 10.77 2.74 0.006
Observations 242
R2 Tjur 0.508

Reported Personal Impact x Gender

s1_glm2 <- glm(VotedHarris ~ man_woman * Impact.c
            + (age.c + edu.c + Demo_Race + Demo_Ethnicity + Country + InfluenceCanada.c),
              data = d1,
              family = binomial)

# summary(s1_glm2)
tab_model(s1_glm2,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.04 0.49 0.41 – 2.65 0.08 0.936
man woman 1.91 0.67 0.96 – 3.85 1.83 0.067
Impact c 1.85 0.25 1.43 – 2.42 4.58 <0.001
age c 0.99 0.01 0.96 – 1.02 -0.91 0.363
edu c 0.95 0.07 0.83 – 1.09 -0.75 0.456
Demo Race [.L] 3.59 4.63 0.30 – 54.90 0.99 0.322
Demo Race [.Q] 1.10 1.30 0.10 – 12.21 0.08 0.936
Demo Race [.C] 2.71 2.47 0.47 – 18.00 1.09 0.274
Demo Race [^4] 1.53 1.01 0.41 – 5.68 0.65 0.518
Demo Race [^5] 1.47 0.69 0.58 – 3.69 0.82 0.413
Demo Ethnicity [1] 2.17 1.56 0.54 – 9.17 1.07 0.283
Country [Canada] 2.72 1.05 1.29 – 5.92 2.58 0.010
InfluenceCanada c 1.33 0.16 1.06 – 1.69 2.46 0.014
man woman × Impact c 1.21 0.27 0.78 – 1.89 0.86 0.390
Observations 242
R2 Tjur 0.356

Full Model: Impact Index & Sexist Attitudes x Gender

s1_glm3 <- glm(VotedHarris ~ man_woman * (Impact.c + HS.c + BS.c + SRQ.c + SDO.c ) 
               + (age.c + edu.c + Demo_Race + Country + InfluenceCanada.c),
              data = d1,
              family = binomial)

# summary(s1_glm3)
tab_model(s1_glm3,
          show.stat = T,
          show.se = T) 
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.31 0.63 0.51 – 3.47 0.57 0.571
man woman 0.69 0.35 0.24 – 1.82 -0.74 0.460
Impact c 2.10 0.37 1.51 – 3.04 4.18 <0.001
HS c 0.43 0.15 0.21 – 0.81 -2.48 0.013
BS c 0.72 0.21 0.41 – 1.27 -1.13 0.260
SRQ c 0.99 0.02 0.95 – 1.03 -0.41 0.680
SDO c 0.56 0.14 0.34 – 0.92 -2.26 0.024
age c 1.00 0.02 0.96 – 1.03 -0.12 0.905
edu c 0.87 0.08 0.73 – 1.03 -1.60 0.109
Demo Race [.L] 0.96 1.29 0.07 – 16.30 -0.03 0.975
Demo Race [.Q] 0.55 0.66 0.05 – 6.55 -0.50 0.617
Demo Race [.C] 2.33 2.32 0.34 – 18.09 0.85 0.397
Demo Race [^4] 0.66 0.51 0.14 – 2.95 -0.53 0.594
Demo Race [^5] 1.37 0.78 0.44 – 4.18 0.54 0.587
Country [Canada] 2.15 0.96 0.90 – 5.27 1.70 0.089
InfluenceCanada c 1.05 0.16 0.77 – 1.42 0.29 0.772
man woman × Impact c 0.92 0.26 0.52 – 1.61 -0.29 0.769
man woman × HS c 0.16 0.11 0.04 – 0.59 -2.61 0.009
man woman × BS c 3.13 1.77 1.06 – 9.89 2.02 0.044
man woman × SRQ c 0.98 0.04 0.91 – 1.06 -0.53 0.594
man woman × SDO c 2.36 1.24 0.86 – 6.86 1.63 0.103
Observations 243
R2 Tjur 0.563

Data Visualization

d1 %>% ggplot(aes(x = Country, 
                  y = VotedHarris_num, 
                  fill = Demo_Gender)) +
  geom_bar(stat = "summary", 
           fun = "mean", 
           position = position_dodge(.9), 
           alpha = .8) +
  stat_summary(geom = "errorbar", 
               color = "black", 
               width = .1, 
               show.legend = F, 
               position = position_dodge(.9)) +
  theme_bw() +
  scale_y_continuous(breaks = seq(0,1,.2)) +
  coord_cartesian(ylim = c(0,1)) +
  scale_fill_manual("Gender", 
                    values = c("#69b3a2", "#404080")) +
  labs(y = "Proportion Voted for Harris",
       title = "Vote Choice by Gender and Country")

d1 %>% ggplot() +
  geom_jitter(aes(x = Impact_index, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender),
              alpha = .7) + 
#  geom_smooth(aes(x = Impact_index, 
#                  y = VotedHarris, 
#                  color = Demo_Gender, 
#                  group = Demo_Gender), 
#              method = "loess", se = F) +
  geom_smooth(aes(x = Impact_index, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender), 
              method = "lm", se = F, fullrange = T) +
  labs(x = "Election Impact on Wellbeing",
       y = "Voted for Harris",
       title = "Likelihood of Voting Harris by Election Impact, Gender, and Country") +
  scale_x_continuous(breaks = seq(1,7,1)) +
  scale_y_discrete(labels = c("No","Yes")) +
  theme_bw() +
  scale_color_manual("Gender", 
                     values = c("#69b3a2", "#404080")) +
  facet_wrap(~Country)

d1 %>% ggplot() +
  geom_jitter(aes(x = HS, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender),
              alpha = .7) + 
#  geom_smooth(aes(x = HS, 
#                  y = VotedHarris, 
#                  color = Demo_Gender, 
#                  group = Demo_Gender), 
#              method = "loess", se = F) +
  geom_smooth(aes(x = HS, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender), 
              method = "lm", se = F, fullrange = T) +
  labs(x = "Hostile Sexism",
       y = "Voted for Harris",
       title = "Likelihood of Voting Harris by Hostile Sexism and Gender") +
  scale_x_continuous(breaks = seq(0,7,1)) +
  scale_y_discrete(labels = c("No","Yes")) +
  coord_cartesian(xlim = c(0,5)) +
  theme_bw() +
  scale_color_manual("Gender", 
                     values = c("#69b3a2", "#404080"))

d1 %>% ggplot() +
  geom_jitter(aes(x = SDO, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender),
              alpha = .7) + 
#  geom_smooth(aes(x = SDO, 
#                  y = VotedHarris, 
#                  color = Demo_Gender, 
#                  group = Demo_Gender), 
#              method = "loess", se = F) +
  geom_smooth(aes(x = SDO, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender), 
              method = "lm", se = F, fullrange = T) +
  labs(x = "Social Dominance Orientation",
       y = "Voted for Harris",
       title = "Likelihood of Voting Harris by SDO, Gender, and Country") +
  scale_x_continuous(breaks = seq(0,7,1)) +
  scale_y_discrete(labels = c("No","Yes")) +
  coord_cartesian(xlim = c(1,7)) +
  theme_bw() +
  scale_color_manual("Gender", 
                     values = c("#69b3a2", "#404080")) 

Study 2

Codes

d2$US_oth <- NA
d2$US_oth[d2$Country == "US"] <- 3
d2$US_oth[d2$Country == "Canada"] <- -1
d2$US_oth[d2$Country == "Aus"] <- -1
d2$US_oth[d2$Country == "UK"] <- -1

d2$Can_AusUK <- NA
d2$Can_AusUK[d2$Country == "US"] <- 0
d2$Can_AusUK[d2$Country == "Canada"] <- 2
d2$Can_AusUK[d2$Country == "Aus"] <- -1
d2$Can_AusUK[d2$Country == "UK"] <- -1

d2$Aus_UK <- NA
d2$Aus_UK[d2$Country == "US"] <- 0
d2$Aus_UK[d2$Country == "Canada"] <- 0
d2$Aus_UK[d2$Country == "Aus"] <- -1
d2$Aus_UK[d2$Country == "UK"] <- 1

d2$Country <- factor(d2$Country, levels = c("US","Canada", "Aus", "UK"))

d2$man_woman <- ifelse(d2$Demo_Gender == "Woman", 1/2, -1/2)

contrasts(d2$Demo_Race) <- contr.poly(nlevels(d2$Demo_Race))

Models

Insufficient power in race, ethnicity variables; omitting for these

Sexist Attitudes x Gender x Country

s2_glm1 <- glm(VotedHarris ~ (HS.c + BS.c + SRQ.c + SDO.c) * (US_oth + Can_AusUK + Aus_UK) * man_woman
            + (age.c + edu.c + InfluenceCanada.c),
              data = d2,
              family = binomial)

# summary(s2_glm1)
tab_model(s2_glm1,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.83 0.25 1.40 – 2.41 4.37 <0.001
HS c 0.61 0.10 0.44 – 0.85 -2.90 0.004
BS c 1.02 0.18 0.72 – 1.45 0.12 0.906
SRQ c 0.96 0.01 0.93 – 0.98 -3.76 <0.001
SDO c 0.93 0.12 0.72 – 1.21 -0.54 0.590
US oth 0.82 0.06 0.71 – 0.95 -2.58 0.010
Can AusUK 0.79 0.09 0.64 – 0.98 -2.13 0.033
Aus UK 0.77 0.16 0.50 – 1.14 -1.29 0.197
man woman 1.36 0.37 0.79 – 2.33 1.13 0.260
age c 1.03 0.01 1.01 – 1.05 2.65 0.008
edu c 0.96 0.04 0.89 – 1.05 -0.89 0.372
InfluenceCanada c 1.92 0.16 1.65 – 2.27 8.10 <0.001
HS c × US oth 0.92 0.08 0.76 – 1.10 -0.93 0.355
HS c × Can AusUK 1.11 0.15 0.84 – 1.46 0.73 0.463
HS c × Aus UK 1.22 0.30 0.75 – 2.00 0.80 0.423
BS c × US oth 0.95 0.09 0.79 – 1.14 -0.57 0.567
BS c × Can AusUK 0.85 0.13 0.63 – 1.13 -1.14 0.256
BS c × Aus UK 0.97 0.26 0.57 – 1.64 -0.11 0.910
SRQ c × US oth 1.01 0.01 1.00 – 1.03 1.79 0.074
SRQ c × Can AusUK 1.02 0.01 1.00 – 1.04 2.23 0.026
SRQ c × Aus UK 1.01 0.02 0.97 – 1.04 0.38 0.706
SDO c × US oth 0.95 0.07 0.82 – 1.09 -0.70 0.482
SDO c × Can AusUK 0.79 0.09 0.64 – 0.97 -2.19 0.029
SDO c × Aus UK 1.24 0.25 0.84 – 1.85 1.08 0.278
HS c × man woman 1.48 0.51 0.76 – 2.93 1.13 0.258
BS c × man woman 1.03 0.37 0.51 – 2.08 0.09 0.928
SRQ c × man woman 1.00 0.02 0.96 – 1.05 0.05 0.958
SDO c × man woman 0.90 0.24 0.53 – 1.50 -0.41 0.679
US oth × man woman 1.00 0.15 0.74 – 1.34 -0.02 0.988
Can AusUK × man woman 1.01 0.22 0.66 – 1.56 0.05 0.960
Aus UK × man woman 1.17 0.49 0.52 – 2.68 0.38 0.706
(HS c × US oth) × man
woman
1.12 0.20 0.79 – 1.62 0.62 0.534
(HS c × Can AusUK) × man
woman
0.68 0.19 0.39 – 1.18 -1.36 0.172
(HS c × Aus UK) × man
woman
2.64 1.32 1.00 – 7.15 1.95 0.052
(BS c × US oth) × man
woman
1.00 0.18 0.70 – 1.42 -0.03 0.979
(BS c × Can AusUK) × man
woman
1.90 0.56 1.08 – 3.47 2.17 0.030
(BS c × Aus UK) × man
woman
0.35 0.19 0.12 – 0.99 -1.96 0.050
(SRQ c × US oth) × man
woman
0.99 0.01 0.96 – 1.01 -0.94 0.346
(SRQ c × Can AusUK) × man
woman
0.96 0.02 0.93 – 1.00 -1.80 0.072
(SRQ c × Aus UK) × man
woman
0.97 0.03 0.90 – 1.03 -0.98 0.326
(SDO c × US oth) × man
woman
1.09 0.16 0.83 – 1.45 0.62 0.537
(SDO c × Can AusUK) × man
woman
1.37 0.30 0.90 – 2.11 1.46 0.143
(SDO c × Aus UK) × man
woman
0.74 0.30 0.33 – 1.61 -0.75 0.452
Observations 581
R2 Tjur 0.472

Personal Impact Index x Gender x Country

s2_glm2 <- glm(VotedHarris ~ (US_oth + Can_AusUK + Aus_UK) * man_woman * Impact.c
            + (age.c + edu.c + InfluenceCanada.c),
              data = d2,
              family = binomial)

# summary(s2_glm2)
tab_model(s2_glm2,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.57 0.18 1.26 – 1.96 3.98 <0.001
US oth 0.78 0.05 0.68 – 0.89 -3.65 <0.001
Can AusUK 0.97 0.09 0.81 – 1.17 -0.38 0.707
Aus UK 0.77 0.12 0.57 – 1.04 -1.68 0.094
man woman 2.11 0.48 1.36 – 3.31 3.31 0.001
Impact c 1.59 0.14 1.35 – 1.89 5.47 <0.001
age c 1.03 0.01 1.01 – 1.05 3.21 0.001
edu c 0.97 0.04 0.90 – 1.04 -0.84 0.403
InfluenceCanada c 1.52 0.12 1.32 – 1.78 5.54 <0.001
US oth × man woman 0.80 0.11 0.62 – 1.04 -1.67 0.094
Can AusUK × man woman 0.99 0.18 0.69 – 1.44 -0.07 0.946
Aus UK × man woman 0.78 0.24 0.42 – 1.42 -0.81 0.420
US oth × Impact c 1.10 0.05 1.01 – 1.21 2.15 0.031
Can AusUK × Impact c 1.06 0.07 0.94 – 1.21 0.96 0.336
Aus UK × Impact c 1.02 0.11 0.83 – 1.25 0.21 0.830
man woman × Impact c 1.20 0.18 0.90 – 1.62 1.22 0.224
(US oth × man woman) ×
Impact c
0.93 0.08 0.78 – 1.11 -0.80 0.422
(Can AusUK × man woman) ×
Impact c
1.07 0.13 0.84 – 1.38 0.56 0.574
(Aus UK × man woman) ×
Impact c
1.00 0.21 0.67 – 1.51 0.01 0.994
Observations 581
R2 Tjur 0.337

Ideology x Gender x Country

s2_glm3 <- glm(VotedHarris ~ man_woman * Ideology.c * (US_oth + Can_AusUK + Aus_UK)
            + (age.c + edu.c + InfluenceCanada.c),
              data = d2,
              family = binomial)

# summary(s2_glm3)
tab_model(s2_glm3,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 0.92 0.12 0.71 – 1.18 -0.68 0.498
man woman 2.12 0.56 1.27 – 3.56 2.87 0.004
Ideology c 0.43 0.05 0.34 – 0.52 -7.88 <0.001
US oth 0.81 0.07 0.68 – 0.94 -2.59 0.009
Can AusUK 0.94 0.09 0.77 – 1.14 -0.62 0.535
Aus UK 0.92 0.16 0.65 – 1.30 -0.48 0.630
age c 1.04 0.01 1.02 – 1.06 3.97 <0.001
edu c 0.96 0.04 0.89 – 1.04 -1.02 0.306
InfluenceCanada c 1.65 0.12 1.43 – 1.90 6.90 <0.001
man woman × Ideology c 0.92 0.20 0.60 – 1.40 -0.39 0.699
man woman × US oth 0.91 0.15 0.66 – 1.26 -0.56 0.575
man woman × Can AusUK 0.90 0.18 0.60 – 1.33 -0.54 0.587
man woman × Aus UK 0.80 0.29 0.40 – 1.61 -0.62 0.534
Ideology c × US oth 0.88 0.05 0.77 – 0.99 -2.00 0.046
Ideology c × Can AusUK 0.98 0.08 0.83 – 1.14 -0.30 0.765
Ideology c × Aus UK 1.21 0.19 0.89 – 1.66 1.19 0.234
(man woman × Ideology c)
× US oth
0.86 0.11 0.67 – 1.09 -1.21 0.228
(man woman × Ideology c)
× Can AusUK
0.95 0.15 0.69 – 1.30 -0.33 0.745
(man woman × Ideology c)
× Aus UK
0.92 0.29 0.50 – 1.73 -0.26 0.794
Observations 581
R2 Tjur 0.425

Impact index, sexist attitudes, ideology x Gender x Country

s2_glm4 <- glm(VotedHarris ~ man_woman * (US_oth + Can_AusUK + Aus_UK) * (Impact.c + HS.c + BS.c + SRQ.c + SDO.c + Ideology.c)
               + age.c + edu.c + InfluenceCanada.c,
              data = d2,
              family = binomial)

# summary(s2_glm4)
tab_model(s2_glm4,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.55 0.29 1.08 – 2.25 2.38 0.018
man woman 1.28 0.48 0.62 – 2.67 0.66 0.506
US oth 0.75 0.08 0.61 – 0.93 -2.65 0.008
Can AusUK 0.78 0.11 0.59 – 1.03 -1.75 0.080
Aus UK 0.92 0.26 0.52 – 1.60 -0.31 0.756
Impact c 1.58 0.19 1.26 – 2.00 3.87 <0.001
HS c 0.70 0.14 0.47 – 1.04 -1.77 0.077
BS c 0.88 0.19 0.58 – 1.33 -0.61 0.544
SRQ c 0.96 0.01 0.94 – 0.99 -2.89 0.004
SDO c 1.08 0.17 0.80 – 1.46 0.48 0.632
Ideology c 0.59 0.08 0.45 – 0.76 -3.94 <0.001
age c 1.03 0.01 1.01 – 1.06 2.94 0.003
edu c 0.95 0.05 0.86 – 1.04 -1.14 0.255
InfluenceCanada c 1.47 0.16 1.20 – 1.83 3.62 <0.001
man woman × US oth 0.91 0.19 0.60 – 1.37 -0.44 0.657
man woman × Can AusUK 0.99 0.28 0.56 – 1.74 -0.03 0.978
man woman × Aus UK 1.03 0.59 0.33 – 3.13 0.05 0.959
man woman × Impact c 1.05 0.22 0.70 – 1.58 0.25 0.805
man woman × HS c 1.73 0.71 0.79 – 3.92 1.34 0.179
man woman × BS c 0.83 0.35 0.36 – 1.88 -0.44 0.661
man woman × SRQ c 1.01 0.03 0.96 – 1.07 0.35 0.726
man woman × SDO c 1.00 0.31 0.54 – 1.82 -0.02 0.988
man woman × Ideology c 0.74 0.20 0.43 – 1.24 -1.15 0.250
US oth × Impact c 1.10 0.07 0.98 – 1.26 1.47 0.143
US oth × HS c 1.01 0.13 0.79 – 1.30 0.12 0.908
US oth × BS c 0.78 0.10 0.59 – 0.99 -1.94 0.052
US oth × SRQ c 1.02 0.01 1.00 – 1.04 2.38 0.017
US oth × SDO c 0.95 0.09 0.79 – 1.13 -0.61 0.539
US oth × Ideology c 0.89 0.07 0.76 – 1.03 -1.54 0.123
Can AusUK × Impact c 1.00 0.08 0.86 – 1.17 0.01 0.988
Can AusUK × HS c 1.14 0.18 0.85 – 1.56 0.87 0.382
Can AusUK × BS c 0.84 0.14 0.62 – 1.16 -1.06 0.290
Can AusUK × SRQ c 1.02 0.01 1.00 – 1.05 2.10 0.036
Can AusUK × SDO c 0.87 0.11 0.68 – 1.10 -1.15 0.251
Can AusUK × Ideology c 0.94 0.10 0.77 – 1.16 -0.55 0.586
Aus UK × Impact c 1.27 0.19 0.95 – 1.72 1.60 0.110
Aus UK × HS c 1.29 0.36 0.75 – 2.25 0.92 0.359
Aus UK × BS c 0.95 0.27 0.55 – 1.67 -0.16 0.871
Aus UK × SRQ c 1.00 0.02 0.96 – 1.03 -0.23 0.817
Aus UK × SDO c 1.29 0.28 0.85 – 1.99 1.17 0.242
Aus UK × Ideology c 1.02 0.21 0.68 – 1.52 0.09 0.925
(man woman × US oth) ×
Impact c
0.95 0.12 0.73 – 1.23 -0.38 0.707
(man woman × US oth) × HS
c
1.15 0.29 0.72 – 1.92 0.57 0.569
(man woman × US oth) × BS
c
0.83 0.21 0.48 – 1.35 -0.74 0.459
(man woman × US oth) ×
SRQ c
1.00 0.02 0.96 – 1.03 -0.24 0.808
(man woman × US oth) ×
SDO c
1.26 0.23 0.89 – 1.82 1.26 0.207
(man woman × US oth) ×
Ideology c
0.74 0.11 0.54 – 0.99 -2.03 0.043
(man woman × Can AusUK) ×
Impact c
1.08 0.17 0.80 – 1.48 0.52 0.604
(man woman × Can AusUK) ×
HS c
0.64 0.20 0.35 – 1.18 -1.44 0.150
(man woman × Can AusUK) ×
BS c
1.82 0.59 0.98 – 3.50 1.86 0.062
(man woman × Can AusUK) ×
SRQ c
0.97 0.02 0.93 – 1.01 -1.43 0.153
(man woman × Can AusUK) ×
SDO c
1.62 0.40 1.01 – 2.66 1.97 0.049
(man woman × Can AusUK) ×
Ideology c
0.93 0.19 0.62 – 1.41 -0.33 0.738
(man woman × Aus UK) ×
Impact c
1.17 0.35 0.65 – 2.12 0.53 0.593
(man woman × Aus UK) × HS
c
3.55 1.98 1.20 – 10.92 2.26 0.024
(man woman × Aus UK) × BS
c
0.37 0.21 0.12 – 1.13 -1.73 0.083
(man woman × Aus UK) ×
SRQ c
0.96 0.04 0.89 – 1.03 -1.10 0.270
(man woman × Aus UK) ×
SDO c
0.73 0.32 0.30 – 1.68 -0.73 0.464
(man woman × Aus UK) ×
Ideology c
0.76 0.31 0.34 – 1.70 -0.67 0.503
Observations 581
R2 Tjur 0.558

Data Visualization

d2 %>% ggplot(aes(x = Country, 
                  y = VotedHarris_num, 
                  fill = Demo_Gender)) +
  geom_bar(stat = "summary", 
           fun = "mean", 
           position = position_dodge(.9), 
           alpha = .8) +
  stat_summary(geom = "errorbar", 
               color = "black", 
               width = .1, 
               show.legend = F, 
               position = position_dodge(.9)) +
  theme_bw() +
  scale_y_continuous(breaks = seq(0,1,.2)) +
  coord_cartesian(ylim = c(0,1)) +
  scale_fill_manual("Gender", 
                    values = c("#69b3a2", "#404080")) +
  labs(y = "Proportion Voted for Harris",
       title = "Vote Choice by Gender and Country")

d2 %>% ggplot(aes(x = Ideology, y = VotedHarris_num, fill = Demo_Gender, color = Demo_Gender)) +
  geom_jitter(size = 1, alpha = .7) +
#  geom_smooth(method = "loess", se = F, fullrange = T) +
  geom_smooth(method = "lm", se = F, fullrange = T, fullrange = T) +
  theme_bw() +
#  facet_wrap(~Country) +
  scale_fill_manual(values = c("#69b3a2", "#404080"))+
  scale_color_manual(values = c("#69b3a2", "#404080")) +
  labs(fill = "Gender", color = "Gender",
       y = "Voted for Harris",
       x = "Ideology (Liberal to Conseravtive)",
       title = "Likelihood of Voting Harris by Ideology and Gender",
       caption= "Interestingly, ideology has narrower range for 'yes' voters than 'no' voters") +
  scale_y_continuous(breaks = c(0,1),
                     labels = c("No", 
                              "Yes")) +
  scale_x_continuous(breaks = seq(1,7,1))

d2 %>% ggplot() +
  geom_jitter(aes(x = Impact_index, 
                  y = VotedHarris_num),
              color = "#69b3a2",
              size = 1) + 
 # geom_smooth(aes(x = Impact_index, 
 #                 y = VotedHarris_num), 
 #             method = "loess", se = F, color = "#404080") +
  geom_smooth(aes(x = Impact_index, 
                  y = VotedHarris_num), 
              method = "lm", se = F, color = "#404080", fullrange = T) +
  labs(x = "Election Impact on Wellbeing",
       y = "Voted for Harris",
       title = "Likelihood of Voting Harris by Impact Index") +
  scale_x_continuous(breaks = seq(1,7,1)) +
  scale_y_continuous(labels = c("No","Yes"),
                    breaks = c(0,1)) +
  theme_bw()

d2 %>% ggplot() +
  geom_jitter(aes(x = SRQ, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender),
              size = 1, alpha = .7) + 
#  geom_smooth(aes(x = SRQ, 
#                  y = VotedHarris, 
#                  color = Demo_Gender, 
#                  group = Demo_Gender), 
#              method = "loess", se = F, fullrange = T) +
  geom_smooth(aes(x = SRQ, 
                  y = VotedHarris, 
                  color = Demo_Gender, 
                  group = Demo_Gender), 
              method = "lm", se = F, fullrange = T, fullrange = T) +
  labs(x = "Social Roles Questionnaire",
       y = "Voted for Harris",
       title = "Likelihood of Voting Harris by SDO, Gender, and Country") +
  scale_x_continuous(breaks = seq(0,80,10)) +
  scale_y_discrete(labels = c("No", 
                              "Yes")) +
  coord_cartesian(xlim = c(0,80)) +
  theme_bw() +
  scale_color_manual("Gender", 
                     values = c("#69b3a2", "#404080")) +
  facet_wrap(~Country)

Study 3

Codes

d3$US_oth <- NA
d3$US_oth[d3$Country == "US"] <- 3
d3$US_oth[d3$Country == "Canada"] <- -1
d3$US_oth[d3$Country == "Aus"] <- -1
d3$US_oth[d3$Country == "UK"] <- -1

d3$Can_AusUK <- NA
d3$Can_AusUK[d3$Country == "US"] <- 0
d3$Can_AusUK[d3$Country == "Canada"] <- 2
d3$Can_AusUK[d3$Country == "Aus"] <- -1
d3$Can_AusUK[d3$Country == "UK"] <- -1

d3$Aus_UK <- NA
d3$Aus_UK[d3$Country == "US"] <- 0
d3$Aus_UK[d3$Country == "Canada"] <- 0
d3$Aus_UK[d3$Country == "Aus"] <- -1
d3$Aus_UK[d3$Country == "UK"] <- 1

d3$Country <- factor(d3$Country, levels = c("US","Canada", "Aus", "UK"))

d3$man_woman <- ifelse(d3$Demo_Gender == "Woman", 1/2, -1/2)

Models

Sexist attitudes x Country x Gender

s3_glm1 <- glm(VotedHarris ~ (HS.c + BS.c + SRQ.c + SDO.c) * (US_oth + Can_AusUK + Aus_UK) * man_woman
            + (age.c + edu.c),
              data = d3,
              family = binomial)

# summary(s3_glm1)
tab_model(s3_glm1,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.54 0.13 1.30 – 1.83 4.95 <0.001
HS c 0.47 0.05 0.37 – 0.58 -6.78 <0.001
BS c 0.95 0.11 0.76 – 1.19 -0.44 0.662
SRQ c 0.98 0.01 0.96 – 0.99 -3.01 0.003
SDO c 0.80 0.06 0.68 – 0.93 -2.79 0.005
US oth 0.97 0.05 0.87 – 1.07 -0.69 0.491
Can AusUK 1.07 0.08 0.92 – 1.24 0.85 0.396
Aus UK 1.09 0.13 0.87 – 1.38 0.75 0.453
man woman 0.91 0.16 0.64 – 1.28 -0.55 0.584
age c 1.00 0.01 0.98 – 1.00 -0.65 0.517
edu c 1.02 0.03 0.97 – 1.09 0.80 0.425
HS c × US oth 1.11 0.07 0.98 – 1.26 1.69 0.091
HS c × Can AusUK 0.93 0.09 0.77 – 1.12 -0.74 0.456
HS c × Aus UK 1.22 0.19 0.90 – 1.66 1.24 0.214
BS c × US oth 1.04 0.07 0.92 – 1.18 0.66 0.507
BS c × Can AusUK 0.95 0.09 0.78 – 1.15 -0.55 0.582
BS c × Aus UK 1.01 0.17 0.73 – 1.40 0.07 0.945
SRQ c × US oth 1.00 0.00 0.99 – 1.01 -0.08 0.940
SRQ c × Can AusUK 0.99 0.01 0.98 – 1.00 -1.33 0.185
SRQ c × Aus UK 0.98 0.01 0.96 – 1.01 -1.37 0.169
SDO c × US oth 0.88 0.04 0.80 – 0.95 -3.02 0.003
SDO c × Can AusUK 1.03 0.07 0.91 – 1.18 0.49 0.622
SDO c × Aus UK 0.96 0.11 0.77 – 1.21 -0.32 0.751
HS c × man woman 1.19 0.27 0.77 – 1.87 0.78 0.433
BS c × man woman 1.26 0.29 0.80 – 1.98 0.99 0.321
SRQ c × man woman 0.97 0.02 0.94 – 1.00 -1.72 0.085
SDO c × man woman 1.06 0.17 0.77 – 1.46 0.37 0.710
US oth × man woman 0.88 0.09 0.72 – 1.07 -1.24 0.213
Can AusUK × man woman 0.94 0.14 0.70 – 1.25 -0.41 0.679
Aus UK × man woman 0.66 0.15 0.42 – 1.04 -1.78 0.075
(HS c × US oth) × man
woman
0.76 0.10 0.59 – 0.97 -2.17 0.030
(HS c × Can AusUK) × man
woman
1.16 0.22 0.80 – 1.70 0.75 0.452
(HS c × Aus UK) × man
woman
0.78 0.25 0.42 – 1.44 -0.78 0.438
(BS c × US oth) × man
woman
1.00 0.13 0.78 – 1.29 0.01 0.992
(BS c × Can AusUK) × man
woman
0.77 0.15 0.52 – 1.11 -1.39 0.165
(BS c × Aus UK) × man
woman
0.56 0.18 0.29 – 1.06 -1.77 0.076
(SRQ c × US oth) × man
woman
1.01 0.01 1.00 – 1.03 1.68 0.094
(SRQ c × Can AusUK) × man
woman
1.01 0.01 0.99 – 1.04 0.84 0.402
(SRQ c × Aus UK) × man
woman
1.04 0.02 0.99 – 1.09 1.67 0.094
(SDO c × US oth) × man
woman
0.87 0.08 0.73 – 1.03 -1.61 0.106
(SDO c × Can AusUK) × man
woman
0.83 0.11 0.63 – 1.08 -1.37 0.169
(SDO c × Aus UK) × man
woman
1.19 0.28 0.75 – 1.89 0.75 0.451
Observations 979
R2 Tjur 0.336

Ideology x Country x Gender

s3_glm2 <- glm(VotedHarris ~ Ideology.c * (US_oth + Can_AusUK + Aus_UK) * man_woman
            + (age.c + edu.c),
              data = d3,
              family = binomial)

# summary(s3_glm2)
tab_model(s3_glm2,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 0.97 0.08 0.82 – 1.14 -0.37 0.711
Ideology c 0.40 0.03 0.35 – 0.46 -13.43 <0.001
US oth 0.98 0.05 0.88 – 1.09 -0.33 0.740
Can AusUK 1.02 0.07 0.90 – 1.16 0.37 0.711
Aus UK 1.15 0.13 0.91 – 1.44 1.16 0.245
man woman 1.40 0.24 1.00 – 1.95 1.97 0.049
age c 1.00 0.00 0.99 – 1.00 -0.59 0.554
edu c 1.06 0.03 1.00 – 1.12 2.01 0.044
Ideology c × US oth 0.89 0.04 0.81 – 0.96 -2.70 0.007
Ideology c × Can AusUK 1.00 0.05 0.90 – 1.10 0.00 0.997
Ideology c × Aus UK 1.01 0.09 0.84 – 1.21 0.05 0.957
Ideology c × man woman 1.14 0.15 0.87 – 1.49 0.96 0.336
US oth × man woman 0.80 0.09 0.65 – 0.99 -2.03 0.043
Can AusUK × man woman 1.09 0.14 0.85 – 1.41 0.69 0.491
Aus UK × man woman 0.84 0.20 0.53 – 1.33 -0.74 0.461
(Ideology c × US oth) ×
man woman
0.82 0.07 0.68 – 0.96 -2.38 0.018
(Ideology c × Can AusUK)
× man woman
0.93 0.10 0.76 – 1.14 -0.68 0.494
(Ideology c × Aus UK) ×
man woman
1.33 0.25 0.92 – 1.94 1.53 0.127
Observations 979
R2 Tjur 0.313

Warmth & competence

s3_glm2 <- glm(VotedHarris ~ (warmth.c + competence.c) * (US_oth + Can_AusUK + Aus_UK) * man_woman
            + (age.c + edu.c),
              data = d3,
              family = binomial)

summary(s3_glm2)
## 
## Call:
## glm(formula = VotedHarris ~ (warmth.c + competence.c) * (US_oth + 
##     Can_AusUK + Aus_UK) * man_woman + (age.c + edu.c), family = binomial, 
##     data = d3)
## 
## Coefficients:
##                                    Estimate Std. Error z value Pr(>|z|)    
## (Intercept)                       0.3389985  0.0692490   4.895 9.81e-07 ***
## warmth.c                          0.2675652  0.0642497   4.164 3.12e-05 ***
## competence.c                     -0.1615851  0.1065984  -1.516 0.129562    
## US_oth                            0.0349062  0.0395440   0.883 0.377388    
## Can_AusUK                        -0.0002646  0.0564642  -0.005 0.996260    
## Aus_UK                           -0.0468534  0.0994646  -0.471 0.637600    
## man_woman                         0.5212620  0.1382042   3.772 0.000162 ***
## age.c                            -0.0024989  0.0033456  -0.747 0.455108    
## edu.c                             0.0615950  0.0248899   2.475 0.013335 *  
## warmth.c:US_oth                   0.0175047  0.0348215   0.503 0.615177    
## warmth.c:Can_AusUK               -0.0886020  0.0530949  -1.669 0.095167 .  
## warmth.c:Aus_UK                  -0.0850188  0.0949122  -0.896 0.370379    
## competence.c:US_oth              -0.0779626  0.0595107  -1.310 0.190176    
## competence.c:Can_AusUK            0.0020077  0.0869716   0.023 0.981583    
## competence.c:Aus_UK              -0.0924168  0.1557753  -0.593 0.553001    
## warmth.c:man_woman               -0.2148910  0.1286956  -1.670 0.094967 .  
## competence.c:man_woman            0.2391955  0.2133682   1.121 0.262268    
## US_oth:man_woman                 -0.0927734  0.0784277  -1.183 0.236843    
## Can_AusUK:man_woman               0.0919844  0.1126824   0.816 0.414320    
## Aus_UK:man_woman                 -0.2447386  0.1986164  -1.232 0.217868    
## warmth.c:US_oth:man_woman        -0.0181967  0.0696461  -0.261 0.793881    
## warmth.c:Can_AusUK:man_woman     -0.0094764  0.1059794  -0.089 0.928750    
## warmth.c:Aus_UK:man_woman        -0.1727156  0.1899533  -0.909 0.363217    
## competence.c:US_oth:man_woman     0.0221224  0.1190105   0.186 0.852534    
## competence.c:Can_AusUK:man_woman -0.2729341  0.1739395  -1.569 0.116617    
## competence.c:Aus_UK:man_woman     0.4704317  0.3115116   1.510 0.131003    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 1326.8  on 977  degrees of freedom
## Residual deviance: 1259.4  on 952  degrees of freedom
##   (17 observations deleted due to missingness)
## AIC: 1311.4
## 
## Number of Fisher Scoring iterations: 5
tab_model(s3_glm2,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.40 0.10 1.23 – 1.61 4.90 <0.001
warmth c 1.31 0.08 1.15 – 1.48 4.16 <0.001
competence c 0.85 0.09 0.69 – 1.05 -1.52 0.130
US oth 1.04 0.04 0.96 – 1.12 0.88 0.377
Can AusUK 1.00 0.06 0.90 – 1.12 -0.00 0.996
Aus UK 0.95 0.09 0.79 – 1.16 -0.47 0.638
man woman 1.68 0.23 1.29 – 2.21 3.77 <0.001
age c 1.00 0.00 0.99 – 1.00 -0.75 0.455
edu c 1.06 0.03 1.01 – 1.12 2.47 0.013
warmth c × US oth 1.02 0.04 0.95 – 1.09 0.50 0.615
warmth c × Can AusUK 0.92 0.05 0.82 – 1.02 -1.67 0.095
warmth c × Aus UK 0.92 0.09 0.76 – 1.11 -0.90 0.370
competence c × US oth 0.92 0.06 0.82 – 1.04 -1.31 0.190
competence c × Can AusUK 1.00 0.09 0.84 – 1.19 0.02 0.982
competence c × Aus UK 0.91 0.14 0.67 – 1.24 -0.59 0.553
warmth c × man woman 0.81 0.10 0.63 – 1.04 -1.67 0.095
competence c × man woman 1.27 0.27 0.84 – 1.93 1.12 0.262
US oth × man woman 0.91 0.07 0.78 – 1.06 -1.18 0.237
Can AusUK × man woman 1.10 0.12 0.88 – 1.37 0.82 0.414
Aus UK × man woman 0.78 0.16 0.53 – 1.16 -1.23 0.218
(warmth c × US oth) × man
woman
0.98 0.07 0.86 – 1.12 -0.26 0.794
(warmth c × Can AusUK) ×
man woman
0.99 0.10 0.80 – 1.22 -0.09 0.929
(warmth c × Aus UK) × man
woman
0.84 0.16 0.58 – 1.22 -0.91 0.363
(competence c × US oth) ×
man woman
1.02 0.12 0.81 – 1.29 0.19 0.853
(competence c × Can
AusUK) × man woman
0.76 0.13 0.54 – 1.07 -1.57 0.117
(competence c × Aus UK) ×
man woman
1.60 0.50 0.87 – 2.96 1.51 0.131
Observations 978
R2 Tjur 0.067

Sexist Attitudes x Gender x Country + Ideology (can’t handle fully interactive model)

s3_glm3 <- glm(VotedHarris ~ (HS.c + BS.c + SRQ.c + SDO.c + Ideology.c) * (warmth.c + competence.c) * man_woman * (US_oth + Can_AusUK + Aus_UK) 
            + (age.c + edu.c),
              data = d3,
              family = binomial)

summary(s3_glm3)
## 
## Call:
## glm(formula = VotedHarris ~ (HS.c + BS.c + SRQ.c + SDO.c + Ideology.c) * 
##     (warmth.c + competence.c) * man_woman * (US_oth + Can_AusUK + 
##     Aus_UK) + (age.c + edu.c), family = binomial, data = d3)
## 
## Coefficients:
##                                              Estimate Std. Error z value
## (Intercept)                                  0.265164   0.142103   1.866
## HS.c                                        -0.959002   0.193753  -4.950
## BS.c                                         0.554932   0.259281   2.140
## SRQ.c                                       -0.060247   0.014301  -4.213
## SDO.c                                        0.244764   0.137135   1.785
## Ideology.c                                  -1.092452   0.132677  -8.234
## warmth.c                                     0.390115   0.145017   2.690
## competence.c                                -0.531197   0.241003  -2.204
## man_woman                                    0.195659   0.284809   0.687
## US_oth                                      -0.001937   0.095561  -0.020
## Can_AusUK                                    0.276661   0.109865   2.518
## Aus_UK                                       0.253666   0.174195   1.456
## age.c                                       -0.001482   0.002211  -0.670
## edu.c                                        0.022312   0.039815   0.560
## HS.c:warmth.c                                0.188844   0.170756   1.106
## HS.c:competence.c                            0.498607   0.317420   1.571
## BS.c:warmth.c                               -0.270589   0.173316  -1.561
## BS.c:competence.c                           -0.708917   0.416966  -1.700
## SRQ.c:warmth.c                               0.007677   0.013101   0.586
## SRQ.c:competence.c                           0.020643   0.020027   1.031
## SDO.c:warmth.c                               0.012918   0.131066   0.099
## SDO.c:competence.c                          -0.239088   0.207592  -1.152
## Ideology.c:warmth.c                          0.144405   0.104885   1.377
## Ideology.c:competence.c                      0.372628   0.184215   2.023
## HS.c:man_woman                               0.601137   0.387126   1.553
## BS.c:man_woman                               0.676628   0.519164   1.303
## SRQ.c:man_woman                             -0.088260   0.028586  -3.088
## SDO.c:man_woman                             -0.073392   0.273399  -0.268
## Ideology.c:man_woman                         0.037064   0.265193   0.140
## warmth.c:man_woman                          -0.720897   0.290361  -2.483
## competence.c:man_woman                      -0.235972   0.482312  -0.489
## HS.c:US_oth                                 -0.149733   0.137004  -1.093
## HS.c:Can_AusUK                               0.023498   0.141117   0.167
## HS.c:Aus_UK                                  0.238263   0.228545   1.043
## BS.c:US_oth                                  0.483384   0.221399   2.183
## BS.c:Can_AusUK                              -0.098818   0.136351  -0.725
## BS.c:Aus_UK                                  0.058860   0.232138   0.254
## SRQ.c:US_oth                                -0.010620   0.010036  -1.058
## SRQ.c:Can_AusUK                             -0.025307   0.011330  -2.234
## SRQ.c:Aus_UK                                -0.013881   0.015386  -0.902
## SDO.c:US_oth                                -0.170401   0.080643  -2.113
## SDO.c:Can_AusUK                              0.144690   0.121404   1.192
## SDO.c:Aus_UK                                -0.156211   0.170523  -0.916
## Ideology.c:US_oth                           -0.177876   0.089041  -1.998
## Ideology.c:Can_AusUK                        -0.031200   0.103377  -0.302
## Ideology.c:Aus_UK                            0.086786   0.161552   0.537
## warmth.c:US_oth                             -0.076198   0.090256  -0.844
## warmth.c:Can_AusUK                          -0.008904   0.118090  -0.075
## warmth.c:Aus_UK                             -0.258134   0.189482  -1.362
## competence.c:US_oth                         -0.344361   0.173678  -1.983
## competence.c:Can_AusUK                      -0.180617   0.172947  -1.044
## competence.c:Aus_UK                         -0.347484   0.279624  -1.243
## man_woman:US_oth                             0.002798   0.191085   0.015
## man_woman:Can_AusUK                         -0.208070   0.219382  -0.948
## man_woman:Aus_UK                            -0.362026   0.348320  -1.039
## HS.c:warmth.c:man_woman                      0.090326   0.341047   0.265
## HS.c:competence.c:man_woman                 -0.358153   0.634713  -0.564
## BS.c:warmth.c:man_woman                     -0.142974   0.345474  -0.414
## BS.c:competence.c:man_woman                 -0.476809   0.834760  -0.571
## SRQ.c:warmth.c:man_woman                     0.048439   0.026195   1.849
## SRQ.c:competence.c:man_woman                 0.084956   0.040195   2.114
## SDO.c:warmth.c:man_woman                    -0.456008   0.262137  -1.740
## SDO.c:competence.c:man_woman                -0.494934   0.414371  -1.194
## Ideology.c:warmth.c:man_woman               -0.224152   0.209665  -1.069
## Ideology.c:competence.c:man_woman           -0.021255   0.367097  -0.058
## HS.c:warmth.c:US_oth                         0.272031   0.102363   2.658
## HS.c:warmth.c:Can_AusUK                     -0.351339   0.147533  -2.381
## HS.c:warmth.c:Aus_UK                        -0.130904   0.215738  -0.607
## HS.c:competence.c:US_oth                     0.491015   0.226397   2.169
## HS.c:competence.c:Can_AusUK                 -0.072288   0.212499  -0.340
## HS.c:competence.c:Aus_UK                    -0.527126   0.404055  -1.305
## BS.c:warmth.c:US_oth                        -0.021463   0.106748  -0.201
## BS.c:warmth.c:Can_AusUK                     -0.140728   0.142505  -0.988
## BS.c:warmth.c:Aus_UK                        -0.078694   0.223664  -0.352
## BS.c:competence.c:US_oth                    -0.501192   0.353837  -1.416
## BS.c:competence.c:Can_AusUK                  0.196793   0.213019   0.924
## BS.c:competence.c:Aus_UK                     0.143078   0.395255   0.362
## SRQ.c:warmth.c:US_oth                       -0.004001   0.007286  -0.549
## SRQ.c:warmth.c:Can_AusUK                     0.021022   0.012378   1.698
## SRQ.c:warmth.c:Aus_UK                        0.010898   0.015926   0.684
## SRQ.c:competence.c:US_oth                    0.012108   0.013580   0.892
## SRQ.c:competence.c:Can_AusUK                 0.025602   0.014609   1.752
## SRQ.c:competence.c:Aus_UK                    0.040010   0.025752   1.554
## SDO.c:warmth.c:US_oth                       -0.053354   0.074918  -0.712
## SDO.c:warmth.c:Can_AusUK                     0.056411   0.116554   0.484
## SDO.c:warmth.c:Aus_UK                       -0.273294   0.169018  -1.617
## SDO.c:competence.c:US_oth                   -0.051231   0.116215  -0.441
## SDO.c:competence.c:Can_AusUK                -0.418453   0.181472  -2.306
## SDO.c:competence.c:Aus_UK                    0.167761   0.279069   0.601
## Ideology.c:warmth.c:US_oth                  -0.003985   0.061601  -0.065
## Ideology.c:warmth.c:Can_AusUK               -0.028270   0.085260  -0.332
## Ideology.c:warmth.c:Aus_UK                   0.052819   0.146266   0.361
## Ideology.c:competence.c:US_oth              -0.009394   0.115326  -0.081
## Ideology.c:competence.c:Can_AusUK            0.253677   0.155049   1.636
## Ideology.c:competence.c:Aus_UK               0.076342   0.223268   0.342
## HS.c:man_woman:US_oth                       -0.328005   0.273914  -1.197
## HS.c:man_woman:Can_AusUK                     0.432134   0.281665   1.534
## HS.c:man_woman:Aus_UK                       -0.383786   0.456073  -0.842
## BS.c:man_woman:US_oth                        0.618835   0.442789   1.398
## BS.c:man_woman:Can_AusUK                    -0.146415   0.272265  -0.538
## BS.c:man_woman:Aus_UK                       -0.904253   0.464611  -1.946
## SRQ.c:man_woman:US_oth                      -0.018464   0.020071  -0.920
## SRQ.c:man_woman:Can_AusUK                   -0.026403   0.022654  -1.165
## SRQ.c:man_woman:Aus_UK                       0.041261   0.030783   1.340
## SDO.c:man_woman:US_oth                      -0.125695   0.161277  -0.779
## SDO.c:man_woman:Can_AusUK                    0.055359   0.242299   0.228
## SDO.c:man_woman:Aus_UK                       0.438142   0.341047   1.285
## Ideology.c:man_woman:US_oth                 -0.225506   0.177927  -1.267
## Ideology.c:man_woman:Can_AusUK              -0.361719   0.206375  -1.753
## Ideology.c:man_woman:Aus_UK                 -0.185331   0.323379  -0.573
## warmth.c:man_woman:US_oth                   -0.216485   0.180670  -1.198
## warmth.c:man_woman:Can_AusUK                 0.195885   0.235493   0.832
## warmth.c:man_woman:Aus_UK                    0.458115   0.378902   1.209
## competence.c:man_woman:US_oth                0.179076   0.348075   0.514
## competence.c:man_woman:Can_AusUK            -0.138624   0.345741  -0.401
## competence.c:man_woman:Aus_UK                0.884653   0.557302   1.587
## HS.c:warmth.c:man_woman:US_oth               0.034205   0.204660   0.167
## HS.c:warmth.c:man_woman:Can_AusUK           -0.190489   0.294934  -0.646
## HS.c:warmth.c:man_woman:Aus_UK               0.946586   0.431367   2.194
## HS.c:competence.c:man_woman:US_oth           0.116581   0.452418   0.258
## HS.c:competence.c:man_woman:Can_AusUK       -0.522314   0.424362  -1.231
## HS.c:competence.c:man_woman:Aus_UK           1.579387   0.807190   1.957
## BS.c:warmth.c:man_woman:US_oth               0.778461   0.213700   3.643
## BS.c:warmth.c:man_woman:Can_AusUK           -0.425588   0.284961  -1.493
## BS.c:warmth.c:man_woman:Aus_UK              -0.258998   0.448291  -0.578
## BS.c:competence.c:man_woman:US_oth          -1.605976   0.707446  -2.270
## BS.c:competence.c:man_woman:Can_AusUK       -0.077058   0.426045  -0.181
## BS.c:competence.c:man_woman:Aus_UK          -0.201938   0.789780  -0.256
## SRQ.c:warmth.c:man_woman:US_oth             -0.052441   0.014530  -3.609
## SRQ.c:warmth.c:man_woman:Can_AusUK           0.074184   0.024748   2.998
## SRQ.c:warmth.c:man_woman:Aus_UK             -0.045054   0.031878  -1.413
## SRQ.c:competence.c:man_woman:US_oth          0.033750   0.027134   1.244
## SRQ.c:competence.c:man_woman:Can_AusUK       0.023643   0.029238   0.809
## SRQ.c:competence.c:man_woman:Aus_UK         -0.060168   0.051530  -1.168
## SDO.c:warmth.c:man_woman:US_oth              0.468748   0.149763   3.130
## SDO.c:warmth.c:man_woman:Can_AusUK          -0.514723   0.233743  -2.202
## SDO.c:warmth.c:man_woman:Aus_UK              0.133657   0.338184   0.395
## SDO.c:competence.c:man_woman:US_oth          0.078889   0.233238   0.338
## SDO.c:competence.c:man_woman:Can_AusUK      -0.230794   0.362576  -0.637
## SDO.c:competence.c:man_woman:Aus_UK         -0.114143   0.558777  -0.204
## Ideology.c:warmth.c:man_woman:US_oth        -0.251796   0.123414  -2.040
## Ideology.c:warmth.c:man_woman:Can_AusUK      0.275568   0.170486   1.616
## Ideology.c:warmth.c:man_woman:Aus_UK         0.674084   0.292522   2.304
## Ideology.c:competence.c:man_woman:US_oth     0.177384   0.231543   0.766
## Ideology.c:competence.c:man_woman:Can_AusUK  0.298249   0.310910   0.959
## Ideology.c:competence.c:man_woman:Aus_UK     0.272896   0.447761   0.609
##                                             Pr(>|z|)    
## (Intercept)                                 0.062042 .  
## HS.c                                        7.44e-07 ***
## BS.c                                        0.032333 *  
## SRQ.c                                       2.52e-05 ***
## SDO.c                                       0.074286 .  
## Ideology.c                                   < 2e-16 ***
## warmth.c                                    0.007143 ** 
## competence.c                                0.027517 *  
## man_woman                                   0.492094    
## US_oth                                      0.983826    
## Can_AusUK                                   0.011796 *  
## Aus_UK                                      0.145332    
## age.c                                       0.502643    
## edu.c                                       0.575213    
## HS.c:warmth.c                               0.268757    
## HS.c:competence.c                           0.116227    
## BS.c:warmth.c                               0.118465    
## BS.c:competence.c                           0.089097 .  
## SRQ.c:warmth.c                              0.557860    
## SRQ.c:competence.c                          0.302634    
## SDO.c:warmth.c                              0.921487    
## SDO.c:competence.c                          0.249435    
## Ideology.c:warmth.c                         0.168578    
## Ideology.c:competence.c                     0.043095 *  
## HS.c:man_woman                              0.120466    
## BS.c:man_woman                              0.192471    
## SRQ.c:man_woman                             0.002018 ** 
## SDO.c:man_woman                             0.788360    
## Ideology.c:man_woman                        0.888849    
## warmth.c:man_woman                          0.013037 *  
## competence.c:man_woman                      0.624662    
## HS.c:US_oth                                 0.274434    
## HS.c:Can_AusUK                              0.867752    
## HS.c:Aus_UK                                 0.297171    
## BS.c:US_oth                                 0.029013 *  
## BS.c:Can_AusUK                              0.468616    
## BS.c:Aus_UK                                 0.799838    
## SRQ.c:US_oth                                0.289940    
## SRQ.c:Can_AusUK                             0.025502 *  
## SRQ.c:Aus_UK                                0.366965    
## SDO.c:US_oth                                0.034599 *  
## SDO.c:Can_AusUK                             0.233338    
## SDO.c:Aus_UK                                0.359628    
## Ideology.c:US_oth                           0.045750 *  
## Ideology.c:Can_AusUK                        0.762801    
## Ideology.c:Aus_UK                           0.591130    
## warmth.c:US_oth                             0.398535    
## warmth.c:Can_AusUK                          0.939895    
## warmth.c:Aus_UK                             0.173098    
## competence.c:US_oth                         0.047395 *  
## competence.c:Can_AusUK                      0.296324    
## competence.c:Aus_UK                         0.213984    
## man_woman:US_oth                            0.988318    
## man_woman:Can_AusUK                         0.342907    
## man_woman:Aus_UK                            0.298641    
## HS.c:warmth.c:man_woman                     0.791126    
## HS.c:competence.c:man_woman                 0.572567    
## BS.c:warmth.c:man_woman                     0.678986    
## BS.c:competence.c:man_woman                 0.567869    
## SRQ.c:warmth.c:man_woman                    0.064437 .  
## SRQ.c:competence.c:man_woman                0.034550 *  
## SDO.c:warmth.c:man_woman                    0.081934 .  
## SDO.c:competence.c:man_woman                0.232312    
## Ideology.c:warmth.c:man_woman               0.285025    
## Ideology.c:competence.c:man_woman           0.953827    
## HS.c:warmth.c:US_oth                        0.007872 ** 
## HS.c:warmth.c:Can_AusUK                     0.017246 *  
## HS.c:warmth.c:Aus_UK                        0.544002    
## HS.c:competence.c:US_oth                    0.030096 *  
## HS.c:competence.c:Can_AusUK                 0.733719    
## HS.c:competence.c:Aus_UK                    0.192033    
## BS.c:warmth.c:US_oth                        0.840651    
## BS.c:warmth.c:Can_AusUK                     0.323382    
## BS.c:warmth.c:Aus_UK                        0.724958    
## BS.c:competence.c:US_oth                    0.156644    
## BS.c:competence.c:Can_AusUK                 0.355575    
## BS.c:competence.c:Aus_UK                    0.717361    
## SRQ.c:warmth.c:US_oth                       0.582953    
## SRQ.c:warmth.c:Can_AusUK                    0.089448 .  
## SRQ.c:warmth.c:Aus_UK                       0.493808    
## SRQ.c:competence.c:US_oth                   0.372596    
## SRQ.c:competence.c:Can_AusUK                0.079696 .  
## SRQ.c:competence.c:Aus_UK                   0.120270    
## SDO.c:warmth.c:US_oth                       0.476362    
## SDO.c:warmth.c:Can_AusUK                    0.628392    
## SDO.c:warmth.c:Aus_UK                       0.105888    
## SDO.c:competence.c:US_oth                   0.659340    
## SDO.c:competence.c:Can_AusUK                0.021117 *  
## SDO.c:competence.c:Aus_UK                   0.547744    
## Ideology.c:warmth.c:US_oth                  0.948424    
## Ideology.c:warmth.c:Can_AusUK               0.740213    
## Ideology.c:warmth.c:Aus_UK                  0.718011    
## Ideology.c:competence.c:US_oth              0.935078    
## Ideology.c:competence.c:Can_AusUK           0.101817    
## Ideology.c:competence.c:Aus_UK              0.732405    
## HS.c:man_woman:US_oth                       0.231123    
## HS.c:man_woman:Can_AusUK                    0.124977    
## HS.c:man_woman:Aus_UK                       0.400067    
## BS.c:man_woman:US_oth                       0.162238    
## BS.c:man_woman:Can_AusUK                    0.590738    
## BS.c:man_woman:Aus_UK                       0.051624 .  
## SRQ.c:man_woman:US_oth                      0.357591    
## SRQ.c:man_woman:Can_AusUK                   0.243830    
## SRQ.c:man_woman:Aus_UK                      0.180117    
## SDO.c:man_woman:US_oth                      0.435762    
## SDO.c:man_woman:Can_AusUK                   0.819277    
## SDO.c:man_woman:Aus_UK                      0.198899    
## Ideology.c:man_woman:US_oth                 0.205011    
## Ideology.c:man_woman:Can_AusUK              0.079649 .  
## Ideology.c:man_woman:Aus_UK                 0.566572    
## warmth.c:man_woman:US_oth                   0.230824    
## warmth.c:man_woman:Can_AusUK                0.405518    
## warmth.c:man_woman:Aus_UK                   0.226640    
## competence.c:man_woman:US_oth               0.606920    
## competence.c:man_woman:Can_AusUK            0.688459    
## competence.c:man_woman:Aus_UK               0.112425    
## HS.c:warmth.c:man_woman:US_oth              0.867268    
## HS.c:warmth.c:man_woman:Can_AusUK           0.518364    
## HS.c:warmth.c:man_woman:Aus_UK              0.028208 *  
## HS.c:competence.c:man_woman:US_oth          0.796651    
## HS.c:competence.c:man_woman:Can_AusUK       0.218390    
## HS.c:competence.c:man_woman:Aus_UK          0.050389 .  
## BS.c:warmth.c:man_woman:US_oth              0.000270 ***
## BS.c:warmth.c:man_woman:Can_AusUK           0.135307    
## BS.c:warmth.c:man_woman:Aus_UK              0.563436    
## BS.c:competence.c:man_woman:US_oth          0.023201 *  
## BS.c:competence.c:man_woman:Can_AusUK       0.856470    
## BS.c:competence.c:man_woman:Aus_UK          0.798191    
## SRQ.c:warmth.c:man_woman:US_oth             0.000307 ***
## SRQ.c:warmth.c:man_woman:Can_AusUK          0.002721 ** 
## SRQ.c:warmth.c:man_woman:Aus_UK             0.157559    
## SRQ.c:competence.c:man_woman:US_oth         0.213563    
## SRQ.c:competence.c:man_woman:Can_AusUK      0.418715    
## SRQ.c:competence.c:man_woman:Aus_UK         0.242960    
## SDO.c:warmth.c:man_woman:US_oth             0.001748 ** 
## SDO.c:warmth.c:man_woman:Can_AusUK          0.027659 *  
## SDO.c:warmth.c:man_woman:Aus_UK             0.692681    
## SDO.c:competence.c:man_woman:US_oth         0.735187    
## SDO.c:competence.c:man_woman:Can_AusUK      0.524425    
## SDO.c:competence.c:man_woman:Aus_UK         0.838140    
## Ideology.c:warmth.c:man_woman:US_oth        0.041325 *  
## Ideology.c:warmth.c:man_woman:Can_AusUK     0.106015    
## Ideology.c:warmth.c:man_woman:Aus_UK        0.021201 *  
## Ideology.c:competence.c:man_woman:US_oth    0.443620    
## Ideology.c:competence.c:man_woman:Can_AusUK 0.337419    
## Ideology.c:competence.c:man_woman:Aus_UK    0.542215    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 1326.79  on 977  degrees of freedom
## Residual deviance:  708.69  on 832  degrees of freedom
##   (17 observations deleted due to missingness)
## AIC: 1000.7
## 
## Number of Fisher Scoring iterations: 8
tab_model(s3_glm3,
          show.stat = T,
          show.se = T)
  Voted Harris
Predictors Odds Ratios std. Error CI Statistic p
(Intercept) 1.30 0.19 0.99 – 1.73 1.87 0.062
HS c 0.38 0.07 0.26 – 0.55 -4.95 <0.001
BS c 1.74 0.45 1.10 – 3.09 2.14 0.032
SRQ c 0.94 0.01 0.91 – 0.97 -4.21 <0.001
SDO c 1.28 0.18 0.98 – 1.68 1.78 0.074
Ideology c 0.34 0.04 0.25 – 0.43 -8.23 <0.001
warmth c 1.48 0.21 1.11 – 1.97 2.69 0.007
competence c 0.59 0.14 0.36 – 0.94 -2.20 0.028
man woman 1.22 0.35 0.69 – 2.13 0.69 0.492
US oth 1.00 0.10 0.83 – 1.21 -0.02 0.984
Can AusUK 1.32 0.14 1.07 – 1.65 2.52 0.012
Aus UK 1.29 0.22 0.92 – 1.83 1.46 0.145
age c 1.00 0.00 NA – 1.00 -0.67 0.503
edu c 1.02 0.04 0.95 – 1.11 0.56 0.575
HS c × warmth c 1.21 0.21 0.87 – 1.70 1.11 0.269
HS c × competence c 1.65 0.52 0.90 – 3.14 1.57 0.116
BS c × warmth c 0.76 0.13 0.54 – 1.07 -1.56 0.118
BS c × competence c 0.49 0.21 0.20 – 1.04 -1.70 0.089
SRQ c × warmth c 1.01 0.01 0.98 – 1.03 0.59 0.558
SRQ c × competence c 1.02 0.02 0.98 – 1.06 1.03 0.303
SDO c × warmth c 1.01 0.13 0.78 – 1.31 0.10 0.921
SDO c × competence c 0.79 0.16 0.52 – 1.19 -1.15 0.249
Ideology c × warmth c 1.16 0.12 0.94 – 1.42 1.38 0.169
Ideology c × competence c 1.45 0.27 1.02 – 2.11 2.02 0.043
HS c × man woman 1.82 0.71 0.84 – 3.87 1.55 0.120
BS c × man woman 1.97 1.02 0.78 – 6.19 1.30 0.192
SRQ c × man woman 0.92 0.03 0.86 – 0.97 -3.09 0.002
SDO c × man woman 0.93 0.25 0.54 – 1.59 -0.27 0.788
Ideology c × man woman 1.04 0.28 0.61 – 1.74 0.14 0.889
warmth c × man woman 0.49 0.14 0.27 – 0.85 -2.48 0.013
competence c × man woman 0.79 0.38 0.30 – 2.02 -0.49 0.625
HS c × US oth 0.86 0.12 0.64 – 1.11 -1.09 0.274
HS c × Can AusUK 1.02 0.14 0.77 – 1.35 0.17 0.868
HS c × Aus UK 1.27 0.29 0.81 – 2.00 1.04 0.297
BS c × US oth 1.62 0.36 1.13 – 2.72 2.18 0.029
BS c × Can AusUK 0.91 0.12 0.69 – 1.19 -0.72 0.469
BS c × Aus UK 1.06 0.25 0.67 – 1.68 0.25 0.800
SRQ c × US oth 0.99 0.01 0.97 – 1.01 -1.06 0.290
SRQ c × Can AusUK 0.98 0.01 0.95 – 1.00 -2.23 0.026
SRQ c × Aus UK 0.99 0.02 0.96 – 1.02 -0.90 0.367
SDO c × US oth 0.84 0.07 0.72 – 0.98 -2.11 0.035
SDO c × Can AusUK 1.16 0.14 0.92 – 1.48 1.19 0.233
SDO c × Aus UK 0.86 0.15 0.61 – 1.19 -0.92 0.360
Ideology c × US oth 0.84 0.07 0.69 – 0.98 -2.00 0.046
Ideology c × Can AusUK 0.97 0.10 0.78 – 1.18 -0.30 0.763
Ideology c × Aus UK 1.09 0.18 0.80 – 1.51 0.54 0.591
warmth c × US oth 0.93 0.08 0.77 – 1.11 -0.84 0.399
warmth c × Can AusUK 0.99 0.12 0.79 – 1.25 -0.08 0.940
warmth c × Aus UK 0.77 0.15 0.53 – 1.11 -1.36 0.173
competence c × US oth 0.71 0.12 0.49 – 0.99 -1.98 0.047
competence c × Can AusUK 0.83 0.14 0.59 – 1.17 -1.04 0.296
competence c × Aus UK 0.71 0.20 0.40 – 1.22 -1.24 0.214
man woman × US oth 1.00 0.19 0.69 – 1.47 0.01 0.988
man woman × Can AusUK 0.81 0.18 0.52 – 1.25 -0.95 0.343
man woman × Aus UK 0.70 0.24 0.35 – 1.37 -1.04 0.299
(HS c × warmth c) × man
woman
1.09 0.37 0.56 – 2.17 0.26 0.791
(HS c × competence c) ×
man woman
0.70 0.44 0.21 – 2.52 -0.56 0.573
(BS c × warmth c) × man
woman
0.87 0.30 0.44 – 1.72 -0.41 0.679
(BS c × competence c) ×
man woman
0.62 0.52 0.10 – 2.80 -0.57 0.568
(SRQ c × warmth c) × man
woman
1.05 0.03 1.00 – 1.11 1.85 0.064
(SRQ c × competence c) ×
man woman
1.09 0.04 1.01 – 1.18 2.11 0.035
(SDO c × warmth c) × man
woman
0.63 0.17 0.38 – 1.05 -1.74 0.082
(SDO c × competence c) ×
man woman
0.61 0.25 0.27 – 1.36 -1.19 0.232
(Ideology c × warmth c) ×
man woman
0.80 0.17 0.53 – 1.20 -1.07 0.285
(Ideology c × competence
c) × man woman
0.98 0.36 0.48 – 2.04 -0.06 0.954
(HS c × warmth c) × US
oth
1.31 0.13 1.09 – 1.64 2.66 0.008
(HS c × warmth c) × Can
AusUK
0.70 0.10 0.52 – 0.93 -2.38 0.017
(HS c × warmth c) × Aus
UK
0.88 0.19 0.57 – 1.34 -0.61 0.544
(HS c × competence c) ×
US oth
1.63 0.37 1.08 – 2.66 2.17 0.030
(HS c × competence c) ×
Can AusUK
0.93 0.20 0.61 – 1.42 -0.34 0.734
(HS c × competence c) ×
Aus UK
0.59 0.24 0.26 – 1.28 -1.30 0.192
(BS c × warmth c) × US
oth
0.98 0.10 0.79 – 1.21 -0.20 0.841
(BS c × warmth c) × Can
AusUK
0.87 0.12 0.65 – 1.15 -0.99 0.323
(BS c × warmth c) × Aus
UK
0.92 0.21 0.59 – 1.42 -0.35 0.725
(BS c × competence c) ×
US oth
0.61 0.21 0.27 – 1.10 -1.42 0.157
(BS c × competence c) ×
Can AusUK
1.22 0.26 0.79 – 1.84 0.92 0.356
(BS c × competence c) ×
Aus UK
1.15 0.46 0.53 – 2.53 0.36 0.717
(SRQ c × warmth c) × US
oth
1.00 0.01 0.98 – 1.01 -0.55 0.583
(SRQ c × warmth c) × Can
AusUK
1.02 0.01 1.00 – 1.05 1.70 0.089
(SRQ c × warmth c) × Aus
UK
1.01 0.02 0.98 – 1.04 0.68 0.494
(SRQ c × competence c) ×
US oth
1.01 0.01 0.99 – 1.04 0.89 0.373
(SRQ c × competence c) ×
Can AusUK
1.03 0.01 1.00 – 1.06 1.75 0.080
(SRQ c × competence c) ×
Aus UK
1.04 0.03 0.99 – 1.10 1.55 0.120
(SDO c × warmth c) × US
oth
0.95 0.07 0.82 – 1.10 -0.71 0.476
(SDO c × warmth c) × Can
AusUK
1.06 0.12 0.84 – 1.33 0.48 0.628
(SDO c × warmth c) × Aus
UK
0.76 0.13 0.54 – 1.05 -1.62 0.106
(SDO c × competence c) ×
US oth
0.95 0.11 0.75 – 1.19 -0.44 0.659
(SDO c × competence c) ×
Can AusUK
0.66 0.12 0.45 – 0.92 -2.31 0.021
(SDO c × competence c) ×
Aus UK
1.18 0.33 0.66 – 2.02 0.60 0.548
(Ideology c × warmth c) ×
US oth
1.00 0.06 0.88 – 1.13 -0.06 0.948
(Ideology c × warmth c) ×
Can AusUK
0.97 0.08 0.82 – 1.15 -0.33 0.740
(Ideology c × warmth c) ×
Aus UK
1.05 0.15 0.79 – 1.41 0.36 0.718
(Ideology c × competence
c) × US oth
0.99 0.11 0.80 – 1.27 -0.08 0.935
(Ideology c × competence
c) × Can AusUK
1.29 0.20 0.95 – 1.76 1.64 0.102
(Ideology c × competence
c) × Aus UK
1.08 0.24 0.70 – 1.68 0.34 0.732
(HS c × man woman) × US
oth
0.72 0.20 0.40 – 1.20 -1.20 0.231
(HS c × man woman) × Can
AusUK
1.54 0.43 0.90 – 2.73 1.53 0.125
(HS c × man woman) × Aus
UK
0.68 0.31 0.28 – 1.67 -0.84 0.400
(BS c × man woman) × US
oth
1.86 0.82 0.88 – 5.22 1.40 0.162
(BS c × man woman) × Can
AusUK
0.86 0.24 0.50 – 1.48 -0.54 0.591
(BS c × man woman) × Aus
UK
0.40 0.19 0.16 – 0.99 -1.95 0.052
(SRQ c × man woman) × US
oth
0.98 0.02 0.94 – 1.02 -0.92 0.358
(SRQ c × man woman) × Can
AusUK
0.97 0.02 0.93 – 1.02 -1.17 0.244
(SRQ c × man woman) × Aus
UK
1.04 0.03 0.98 – 1.11 1.34 0.180
(SDO c × man woman) × US
oth
0.88 0.14 0.64 – 1.21 -0.78 0.436
(SDO c × man woman) × Can
AusUK
1.06 0.26 0.65 – 1.71 0.23 0.819
(SDO c × man woman) × Aus
UK
1.55 0.53 0.80 – 3.06 1.28 0.199
(Ideology c × man woman)
× US oth
0.80 0.14 0.54 – 1.11 -1.27 0.205
(Ideology c × man woman)
× Can AusUK
0.70 0.14 0.46 – 1.04 -1.75 0.080
(Ideology c × man woman)
× Aus UK
0.83 0.27 0.43 – 1.55 -0.57 0.567
(warmth c × man woman) ×
US oth
0.81 0.15 0.56 – 1.13 -1.20 0.231
(warmth c × man woman) ×
Can AusUK
1.22 0.29 0.76 – 1.93 0.83 0.406
(warmth c × man woman) ×
Aus UK
1.58 0.60 0.76 – 3.41 1.21 0.227
(competence c × man
woman) × US oth
1.20 0.42 0.58 – 2.35 0.51 0.607
(competence c × man
woman) × Can AusUK
0.87 0.30 0.44 – 1.71 -0.40 0.688
(competence c × man
woman) × Aus UK
2.42 1.35 0.82 – 7.38 1.59 0.112
(HS c × warmth c × man
woman) × US oth
1.03 0.21 0.71 – 1.61 0.17 0.867
(HS c × warmth c × man
woman) × Can AusUK
0.83 0.24 0.46 – 1.47 -0.65 0.518
(HS c × warmth c × man
woman) × Aus UK
2.58 1.11 1.13 – 6.17 2.19 0.028
(HS c × competence c ×
man woman) × US oth
1.12 0.51 0.48 – 2.95 0.26 0.797
(HS c × competence c ×
man woman) × Can AusUK
0.59 0.25 0.25 – 1.35 -1.23 0.218
(HS c × competence c ×
man woman) × Aus UK
4.85 3.92 1.04 – 25.03 1.96 0.050
(BS c × warmth c × man
woman) × US oth
2.18 0.47 1.46 – 3.43 3.64 <0.001
(BS c × warmth c × man
woman) × Can AusUK
0.65 0.19 0.37 – 1.13 -1.49 0.135
(BS c × warmth c × man
woman) × Aus UK
0.77 0.35 0.32 – 1.88 -0.58 0.563
(BS c × competence c ×
man woman) × US oth
0.20 0.14 0.04 – 0.66 -2.27 0.023
(BS c × competence c ×
man woman) × Can AusUK
0.93 0.39 0.40 – 2.17 -0.18 0.856
(BS c × competence c ×
man woman) × Aus UK
0.82 0.65 0.17 – 3.87 -0.26 0.798
(SRQ c × warmth c × man
woman) × US oth
0.95 0.01 0.92 – 0.98 -3.61 <0.001
(SRQ c × warmth c × man
woman) × Can AusUK
1.08 0.03 1.03 – 1.13 3.00 0.003
(SRQ c × warmth c × man
woman) × Aus UK
0.96 0.03 0.90 – 1.02 -1.41 0.158
(SRQ c × competence c ×
man woman) × US oth
1.03 0.03 0.98 – 1.10 1.24 0.214
(SRQ c × competence c ×
man woman) × Can AusUK
1.02 0.03 0.97 – 1.08 0.81 0.419
(SRQ c × competence c ×
man woman) × Aus UK
0.94 0.05 0.85 – 1.04 -1.17 0.243
(SDO c × warmth c × man
woman) × US oth
1.60 0.24 1.20 – 2.17 3.13 0.002
(SDO c × warmth c × man
woman) × Can AusUK
0.60 0.14 0.37 – 0.93 -2.20 0.028
(SDO c × warmth c × man
woman) × Aus UK
1.14 0.39 0.59 – 2.26 0.40 0.693
(SDO c × competence c ×
man woman) × US oth
1.08 0.25 0.68 – 1.71 0.34 0.735
(SDO c × competence c ×
man woman) × Can AusUK
0.79 0.29 0.40 – 1.65 -0.64 0.524
(SDO c × competence c ×
man woman) × Aus UK
0.89 0.50 0.30 – 2.85 -0.20 0.838
(Ideology c × warmth c ×
man woman) × US oth
0.78 0.10 0.60 – 0.98 -2.04 0.041
(Ideology c × warmth c ×
man woman) × Can AusUK
1.32 0.22 0.95 – 1.85 1.62 0.106
(Ideology c × warmth c ×
man woman) × Aus UK
1.96 0.57 1.11 – 3.52 2.30 0.021
(Ideology c × competence
c × man woman) × US oth
1.19 0.28 0.77 – 1.95 0.77 0.444
(Ideology c × competence
c × man woman) × Can
AusUK
1.35 0.42 0.74 – 2.52 0.96 0.337
(Ideology c × competence
c × man woman) × Aus UK
1.31 0.59 0.55 – 3.19 0.61 0.542
Observations 978
R2 Tjur 0.518

Data Visualization

d3 %>% ggplot(aes(x = Country, 
                  y = VotedHarris_num, 
                  fill = Demo_Gender)) +
  geom_bar(stat = "summary", 
           fun = "mean", 
           position = position_dodge(.9), 
           alpha = .8) +
  stat_summary(geom = "errorbar", 
               color = "black", 
               width = .1, 
               show.legend = F, 
               position = position_dodge(.9)) +
  theme_bw() +
  scale_y_continuous(breaks = seq(0,1,.2)) +
  coord_cartesian(ylim = c(0,1)) +
  scale_fill_manual("Gender", 
                    values = c("#69b3a2", "#404080")) +
  labs(y = "Proportion Voted for Harris",
       title = "Vote Choice by Gender and Country")

d3 %>% ggplot(aes(x = HS, y = VotedHarris_num, group = Demo_Gender, color = Demo_Gender)) +
  geom_jitter(size = 1, alpha = .7) +
#  geom_smooth(method = "loess", se = F, fullrange = T) +
  geom_smooth(method = "lm", se = F, fullrange = T, fullrange = T) +
  theme_bw() +
#  facet_wrap(~Country) +
  scale_fill_manual(values = c("#69b3a2", "#404080"))+
  scale_color_manual(values = c("#69b3a2", "#404080")) +
  labs(fill = "Gender", color = "Gender",
       y = "Voted for Harris",
       x = "Hostile Sexism",
       title = "Likelihood of Voting Harris by Hostile Sexism, Gender, and Country") +
  scale_y_continuous(breaks = c(0,1),
                     labels = c("No", 
                              "Yes")) +
  scale_x_continuous(breaks = seq(0,7,1)) +
  coord_cartesian(xlim = c(0,5))

d3 %>% ggplot(aes(x = SRQ, y = VotedHarris_num, group = Demo_Gender, color = Demo_Gender)) +
  geom_jitter(size = 1, alpha = .7) +
#  geom_smooth(method = "loess", se = F, fullrange = T) +
  geom_smooth(method = "lm", se = F, fullrange = T,  fullrange = T) +
  theme_bw() +
#  facet_wrap(~Country) +
  scale_fill_manual(values = c("#69b3a2", "#404080"))+
  scale_color_manual(values = c("#69b3a2", "#404080")) +
  labs(fill = "Gender", color = "Gender",
       y = "Voted for Harris",
       x = "Social Roles Questionnaire",
       title = "Likelihood of Voting Harris by SRQ, Gender, and Country") +
  scale_y_continuous(breaks = c(0,1),
                     labels = c("No", 
                              "Yes")) +
  scale_x_continuous(breaks = seq(0,90,10)) +
  coord_cartesian(xlim = c(0,90))

d3 %>% ggplot(aes(x = Ideology, y = VotedHarris_num, group = Demo_Gender, color = Demo_Gender)) +
  geom_jitter(size = 1, alpha = .7) +
#  geom_smooth(method = "loess", se = F, fullrange = T) +
  geom_smooth(method = "lm", se = F, fullrange = T,  fullrange = T) +
  theme_bw() +
  facet_wrap(~Country) +
  scale_fill_manual(values = c("#69b3a2", "#404080"))+
  scale_color_manual(values = c("#69b3a2", "#404080")) +
  labs(fill = "Gender", color = "Gender",
       y = "Voted for Harris",
       x = "Ideology (Liberal to Conseravtive)",
       title = "Likelihood of Voting Harris by Ideology, Gender, and Country") +
  scale_y_continuous(breaks = c(0,1),
                     labels = c("No", 
                              "Yes")) +
  scale_x_continuous(breaks = seq(1,7,1)) +
  coord_cartesian(xlim = c(1,7))

Analyses: Candidate Experiment

  • Note on analyses: These models are ordinary least squares logistic regressions, evaluating how predictors affect likelihood of moving from one candidate category to the next
    • Categories are ordered as: 1) Masculine man; 2) Feminine man; 3) Masculine woman; 4) Feminine woman
  • Model results will be presented as odds ratios
    • ORs greater than 1 indicate an association with increasing along the candidate category spectrum, toward “feminine woman”
    • ORs less than 1 indicate an association with decreasing along the candidate category spectrum, toward “masculine man”
    • (the “Statistic” column in each model table is the Z-score assoc. with each OR)
  • Models control for age & education
  • Continuous predictors are mean-centered

Codes

d3$US_oth <- NA
d3$US_oth[d3$Country == "US"] <- 3
d3$US_oth[d3$Country == "Canada"] <- -1
d3$US_oth[d3$Country == "Aus"] <- -1
d3$US_oth[d3$Country == "UK"] <- -1

d3$Can_AusUK <- NA
d3$Can_AusUK[d3$Country == "US"] <- 0
d3$Can_AusUK[d3$Country == "Canada"] <- 2
d3$Can_AusUK[d3$Country == "Aus"] <- -1
d3$Can_AusUK[d3$Country == "UK"] <- -1

d3$Aus_UK <- NA
d3$Aus_UK[d3$Country == "US"] <- 0
d3$Aus_UK[d3$Country == "Canada"] <- 0
d3$Aus_UK[d3$Country == "Aus"] <- -1
d3$Aus_UK[d3$Country == "UK"] <- 1

d3$Country <- factor(d3$Country, levels = c("US","Canada", "Aus", "UK"))
contrasts(d3$Country) <- contr.helmert(nlevels(d3$Country))

d3$man_woman <- ifelse(d3$Demo_Gender == "Woman", 1/2, -1/2)


contrasts(d3$Demo_Race) <- contr.poly(nlevels(d3$Demo_Race))

d3$LikelihoodofVoting <- as.factor(d3$LikelihoodofVoting)

Models

Interactive models

s3_olr1 <- clm(LikelihoodofVoting ~ man_woman * (US_oth + Can_AusUK + Aus_UK)
               + age.c + edu.c,
               data = d3,
               link = "logit"
)

# summary(s3_olr1)
tab_model(s3_olr1,
          show.se = T,
          show.stat = T)
  Likelihoodof Voting
Predictors Odds Ratios std. Error CI Statistic p
1|2 0.23 0.02 0.19 – 0.27 -17.87 <0.001
2|3 0.58 0.04 0.51 – 0.66 -8.11 <0.001
3|4 1.75 0.12 1.54 – 2.00 8.26 <0.001
man woman 1.81 0.21 1.44 – 2.29 5.02 <0.001
US oth 0.99 0.03 0.93 – 1.06 -0.18 0.858
Can AusUK 0.87 0.04 0.79 – 0.95 -3.04 0.002
Aus UK 0.80 0.07 0.68 – 0.95 -2.56 0.010
age c 1.00 0.00 1.00 – 1.00 -1.17 0.242
edu c 1.09 0.02 1.05 – 1.14 4.06 <0.001
man woman × US oth 0.91 0.06 0.80 – 1.04 -1.35 0.177
man woman × Can AusUK 0.99 0.09 0.82 – 1.20 -0.06 0.949
man woman × Aus UK 1.09 0.19 0.78 – 1.52 0.50 0.617
Observations 980
R2 Nagelkerke 0.067
s3_olr2 <- clm(LikelihoodofVoting ~ man_woman * (US_oth + Can_AusUK + Aus_UK) * (HS.c + BS.c + SDO.c + SRQ.c)
               + age.c + edu.c,
               data = d3,
               link = "logit"
)

# summary(s3_olr2)
tab_model(s3_olr2,
          show.stat = T,
          show.se=T)
  Likelihoodof Voting
Predictors Odds Ratios std. Error CI Statistic p
1|2 0.17 0.02 0.14 – 0.21 -18.44 <0.001
2|3 0.53 0.04 0.46 – 0.62 -8.21 <0.001
3|4 2.02 0.16 1.74 – 2.36 9.15 <0.001
man woman 1.21 0.16 0.94 – 1.55 1.44 0.149
US oth 0.95 0.04 0.88 – 1.02 -1.38 0.167
Can AusUK 0.87 0.05 0.78 – 0.96 -2.72 0.007
Aus UK 0.88 0.08 0.73 – 1.05 -1.40 0.161
HS c 0.62 0.05 0.53 – 0.73 -5.83 <0.001
BS c 1.14 0.10 0.96 – 1.35 1.47 0.140
SDO c 0.73 0.05 0.64 – 0.82 -5.00 <0.001
SRQ c 0.98 0.01 0.97 – 1.00 -2.63 0.009
age c 1.00 0.00 1.00 – 1.00 0.00 0.999
edu c 1.07 0.03 1.02 – 1.12 2.84 0.004
man woman × US oth 0.89 0.07 0.77 – 1.03 -1.57 0.117
man woman × Can AusUK 1.05 0.11 0.85 – 1.29 0.43 0.665
man woman × Aus UK 1.16 0.21 0.81 – 1.67 0.80 0.426
man woman × HS c 0.73 0.12 0.53 – 1.01 -1.91 0.056
man woman × BS c 0.87 0.15 0.62 – 1.22 -0.81 0.420
man woman × SDO c 1.13 0.14 0.88 – 1.44 0.95 0.342
man woman × SRQ c 1.03 0.01 1.00 – 1.05 2.25 0.024
US oth × HS c 1.01 0.05 0.92 – 1.11 0.24 0.813
US oth × BS c 0.98 0.05 0.90 – 1.08 -0.37 0.715
US oth × SDO c 0.99 0.03 0.93 – 1.06 -0.21 0.837
US oth × SRQ c 1.00 0.00 1.00 – 1.01 0.88 0.381
Can AusUK × HS c 1.03 0.07 0.90 – 1.16 0.40 0.688
Can AusUK × BS c 0.98 0.07 0.85 – 1.13 -0.25 0.804
Can AusUK × SDO c 0.96 0.05 0.87 – 1.06 -0.87 0.384
Can AusUK × SRQ c 1.01 0.00 1.00 – 1.02 1.59 0.111
Aus UK × HS c 1.03 0.12 0.82 – 1.29 0.25 0.804
Aus UK × BS c 0.82 0.11 0.64 – 1.06 -1.48 0.138
Aus UK × SDO c 0.98 0.09 0.81 – 1.18 -0.21 0.831
Aus UK × SRQ c 1.01 0.01 0.99 – 1.03 0.90 0.370
(man woman × US oth) × HS
c
0.89 0.08 0.74 – 1.06 -1.30 0.194
(man woman × US oth) × BS
c
0.98 0.09 0.81 – 1.18 -0.22 0.827
(man woman × US oth) ×
SDO c
1.00 0.07 0.87 – 1.14 -0.03 0.974
(man woman × US oth) ×
SRQ c
1.01 0.01 1.00 – 1.02 1.25 0.211
(man woman × Can AusUK) ×
HS c
1.03 0.13 0.80 – 1.33 0.26 0.798
(man woman × Can AusUK) ×
BS c
0.82 0.12 0.62 – 1.09 -1.37 0.170
(man woman × Can AusUK) ×
SDO c
1.15 0.12 0.94 – 1.40 1.35 0.176
(man woman × Can AusUK) ×
SRQ c
1.01 0.01 0.99 – 1.03 0.75 0.455
(man woman × Aus UK) × HS
c
1.00 0.23 0.63 – 1.57 -0.01 0.996
(man woman × Aus UK) × BS
c
0.93 0.24 0.56 – 1.55 -0.28 0.779
(man woman × Aus UK) ×
SDO c
1.09 0.20 0.75 – 1.57 0.44 0.662
(man woman × Aus UK) ×
SRQ c
1.00 0.02 0.96 – 1.04 -0.10 0.917
Observations 980
R2 Nagelkerke 0.305
s3_olr3 <- clm(LikelihoodofVoting ~ (US_oth + Can_AusUK + Aus_UK) * man_woman + (warmth.c + competence.c)
               + age.c + edu.c,
               data = d3,
               link = "logit"
)

# summary(s3_olr3)
tab_model(s3_olr3,
          show.stat = T,
          show.se = T)
  Likelihoodof Voting
Predictors Odds Ratios std. Error CI Statistic p
1|2 0.21 0.02 0.18 – 0.25 -18.26 <0.001
2|3 0.54 0.04 0.47 – 0.62 -8.68 <0.001
3|4 1.77 0.12 1.54 – 2.03 8.14 <0.001
US oth 1.00 0.03 0.93 – 1.07 -0.03 0.977
Can AusUK 0.87 0.04 0.79 – 0.96 -2.87 0.004
Aus UK 0.86 0.07 0.72 – 1.02 -1.75 0.081
man woman 1.86 0.23 1.46 – 2.37 5.01 <0.001
warmth c 1.68 0.10 1.50 – 1.88 8.98 <0.001
competence c 0.60 0.06 0.49 – 0.72 -5.28 <0.001
age c 1.00 0.00 1.00 – 1.00 -1.57 0.117
edu c 1.10 0.02 1.06 – 1.15 4.40 <0.001
US oth × man woman 0.93 0.06 0.81 – 1.06 -1.10 0.273
Can AusUK × man woman 1.01 0.10 0.83 – 1.22 0.09 0.931
Aus UK × man woman 1.10 0.19 0.78 – 1.55 0.56 0.577
Observations 979
R2 Nagelkerke 0.165
s3_olr4 <- clm(LikelihoodofVoting ~ man_woman * (US_oth + Can_AusUK + Aus_UK) * Ideology.c
               + age.c + edu.c,
               data = d3,
               link = "logit"
)

# summary(s3_olr4)
tab_model(s3_olr4,
          show.stat = T,
          show.se = T)
  Likelihoodof Voting
Predictors Odds Ratios std. Error CI Statistic p
1|2 0.25 0.02 0.21 – 0.30 -15.84 <0.001
2|3 0.70 0.05 0.61 – 0.81 -4.82 <0.001
3|4 2.39 0.18 2.05 – 2.78 11.28 <0.001
man woman 1.47 0.19 1.14 – 1.89 3.00 0.003
US oth 0.99 0.04 0.92 – 1.06 -0.28 0.777
Can AusUK 0.89 0.05 0.81 – 0.99 -2.22 0.026
Aus UK 0.88 0.08 0.73 – 1.06 -1.32 0.187
Ideology c 0.62 0.03 0.57 – 0.68 -10.91 <0.001
age c 1.00 0.00 1.00 – 1.00 -0.33 0.742
edu c 1.09 0.02 1.04 – 1.14 3.67 <0.001
man woman × US oth 0.92 0.07 0.80 – 1.06 -1.17 0.240
man woman × Can AusUK 1.08 0.11 0.88 – 1.32 0.76 0.448
man woman × Aus UK 1.23 0.23 0.85 – 1.78 1.11 0.268
man woman × Ideology c 0.99 0.08 0.84 – 1.17 -0.09 0.925
US oth × Ideology c 1.02 0.02 0.98 – 1.06 0.99 0.320
Can AusUK × Ideology c 1.06 0.04 0.99 – 1.14 1.75 0.081
Aus UK × Ideology c 0.94 0.06 0.82 – 1.07 -1.00 0.320
(man woman × US oth) ×
Ideology c
0.98 0.04 0.90 – 1.07 -0.46 0.644
(man woman × Can AusUK) ×
Ideology c
1.03 0.07 0.90 – 1.18 0.39 0.698
(man woman × Aus UK) ×
Ideology c
1.24 0.17 0.96 – 1.62 1.63 0.103
Observations 980
R2 Nagelkerke 0.202
s3_olr5 <- clm(LikelihoodofVoting ~ man_woman * (US_oth + Can_AusUK + Aus_UK) * (Ideology.c + warmth.c + competence.c)
               + age.c + edu.c,
               data = d3,
               link = "logit"
)

# summary(s3_olr5)
tab_model(s3_olr5,
          show.stat = T,
          show.se = T)
  Likelihoodof Voting
Predictors Odds Ratios std. Error CI Statistic p
1|2 0.22 0.02 0.19 – 0.27 -16.22 <0.001
2|3 0.65 0.05 0.56 – 0.76 -5.58 <0.001
3|4 2.37 0.19 2.02 – 2.77 10.65 <0.001
man woman 1.54 0.21 1.18 – 2.01 3.21 0.001
US oth 1.00 0.04 0.93 – 1.08 0.13 0.899
Can AusUK 0.92 0.05 0.83 – 1.02 -1.57 0.117
Aus UK 0.90 0.09 0.74 – 1.09 -1.09 0.277
Ideology c 0.64 0.03 0.58 – 0.69 -9.93 <0.001
warmth c 1.61 0.10 1.43 – 1.82 7.74 <0.001
competence c 0.62 0.06 0.51 – 0.76 -4.57 <0.001
age c 1.00 0.00 1.00 – 1.00 -0.85 0.396
edu c 1.09 0.03 1.04 – 1.15 3.75 <0.001
man woman × US oth 0.94 0.07 0.81 – 1.09 -0.85 0.397
man woman × Can AusUK 1.11 0.12 0.90 – 1.37 0.96 0.338
man woman × Aus UK 1.32 0.26 0.90 – 1.96 1.41 0.157
man woman × Ideology c 0.99 0.09 0.83 – 1.18 -0.07 0.944
man woman × warmth c 0.87 0.11 0.68 – 1.10 -1.15 0.250
man woman × competence c 0.76 0.16 0.50 – 1.13 -1.36 0.175
US oth × Ideology c 1.02 0.02 0.98 – 1.07 0.94 0.349
US oth × warmth c 0.97 0.03 0.91 – 1.03 -0.97 0.332
US oth × competence c 0.96 0.06 0.85 – 1.07 -0.77 0.441
Can AusUK × Ideology c 1.07 0.04 1.00 – 1.15 1.98 0.048
Can AusUK × warmth c 0.95 0.05 0.87 – 1.05 -0.95 0.343
Can AusUK × competence c 1.03 0.08 0.87 – 1.21 0.33 0.739
Aus UK × Ideology c 0.96 0.07 0.84 – 1.11 -0.51 0.608
Aus UK × warmth c 0.89 0.08 0.74 – 1.06 -1.31 0.192
Aus UK × competence c 0.92 0.14 0.68 – 1.23 -0.57 0.568
(man woman × US oth) ×
Ideology c
0.97 0.04 0.89 – 1.06 -0.71 0.475
(man woman × US oth) ×
warmth c
0.89 0.06 0.78 – 1.01 -1.80 0.072
(man woman × US oth) ×
competence c
0.95 0.11 0.76 – 1.20 -0.42 0.675
(man woman × Can AusUK) ×
Ideology c
1.03 0.07 0.89 – 1.18 0.35 0.728
(man woman × Can AusUK) ×
warmth c
0.99 0.10 0.82 – 1.21 -0.06 0.954
(man woman × Can AusUK) ×
competence c
0.82 0.14 0.59 – 1.13 -1.20 0.231
(man woman × Aus UK) ×
Ideology c
1.16 0.16 0.88 – 1.53 1.05 0.295
(man woman × Aus UK) ×
warmth c
1.33 0.25 0.92 – 1.91 1.52 0.127
(man woman × Aus UK) ×
competence c
1.65 0.50 0.92 – 2.97 1.67 0.094
Observations 979
R2 Nagelkerke 0.282

Full model

gender main effects are pretty unstable across these models

s3_olr6 <- clm(LikelihoodofVoting ~ man_woman 
               * (warmth.c + competence.c + Ideology.c + HS.c + BS.c + SDO.c + SRQ.c)
               * (US_oth + Can_AusUK + Aus_UK)
              + age.c + edu.c,
               data = d3,
               link = "logit"
)

# summary(s3_olr6)
tab_model(s3_olr6,
          show.stat = T,
          show.se = T)
  Likelihoodof Voting
Predictors Odds Ratios std. Error CI Statistic p
1|2 0.16 0.02 0.13 – 0.20 -17.18 <0.001
2|3 0.54 0.05 0.45 – 0.64 -7.12 <0.001
3|4 2.27 0.20 1.91 – 2.70 9.28 <0.001
man woman 1.21 0.18 0.90 – 1.63 1.28 0.200
warmth c 1.52 0.10 1.34 – 1.73 6.30 <0.001
competence c 0.52 0.06 0.42 – 0.64 -6.01 <0.001
Ideology c 0.85 0.05 0.76 – 0.95 -2.80 0.005
HS c 0.64 0.06 0.54 – 0.77 -4.96 <0.001
BS c 1.12 0.10 0.93 – 1.34 1.20 0.231
SDO c 0.85 0.06 0.74 – 0.97 -2.35 0.019
SRQ c 0.98 0.01 0.97 – 0.99 -3.30 0.001
US oth 0.98 0.04 0.91 – 1.07 -0.41 0.681
Can AusUK 0.92 0.06 0.81 – 1.04 -1.38 0.167
Aus UK 0.86 0.10 0.70 – 1.07 -1.31 0.189
age c 1.00 0.00 1.00 – 1.00 -0.57 0.569
edu c 1.07 0.03 1.02 – 1.13 2.74 0.006
man woman × warmth c 0.95 0.13 0.73 – 1.23 -0.41 0.684
man woman × competence c 0.76 0.17 0.50 – 1.17 -1.24 0.215
man woman × Ideology c 0.87 0.10 0.70 – 1.09 -1.21 0.225
man woman × HS c 0.77 0.14 0.54 – 1.08 -1.50 0.134
man woman × BS c 0.85 0.16 0.59 – 1.23 -0.87 0.385
man woman × SDO c 1.09 0.16 0.83 – 1.44 0.63 0.526
man woman × SRQ c 1.03 0.01 1.00 – 1.05 2.23 0.026
man woman × US oth 0.89 0.07 0.76 – 1.05 -1.34 0.179
man woman × Can AusUK 1.04 0.13 0.82 – 1.32 0.30 0.765
man woman × Aus UK 1.48 0.33 0.96 – 2.28 1.76 0.078
warmth c × US oth 0.95 0.03 0.89 – 1.02 -1.38 0.167
warmth c × Can AusUK 0.94 0.05 0.85 – 1.05 -1.04 0.298
warmth c × Aus UK 0.96 0.10 0.79 – 1.17 -0.40 0.691
competence c × US oth 0.95 0.06 0.84 – 1.08 -0.75 0.453
competence c × Can AusUK 1.09 0.09 0.92 – 1.29 0.95 0.342
competence c × Aus UK 0.91 0.14 0.67 – 1.24 -0.60 0.546
Ideology c × US oth 1.03 0.03 0.97 – 1.09 1.05 0.293
Ideology c × Can AusUK 1.05 0.05 0.96 – 1.14 0.97 0.333
Ideology c × Aus UK 0.90 0.08 0.76 – 1.07 -1.18 0.238
HS c × US oth 0.99 0.05 0.90 – 1.09 -0.24 0.807
HS c × Can AusUK 1.04 0.07 0.90 – 1.19 0.55 0.584
HS c × Aus UK 1.01 0.13 0.79 – 1.30 0.10 0.918
BS c × US oth 1.01 0.05 0.91 – 1.11 0.18 0.857
BS c × Can AusUK 0.97 0.07 0.83 – 1.12 -0.46 0.648
BS c × Aus UK 0.88 0.12 0.67 – 1.16 -0.88 0.380
SDO c × US oth 0.97 0.04 0.90 – 1.04 -0.93 0.352
SDO c × Can AusUK 0.95 0.06 0.85 – 1.07 -0.88 0.379
SDO c × Aus UK 1.03 0.11 0.84 – 1.26 0.26 0.793
SRQ c × US oth 1.00 0.00 1.00 – 1.01 0.82 0.410
SRQ c × Can AusUK 1.01 0.01 1.00 – 1.02 1.60 0.109
SRQ c × Aus UK 1.01 0.01 0.99 – 1.03 0.77 0.442
(man woman × warmth c) ×
US oth
0.85 0.06 0.74 – 0.97 -2.36 0.018
(man woman × warmth c) ×
Can AusUK
1.04 0.11 0.84 – 1.29 0.36 0.721
(man woman × warmth c) ×
Aus UK
1.35 0.27 0.91 – 1.99 1.50 0.135
(man woman × competence
c) × US oth
0.93 0.12 0.73 – 1.18 -0.61 0.543
(man woman × competence
c) × Can AusUK
0.91 0.16 0.65 – 1.28 -0.54 0.592
(man woman × competence
c) × Aus UK
1.39 0.44 0.75 – 2.59 1.05 0.293
(man woman × Ideology c)
× US oth
0.96 0.06 0.85 – 1.07 -0.78 0.433
(man woman × Ideology c)
× Can AusUK
0.91 0.08 0.76 – 1.09 -0.98 0.327
(man woman × Ideology c)
× Aus UK
1.33 0.24 0.94 – 1.89 1.60 0.110
(man woman × HS c) × US
oth
0.92 0.09 0.76 – 1.12 -0.82 0.414
(man woman × HS c) × Can
AusUK
1.18 0.17 0.89 – 1.56 1.16 0.246
(man woman × HS c) × Aus
UK
0.79 0.20 0.48 – 1.29 -0.93 0.350
(man woman × BS c) × US
oth
1.01 0.10 0.83 – 1.23 0.13 0.896
(man woman × BS c) × Can
AusUK
0.80 0.12 0.60 – 1.08 -1.44 0.149
(man woman × BS c) × Aus
UK
0.94 0.26 0.54 – 1.62 -0.23 0.818
(man woman × SDO c) × US
oth
1.00 0.08 0.86 – 1.16 0.04 0.971
(man woman × SDO c) × Can
AusUK
1.26 0.15 1.00 – 1.59 1.96 0.049
(man woman × SDO c) × Aus
UK
1.06 0.22 0.70 – 1.60 0.28 0.780
(man woman × SRQ c) × US
oth
1.01 0.01 0.99 – 1.02 0.77 0.440
(man woman × SRQ c) × Can
AusUK
1.01 0.01 0.99 – 1.03 0.53 0.599
(man woman × SRQ c) × Aus
UK
1.00 0.02 0.96 – 1.03 -0.23 0.817
Observations 979
R2 Nagelkerke 0.391

Assumption

H0: Parallel Regression Assumption holds

Data Visualization

Experiment Vote Choice

d3 %>% ggplot(aes(x = as.numeric(LikelihoodofVoting), fill = Demo_Gender)) +
  geom_histogram( color="#e9ecef", alpha=0.8, position = 'dodge',
                  binwidth = .5) +
  scale_fill_manual(values=c("#69b3a2", "#404080")) +
  theme_bw() +
  scale_x_continuous( labels = c("Masc Man",
                                 "Femme Man",
                                 "Masc Woman",
                                 "Femme Woman"),
                       breaks = seq(1,4,1)) +
  scale_y_continuous(breaks = seq(0,80,10)) +
  labs(x = "Candidate Selected",
       y = "Frequency",
       fill = "Participant
Gender",
title = "Number Voting for Each Candidate by Gender and Country") +
  facet_wrap(~Country)

Predicting experimental vote choice

s3_olr.plot <- clm(LikelihoodofVoting ~ Demo_Gender * (Ideology + WomanTraits_warmth + WomanTraits_competence + HS + BS + SDO + SRQ) * Country,
               data = d3,
               link = "logit"
)


## Plot 1: Female candidate traits: Warmth

preds_olr.plot <- predictions(
  s3_olr.plot,
  newdata = datagrid(
    
# Plotting variables
    WomanTraits_warmth = seq(
      min(d3$WomanTraits_warmth, na.rm = TRUE), 
      max(d3$WomanTraits_warmth, na.rm = TRUE), 
      length.out = 100
    ),
    Demo_Gender = unique(na.omit(d3$Demo_Gender)),
    Country = unique(na.omit(d3$Country)),
    
# Variables to hold constant
    Ideology = mean(d3$Ideology, na.rm = TRUE),
    WomanTraits_competence = mean(d3$WomanTraits_competence, na.rm = TRUE),
    HS = mean(d3$HS, na.rm = TRUE), 
    BS = mean(d3$BS, na.rm = TRUE),
    SDO = mean(d3$SDO, na.rm = TRUE),
    SRQ = mean(d3$SRQ, na.rm = TRUE)
  )
)

# Ensure prediction data is sorted cleanly for smooth lines
preds <- preds_olr.plot[order(preds_olr.plot$Demo_Gender, preds_olr.plot$group, preds_olr.plot$WomanTraits_warmth), ]

preds$group_labeled <- factor(
  preds$group, 
  labels = c("Masc Man",
             "Femme Man",
             "Masc Woman",
             "Femme Woman")
)

# Plot
preds %>% ggplot(aes(x = WomanTraits_warmth, y = estimate, color = factor(group_labeled), group = factor(group_labeled))) +
  geom_line(linewidth = 1) +
  geom_ribbon(aes(ymin = conf.low, ymax = conf.high, fill = factor(group_labeled)), alpha = 0.15, color = NA) +
  facet_grid(Demo_Gender ~ Country) +
  labs(
    x = "Woman Traits Warmth",
    y = "Predicted Probability of Voting",
    color = "Candidate",
    fill = "Candidate",
    caption = "Note: Model-predicted values. Other continuous predictors held at mean.",
    title = "Likelihood of voting for each candidate by warmth importance, shown by gender and country"
  ) +
  scale_x_continuous(breaks = seq(1,7,1)) +
  theme_minimal()

## Plot 2: SRQ

preds_olr.plot2 <- predictions(
  s3_olr.plot,
  newdata = datagrid(
    
# Plotting variables
    SRQ = seq(
      min(d3$SRQ, na.rm = TRUE), 
      max(d3$SRQ, na.rm = TRUE), 
      length.out = 100
    ),
    Demo_Gender = unique(na.omit(d3$Demo_Gender)),
    Country = unique(na.omit(d3$Country)),
    
# Variables to hold constant
    Ideology = mean(d3$Ideology, na.rm = TRUE), 
    WomanTraits_competence = mean(d3$WomanTraits_competence, na.rm = TRUE),
    WomanTraits_warmth = mean(d3$WomanTraits_warmth, na.rm = TRUE),
    HS = mean(d3$HS, na.rm = TRUE),
    BS = mean(d3$BS, na.rm = TRUE),
    SDO = mean(d3$SDO, na.rm = TRUE)
  )
)

# Ensure prediction data is sorted cleanly for smooth lines
preds2 <- preds_olr.plot2[order(preds_olr.plot2$Demo_Gender, preds_olr.plot2$group, preds_olr.plot2$WomanTraits_warmth), ]

preds2$group_labeled <- factor(
  preds2$group, 
  labels = c("Masc Man",
             "Femme Man",
             "Masc Woman",
             "Femme Woman")
)

# Plot
preds2 %>% ggplot(aes(x = SRQ, y = estimate, color = factor(group_labeled), group = factor(group_labeled))) +
  geom_line(linewidth = 1) +
  geom_ribbon(aes(ymin = conf.low, ymax = conf.high, fill = factor(group_labeled)), alpha = 0.15, color = NA) +
  facet_grid(Demo_Gender ~ Country) +
  labs(
    x = "Social Roles Questionnaire (0 to 100)",
    y = "Predicted Probability of Voting",
    color = "Candidate",
    fill = "Candidate",
    caption = "Note: Model-predicted values. Other continuous predictors held at mean.",
    title = "Likelihood of voting for each candidate by SRQ, shown by gender and country"
  ) +
  scale_x_continuous(breaks = seq(0,100,25)) +
  theme_minimal()

Summary of Findings

DV: Voting for Harris

Conditional Effects

Across all three studies, there is no effect of gender in likelihood of voting for Harris (compared with any other choice) once controlling for other variables.

Hostile sexism consistently negatively predicts likelihood of voting Harris across all three studies, and social roles questionnaire negatively predicts it in studies 2 and 3. Impact index positively predicts and ideology (increasing conservatism) negatively predicts likelihood of voting Harris when measured (studies 1 and 2 for impact index, studies 2 and 3 for ideology).

Interaction Effects

Across all three studies, there were only two cross-country higher-order interactions (SRQ in Study 2, Ideology in Study 3). Higher-order gender interactions were only observed in study 1 (with hostile sexism and benevolent sexism).

DV: Experimental Candidate Selection

Conditional Effects

Participants’ rated importance of a female candidate’s warmth positively predicted voting for the feminine woman candidate, whereas rated importance of a female candidate’s competence negatively predicted voting for the feminine woman. Increasing conservatism, hostile sexism, social dominance orientation, and social roles questionnaire also negatively predicted voting for the feminine woman candidate. Education positively predicted voting for the feminine woman.

Interaction Effects

There was an interaction between participant gender (man vs. woman) and social roles questionnaire on voting for the feminine woman candidate. There was also a three-way interaction between gender, country, and rated importance of warmth in likelihood of voting for the feminine woman candidate.