Media Censorship

Media Censorship Data

media_censorship_data <- final_data |> 
  filter(condition %in% c(2, 5, 8, 11))


#support for candidate
media_censorship_data |> 
  ggplot(aes(x = support_election)) +
  geom_bar() +
  facet_wrap(~politician_gender + politician_party)

media_censorship_data |> 
  ggplot(aes(x = support_election_lumped)) +
  geom_bar() +
  facet_wrap(~politician_gender + politician_party)

#recognition of antidemocratic behavior

media_censorship_data |> 
  ggplot(aes(x = threaten_country)) +
  geom_bar() +
  facet_wrap(~politician_gender + politician_party)

media_censorship_data |> 
  ggplot(aes(x = threaten_country_lumped)) +
  geom_bar() +
  facet_wrap(~politician_gender + politician_party)

#Opinion on bill

media_censorship_data |> 
  ggplot(aes(x = opinion_bill)) +
  geom_bar() +
  facet_wrap(~politician_gender + politician_party)

media_censorship_data |> 
  ggplot(aes(x = opinion_bill_lumped)) +
  geom_bar() +
  facet_wrap(~politician_gender + politician_party)

media_difference_in_means <- media_censorship_data |> 
  group_by(politician_party, politician_gender) |> 
  summarize(mean_support_for_politician = mean(support_election, na.rm = TRUE),
            sd_support_for_politician = sd(support_election, na.rm = TRUE),
            mean_antidemocratic_recognition = mean(threaten_country, na.rm = TRUE),
            sd_antidemocratic_recognition = sd(threaten_country, na.rm = TRUE),
            mean_support_for_bill = mean(opinion_bill, na.rm = TRUE),
            sd_support_for_bill = sd(opinion_bill, na.rm = TRUE),
            n = n(),
            .groups = "drop") |> 
  print()
# A tibble: 4 × 9
  politician_party politician_gender mean_support_for_politician
  <fct>            <fct>                                   <dbl>
1 Democrat         Male Politician                          3.79
2 Democrat         Female Politician                        4.05
3 Republican       Male Politician                          3.54
4 Republican       Female Politician                        3.37
# ℹ 6 more variables: sd_support_for_politician <dbl>,
#   mean_antidemocratic_recognition <dbl>, sd_antidemocratic_recognition <dbl>,
#   mean_support_for_bill <dbl>, sd_support_for_bill <dbl>, n <int>

Hypothesis 1: support for a politician increases under female candidates, regardless of party

#simple model

media_h1 <- lm(support_election ~ politician_gender * politician_party, data = media_censorship_data)

summary(media_h1)

