↖︎ Vishal Singh

Case 4 of 5 · Year-over-year comparisons and omitted-variable bias

Case 4 of 5 · Year-over-year comparisons and omitted-variable bias

Coca-Cola, a rumor and a different constituency

In the four weeks from February 3, 2025, Coca-Cola's share of soda units fell 4.25 points in ZIPs more than half Hispanic and 0.63 points in ZIPs under 10% Hispanic; a partisan reading of the same data disappears once Hispanic share is in the regression.

Author
AffiliationNYU Stern
Published

In the four weeks from February 3, 2025, Coca-Cola's share of soda units fell 4.25 points in ZIPs more than half Hispanic and 0.63 points in ZIPs under 10% Hispanic; a partisan reading of the same data disappears once Hispanic share is in the regression.

Author

Vishal Singh

Affiliation

NYU Stern School of Business

Published

14 September 2026

Updated

2 October 2026

In the four weeks from February 3, 2025, Coca-Cola's share of soda units fell 4.25 percentage points in ZIP codes more than half Hispanic and 0.63 points in ZIPs under 10% Hispanic. Both figures are measured against the same four weeks of 2024 and against each ZIP's own gap to 2024 in the three weeks before. The data are convenience-store sales in 2,469 ZIPs. A regression of the same change on Republican vote share alone gives a slope of +0.24 points per 10 points of Republican share (standard error 0.04), which reads as Democratic-leaning ZIPs cutting back more. With Hispanic share in the regression the slope is +0.07 (0.04), and comparing ZIPs only within the same county it is −0.06 (0.07). This case builds the comparison one step at a time and shows how an omitted variable produced the first number.

A boycott list, a day of action and a false rumor

In late January 2025 a "Latino Freeze" spending boycott circulated on TikTok. Its do-not-shop list named Coca-Cola alongside Target, Walmart and Amazon. Monday, February 3 was a national "Day Without Immigrants", a call to skip work, school and shopping. Between February 5 and 8, posts in English and Spanish claimed that Coca-Cola had called immigration agents on Latino workers at a Texas bottling plant. Fact-checkers at Snopes, PolitiFact and Lead Stories found no evidence, and the plant named in the posts does not exist.

Case 3 found Bud Light's loss concentrated in Republican-leaning ZIPs after the brand's own marketing. Here Coca-Cola had changed nothing about its products, prices or marketing, and the people called on to act were defined by ethnicity rather than party. This case asks whether the method from case 3 still works.

Timeline, with the week numbers used in the data (week 0 starts on Monday, February 3, 2025)
DateWeekEvent
Dec 30 to Jan 19−5 to −3Baseline weeks
Jan 20−2Presidential inauguration; the week is left out of the baseline
Late January−2 to −1The "Latino Freeze" list, naming Coca-Cola, circulates on TikTok
Jan 31 to Feb 2−1Weekend on which sales begin to fall in ZIPs over 30% Hispanic
Feb 30"Day Without Immigrants"
Feb 5 to 80False posts claim Coca-Cola called immigration agents on Latino workers at a Texas plant
Feb 8 to 90U.S. searches for a Coca-Cola boycott peak
AfterwardsFact-checkers publish findings that the claim has no evidence behind it

The data: 2,469 ZIP codes and 79 weeks of soda sales

The unit is a ZIP code and a week. The sales are point-of-sale records from convenience stores, with at least two stores pooled in every ZIP in every week.Data. Weekly share of carbonated soft drink units sold, January 1, 2024 to June 30, 2025, from convenience-store point-of-sale data (Dewey). ZIP demographics are Census and election-return estimates supplied with the case files. No survey weights: each ZIP counts once.

