Mayor 2023: Back of the Envelope, Part 2

As the Mayor’s race heats up, I’m doing a series establishing some baseline numbers. What follows are simplistic calculations using reasonable assumptions. Welcome to the Back of the Envelope. See Part 1 here.

Will candidates split the vote?

In Part 1, I argued that Philadelphia’s votes will probably be constituted by 35% from the Black Voter wards, 30% from Wealthy Progressive wards, 25% from White Moderate wards, and 10% from Hispanic Voter wards.

This means no single bloc can win the election on its own, and the winner will need to dominate at least one bloc, and do well enough in the others.

The question on my mind is whether candidates could run for the same bloc, and split the vote. With 10 candidates announced, 20% of the vote could be enough to win. If three candidates each try to consolidate the Black Voter divisions, or the Progressive divisions, could a single candidate optimize for the White Moderate divisions and win? And what neighborhoods will candidates presumably do best in?

Below, I consider candidates’ past elections to guess which blocs they’ll be targetting.

View code
library(tidyverse)
library(sf)

source("../../admin_scripts/util.R")


setwd("C:/Users/Jonathan Tannen/Dropbox/sixty_six/posts/council_ballot_position_23/")

df_major_type <- readRDS("../../data/processed_data/df_major_type_20230116.Rds")
df_major <- df_major_type %>%
  group_by(office, candidate, party, warddiv, year, election_type, district, ward, is_topline_office) %>%
  summarise(votes = sum(votes))

df_major <- df_major %>% 
  group_by(year, election_type, office, district, warddiv) %>%
  mutate(pvote = votes / sum(votes)) %>%
  ungroup()

topline_votes <- df_major %>% 
  filter(is_topline_office) %>%
  group_by(election_type, year) %>%
  summarise(votes = sum(votes)) %>%
  mutate(
    year = asnum(year),
    cycle = case_when(
      year %% 4 == 0 ~ "President",
      year %% 4 == 1 ~ "District Attorney",
      year %% 4 == 2 ~ "Governor",
      year %% 4 == 3 ~ "Mayor"
    )
  )

cycle_colors <- c("President" = strong_red, "Mayor" = strong_blue, "District Attorney" = strong_green, "Governor" = strong_orange)

office_votes = df_major %>%
  group_by(
    year, warddiv, election_type, office, district
  ) %>%
  summarise(total_votes = sum(votes))

format_name <- function(x){
  x <- stringr::str_to_lower(x)
  x <- gsub("(^|\\s)([a-z])", "\\1\\U\\2", x, perl = TRUE)
  x
}
divs <- st_read("../../data/gis/warddivs/202011/Political_Divisions.shp") %>%
  mutate(warddiv = pretty_div(DIVISION_N))source("../../data/prep_data/div_svd_time_util.R")
div_cat_fn <- readRDS("../../data/processed_data/svd_time_20230116.RDS")

div_cats <- div_cat_fn %>% get_row_cats(2017) %>% rename(warddiv = row_id)

cats <- c(
  "Black Voters",
  "Wealthy Progressives",
  "Hispanic Voters",
  "White Moderates"
)

cat_colors <- c(light_blue, light_red, light_orange, light_green)
names(cat_colors) <- cats

What blocs will each candidate target?

For each candidate, let’s look at some relevant elections. First, the At Large Council candidates, for whom we have the best, city-wide, reference elections.

View code
relevant_elections <- tribble(
  ~candidate, ~year, ~election_type, ~office,
  "ALLAN DOMB", 2015, "primary", "COUNCIL AT LARGE",
  "ALLAN DOMB", 2019, "primary", "COUNCIL AT LARGE",
  "DEREK S GREEN", 2015, "primary", "COUNCIL AT LARGE",
  "DEREK S GREEN", 2019, "primary", "COUNCIL AT LARGE",
  "HELEN GYM", 2015, "primary", "COUNCIL AT LARGE",
  "HELEN GYM", 2019, "primary", "COUNCIL AT LARGE",
  "REBECCA RHYNHART", 2017, "primary", "CITY CONTROLLER",
  # "REBECCA RHYNHART", 2021, "primary", "CITY CONTROLLER",
) %>% mutate(year = as.character(year))

