2024 π Daylatest newsbuy art
Trance opera—Spente le Stellebe dramaticmore quotes
very clickable
visualization + math

π Day 2024 Art Posters - A community garden of digits of π
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca buy artwork
2024 π DAY | 768 digits of `\pi` as a garden at night. Explore the gardens (BUY ARTWORK)

`\pi` Day 2021 Art Posters - A forest of `\pi` (a Lindenmayer system)


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2021 π DAY | Good things grow for those who wait.' edition.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2019 π DAY | Hundreds of digits, hundreds of languages and a special kids' edition.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2018 π DAY | Street maps to new destinations.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2017 π DAY | Imagine the sky in a new way.


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2016 π APPROXIMATION DAY | What would happen if about right was right.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2016 π DAY | These digits really fall for each other.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2015 π DAY | A transcendental experience.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2014 π APPROXIMATION DAY | Spirals into roughness.


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2014 π DAY | Hypnotizes you into looking.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2014 π DAY | Come into the fold.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2013 π DAY | Where it started.

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
CIRCULAR π ART | And other distractions.

On March 14th celebrate `\pi` Day. Hug `\pi`—find a way to do it.

For those who favour `\tau=2\pi` will have to postpone celebrations until July 26th. That's what you get for thinking that `\pi` is wrong. I sympathize with this position and have `\tau` day art too!

If you're not into details, you may opt to party on July 22nd, which is `\pi` approximation day (`\pi` ≈ 22/7). It's 20% more accurate that the official `\pi` day!

Finally, if you believe that `\pi = 3`, you should read why `\pi` is not equal to 3.

Most of the art is available for purchase as framed prints and, yes, even pillows. Sleep's never been more important — I take custom requests.

The trees along this city street,
Save for the traffic and the trains,
Would make a sound as thin and sweet
As trees in country lanes.
—Edna St. Vincent Millay (City Trees)


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca buy artwork
768 digits of `\pi` as a forest of trees. Underwater technicolor edition. (BUY ARTWORK)

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca buy artwork
768 digits of `\pi` as a forest of trees. Scorched technicolor edition. (BUY ARTWORK)