The data
UnitZIP code by week
ZIP codes2,469, in 47 states
Weeks79, from January 1, 2024 to June 30, 2025 (weeks −57 to +21 around February 3, 2025)
OutcomeCoca-Cola (regular, Diet and Zero together) as a percent of the ZIP's soda units
Other brandsPepsi, Sprite, Fanta, Mountain Dew, Dr Pepper, Mexican-import sodas such as Jarritos, all other sodas; the shares add to 100
ZIP traitsHispanic share, foreign-born share, Black share, adults with some college, median income, density, Republican share of the two-party vote in 2020

Coca-Cola's level differs by neighborhood before anything happens. Over the 57 weeks before February 3, 2025 it averaged 24.0% of soda units in ZIPs under 10% Hispanic and 39.1% in ZIPs more than half Hispanic. A comparison of levels would mix the boycott with that permanent gap, so every comparison below is a change.

ZIP groups by Hispanic share
Hispanic share of the ZIPZIPsCoca-Cola share before the event, % of soda units
Under 10%1,70124.0
10 to 30%52929.0
30 to 50%13733.5
50% or more10239.1

Step 1: compare each week with the same week of 2024

Soda sales may have calendar patterns, and the Super Bowl falls in the first days of February. The cleanest correction subtracts, for each ZIP, its share in the same week 52 weeks earlier:

Δi,t = Sharei,t − Sharei,t−52

Anything that recurs in the same week every year cancels. The effect for ZIP i is the average of Δ over weeks 0 to 3 minus its average over the three baseline weeks, −5 to −3 (December 30 to January 19):

Effecti = 14 ∑3k=0 Δi,k − 13 ∑−3k=−5 Δi,k

Weeks −2 and −1 are in neither period. The inauguration fell in week −2, and a dip had begun in week −1.

Interactive

Source: PDI convenience-store soda sales (Dewey), 2,469 ZIPs, January 2024 to June 2025. Unweighted means across ZIPs in each Hispanic-share group: the raw view shows ZIPs under 10% Hispanic (1,701) and 50% or more (102); the change view shows all four groups (1,701, 529, 137 and 102). The shaded band is February 5 to 8.

In these data the calendar correction is small. In every week from January 1 to July 1, 2024 the share in ZIPs 50% or more Hispanic stayed between 38.4% and 39.8%, and subtracting the same calendar weeks of 2024 moves a group's estimate by at most 0.2 points.

Change from the three baseline weeks to weeks 0 to 3, percentage points of soda units
ZIP group2025 onlySame calendar weeks of 20242025 minus 2024 (the estimate)
Under 10% Hispanic−0.85−0.22−0.63
10 to 30%−2.37+0.07−2.44
30 to 50%−3.48−0.03−3.45
50% or more−4.40−0.15−4.25
Slope on Hispanic share, per 10 points (s.e.)−0.68 (0.04)+0.03 (0.03)−0.71 (0.05)

The year-over-year version costs some precision: the standard error of the slope rises from 0.04 to 0.05 because the 2024 weeks add their own noise. In return it gives a placebo. Run on 2024, when nothing happened, the same arithmetic gives a slope of +0.03 (0.03) per 10 points of Hispanic share, so the 2025 gradient is not a calendar pattern that recurs every February.3

Step 2: the loss rises with Hispanic share and the first dip comes before the rumor

Figure 2 plots the weekly change for four groups of ZIPs, each relative to its own baseline weeks. Over weeks 0 to 3 the mean change is −0.63 points in ZIPs under 10% Hispanic, −2.44 in ZIPs 10 to 30% Hispanic, −3.45 in ZIPs 30 to 50% and −4.25 in ZIPs 50% or more. The 50%-or-more group is 4.5 points below its baseline in week 0 and 5.0 in week 1, its lowest point.

Interactive

Coca-Cola's share fell in order of Hispanic share from the week of February 3, and the 50%-or-more group had recovered by late June

The groups separate in the week of February 3 and converge again over five months. Mean weekly change in Coca-Cola's share of soda units against the same week of 2024, minus each ZIP's average over weeks −5 to −3, in percentage points. Hover for the value and 95% interval in each week.