comp_results <- df_major %>% 
  inner_join(relevant_elections) %>%
  left_join(office_votes) %>%
  group_by(candidate, year, election_type, office) %>%
  mutate(
    pvote_total = sum(votes, na.rm=TRUE) / sum(total_votes, na.rm=TRUE)
  ) %>%
  ungroup()

cat_results <- comp_results %>% 
  left_join(div_cats %>% select(-year)) %>%
  group_by(year, cat, candidate, election_type, office) %>%
  summarise(pvote = sum(votes) / sum(total_votes))


ggplot(
  divs %>% 
    left_join(comp_results) %>%
    filter(office == "COUNCIL AT LARGE")
) +
  geom_sf(aes(fill=100 * pmin(pvote, 0.20)), color=NA) +
  facet_grid(format_name(candidate) ~ year) +
  scale_fill_viridis_c("% of Votes") +
  theme_map_sixtysix() %+replace%
  theme(legend.position = "right") +
  labs(
    title="City Council At Large Results"
  )

In 2015, Domb did best in the River Wards, the Northeast, South Philly, and immediate Center City. Derek Green did best among Black divisions in West, and North Philly. Helen Gym did best in the ring around Center City and Chestnut Hill / Mount Airy.

In 2019, we see the effects of incumbency. Gym led the pack by dominating the Wealthy Progressive and Black Wards from the first column, Green did well enough everywhere despite losing his top ballot position, and Domb improved across Black Wards thanks to incumbency and the party endorsement.

View code
ggplot(
  cat_results %>% filter(office == "COUNCIL AT LARGE"),
  aes(x=cat, y=100*pvote)
) +
  geom_bar(aes(fill=cat), color=NA, stat="identity") +
  facet_grid(year ~ format_name(candidate)) +
  scale_fill_manual(values=cat_colors, guide=FALSE) +
  theme_sixtysix() %+replace%
  theme(axis.text.x = element_text(angle=45, hjust = 1, vjust=1)) +
  labs(
    title="City Council At Large by Bloc",
    x=NULL,
    y = "% of Votes"
  )

Rhynhart won a competitive primary in 2017, albeit in the lower profile City Controller race. She did so by cleaning up in the Wealthy Progressive wards around Center City, though did beat 50% in the Black and Hispanic Voter blocs.

View code
ggplot(
  divs %>% 
    left_join(comp_results) %>%
    filter(office == "CITY CONTROLLER")
) +
  geom_sf(aes(fill=100 * pvote), color=NA) +
  # facet_grid(format_name(candidate) ~ year) +
  scale_fill_viridis_c("% of Votes") +
  theme_map_sixtysix() %+replace%
  theme(legend.position = "right") +
  labs(
    title="Rebecca Rhynhart in the 2017 Primary"
  )

View code
ggplot(
  cat_results %>% filter(office == "CITY CONTROLLER"),
  aes(x=cat, y=100*pvote)
) +
  geom_bar(aes(fill=cat), color=NA, stat="identity") +
  # facet_grid(year ~ format_name(candidate)) +
  scale_fill_manual(values=cat_colors, guide=FALSE) +
  theme_sixtysix() %+replace%
  theme(axis.text.x = element_text(angle=45, hjust = 1, vjust=1)) +
  labs(
    title="Rebecca Rhynhart's 2017 Primary by Bloc",
    x=NULL,
    y = "% of Votes"
  )

Three candidates come from districted seats: Cherelle Parker, Maria Quiñones-Sánchez, and Amen Brown. We don’t have city-wide elections for them, but can reasonably guess that they’ll do well in the blocs that their districts constitute.

