Author Archive


November 22, 2015 Leave a comment



Rational ratios of angular speed…mesmerizing….

Categories: Uncategorized


October 26, 2015 Leave a comment


Categories: Uncategorized

Balloons and Bubbles

September 19, 2015 Leave a comment

I have seen a number of enjoyable presentations with animated bubble charts or plots and just wanted to play with this idea. This toy and meaningless example is a start.

The data:

ausb = WolframAlpha[
"birth rate Australia 1970 to 2010", \
{{"History:AnnualBirths:CountryData", 1}, "TimeSeriesData"}];
ausd = WolframAlpha[
"death rate Australia 1970 to 2010", \
{{"History:AnnualDeaths:CountryData", 1}, "TimeSeriesData"}];
ausp = WolframAlpha[
"population Australia 1970 to 2010", \
{{"SensiblePlot:Population:CountryData", 1}, "TimeSeriesData"}];
germb = WolframAlpha[
"birth rate Germany 1970 to 2010", \
{{"History:AnnualBirths:CountryData", 1}, "TimeSeriesData"}];
germd = WolframAlpha[
"death rate Germany 1970 to 2010", \
{{"History:AnnualDeaths:CountryData", 1}, "TimeSeriesData"}];
germp = WolframAlpha[
"population Germany 1970 to 2010", \
{{"SensiblePlot:Population:CountryData", 1}, "TimeSeriesData"}];

The plot:

all = Select[
GatherBy[ Join[ausb, ausd, ausp[[1]], germb, germd, germp[[1]]],
First], Length[#] == 6 &];
yrs = all[[All, 1, 1, 1]];
BubbleChart[List /@ #1, PlotRange -> {{0, 2000000}, {0, 1000000}},
ChartLegends -> {"Australia", "Germany"},
FrameLabel -> {"Annual births", "Annual deaths"}, PlotLabel -> #2,
BaseStyle -> 12,
Epilog -> {Dashed,
Line[{{0, 0}, {2000000,
2000000}}]}] &, {({{#1[[2, 1]], #2[[2, 1]], #3[[2,
1]]}, {#4[[2, 1]], #5[[2, 1]], #6[[2, 1]]}} & @@@ all),

The result:


I derived some amusement that the “rapid drop” in annual births below “annual deaths” for Germany was associated with a slow “deflation” of the population “balloon”/bubble.


Categories: Uncategorized

Statistics Done Wrong

August 17, 2015 Leave a comment


This is an excellent and important book.  The author covers important errors, biases and misconceptions in scientific studies and their interpretation. The writing  is clear, engaging and entertaining style. Each chapter ends with tips to prevent or address the issues raised in the chapter. The final chapter is an exhortation to all stakeholders, scientists, publishers, students and the general audience for culture change to improve the quality of our scientific debate and development.

The book highlights the prevalence of these misconceptions and as a reader I am glad to be woken from my self-satisfied but delusional slumber to think “slow” as well as fast (ala Kahneman).


Categories: books, Mathematics

More Tupperware

July 11, 2015 Leave a comment

This is post is motivated by a question on Mathematica Stackexchange and the interesting link posted in the question. My previous post had subtle errors. This allowed me to play with the higher resolution self referential formula:

g[x_, y_] :=
Boole[Mod[Floor[Floor[y/61] 2^(-61 x - Mod[y, 61])], 2] == 1]
w[nu_] := ArrayPlot[Table[g[j, k], {k, nu + 60, nu, -1}, {j, 0, 375}]]
btupf[s_] := Module[{i, m, r},
i = Rasterize[s, ImageSize -> {376, 61}];
m = Map[Boole[Max@# < 1] &, ImageData[i], {2}]; r = 61 FromDigits[Flatten[Reverse@Transpose[m]], 2]; ArrayPlot[Table[g[j, k], {k, r + 60, r, -1}, {j, 0, 375}]]]

btupn[s_] := Module[{i, m, r}, i = Rasterize[s, ImageSize -> {376, 61}];
m = Map[Boole[Max@# < 1] &, ImageData[i], {2}];
r = 61 FromDigits[Flatten[Reverse@Transpose[m]], 2]]

Here g is the function,
plots a given number, btupf allows you to put in text, convert to number and array plot, btupn gives you the number.

So w[nn] where nn is this number yields:


Categories: Mathematica

Why Some Say The Moon

June 22, 2015 Leave a comment

It has been a challenging time recently and I have been plagued by ill health…in spare time I have been musing with TimelinePlot,










I still vividly remember watching the lunar landing on a black and white television…despite the grainy image and the staccato noisy audio I was transfixed.

Categories: Mathematica


May 30, 2015 1 comment

I particularly enjoyed a Numberphile video on the”everything formula”. Well here is a version of my tupper number:

Here is my first attempt at coding (unfortunately wordpress does not correctly process my code tag, hence the image):  
tupf produces the array plot and tupn the number. The above number was produced using tupn[Style["u b p d q n",20]]…not perfect but fun.

Categories: Mathematica

Get every new post delivered to your Inbox.

Join 41 other followers