Source: PDI convenience-store soda sales (Dewey), 2,469 ZIPs, week of December 30, 2024 to June 30, 2025; unweighted ZIP means. Intervals are the mean plus or minus 1.96 standard errors across ZIPs.

The groups begin to separate before the rumor. In the week of January 27, which contains the weekend of January 31 to February 2, the 30-to-50% group was 1.5 points below its baseline (95% interval 0.8 to 2.1) and the 10-to-30% group 0.7 below (0.3 to 1.0). The weekly series is noisy for the smaller groups. The 30-to-50% group, with 137 ZIPs, was already 1.0 below in week −2, and the 50%-or-more group, with 102, was 0.5 below in week −1 with an interval of −1.2 to +0.3. The posts of February 5 to 8 fall inside week 0, together with the Day Without Immigrants, so the weekly data cannot give the rumor's own contribution to the 4.5-point drop in that week.

Step 3: one dot per ZIP

As in case 3, each ZIP is reduced to one number, its effect over weeks 0 to 3, and plotted against its Hispanic share. One number per ZIP keeps a ZIP's seven weeks from counting as seven independent observations, which would make the standard errors too small.2 The fitted line is

Effecti = a + bH × Hispanic sharei + ei,   bH = −0.071 per point

or −0.71 points for each 10 points of Hispanic share (standard error 0.05). A ZIP that is 50% Hispanic lost about 3.6 points more than a ZIP with almost no Hispanic residents. Hispanic share accounts for 10.7% of the variation in ZIP changes, which have a standard deviation of 3.3 points, so the line is a tendency with wide scatter around it.

The menu in Figure 3 refits the line within bands of Republican vote share. The slope is between −0.58 and −0.82 per 10 points in each of the four bands, which is the first sign that the Hispanic gradient does not depend on partisanship.

Interactive
ZIPs in the fit
Slope per 10 points of Hispanic share
95% interval
Mean change in the fit sample

Source: PDI convenience-store soda sales (Dewey), 2,469 ZIPs; ZIP traits from the Census (Hispanic share) and 2020 presidential returns (Republican share). OLS with heteroskedasticity-robust standard errors, unweighted. Change is the average of weeks 0 to 3 minus the average of weeks −5 to −3, each against the same week 52 weeks earlier.

Step 4: the partisan slope is Hispanic share in disguise

An analyst fresh from case 3 regresses the same effect on Republican vote share instead:

Effecti = a + bR × Republican sharei + ei,   bR = +0.024 per point

That is +0.24 per 10 points, with a standard error of 0.04. ZIPs that voted more Republican in 2020 lost less, which seems to support the headline "Democratic-leaning ZIPs boycotted Coca-Cola". Put Hispanic share into the regression as well and the Republican slope falls to +0.07 (0.04), 72% smaller and within two standard errors of zero. Comparing each ZIP only with ZIPs in the same county gives −0.06 (0.07) for Republican share and −0.43 (0.11) for Hispanic share.

Change in Coca-Cola's share, weeks 0 to 3, regressed on ZIP traits (points per 10 points of the trait; standard errors in parentheses)
Regressor(1) Republican only(2) Hispanic only(3) Both(4) Both, same county
Hispanic share−0.71 (0.05)−0.69 (0.05)−0.43 (0.11)
Republican share+0.24 (0.04)+0.07 (0.04)−0.06 (0.07)
ZIPs2,4692,4692,4691,819

Columns 1 to 3: ZIP-level OLS with robust standard errors. Column 4: ZIP-by-week panel (weeks −5 to −3 and 0 to 3) with ZIP and county-by-week fixed effects and standard errors clustered by ZIP; only counties with at least two ZIPs enter.