View code
district_elections <- tribble(
  ~candidate, ~year, ~election_type, ~office,
  "CHERELLE L PARKER", 2019, "primary", "DISTRICT COUNCIL",
  "MARIA QUIÑONES SÁNCHEZ", 2019, "primary", "DISTRICT COUNCIL",
  "AMEN BROWN", 2022, "primary", "REPRESENTATIVE IN THE GENERAL ASSEMBLY"
) %>% mutate(year = as.character(year))

districts <- 
  divs %>% inner_join(
    df_major %>% 
    inner_join(district_elections) %>%
    select(office, district, warddiv, candidate)
  ) %>% group_by(office, district, candidate) %>% 
  summarise(
    geometry=st_union(geometry)
  )
districts <- nngeo::st_remove_holes(districts)


phila_whole <- st_read("../../data/gis/warddivs/201911/Political_Wards.shp") %>% 
  st_union()ggplot(phila_whole) + 
  geom_sf(color=NA) + 
  geom_sf(data=districts, aes(fill = candidate), color=NA) +
  geom_sf_text(
    data=districts %>% mutate(geom = st_centroid(geom)),
    aes(label = sprintf(
      "%s\n(%s-%s)", 
      format_name(candidate), 
      case_when(office == "DISTRICT COUNCIL" ~ "Council", TRUE ~ "PA"),
      district
    )),
    fontface="bold"
  ) +
  scale_fill_discrete(guide=F) +
  theme_map_sixtysix() +
  labs(
    title="Mayoral Candidates' Districts"
  )

Parker and Brown represent predominantly-Black districts, Quiñones-Sánchez a predominantly Hispanic one.

That leaves Warren Bloom, Jeff Brown, and Jimmy DeLeon as candidates without instructive elections. Brown is spending a lot of money, and running on a platform that I expect to do better in the White Moderate bloc. Bloom has run six times before without making any ripples, and DeLeon is a… longshot.

The crowded lanes

Where does that leave us? We have three candidates who have done best in Black wards, Parker, Green, and Amen Brown; two candidates who have done best in Wealthy Progressive wards, Gym and Rhynhart; two we can expect to do best in White Moderate wards, Domb and Jeff Brown; one who’ll do best in the Hispanic bloc, Quiñones-Sánchez; and two without an electoral history, Bloom and Deleon.

Obviously these candidates won’t split the blocs evenly, so we can’t just divide the vote by the number of candidates. And winning will require doing well enough in other blocs to hit a winning profile. The candidates from At Large council–Green, Gym, and Domb–probably enter with the best name recognition, and might have higher baselines even in their worse blocs.

That was clear in the last two competitive mayoral primaries of 2007 and 2015. In 2007, Nutter dominated the Wealthy Progressive bloc, but also won the Black Voter bloc and broke 20% in White Moderate and Hispanic Voter blocs. In 2015, Kenney did best in the White Moderate and Wealthy Progressive blocs, but actually won in all of them. Both Kenney and Nutter were City Councilmembers.

View code
df_major %>% 
  filter(
    office == "MAYOR",
    election_type == "primary",
    party == "DEMOCRATIC",
    year %in% c(2015, 2007)
  ) %>% 
  left_join(div_cats %>% select(-year)) %>%
  group_by(cat, year, candidate) %>%
  summarise(votes=sum(votes, na.rm=TRUE)) %>%
  group_by(cat, year) %>%
  mutate(pvote = votes / sum(votes)) %>%
  ggplot(
    aes(x = cat, y=100*pvote) 
  ) +
  geom_text(aes(label = format_name(candidate))) +
  facet_grid(year ~.) +
  geom_hline(yintercept=0)+
  theme_bw() +
  labs(
    x=NULL,
    y="% of Votes",
    title="Results in each bloc",
    subtitle="2007 and 2015 Mayoral Primaries"
  )

With petitions due soon and the race maturing, we will probably see several candidates drop out ahead of the election. That could help a candidate consolidate their bloc. Without drop-outs, though, this election would probably prove more fractured than 2007 and 2015. But the broader point holds: the winner will probably do well in two blocs, and well-enough in the rest.