Call:
lm(formula = support_election ~ politician_gender * politician_party, 
    data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.0506 -1.7914  0.2086  1.4581  3.6289 

Coefficients:
                                                              Estimate
(Intercept)                                                     3.7914
politician_genderFemale Politician                              0.2592
politician_partyRepublican                                     -0.2495
politician_genderFemale Politician:politician_partyRepublican  -0.4301
                                                              Std. Error
(Intercept)                                                       0.1454
politician_genderFemale Politician                                0.2072
politician_partyRepublican                                        0.2082
politician_genderFemale Politician:politician_partyRepublican     0.2947
                                                              t value Pr(>|t|)
(Intercept)                                                    26.078   <2e-16
politician_genderFemale Politician                              1.251    0.211
politician_partyRepublican                                     -1.198    0.231
politician_genderFemale Politician:politician_partyRepublican  -1.459    0.145
                                                                 
(Intercept)                                                   ***
politician_genderFemale Politician                               
politician_partyRepublican                                       
politician_genderFemale Politician:politician_partyRepublican    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.856 on 631 degrees of freedom
Multiple R-squared:  0.01886,   Adjusted R-squared:  0.0142 
F-statistic: 4.044 on 3 and 631 DF,  p-value: 0.007297
plot_model(media_h1, terms = c("politician_party", "politician_gender"),
           type= "pred")

#include covariates

media_h1_covariates <- lm(support_election ~ politician_gender * politician_party + survey_partyid + respondent_gender, data = media_censorship_data)

summary(media_h1_covariates)

Call:
lm(formula = support_election ~ politician_gender * politician_party + 
    survey_partyid + respondent_gender, data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.4845 -1.7168  0.2729  1.3357  3.7377 

Coefficients:
                                                              Estimate
(Intercept)                                                    4.20861
politician_genderFemale Politician                             0.24901
politician_partyRepublican                                    -0.24395
survey_partyid                                                -0.13222
respondent_genderFemale Respondent                             0.29129
politician_genderFemale Politician:politician_partyRepublican -0.42253
                                                              Std. Error
(Intercept)                                                      0.23723
politician_genderFemale Politician                               0.20523
politician_partyRepublican                                       0.20629
survey_partyid                                                   0.03995
respondent_genderFemale Respondent                               0.14601
politician_genderFemale Politician:politician_partyRepublican    0.29180
                                                              t value Pr(>|t|)
(Intercept)                                                    17.741  < 2e-16
politician_genderFemale Politician                              1.213 0.225472
politician_partyRepublican                                     -1.183 0.237424
survey_partyid                                                 -3.310 0.000988
respondent_genderFemale Respondent                              1.995 0.046470
politician_genderFemale Politician:politician_partyRepublican  -1.448 0.148107
                                                                 
(Intercept)                                                   ***
politician_genderFemale Politician                               
politician_partyRepublican                                       
survey_partyid                                                ***
respondent_genderFemale Respondent                            *  
politician_genderFemale Politician:politician_partyRepublican    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.838 on 629 degrees of freedom
Multiple R-squared:  0.04129,   Adjusted R-squared:  0.03367 
F-statistic: 5.418 on 5 and 629 DF,  p-value: 6.869e-05
plot_model(media_h1_covariates, terms = c("politician_party", "politician_gender"),
           type= "pred")

statistically significant difference between the effects of female voters under dems versus republicans

but cannot say that there is a significant difference between male and female under each condition.

still, can comment on the direction?

Hypothesis 2: this effect will be more pronounced among outparty voters than inparty voters

H2a: Overall effect

#simple model

media_h2 <- lm(support_election ~ politician_gender * inparty_outparty, data = media_censorship_data)

summary(media_h2)

Call:
lm(formula = support_election ~ politician_gender * inparty_outparty, 
    data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.8983 -1.0050  0.1017  1.1017  4.0000 

Coefficients:
                                                           Estimate Std. Error
(Intercept)                                                3.000000   0.119181
politician_genderFemale Politician                         0.005025   0.167272
inparty_outpartyInparty                                    1.704000   0.190093
politician_genderFemale Politician:inparty_outpartyInparty 0.189280   0.270450
                                                           t value Pr(>|t|)    
(Intercept)                                                 25.172   <2e-16 ***
politician_genderFemale Politician                           0.030    0.976    
inparty_outpartyInparty                                      8.964   <2e-16 ***
politician_genderFemale Politician:inparty_outpartyInparty   0.700    0.484    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.656 on 631 degrees of freedom
Multiple R-squared:  0.2194,    Adjusted R-squared:  0.2156 
F-statistic:  59.1 on 3 and 631 DF,  p-value: < 2.2e-16
plot_model(media_h2, terms = c("inparty_outparty", "politician_gender"),
           type= "pred")

#include covariates

media_h2_covariates <- lm(support_election ~ politician_gender * inparty_outparty + survey_partyid + respondent_gender, data = media_censorship_data)

summary(media_h2_covariates)

Call:
lm(formula = support_election ~ politician_gender * inparty_outparty + 
    survey_partyid + respondent_gender, data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.9188 -1.3777  0.1189  1.1290  4.1290 

Coefficients:
                                                           Estimate Std. Error
(Intercept)                                                 3.44936    0.20061
politician_genderFemale Politician                          0.01012    0.16466
inparty_outpartyInparty                                     1.73455    0.18727
survey_partyid                                             -0.14459    0.03545
respondent_genderFemale Respondent                          0.31427    0.12940
politician_genderFemale Politician:inparty_outpartyInparty  0.15858    0.26635
                                                           t value Pr(>|t|)    
(Intercept)                                                 17.194  < 2e-16 ***
politician_genderFemale Politician                           0.061   0.9510    
inparty_outpartyInparty                                      9.262  < 2e-16 ***
survey_partyid                                              -4.078 5.12e-05 ***
respondent_genderFemale Respondent                           2.429   0.0154 *  
politician_genderFemale Politician:inparty_outpartyInparty   0.595   0.5518    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.63 on 629 degrees of freedom
Multiple R-squared:  0.246, Adjusted R-squared:   0.24 
F-statistic: 41.04 on 5 and 629 DF,  p-value: < 2.2e-16
plot_model(media_h2_covariates, terms = c("inparty_outparty", "politician_gender"),
           type= "pred")

#compare this to next chunk - is it needed?

H2b: Split to look at just Democratic candidates, then just Republican

#################repeated for just democratic candidate ###############

dem_politician_media_data <- media_censorship_data |> 
  filter(politician_party == "Democrat")

media_h2a <- lm(support_election ~ politician_gender * inparty_outparty, data = dem_politician_media_data)

summary(media_h2a)

Call:
lm(formula = support_election ~ politician_gender * inparty_outparty, 
    data = dem_politician_media_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.8529 -1.4444  0.3699  1.1471  3.8889 

Coefficients:
                                                           Estimate Std. Error
(Intercept)                                                  3.1111     0.1767
politician_genderFemale Politician                           0.3333     0.2499
inparty_outpartyInparty                                      1.5190     0.2640
politician_genderFemale Politician:inparty_outpartyInparty  -0.1105     0.3772
                                                           t value Pr(>|t|)    
(Intercept)                                                 17.608  < 2e-16 ***
politician_genderFemale Politician                           1.334    0.183    
inparty_outpartyInparty                                      5.753 2.06e-08 ***
politician_genderFemale Politician:inparty_outpartyInparty  -0.293    0.770    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.676 on 317 degrees of freedom
Multiple R-squared:  0.1644,    Adjusted R-squared:  0.1565 
F-statistic: 20.79 on 3 and 317 DF,  p-value: 2.538e-12
plot_model(media_h2a, terms = c("inparty_outparty", "politician_gender"),
           type= "pred")

#####################just republican candidate ##################
rep_politician_media_data <- media_censorship_data |> 
  filter(politician_party == "Republican")

media_h2b <- lm(support_election ~ politician_gender * inparty_outparty, data = rep_politician_media_data)

summary(media_h2b)

Call:
lm(formula = support_election ~ politician_gender * inparty_outparty, 
    data = rep_politician_media_data)

Residuals:
   Min     1Q Median     3Q    Max 
-3.960 -1.642  0.040  1.192  4.358 

Coefficients:
                                                           Estimate Std. Error
(Intercept)                                                  2.9029     0.1587
politician_genderFemale Politician                          -0.2607     0.2213
inparty_outpartyInparty                                      1.9048     0.2740
politician_genderFemale Politician:inparty_outpartyInparty   0.4130     0.3882
                                                           t value Pr(>|t|)    
(Intercept)                                                 18.293  < 2e-16 ***
politician_genderFemale Politician                          -1.178    0.240    
inparty_outpartyInparty                                      6.952 2.13e-11 ***
politician_genderFemale Politician:inparty_outpartyInparty   1.064    0.288    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.611 on 310 degrees of freedom
Multiple R-squared:  0.2794,    Adjusted R-squared:  0.2725 
F-statistic: 40.07 on 3 and 310 DF,  p-value: < 2.2e-16
plot_model(media_h2b, terms = c("inparty_outparty", "politician_gender"),
           type= "pred")

For dems, both out and inparty voters see higher effects for female candidates

for Reps, outparty decreases, but inparty increases (albeit not significantly)

Comparing the effects of party

#Just inparty/outparty

media_party1 <- lm(support_election ~ politician_gender + inparty_outparty, data = media_censorship_data)

summary(media_party1)

Call:
lm(formula = support_election ~ politician_gender + inparty_outparty, 
    data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.8382 -1.0407  0.1618  1.1618  4.0368 

Coefficients:
                                   Estimate Std. Error t value Pr(>|t|)    
(Intercept)                         2.96324    0.10694  27.709   <2e-16 ***
politician_genderFemale Politician  0.07743    0.13139   0.589    0.556    
inparty_outpartyInparty             1.79751    0.13516  13.299   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.655 on 632 degrees of freedom
Multiple R-squared:  0.2187,    Adjusted R-squared:  0.2163 
F-statistic: 88.48 on 2 and 632 DF,  p-value: < 2.2e-16
#Just survey_party


media_party2 <- lm(support_election ~ politician_gender + survey_partyid, data = media_censorship_data)

summary(media_party2)

Call:
lm(formula = support_election ~ politician_gender + survey_partyid, 
    data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.1350 -1.7412  0.2588  1.4214  3.6839 

Coefficients:
                                   Estimate Std. Error t value Pr(>|t|)    
(Intercept)                         4.23491    0.20255  20.908   <2e-16 ***
politician_genderFemale Politician  0.03140    0.14740   0.213   0.8314    
survey_partyid                     -0.13127    0.04036  -3.253   0.0012 ** 
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.857 on 632 degrees of freedom
Multiple R-squared:  0.01658,   Adjusted R-squared:  0.01347 
F-statistic: 5.327 on 2 and 632 DF,  p-value: 0.00508
#both

media_party3 <- lm(support_election ~ politician_gender + inparty_outparty * survey_partyid, data = media_censorship_data)

summary(media_party3)

Call:
lm(formula = support_election ~ politician_gender + inparty_outparty * 
    survey_partyid, data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.9789 -1.4391  0.0996  1.2039  3.9917 

Coefficients:
                                       Estimate Std. Error t value Pr(>|t|)    
(Intercept)                             3.72653    0.23714  15.714  < 2e-16 ***
politician_genderFemale Politician      0.07165    0.12988   0.552 0.581350    
inparty_outpartyInparty                 1.50394    0.33260   4.522 7.33e-06 ***
survey_partyid                         -0.17955    0.05022  -3.575 0.000376 ***
inparty_outpartyInparty:survey_partyid  0.07180    0.07112   1.010 0.313037    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.635 on 630 degrees of freedom
Multiple R-squared:  0.2397,    Adjusted R-squared:  0.2349 
F-statistic: 49.66 on 4 and 630 DF,  p-value: < 2.2e-16
vif_media <- vif(media_party3, type = "predictor")
GVIFs computed for predictors
print(vif_media)
                    GVIF Df GVIF^(1/(2*Df))   Interacts With
politician_gender 1.0014  1        1.000700             --  
inparty_outparty  1.0014  3        1.000233   survey_partyid
survey_partyid    1.0014  3        1.000233 inparty_outparty
                                  Other Predictors
politician_gender inparty_outparty, survey_partyid
inparty_outparty                 politician_gender
survey_partyid                   politician_gender

Hypothesis 3: These above effects will be more pronounced among female voters

#hypothesis 3: 

#simple model

media_h3 <- lm(support_election ~ politician_gender * respondent_gender, data = media_censorship_data)

summary(media_h3)

Call:
lm(formula = support_election ~ politician_gender * respondent_gender, 
    data = media_censorship_data)

Residuals:
   Min     1Q Median     3Q    Max 
-2.847 -1.825  0.175  1.408  3.516 

Coefficients:
                                                                      Estimate
(Intercept)                                                             3.4839
politician_genderFemale Politician                                      0.1085
respondent_genderFemale Respondent                                      0.3628
politician_genderFemale Politician:respondent_genderFemale Respondent  -0.1301
                                                                      Std. Error
(Intercept)                                                               0.1500
politician_genderFemale Politician                                        0.2115
respondent_genderFemale Respondent                                        0.2095
politician_genderFemale Politician:respondent_genderFemale Respondent     0.2965
                                                                      t value
(Intercept)                                                            23.224
politician_genderFemale Politician                                      0.513
respondent_genderFemale Respondent                                      1.731
politician_genderFemale Politician:respondent_genderFemale Respondent  -0.439
                                                                      Pr(>|t|)
(Intercept)                                                             <2e-16
politician_genderFemale Politician                                      0.6081
respondent_genderFemale Respondent                                      0.0839
politician_genderFemale Politician:respondent_genderFemale Respondent   0.6609
                                                                         
(Intercept)                                                           ***
politician_genderFemale Politician                                       
respondent_genderFemale Respondent                                    .  
politician_genderFemale Politician:respondent_genderFemale Respondent    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.868 on 631 degrees of freedom
Multiple R-squared:  0.006768,  Adjusted R-squared:  0.002046 
F-statistic: 1.433 on 3 and 631 DF,  p-value: 0.2319
plot_model(media_h3, terms = c("respondent_gender", "politician_gender"),
           type= "pred")

#include covariates


media_h3_covariates <- lm(support_election ~ politician_gender * respondent_gender + politician_party + survey_partyid, data = media_censorship_data)

summary(media_h3_covariates)

Call:
lm(formula = support_election ~ politician_gender * respondent_gender + 
    politician_party + survey_partyid, data = media_censorship_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.4955 -1.7377  0.2425  1.3570  3.6971 

Coefficients:
                                                                      Estimate
(Intercept)                                                            4.28793
politician_genderFemale Politician                                     0.09137
respondent_genderFemale Respondent                                     0.34014
politician_partyRepublican                                            -0.45463
survey_partyid                                                        -0.13261
politician_genderFemale Politician:respondent_genderFemale Respondent -0.10094
                                                                      Std. Error
(Intercept)                                                              0.23937
politician_genderFemale Politician                                       0.20851
respondent_genderFemale Respondent                                       0.20660
politician_partyRepublican                                               0.14621
survey_partyid                                                           0.04002
politician_genderFemale Politician:respondent_genderFemale Respondent    0.29232
                                                                      t value
(Intercept)                                                            17.913
politician_genderFemale Politician                                      0.438
respondent_genderFemale Respondent                                      1.646
politician_partyRepublican                                             -3.109
survey_partyid                                                         -3.314
politician_genderFemale Politician:respondent_genderFemale Respondent  -0.345
                                                                      Pr(>|t|)
(Intercept)                                                            < 2e-16
politician_genderFemale Politician                                    0.661401
respondent_genderFemale Respondent                                    0.100188
politician_partyRepublican                                            0.001959
survey_partyid                                                        0.000974
politician_genderFemale Politician:respondent_genderFemale Respondent 0.729967
                                                                         
(Intercept)                                                           ***
politician_genderFemale Politician                                       
respondent_genderFemale Respondent                                       
politician_partyRepublican                                            ** 
survey_partyid                                                        ***
politician_genderFemale Politician:respondent_genderFemale Respondent    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.841 on 629 degrees of freedom
Multiple R-squared:  0.03828,   Adjusted R-squared:  0.03063 
F-statistic: 5.007 on 5 and 629 DF,  p-value: 0.0001659
plot_model(media_h3_covariates, terms = c("respondent_gender", "politician_gender"),
           type= "pred")

women respondents see no real difference.

Filter by Candidate Gender

male_media_data <- media_censorship_data |> 
  filter(politician_gender == "Male Politician")

lm_male_media <- lm(support_election ~ politician_party * respondent_gender, data = male_media_data)

summary(lm_male_media)

Call:
lm(formula = support_election ~ politician_party * respondent_gender, 
    data = male_media_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.0000 -1.6711  0.3289  1.4474  3.5823 

Coefficients:
                                                              Estimate
(Intercept)                                                     3.5526
politician_partyRepublican                                     -0.1349
respondent_genderFemale Respondent                              0.4474
politician_partyRepublican:respondent_genderFemale Respondent  -0.1940
                                                              Std. Error
(Intercept)                                                       0.2144
politician_partyRepublican                                        0.3003
respondent_genderFemale Respondent                                0.2934
politician_partyRepublican:respondent_genderFemale Respondent     0.4199
                                                              t value Pr(>|t|)
(Intercept)                                                    16.572   <2e-16
politician_partyRepublican                                     -0.449    0.654
respondent_genderFemale Respondent                              1.525    0.128
politician_partyRepublican:respondent_genderFemale Respondent  -0.462    0.644
                                                                 
(Intercept)                                                   ***
politician_partyRepublican                                       
respondent_genderFemale Respondent                               
politician_partyRepublican:respondent_genderFemale Respondent    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.869 on 314 degrees of freedom
Multiple R-squared:  0.01398,   Adjusted R-squared:  0.004558 
F-statistic: 1.484 on 3 and 314 DF,  p-value: 0.2189
plot_model(lm_male_media, terms = c("respondent_gender", "politician_party"),
           type= "pred")

female_media_data <- media_censorship_data |> 
  filter(politician_gender == "Female Politician")

lm_female_media <- lm(support_election ~ politician_party * respondent_gender, data = female_media_data)

summary(lm_female_media)

Call:
lm(formula = support_election ~ politician_party * respondent_gender, 
    data = female_media_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.2073 -1.4231  0.1184  1.5769  3.6790 

Coefficients:
                                                              Estimate
(Intercept)                                                     3.8816
politician_partyRepublican                                     -0.5606
respondent_genderFemale Respondent                              0.3257
politician_partyRepublican:respondent_genderFemale Respondent  -0.2236
                                                              Std. Error
(Intercept)                                                       0.2113
politician_partyRepublican                                        0.2942
respondent_genderFemale Respondent                                0.2933
politician_partyRepublican:respondent_genderFemale Respondent     0.4140
                                                              t value Pr(>|t|)
(Intercept)                                                    18.370   <2e-16
politician_partyRepublican                                     -1.906   0.0576
respondent_genderFemale Respondent                              1.111   0.2676
politician_partyRepublican:respondent_genderFemale Respondent  -0.540   0.5895
                                                                 
(Intercept)                                                   ***
politician_partyRepublican                                    .  
respondent_genderFemale Respondent                               
politician_partyRepublican:respondent_genderFemale Respondent    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.842 on 313 degrees of freedom
Multiple R-squared:  0.03734,   Adjusted R-squared:  0.02811 
F-statistic: 4.047 on 3 and 313 DF,  p-value: 0.007626
plot_model(lm_female_media, terms = c("respondent_gender", "politician_party"),
           type= "pred")

#standardize so can compare side by side?

Filtering by respondent party: Democrats

dem_media_data <- media_censorship_data |> 
  filter(survey_partyid == c(5,6,7))
Warning: There was 1 warning in `filter()`.
ℹ In argument: `survey_partyid == c(5, 6, 7)`.
Caused by warning in `survey_partyid == c(5, 6, 7)`:
! longer object length is not a multiple of shorter object length
lm_dem_media <- lm(support_election ~ politician_gender * politician_party, data = dem_media_data)

summary(lm_dem_media)

Call:
lm(formula = support_election ~ politician_gender * politician_party, 
    data = dem_media_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.7308 -1.0313 -0.1304  1.1190  2.9687 

Coefficients:
                                                              Estimate
(Intercept)                                                    4.54545
politician_genderFemale Politician                             0.18531
politician_partyRepublican                                    -2.51420
politician_genderFemale Politician:politician_partyRepublican -0.08613
                                                              Std. Error
(Intercept)                                                      0.28919
politician_genderFemale Politician                               0.39293
politician_partyRepublican                                       0.37567
politician_genderFemale Politician:politician_partyRepublican    0.54026
                                                              t value Pr(>|t|)
(Intercept)                                                    15.718  < 2e-16
politician_genderFemale Politician                              0.472    0.638
politician_partyRepublican                                     -6.693 1.33e-09
politician_genderFemale Politician:politician_partyRepublican  -0.159    0.874
                                                                 
(Intercept)                                                   ***
politician_genderFemale Politician                               
politician_partyRepublican                                    ***
politician_genderFemale Politician:politician_partyRepublican    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.356 on 99 degrees of freedom
Multiple R-squared:  0.4831,    Adjusted R-squared:  0.4674 
F-statistic: 30.84 on 3 and 99 DF,  p-value: 3.655e-14
plot_model(lm_dem_media, terms = c("politician_party", "politician_gender"),
           type= "pred")

#with covariates


lm_dem_media2 <- lm(support_election ~ politician_gender * politician_party + respondent_gender, data = dem_media_data)

summary(lm_dem_media2)

Call:
lm(formula = support_election ~ politician_gender * politician_party + 
    respondent_gender, data = dem_media_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.2012 -0.7937 -0.1151  0.9382  3.4260 

Coefficients:
                                                              Estimate
(Intercept)                                                    4.11514
politician_genderFemale Politician                             0.08601
politician_partyRepublican                                    -2.54110
respondent_genderFemale Respondent                             0.86063
politician_genderFemale Politician:politician_partyRepublican  0.13361
                                                              Std. Error
(Intercept)                                                      0.30411
politician_genderFemale Politician                               0.37541
politician_partyRepublican                                       0.35788
respondent_genderFemale Respondent                               0.25785
politician_genderFemale Politician:politician_partyRepublican    0.51875
                                                              t value Pr(>|t|)
(Intercept)                                                    13.532  < 2e-16
politician_genderFemale Politician                              0.229   0.8193
politician_partyRepublican                                     -7.100 1.99e-10
respondent_genderFemale Respondent                              3.338   0.0012
politician_genderFemale Politician:politician_partyRepublican   0.258   0.7973
                                                                 
(Intercept)                                                   ***
politician_genderFemale Politician                               
politician_partyRepublican                                    ***
respondent_genderFemale Respondent                            ** 
politician_genderFemale Politician:politician_partyRepublican    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.292 on 98 degrees of freedom
Multiple R-squared:  0.5359,    Adjusted R-squared:  0.5169 
F-statistic: 28.29 on 4 and 98 DF,  p-value: 1.262e-15
plot_model(lm_dem_media2, terms = c("politician_party", "politician_gender"),
           type= "pred")

#need to reconcile partyid from survey, prolific

Filtering by respondent party: Republicans

rep_media_data <- media_censorship_data |> 
  filter(survey_partyid == c(1,2,3))
Warning: There was 1 warning in `filter()`.
ℹ In argument: `survey_partyid == c(1, 2, 3)`.
Caused by warning in `survey_partyid == c(1, 2, 3)`:
! longer object length is not a multiple of shorter object length
lm_rep_media <- lm(support_election ~ politician_gender * politician_party, data = rep_media_data)

summary(lm_rep_media)

Call:
lm(formula = support_election ~ politician_gender * politician_party, 
    data = rep_media_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.2500 -1.0769 -0.0769  1.3684  3.3684 

Coefficients:
                                                              Estimate
(Intercept)                                                     2.0769
politician_genderFemale Politician                              1.5547
politician_partyRepublican                                      2.1731
politician_genderFemale Politician:politician_partyRepublican  -0.9297
                                                              Std. Error
(Intercept)                                                       0.4688
politician_genderFemale Politician                                0.6085
politician_partyRepublican                                        0.7596
politician_genderFemale Politician:politician_partyRepublican     0.9518
                                                              t value Pr(>|t|)
(Intercept)                                                     4.430 4.88e-05
politician_genderFemale Politician                              2.555  0.01358
politician_partyRepublican                                      2.861  0.00607
politician_genderFemale Politician:politician_partyRepublican  -0.977  0.33325
                                                                 
(Intercept)                                                   ***
politician_genderFemale Politician                            *  
politician_partyRepublican                                    ** 
politician_genderFemale Politician:politician_partyRepublican    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.69 on 52 degrees of freedom
Multiple R-squared:  0.2836,    Adjusted R-squared:  0.2423 
F-statistic: 6.863 on 3 and 52 DF,  p-value: 0.000556
plot_model(lm_rep_media, terms = c("politician_party", "politician_gender"),
           type= "pred")

#include covariates


lm_rep_media2 <- lm(support_election ~ politician_gender * politician_party + respondent_gender, data = rep_media_data)

summary(lm_rep_media2)

Call:
lm(formula = support_election ~ politician_gender * politician_party + 
    respondent_gender, data = rep_media_data)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.0645 -1.1532  0.0088  1.5767  3.5767 

Coefficients:
                                                              Estimate
(Intercept)                                                     1.8867
politician_genderFemale Politician                              1.5366
politician_partyRepublican                                      2.1778
respondent_genderFemale Respondent                              0.4946
politician_genderFemale Politician:politician_partyRepublican  -0.9425
                                                              Std. Error
(Intercept)                                                       0.5003
politician_genderFemale Politician                                0.6077
politician_partyRepublican                                        0.7585
respondent_genderFemale Respondent                                0.4590
politician_genderFemale Politician:politician_partyRepublican     0.9504
                                                              t value Pr(>|t|)
(Intercept)                                                     3.771 0.000424
politician_genderFemale Politician                              2.528 0.014595
politician_partyRepublican                                      2.871 0.005938
respondent_genderFemale Respondent                              1.078 0.286293
politician_genderFemale Politician:politician_partyRepublican  -0.992 0.326032
                                                                 
(Intercept)                                                   ***
politician_genderFemale Politician                            *  
politician_partyRepublican                                    ** 
respondent_genderFemale Respondent                               
politician_genderFemale Politician:politician_partyRepublican    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.688 on 51 degrees of freedom
Multiple R-squared:  0.2996,    Adjusted R-squared:  0.2446 
F-statistic: 5.453 on 4 and 51 DF,  p-value: 0.0009841
plot_model(lm_rep_media2, terms = c("politician_party", "politician_gender"),
           type= "pred")

#need to reconcile partyid from survey, prolific

Interesting effect here with female politicians from outparty for reps (this mirrors outparty plot findings but is more significant because does not include moderates).