r = 0.98 between seats and seats_now

Measured on saas-subscriptions: the Pearson correlation between seats and seats_now over 2,500 rows is 0.98.

That is strong and positive. Squared, it says the linear relationship accounts for 96.9% of the variance in either column. That is the number worth quoting rather than r itself: 0.98 sounds like a much stronger claim than 96.9% does.

df[["seats", "seats_now"]].corr(method="pearson")
# also worth running:
df[["seats", "seats_now"]].corr(method="spearman")

Run Spearman next to Pearson every time. They agree when the relationship is roughly linear and diverge when it is monotone but curved — and the gap between them is a free diagnostic that costs one line.

A correlation this size will get read as a mechanism. It is not one. Both columns here are generated from a shared underlying process, which is exactly the situation a confounder creates in real data: two things that move together because a third thing moves both.

Look at the scatter before quoting r. Anscombe's quartet is four datasets with identical correlation coefficients and nothing else in common, and the correlation explorer will draw this pair so you can see which case you are in.

Longer treatment in Correlation, confounding, and Simpson's paradox.