Common Pleas Deep Dive, 2021

Belatedly, I’ve had time to sit down with the 2021 Primary results. Here are some observations.

In November, Philadelphia will elect eight new judges on the Court of Common Pleas. After the May Primary, we know almost certainly who those judges will be; the Democratic nominees will all win.

All eight Democratic nominees are Recommended by the Bar, three Highly. Surprisingly, they don’t include the person in the number one ballot position. And they won with a wide diversity of maps.

View code
library(dplyr)
library(tidyr)
library(ggplot2)
devtools::load_all("../../admin_scripts/sixtysix/")

ballot <- read.csv("../../data/common_pleas/judicial_ballot_position.csv")
res <- readxl::read_xlsx("C:/Users/Jonathan Tannen/Downloads/2021_primary (1).xlsx")
res <- res %>%
  pivot_longer(
    cols=`JUSTICE OF THE\r\nSUPREME COURT DEM\r\nMARIA MCLAUGHLIN`:`QUESTION #5\r\nNO`,
    names_to="candidate",
    values_to="votes"
  )
names(res) <- gsub("(\\r|\\n)+", " ", names(res))
names(res) <- gsub("\\s", "_", tolower(names(res)))

res$vote_type <- case_when(
  res$vote_type == "E" ~ "Election Day",
  res$vote_type == "M" ~ "Mail",
  res$vote_type == "P" ~ "Provisional"
)

res_cp <- res %>%
  filter(
    grepl("^JUDGE OF THE\r\nCOURT OF COMMON PLEAS DEM\r\n", candidate)
  ) %>%
  mutate(
    candidate = gsub("^JUDGE OF THE\r\nCOURT OF COMMON PLEAS DEM\r\n","", candidate)
  )

res_cp <- res_cp %>%
  left_join(ballot %>% filter(year == 2021) %>% mutate(candidate = toupper(name))) %>%
  mutate(name=format_name(name))

assertthat::assert_that(
  res_cp %>% filter(is.na(name)) %>% with(all(candidate == "Write-in"))
)

res_type <- res_cp %>%
  filter(!is.na(name)) %>%
  group_by(name, vote_type, rownumber, colnumber, philacommrec, dcc, inq) %>%
  summarise(votes=sum(votes)) %>%
  group_by(vote_type) %>%
  mutate(pvote = votes/sum(votes))
View code
res_total <- res_type %>% 
  group_by(name, rownumber, colnumber, philacommrec, dcc, inq) %>%
  summarise(votes=sum(votes), .groups="drop") %>%
  mutate(pvote = votes/sum(votes))

ggplot(
  res_total %>% arrange(votes) %>% mutate(winner = rank(-votes) <= 8),
  aes(y=rownumber, x=colnumber)
) +
  geom_tile(
    aes(fill=pvote*100, color=winner),
    size=2
  ) +
  geom_text(
    aes(
      label = ifelse(philacommrec==1, "R", ifelse(philacommrec==2,"HR","")),
      x=colnumber+0.45,
      y=rownumber+0.45
    ),
    color="grey70",
    hjust=1, vjust=0
  ) +
  geom_text(
    aes(
      label = ifelse(dcc==1, "D", ""),
      x=colnumber-0.45,
      y=rownumber+0.45
    ),
    color="grey70",
    hjust=0, vjust=0
  ) +
  geom_text(
    aes(label = sprintf("%s\n%0.1f%%", name, 100*pvote)),
    color="black"
    # fontface="bold"
  ) +
  scale_y_reverse(NULL) +
  scale_x_continuous(NULL)+
  scale_fill_viridis_c(guide=FALSE) +
  scale_color_manual(values=c(`FALSE`=NA, `TRUE`="yellow"), guide=FALSE) +
  annotate(
    "text",
    label="R = Recommended\nHR = Highly Recommended\nD = DCC Endorsed",
    x = 1.6,
    y = 6,
    hjust=0,
    vjust=0.5,
    color="grey70"
  ) +
  theme_sixtysix() %+replace% 
  theme(
    panel.grid.major=element_blank(),
    axis.text=element_blank()
  ) +
  ggtitle(
    "Common Pleas Results",
    "2021 Democratic Primary, arranged as on the ballot. Winners are outlined."
  )