The workbench below runs these regressions in your browser on the 2,469 ZIP changes. Start with Republican share alone, then add Hispanic share, then everything. With all seven regressors the Republican coefficient is −0.12 (0.07) and the Hispanic coefficient −0.63 (0.08). Foreign-born share alone gives −1.16 (0.10) per 10 points because it moves with Hispanic share across ZIPs (correlation 0.73); in the full model it falls to −0.48 (0.17), and the split between the two is imprecise.

Interactive
ZIPs in the regression
R²
Hispanic coefficient, per 10 points
Republican coefficient, per 10 points

Where the short slope comes from

Decompose the slope on
Adding Hispanic share to a regression on Republican share cuts the Republican coefficient from +0.24 to +0.07. Each regression uses the 2,469 ZIP-level changes in Coca-Cola's share over weeks 0 to 3 as the outcome; dots are coefficients and lines 95% intervals. Percentage-point regressors are in units of 10 points, income in $10,000 and density in 1,000 people per square mile. The decomposition shows how much of the short regression's slope is borrowed from each other selected variable. Under "same county" every variable is measured as a deviation from its county mean, using counties with at least two ZIPs, with robust standard errors that are a little narrower than the clustered ones in the table above.

Source: PDI convenience-store soda sales (Dewey), 2,469 ZIPs; ZIP traits from the Census and 2020 presidential returns. OLS with heteroskedasticity-robust (HC1) standard errors, computed in your browser; unweighted.

Why the short slope is too large

The omitted-variable formula says by how much.1 Leave Hispanic share out and the slope on Republican share absorbs part of Hispanic share's relation to the outcome:

bshortR = blongR + blongH × δ

where δ is the slope from regressing Hispanic share on Republican share. Here δ = −0.25: ZIPs 10 points more Republican have, on average, 2.5 points less Hispanic population. In units of 10 points the formula reads

+0.240 = +0.067 + (−0.687) × (−0.252)

The product, +0.173, is 72% of the short slope. The bias needs two conditions: Hispanic share must move with the outcome (bH is not zero) and with the included variable (δ is not zero). Both are negative here, so their product is positive and the short slope is too high. In per-point units the same check is +0.0067 + (−0.0687) × (−0.2519) = +0.0240. The decomposition panel in Figure 4 recomputes it for any set of regressors.

Step 5: Pepsi's gain was about three-quarters of Coca-Cola's loss

Shares add to 100, so every point Coca-Cola lost went to another brand. Per 10 points of Hispanic share, Pepsi's share rose 0.54 points (standard error 0.03) over weeks 0 to 3 while Coca-Cola's fell 0.71, so Pepsi's gain was 76% of the loss. With Republican share also in the regression the figures are +0.49 and −0.69, or 72%. All other sodas together gained 0.17 per 10 points (0.03). The Mexican-import sodas, the brands most closely tied to Hispanic buyers, gained 0.02 (0.01), and none of Sprite, Fanta, Mountain Dew or Dr Pepper moved by more than 0.03.

Interactive

Source: PDI convenience-store soda sales (Dewey), 2,469 ZIPs, week of December 30, 2024 to June 30, 2025; unweighted. Each brand's change is its share minus its share 52 weeks earlier, relative to weeks −5 to −3. "All other sodas" is Mountain Dew, Dr Pepper, Sprite, Fanta and unbranded other sodas.

Step 6: the Hispanic gradient shrank in every block of weeks and was mostly gone by June

Averaged over the 22 weeks from February 3 to June 30, the Hispanic gradient is −0.35 per 10 points (standard error 0.04), half the −0.71 of the first four weeks. Figure 6 splits the period into four-week blocks. The gradient is −0.47 in weeks 4 to 7, −0.31 in weeks 8 to 11 and −0.12 in weeks 20 to 21 (95% interval −0.23 to −0.01). By week 20 the 50%-or-more group, with a mean change of +0.2 points (standard error 0.5), was back to its 2024 pattern, while the 30-to-50% group was still 0.8 points down (0.4). Unlike Bud Light, which case 3 finds had not recovered through 2024, the Coca-Cola gap closed in about five months.