`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca buy artwork
768 digits of `\pi` as a forest of leafless trees. Desolate edition. (BUY ARTWORK)

Welcome to this year's celebration of `\pi` and mathematics.

The theme this year is flower and flowers—in contrast to last year's understandable downturn in mood.

This year's `\pi` poem City Trees by Edna St. Vincent Millay.

This year's `\pi` day song is Sway by Laleh.

The 2015 π Day art takes a Mondrain perspective on π. The art was used in a collaboration with Max Cooper for his track Transcendental Tree Map from the album Yearning for the Infinite. Animation by Nick Cobby and myself. Watch the full show at the Barbican Centre.

The forest was generated using an Lindenmayer system — this part of the process is relatively easy.

Integrating variation into how the trees look while keeping a clean look-and-feel took a lot of experimentation — there are many free parameters in this section.

Lindenmayer system

A Lindenmayer system is a kind of parallel rewriting system. The system starts with some initial state, typically represented by a string. For example FX.

On each iteration, rules are applied to each characters in the state, which replaces the character with some other character. For example, suppose we replace X with [-FX][+FX]. Starting with the axiom FX The first three rounds of rewriting gives each give us

FX
F[-FX][+FX]
F[-F[-FX][+FX]][+F[-FX][+FX]]
F[-F[-F[-FX][+FX]][+F[-FX][+FX]]][+F[-F[-FX][+FX]][+F[-FX][+FX]]]

It doesn't look like much, I know, but the magic happens when you interpret these characters in the context of turtle graphics. Specifically, let

F move forward distance d
- turn left &theta degrees;
+ turn right &theta degrees;
[ spawn a new turtle at this position facing at current angle
] restore previous turtle

If we set `\theta = 15\deg` then the shapes for each iteration are


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
The first three rounds of rewriting of a simple L-system with axiom FX and rewrite rule X = [-FX][+FX].

We can interpret F as a branch and X a point were new branches grow. You can already see that by level 3 two of the branches — there's nothing in the system that prevents this, since each turtle operates independently.

Rules can be as simple or as complicated as you want and new characters to the string can be added to correspond to different turtle actions.

The `\pi` forest Lindenmayer system

The rules for the `\pi` forest are quite simple. The system starts with the axiom fX (grow a trunk and then branch off) and each digit corresponds to a different branching rule — the digit gives the number of branches with zero terminating a branch.

X 0 []
X 1 [---FX][++h]
X 2 [----FX][+++FX]
X 3 [-----FX][-FX][++++FX]
X 4 [------FX][---FX][++FX][+++++FX]
X 5 [-------FX][----FX][-FX][+++FX][+++++++FX]
X 6 [--------FX][------FX][---FX][+FX][+++++FX][++++++++FX]
X 7 [--------FX][-----FX][--FX][FX][+++FX][++++++FX][+++++++++FX]
X 8 [---------FX][------FX][---FX][-FX][++FX][++++FX][+++++++FX][+++++++++FX]

where each – or + correspond to a left or right turn of 8 degrees. Any turns

The digit 9 has a special rule that sprouts a flower at the previous branching point and starts a new tree.

Also, at least step branches grow

F FF
f ff
g gg

with the f step being a fraction `(1-1/\pi)^2 \approx 46%` of the length of F and g being `(1-1/\pi)^4 \approx 22%` of the length of F. These shorter branch allow me to have a shorter trunk (axiom is fX and not FX) and a very short branch used in the rule for the digit 1, which branches off to one side while growing a little stubby branch to the other side.

There is one additional rule in which the turtle is mirrored at a branching point when a digit and the next digit have different parity (one odd and one even).

X = MX

The mirror operator makes the current turtle interpret any left turns as right turns, and vice versa.

The reason for the multiple turns in the rules is for convenience &mdsah; it allows me to define the angles of the branches as a multiple of a fixed value.

The system is grown for 6 iterations.

The first tree

Here is the first tree of the forest – for digits 314159. It ends in a 9 because the rule for 9 is to grow a flower, terminate the current tree and start a new tree with the next set of digits (up to the next 9, and so on).


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
The first tree in the `\pi` forest. Its digits are 314159 -- they end at the first 9 in `\pi`.

You can trace the digits for the tree left-to-right along the leaves. The flower of the 9 will appear on the branch of the digit immediately before it — in this case, a 5. The color of the flower petal is green (5) and the center is dark orange (first digit of the next tree).

The L-system string for this tree is below. See if you can discover what symbol is used for the flower.

ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff
[-----FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
 [---FFFFFFFFFFFFFFFF
  [-------FFFFFFFFO][----FFFFFFFF][-FFFFFFFF][+++FFFFFFFF][+++++++FFFFFFFF]
 ]
 [++gggggggggggggggg]
]
[-FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFM
 [------FFFFFFFFFFFFFFFF][---FFFFFFFFFFFFFFFF][++FFFFFFFFFFFFFFFF][+++++FFFFFFFFFFFFFFFF]
]
[++++FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
 [---FFFFFFFFFFFFFFFF][++gggggggggggggggg]
]

adding variation

To add variation to the tree, I take the tree's digits and create linear congruential generators for each of the following: branch turn angle, branch growth rate, thickness tapering and turtle angle when growing a branch.

Each generator is constructed as follows. The digits of the tree `d_1d_2d_3...d_n` are split into `a = d_1d_2d_3...d_{n-1}`, `c = d_n` and `m = d_1d_2d_3...d_n`. The generator is seeded with `X_0 = d_1` and then the next digit is `X_1 = (a X_0 + c) \; \textrm{mod} \; m`. Each output `X_i` is transformed to `y \times (2X_i/m - 1)`. Because some of the trees have a large number of digits, the operations are performed with an arbitrary precision library.


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
The output of the linear congruential generator (LCG) for the first tree. The first 1,000 samples are shown below a lag plot of two successive samples. As with all LCGs, you get a banding pattern.

The output will be distributed in the range `[-y,y)`. The distribution will depend on the values of `a`, `c` and `m`. It might be reasonably uniform for some combinations and absolutely non-uniform for others. But here, non-uniform doesn't mean horrible — just another rule in the digit forest.

Each type of variation starts with an identical generator.


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
The effect of successively adding different types of variation. Shown are the first 5 trees in the forest with no leaves or flowers. A row contains all the variation in the previous row plus the type noted on the row.

The parameters that influence the amount of variation are irrelevant — I tweaked them until I had something that looked "pretty good".

winds of change are blowing

You might have noticed that because a branch growth step is doubled at each iteration, a branch is composed of 64 turtle steps. This was done so that there's lots of opportunities to sample variation along a branch. And also for boring technical reasons: SVG does not allow a path with varying thickness.

This means that variation can build up along a branch, especially if we allow the turtle to veer quite a bit as the branch grows. Remember, a branch begins to grows in a particular direction but can change course along its growth.

Below I show what happens when I crank up every variation type by 200% and by 300% (e.g. 300% moar branch growth rate jitter!). Things add up pretty quickly — I'm looking at you 79! What are you doing? Who knew the digit forest could be so complicated?


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
The effect of increasing variation of all types to 200% and 300%.

once more with ecosystem flavour

Exploring the space of variation is fun — even at 1am as I'm trying to finish this.

Below are some scenarios to get the imagination going.

The underwater scenario uses an effect in which turtle angle during branch growth is increased and constrained to be mostly up and no branch thickness change.

The prickly forest has a larger branch angle and larger branch length decay as a function of depth.

The dry forest has extreme thickness tapering and branch growth angle variation increases with each branch depth.

Finally, the bat cave is similar to the underwater case except that it has a much narrower branching angle and branch length now increases with branch depth. And one more change that's easy to spot.


`\pi` Day 2021 Art Posters - A forest of Pi (a Lindenmayer system)
 / Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Different combinations of variation can create some funky effects. Here I show the digit forest as imagined underwater, in a desert, in a drought and just hanging out in a bat cave.

So as you can see, the canonical digit forest is just one of a wide range of possible numerical environments.

news + thoughts

Propensity score matching

Mon 16-09-2024

I don’t have good luck in the match points. —Rafael Nadal, Spanish tennis player

In many experimental designs, we need to keep in mind the possibility of confounding variables, which may give rise to bias in the estimate of the treatment effect.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Nature Methods Points of Significance column: Propensity score matching. (read)

If the control and experimental groups aren't matched (or, roughly, similar enough), this bias can arise.

Sometimes this can be dealt with by randomizing, which on average can balance this effect out. When randomization is not possible, propensity score matching is an excellent strategy to match control and experimental groups.

Kurz, C.F., Krzywinski, M. & Altman, N. (2024) Points of significance: Propensity score matching. Nat. Methods 21:1770–1772.

Nasa to send our human genome discs to the Moon

Sat 23-03-2024

We'd like to say a ‘cosmic hello’: mathematics, culture, palaeontology, art and science, and ... human genomes.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
SANCTUARY PROJECT | A cosmic hello of art, science, and genomes. (details)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
SANCTUARY PROJECT | Benoit Faiveley, founder of the Sanctuary project gives the Sanctuary disc a visual check at CEA LeQ Grenoble (image: Vincent Thomas). (details)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
SANCTUARY PROJECT | Sanctuary team examines the Life disc at INRIA Paris Saclay (image: Benedict Redgrove) (details)

Comparing classifier performance with baselines

Fri 22-03-2024

All animals are equal, but some animals are more equal than others. —George Orwell

This month, we will illustrate the importance of establishing a baseline performance level.

Baselines are typically generated independently for each dataset using very simple models. Their role is to set the minimum level of acceptable performance and help with comparing relative improvements in performance of other models.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Nature Methods Points of Significance column: Comparing classifier performance with baselines. (read)

Unfortunately, baselines are often overlooked and, in the presence of a class imbalance, must be established with care.

Megahed, F.M, Chen, Y-J., Jones-Farmer, A., Rigdon, S.E., Krzywinski, M. & Altman, N. (2024) Points of significance: Comparing classifier performance with baselines. Nat. Methods 21:546–548.

Happy 2024 π Day—
sunflowers ho!

Sat 09-03-2024

Celebrate π Day (March 14th) and dig into the digit garden. Let's grow something.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
2024 π DAY | A garden of 1,000 digits of π. (details)

How Analyzing Cosmic Nothing Might Explain Everything

Thu 18-01-2024

Huge empty areas of the universe called voids could help solve the greatest mysteries in the cosmos.

My graphic accompanying How Analyzing Cosmic Nothing Might Explain Everything in the January 2024 issue of Scientific American depicts the entire Universe in a two-page spread — full of nothing.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
How Analyzing Cosmic Nothing Might Explain Everything. Text by Michael Lemonick (editor), art direction by Jen Christiansen (Senior Graphics Editor), source: SDSS

The graphic uses the latest data from SDSS 12 and is an update to my Superclusters and Voids poster.

Michael Lemonick (editor) explains on the graphic:

“Regions of relatively empty space called cosmic voids are everywhere in the universe, and scientists believe studying their size, shape and spread across the cosmos could help them understand dark matter, dark energy and other big mysteries.

To use voids in this way, astronomers must map these regions in detail—a project that is just beginning.

Shown here are voids discovered by the Sloan Digital Sky Survey (SDSS), along with a selection of 16 previously named voids. Scientists expect voids to be evenly distributed throughout space—the lack of voids in some regions on the globe simply reflects SDSS’s sky coverage.”

voids

Sofia Contarini, Alice Pisani, Nico Hamaus, Federico Marulli Lauro Moscardini & Marco Baldi (2023) Cosmological Constraints from the BOSS DR12 Void Size Function Astrophysical Journal 953:46.

Nico Hamaus, Alice Pisani, Jin-Ah Choi, Guilhem Lavaux, Benjamin D. Wandelt & Jochen Weller (2020) Journal of Cosmology and Astroparticle Physics 2020:023.

Sloan Digital Sky Survey Data Release 12

constellation figures

Alan MacRobert (Sky & Telescope), Paulina Rowicka/Martin Krzywinski (revisions & Microscopium)

stars

Hoffleit & Warren Jr. (1991) The Bright Star Catalog, 5th Revised Edition (Preliminary Version).

cosmology

H0 = 67.4 km/(Mpc·s), Ωm = 0.315, Ωv = 0.685. Planck collaboration Planck 2018 results. VI. Cosmological parameters (2018).

Error in predictor variables

Tue 02-01-2024

It is the mark of an educated mind to rest satisfied with the degree of precision that the nature of the subject admits and not to seek exactness where only an approximation is possible. —Aristotle

In regression, the predictors are (typically) assumed to have known values that are measured without error.

Practically, however, predictors are often measured with error. This has a profound (but predictable) effect on the estimates of relationships among variables – the so-called “error in variables” problem.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Nature Methods Points of Significance column: Error in predictor variables. (read)

Error in measuring the predictors is often ignored. In this column, we discuss when ignoring this error is harmless and when it can lead to large bias that can leads us to miss important effects.

Altman, N. & Krzywinski, M. (2024) Points of significance: Error in predictor variables. Nat. Methods 21:4–6.

Background reading

Altman, N. & Krzywinski, M. (2015) Points of significance: Simple linear regression. Nat. Methods 12:999–1000.

Lever, J., Krzywinski, M. & Altman, N. (2016) Points of significance: Logistic regression. Nat. Methods 13:541–542 (2016).

Das, K., Krzywinski, M. & Altman, N. (2019) Points of significance: Quantile regression. Nat. Methods 16:451–452.

Martin Krzywinski | contact | Canada's Michael Smith Genome Sciences CentreBC Cancer Research CenterBC CancerPHSA
Google whack “vicissitudinal corporealization”
{ 10.9.234.152 }