Four candidates won in the first column, three in the second, and one in the third. Three winners were Highly Recommended by the Bar, including Michele Hangley in the second column and Chris Hall in the third, but by itself that rating wasn’t sufficient: John Padova and Mark Moore failed to capitalize on it. There’s some additional work needed to use it to your advantage.

The candidates’ maps are diverse. Nick Kamau and Cateria McCabe won everywhere, though slightly stronger in the Black wards of West and North Philly (and decidedly not the Northeast). Wendi Barish also won everywhere, slightly stronger in Center City and its ring. Betsy Wahl, Chris Hall, and Michele Hangley all won thanks to their strength in the Wealthy Progressive ring around Center City and in Chestnut Hill and Mount Airy. Craig Levin did the opposite, winning the Northeast and West and North Philly, presumably on the strength of his DCC endorsement. And Dan Sulman was the eighth and final winner, with the bright yellow 53rd ward just enough to push him through, where his sister is the Ward Leader.

View code
library(sf)

divs <- st_read("../../data/gis/warddivs/202011/Political_Divisions.shp") %>%
  mutate(warddiv = pretty_div(DIVISION_N))wards <- st_read("../../data/gis/warddivs/201911/Political_Wards.shp") %>%
  mutate(ward=sprintf("%02d", asnum(WARD_NUM)))res_ward_type <- res_cp %>%
  mutate(ward = substr(division, 1, 2)) %>%
  group_by(ward, name, vote_type) %>%
  summarise(votes=sum(votes)) %>%
  group_by(vote_type) %>%
  mutate(pvote=votes/sum(votes))

res_ward <- res_ward_type %>%
  group_by(ward, name) %>%
  summarise(votes=sum(votes)) %>%
  group_by(ward) %>%
  mutate(pvote=votes/sum(votes))

res_ward <- left_join(wards, res_ward)

candidate_order <- res_total %>% arrange(desc(votes)) %>% with(name)

ggplot(
  res_ward %>% 
    filter(!is.na(name)) %>%
    mutate(name=factor(name, levels=candidate_order))
) +
  geom_sf(aes(fill=100*pvote), color=NA) +
  scale_fill_viridis_c("Vote %") +
  facet_wrap(~ name) +
  theme_map_sixtysix() %+replace% theme(legend.position="right") +
  ggtitle("Common Pleas Results", "2021 Democratic Primary")

Caroline Turner was the first runner up, and the first candidate to fail to win from the top ballot position since at least 2007 (which is all the ballot layouts I can find). But she did really well in the 1st and 2nd Wards, which now deserve a name.

The Reclaim Wards

When clicking through the results online, I saw a cut that made me laugh out loud.

Results from Division 01-01

The top eight winners in 01-01 each received more than 9.99% of the vote. Ninth place? Only 2.5%. This is the kind of electoral coordination party bosses dream of.

In fact, that consolidation is true of the entire first ward (covering East Passyunk in South Philly).