Interactive

The Hispanic-share gradient fell from −0.71 per 10 points in weeks 0 to 3 to −0.12 in weeks 20 to 21

The gradient shrinks in every block, and the Republican slope stays near zero once Hispanic share is in the model. Slope of the change in Coca-Cola's share on ZIP characteristics, by four-week block after February 3, 2025, in points per 10 points of the characteristic. Choose whether the two variables enter separate regressions or one.

Source: PDI convenience-store soda sales (Dewey), 2,469 ZIPs, weeks 0 to 21 (February 3 to June 30, 2025). Each block is the average of the weekly changes against the same week of 2024, relative to weeks −5 to −3; the last block has two weeks. Unweighted OLS with robust standard errors; ZIPs treated as independent.

Two boycotts compared

Bud Light, 2023, and Coca-Cola, 2025
Bud Light, 2023Coca-Cola, 2025
TriggerThe brand's own marketingA boycott list and a false rumor; the brand changed nothing
Where the loss concentratedRepublican ZIPsHispanic ZIPs; no partisan gradient once Hispanic share is in the model
SpeedBreak in the first week, most of the drop within a monthBreak within days, deepest in the second week
PersistenceNo recovery through 2024Largely recovered in about five months
Where the share wentRival light lagers and other beerPepsi

The Bud Light loss followed the Republican vote and had not recovered through 2024. The Coca-Cola loss followed Hispanic share, and its gradient was −0.12 per 10 points by weeks 20 to 21.

What these data cannot show

ZIP characteristics describe places, not people. A ZIP that is 50% Hispanic has residents who are not Hispanic, and the data do not say who bought less Coca-Cola there. A slope across ZIPs compares places with different Hispanic shares; reading it as the behavior of Hispanic shoppers infers individual behavior from group data and can be wrong in either direction.

The stores are convenience stores, and the ZIPs lie mostly in the South (1,309) and Midwest (753), with 227 in the West and 180 in the Northeast. Grocery chains, restaurants and fountain sales are not covered.

Sales in these stores cannot separate a boycott from anything else that changed in the same weeks in the same places. Tariff news, a promotion, a distribution change or a shift in which stores report would appear as a 2025 change in Coca-Cola's share if it moved with Hispanic share. The 2024 placebo rules out a calendar pattern that repeats each year, not an event specific to 2025. The year-over-year comparison also assumes that, without the boycott, each ZIP's gap to 2024 would have stayed near its baseline level, a form of the parallel-trends assumption that the baseline weeks can check only in part.3

The first week of the loss partly pre-dates the rumor, and week 0 contains both the Day Without Immigrants and the posts of February 5 to 8, so these data give no separate estimate of the rumor. The groups of 102 and 137 ZIPs give single-week estimates with standard errors of 0.3 to 0.6 points. Standard errors in the workbench treat ZIPs as independent, but ZIPs in the same county are correlated: clustering by county instead of by ZIP raises the within-county Hispanic standard error from 0.11 to 0.12 per 10 points. Finally, a fall in Coca-Cola's share and a rise in Pepsi's describe the mix of units sold, and say nothing here about whether total soda sales changed.

Questions for discussion

  1. The rumor was false and the plant does not exist, yet sales fell within days. What should Coca-Cola have done in the first 72 hours, who should have heard it from the company, and what would a public denial risk?
  2. Using only weekly store sales by ZIP, design a monitoring report that would have flagged the break in the week of February 3 or earlier. Which neighborhoods would it group, what threshold would trigger a call, and how many false alarms would you accept?
  3. A colleague sees column 1 of the regression table and proposes a campaign in Democratic-leaning ZIPs. What do you ask before approving it, and which column do you show in reply?
  4. Suppose Hispanic share were positively correlated with Republican share (δ = +0.25) and the other numbers stayed the same. Compute the short slope on Republican share from the formula and say whether the short regression would then understate or overstate the long one.
  5. The weekly data show a dip in the week of January 27, before the rumor. List two explanations that do not involve the rumor and say what additional data, such as daily sales, would distinguish them.
  6. Bud Light's loss had not recovered through 2024 and Coca-Cola's gradient was −0.12 per 10 points by weeks 20 to 21. Which features of the two episodes might explain the difference, and which of them can these data speak to? Does Pepsi's gain in the same ZIPs change how you would measure the recovery?