View code
ward_bar <- function(ward, endorsements){
  df <- res_ward %>%
  as.data.frame() %>%
  filter(ward==!!ward, !is.na(name)) %>%
  arrange(desc(votes)) %>%
  # filter(1:n() <= 10) %>%
  mutate(name=factor(name, levels=name)) %>%
  mutate(
    last_name = gsub(".* ([A-Za-z]+)$", "\\1", name),
    endorsed= last_name %in% endorsements,
    reclaim = last_name %in% c("Hall", "Hangley", "Kamau", "Barish","Sulman", "McCabe", "Turner", "Wahl")
  )
  
  ggplot(df, aes(x=name, y=votes)) +
  geom_bar(stat="identity", aes(color=endorsed, fill=reclaim), size=1.2) +
      geom_vline(xintercept=8.5, linetype="dashed") +
  scale_color_manual(
    NULL,
    values=c(`TRUE` = "goldenrod", `FALSE`=NA),
    labels=sprintf(c("Not %s Endorsed", "%s Endorsed"), ward)
  ) +
  scale_fill_manual(
    NULL,
    values=c(`TRUE`="grey30", `FALSE`="grey60"),
    labels=c("Not Reclaim", "Reclaim")
  ) +
  theme_sixtysix() %+replace% 
  theme(axis.text.x = element_text(angle=45, hjust=1, vjust=1))+
  labs(
    title=sprintf("Common Pleas Results in Ward %s", ward),
    subtitle="2021 Democratic Primary",
    x=NULL,
    y="Votes"
  ) 

}

ward_bar(
  "01", 
  c("Hall", "Hangley", "Kamau", "Barish","Sulman", "McCabe", "Turner", "Wahl")
)  

It’s usually impossible to separate the many overlapping endorsements. Was it the Bar that brought the win, the DCC, or the Ward? But these eight winners are exactly the ward’s endorsed candidates. They were also the full Reclaim slate, so it’s impossible to separate the Ward’s power from Reclaim. But the deciding factors were almost certainly these two.

The four wards with the biggest gap between eighth and ninth place–suggesting the strongest slate power–were South Philly’s 1st and 2nd and West Philly’s 27th and 46th.

The 2nd Ward, just to the 1st’s North, had slightly different endorsements than Reclaim. The six candidates who had both a Reclaim and 2nd Ward endorsement did best. Barish came in seventh with only a Reclaim endorsement, Levin in eighth with only the 2nd, and Sulman in ninth with only Reclaim.

View code
ward_bar(
  "02", 
  c("Hall", "Hangley", "Kamau", "Levin", "McCabe", "Turner", "Wahl")
)  

In the 27th (where, full disclosure, I’m a committeeperson), the Reclaim endorsements appear to have carried the day: Turner and Hall won, while 27-endorsed Moore and Levin didn’t.

View code
ward_bar(
  "27",
      c("Barish", "Moore", "Hangley", "Kamau", "Levin", "McCabe", "Sulman", "Wahl")
)

I haven’t found the 46th ward endorsements, but the Reclaim slate cleaned up.

View code
ward_bar(
  "46",
      c()
) + labs(subtitle="2021 Primary. Endorsements not available.")

Not only was the Reclaim slate particularly strong in these wards, but the gap between eighth and ninth position make it clear it was the Reclaim endorsement itself, and not one of the other pregressive slates that drove voters.

But the Reclaim endorsement wasn’t itself enough to win across the city. Caroline Turner came in ninth despite it and first ballot position, mostly due to poor results in the Black wards.

Wealthy Progressive divisions did consolidate their votes in a way other divisions didn’t. Grouping divisions by my Voting Blocs shows that the Wealthy Progressive divisions’ preferred candidates did better there than the preferred candidates of other blocs’.

View code
devtools::load_all("C:/Users/Jonathan Tannen/Dropbox/sixty_six/posts/svdcov/")
svd_time <- readRDS("../../data/processed_data/svd_time_20210813.RDS")
div_cats <- get_row_cats(svd_time, 2020)

# ggplot(divs %>% left_join(div_cats, by=c("warddiv"="row_id"))) +
#   geom_sf(aes(fill=cat), color=NA)

res_cat <- res_cp %>% left_join(div_cats, by=c("division"="row_id")) %>%
  filter(!is.na(name)) %>%
  group_by(name, cat) %>%
  summarise(votes=sum(votes)) %>%
  group_by(cat) %>%
  mutate(
    rank = rank(desc(votes)),
    pvote=votes/sum(votes)
  )

ggplot(res_cat, aes(x=rank, y=100*pvote)) +
  geom_bar(stat="identity", fill="grey50") +
  facet_grid(cat ~ .) +
  geom_text(
    aes(label=gsub(".* ([A-Za-z]+)$", "\\1", name)),
    y=0.3,
    angle=90,
    hjust=0
  )+
  theme_sixtysix() +
  geom_vline(xintercept=8.5, linetype="dashed") +
  labs(x="Rank", y="% of Vote", title="Results by Division Bloc")

The top eight candidates in the Wealthy Progressive divisions received on average 9.1% of the vote, compared to 8.2% in the Black Voter divisions, 7.9% in the Hispanic Voter, and 7.5% in the White Moderate. To measure another way, Common Pleas candidates had a gini coefficient of votes of 0.29 in Wealthy Progressive divisions, compared to 0.22, 0.19, and 0.14 in the other three blocs, indicating greater inequality in the votes candidates received, and thus more separation.

This is probably due to voters there being more likely to look up recommendations either on or before election day, and having consolidated preferences.

View code
gini <- function(x){
  outer_sum <- outer(x, x, FUN="-")
  gini <- sum(abs(outer_sum)) / (2 * length(x)^2 * mean(x))
  return(gini)
}

# res_cat %>%
#   group_by(cat) %>%
#   summarise(
#     gini = gini(pvote),
#     mean=mean(pvote[rank <= 8])
#   )

Mail-In Votes

Entering this election, I was especially interested if mail-in ballots would have different dynamics than in-person voting. When people vote at the kitchen table, likely over days, will ballot position matter less? Will endorsements matter more?

In total, 33% of the votes for CP came by Mail, vs 66% on Election Day (and 1% Provisionals). The Wealthy Progressive divisions were more likely to use mail: 38% of their votes were by mail, compared to 34% in the White Moderate divisions, and 29 and 28% in the Black and Hispanic Voter divisions.

View code
# res_cp %>% group_by(vote_type) %>%
#   summarise(votes=sum(votes)) %>%
#   mutate(pct=votes/sum(votes))
# 
# res_cp %>% 
#   left_join(div_cats, by=c("division"="row_id")) %>%
#   group_by(cat, vote_type) %>%
#   summarise(votes=sum(votes)) %>%
#   group_by(cat) %>%
#   mutate(pct=votes/sum(votes))

The candidates who did better by mail were all in the bottom right of the ballot. The top three were also Highly Recommended, suggesting that endorsements were more likely to be looked up by people voting by mail.

View code
res_votetype <- res_cp %>%
  filter(!is.na(name)) %>%
  group_by(vote_type, name, rownumber, colnumber, philacommrec, dcc) %>%
  summarise(votes=sum(votes)) %>%
  group_by(vote_type) %>%
  mutate(pvote=votes/sum(votes)) %>%
  group_by(name) %>%
  mutate(votes_total=sum(votes)) %>%
  ungroup() %>%
  pivot_wider(names_from=vote_type, values_from=c(votes, pvote)) %>%
  mutate(pvote_total=votes_total / sum(votes_total)) %>%
  arrange(desc(pvote_Mail - `pvote_Election Day`))