Replicate this analysis

The folder 03_bud_light_boycott/extension_coke_2025 has the weekly ZIP panel (coke_zip_week.csv, 195,051 rows), the ZIP demographics and a README. The steps: pivot to one row per ZIP and one column per week; subtract each week's share from the share 52 weeks earlier; average weeks 0 to 3 and weeks −5 to −3 and take the difference; regress it on pct_hispanic, then add rep_share_2020; regress pct_hispanic on rep_share_2020 and check that the short slope equals the long slope plus bH times δ. In Excel the regressions are LINEST on the ZIP-level changes, and a pivot table by Hispanic-share bin reproduces the group means in Figure 2.

Data sources

Soda sales: PDI Technologies point-of-sale data from convenience stores, obtained through Dewey under a research licence. ZIP demographics: U.S. Census Bureau American Community Survey estimates and Republican share of the two-party vote in the 2020 presidential election, as supplied in zip_demographics.csv. The false claim about a Texas plant was checked by Snopes, PolitiFact and Lead Stories; the README names them and this page does not link to individual articles. The files are derived from licensed data; check the Dewey and PDI licence terms before posting them outside the course.

Methods note

Sample: 2,469 ZIPs with at least two stores reporting in each of 79 weeks (January 1, 2024 to June 30, 2025); week 0 is the week of Monday, February 3, 2025 and the data end with week 21. The year-over-year change Δ is a ZIP's share minus its share 52 weeks earlier, so the earliest usable week is −5. A ZIP's effect over weeks a to b is the mean of Δ over those weeks minus its mean over weeks −5 to −3. The pre-event levels (24.0%, 29.0%, 33.5% and 39.1%) are ZIP means over weeks −57 to −1. ZIP groups use Hispanic share cutoffs of 10, 30 and 50%. Regressions are unweighted OLS on the ZIP-level effects with HC1 standard errors; regressors in percent are scaled to 10-point units, income to $10,000 and density to 1,000 people per square mile. Column 4 of the regression table is a ZIP-by-week panel with ZIP and county-by-week fixed effects and standard errors clustered by ZIP (1,819 ZIPs in 447 counties with at least two ZIPs); the workbench's "same county" option uses the equivalent ZIP-level regression on county-demeaned variables, where the coefficients match and robust standard errors are 0.09 and 0.06. The Pepsi figures use the same effect for Pepsi's share. Group means in Figure 2 carry intervals of 1.96 times the standard deviation across ZIPs divided by the square root of the group size. Figure 3 clips the vertical axis to −12 to +8 points; clipped ZIPs are drawn at the edge, and the fitted lines use all ZIPs.

How to cite

@misc{singh2026coke,
  author = {Singh, Vishal},
  title  = {Coca-Cola, a rumor and a different constituency},
  year   = {2026},
  note   = {Teaching case, NYU Stern School of Business},
  url    = {https://vishalsingh.org}
}

References

  1. Angrist, J. D. and Pischke, J.-S. (2009). Mostly Harmless Econometrics: An Empiricist's Companion. Princeton University Press.
  2. Bertrand, M., Duflo, E. and Mullainathan, S. (2004). How much should we trust differences-in-differences estimates? Quarterly Journal of Economics 119(1), 249–275.
  3. Roth, J., Sant'Anna, P. H. C., Bilinski, A. and Poe, J. (2023). What's trending in difference-in-differences? A synthesis of the recent econometrics literature. Journal of Econometrics 235(2), 2218–2244.