ggplot(
  res_votetype %>% arrange(votes_total) %>% mutate(winner = rank(-votes_total) <= 8) %>%
    mutate(diff=pvote_Mail - `pvote_Election Day`),
  aes(y=rownumber, x=colnumber)
) +
  geom_tile(
    aes(fill=100*diff),
    size=2
  ) +
  geom_text(
    aes(
      label = ifelse(philacommrec==1, "R", ifelse(philacommrec==2,"HR","")),
      x=colnumber+0.45,
      y=rownumber+0.45
    ),
    color="grey70",
    hjust=1, vjust=0
  ) +
  geom_text(
    aes(
      label = ifelse(dcc==1, "D", ""),
      x=colnumber-0.45,
      y=rownumber+0.45
    ),
    color="grey70",
    hjust=0, vjust=0
  ) +
  geom_text(
    aes(
      label = sprintf("%s\n%s%0.1f%%, (%0.1f%% - %0.1f%%)", name, ifelse(diff>0, "+", "-"),100*abs(diff), 100*pvote_Mail,100* `pvote_Election Day`)
    ),
    color="black"
    # fontface="bold"
  ) +
  scale_y_reverse(NULL) +
  scale_x_continuous(NULL)+
  scale_fill_viridis_c(guide=FALSE) +
  annotate(
    "text",
    label="R = Recommended\nHR = Highly Recommended\nD = DCC Endorsed",
    x = 1.6,
    y = 6,
    hjust=0,
    vjust=0.5,
    color="grey70"
  ) +
  theme_sixtysix() %+replace% 
  theme(
    panel.grid.major=element_blank(),
    axis.text=element_blank()
  ) +
  ggtitle(
    "Common Pleas: Mail minus In Person.",
    "2021 Democratic Primary, arranged as on the ballot."
  )

Some of these differences are due to differential rates of mail-in usage: Moore, Hangley, and Wahl all did better in Wealthy Progressive wards, and they mailed in more often. We can adjust that by taking the within-Division difference for each candidate, and then taking a weighted average across divisions weighted by total votes. This decomposition leads to basically the same finding but slightly smaller effects: Moore did 1.9 percentage points better by mail within a given division, Hangley 1.2, and Hall 0.9.

Some of this may still be because the mail voters were systematically different from in-person voters, but this within-division comparison is the closest we can get with the available data.

Doing the decomposition by Voting Bloc is interesting. Moore does better everywhere by mail than in person, probably reflecting the important DCC and Bar effects but a lack of endorsements with Election Day presence. In every bloc, the highest differences are the poor-ballot-position, Bar-recommended candidates. Interestingly, the largest differences are in the White Moderate divisions (South Philly and the Northeast), probably reflecting the politicized nature of mail-in voting, and that in those divisions there were political differences between those who voted by mail and in person.

View code
res_votetype_weighted <- res_cp %>%
  select(division, name, vote_type, votes) %>%
  filter(!is.na(name)) %>%
  group_by(division, vote_type) %>%
  mutate(pvote=votes/sum(votes)) %>%
  group_by(division) %>%
  mutate(total_votes=sum(votes)) %>%
  select(-votes) %>%
  pivot_wider(names_from=vote_type, values_from=pvote) %>%
  mutate(diff = `Mail` - `Election Day`) %>%
  group_by(name) %>%
  summarise(
    diff=weighted.mean(diff, w=total_votes, na.rm = T)
  )

res_votetype_weighted_cat <- res_cp %>%
  select(division, name, vote_type, votes) %>%
  filter(!is.na(name)) %>%
  group_by(division, vote_type) %>%
  mutate(pvote=votes/sum(votes)) %>%
  group_by(division) %>%
  mutate(total_votes=sum(votes)) %>%
  select(-votes) %>%
  left_join(div_cats, by=c("division"="row_id")) %>%
  pivot_wider(names_from=vote_type, values_from=pvote) %>%
  mutate(diff = `Mail` - `Election Day`) %>%
  group_by(name, cat) %>%
  summarise(
    diff=weighted.mean(diff, w=total_votes, na.rm = T)
  )

ggplot(
  res_votetype_weighted_cat %>% group_by(cat) %>%
    mutate(rank=rank(desc(diff))),
  aes(x=rank, y=100*diff)
) +
  geom_bar(stat="identity") +
  facet_grid(cat~.) +
  theme_sixtysix() +
  geom_text(aes(label= gsub(".* ([A-Za-z]+)$", "\\1", name)), angle=90, hjust=0, y=0.1)+
  labs(
    title="Mail minus In-Person Results", subtitle="By Voting Bloc",
    y="Mail % minus In-Person %",
    x=NULL
  )

Next: The Effect of the Bar Recommendations!

Coming soon.