Andrew Cooke | Contents | Latest | RSS | Previous | Next


Welcome to my blog, which was once a mailing list of the same name and is still generated by mail. Please reply via the "comment" links.

Always interested in offers/projects/new ideas. Eclectic experience in fields like: numerical computing; Python web; Java enterprise; functional languages; GPGPU; SQL databases; etc. Based in Santiago, Chile; telecommute worldwide. CV; email.

Personal Projects

Choochoo Training Diary

Last 100 entries

[Programming] React Leaflet; AliExpress Independent Sellers; Applebaum - Twilight of Democracy; [Politics] Back + US Elections; [Programming,Exercise] Simple Timer Script; [News] 2019: The year revolt went global; [Politics] The world's most-surveilled cities; [Bike] Hope Freehub; [Restaurant] Mama Chau's (Chinese, Providencia); [Politics] Brexit Podcast; [Diary] Pneumonia; [Politics] Britain's Reichstag Fire moment; install cairo; [Programming] GCC Sanitizer Flags; [GPU, Programming] Per-Thread Program Counters; My Bike Accident - Looking Back One Year; [Python] Geographic heights are incredibly easy!; [Cooking] Cookie Recipe; Efficient, Simple, Directed Maximisation of Noisy Function; And for argparse; Bash Completion in Python; [Computing] Configuring Github Jekyll Locally; [Maths, Link] The Napkin Project; You can Masquerade in Firewalld; [Bike] Servicing Budget (Spring) Forks; [Crypto] CIA Internet Comms Failure; [Python] Cute Rate Limiting API; [Causality] Judea Pearl Lecture; [Security, Computing] Chinese Hardware Hack Of Supermicro Boards; SQLAlchemy Joined Table Inheritance and Delete Cascade; [Translation] The Club; [Computing] Super Potato Bruh; [Computing] Extending Jupyter; Further HRM Details; [Computing, Bike] Activities in ch2; [Books, Link] Modern Japanese Lit; What ended up there; [Link, Book] Logic Book; Update - Garmin Express / Connect; Garmin Forerunner 35 v 230; [Link, Politics, Internet] Government Trolls; [Link, Politics] Why identity politics benefits the right more than the left; SSH Forwarding; A Specification For Repeating Events; A Fight for the Soul of Science; [Science, Book, Link] Lost In Math; OpenSuse Leap 15 Network Fixes; Update; [Book] Galileo's Middle Finger; [Bike] Chinese Carbon Rims; [Bike] Servicing Shimano XT Front Hub HB-M8010; [Bike] Aliexpress Cycling Tops; [Computing] Change to ssh handling of multiple identities?; [Bike] Endura Hummvee Lite II; [Computing] Marble Based Logic; [Link, Politics] Sanity Check For Nuclear Launch; [Link, Science] Entropy and Life; [Link, Bike] Cheap Cycling Jerseys; [Link, Music] Music To Steal 2017; [Link, Future] Simulated Brain Drives Robot; [Link, Computing] Learned Index Structures; Solo Air Equalization; Update: Higher Pressures; Psychology; [Bike] Exercise And Fuel; Continental Race King 2.2; Removing Lowers; Mnesiacs; [Maths, Link] Dividing By Zero; [Book, Review] Ray Monk - Ludwig Wittgenstein: The Duty Of Genius; [Link, Bike, Computing] Evolving Lacing Patterns; [Jam] Strawberry and Orange Jam; [Chile, Privacy] Biometric Check During Mail Delivery; [Link, Chile, Spanish] Article on the Chilean Drought; [Bike] Extended Gear Ratios, Shimano XT M8000 (24/36 Chainring); [Link, Politics, USA] The Future Of American Democracy; Mass Hysteria; [Review, Books, Links] Kazuo Ishiguro - Never Let Me Go; [Link, Books] David Mitchell's Favourite Japanese Fiction; [Link, Bike] Rear Suspension Geometry; [Link, Cycling, Art] Strava Artwork; [Link, Computing] Useful gcc flags; [Link] Voynich Manuscript Decoded; [Bike] Notes on Servicing Suspension Forks; [Links, Computing] Snap, Flatpack, Appimage; [Link, Computing] Oracle is leaving Java (to die); [Link, Politics] Cubans + Ultrasonics; [Book, Link] Laurent Binet; VirtualBox; [Book, Link] No One's Ways; [Link] The Biggest Problem For Cyclists Is Bad Driving; [Computing] Doxygen, Sphinx, Breathe; [Admin] Brokw Recent Permalinks; [Bike, Chile] Buying Bearings in Santiago; [Computing, Opensuse] Upgrading to 42.3; [Link, Physics] First Support for a Physics Theory of Life; [Link, Bike] Peruvian Frame Maker; [Link] Awesome Game Theory Tit-For-Tat Thing; [Food, Review] La Fabbrica - Good Italian Food In Santiago; [Link, Programming] MySQL UTF8 Broken; [Link, Books] Latin American Authors

© 2006-2017 Andrew Cooke (site) / post authors (content).

Uniform Fences Don't Give Uniform Data

From: andrew cooke <andrew@...>

Date: Thu, 29 Mar 2012 10:10:39 -0300

from random import random
from itertools import chain
from collections import Counter

def split(k):
    fences = list(sorted(chain([0], [random() for i in range(k-1)], [1])))
    return [j - i for i,j in zip(fences, fences[1:])]

def bin(n, data):
    return Counter(int(x * n) for x in chain(*data))

def plot(bins, n, w):
    m = max(bins.values())
    for b in sorted(bins):
        print '%4.2f %s' % (b / float(n), '*' * int(bins[b] * w / m))

k = 4
n = 50
print split(k)
b = bin(n, (split(k) for i in range(n*1000/k)))
print b
plot(b, n, 70)

0.00 **********************************************************************
0.02 ******************************************************************
0.04 ******************************************************************
0.06 ***********************************************************
0.08 ************************************************************
0.10 ********************************************************
0.12 ****************************************************
0.14 *************************************************
0.16 ************************************************
0.18 *********************************************
0.20 ********************************************
0.22 *******************************************
0.24 **************************************
0.26 **************************************
0.28 **********************************
0.30 *********************************
0.32 ******************************
0.34 *****************************
0.36 ***************************
0.38 ***************************
0.40 *************************
0.42 *********************
0.44 *********************
0.46 *******************
0.48 *****************
0.50 *****************
0.52 ***************
0.54 **************
0.56 *************
0.58 ***********
0.60 **********
0.62 *********
0.64 ********
0.66 ********
0.68 *******
0.70 ******
0.72 *****
0.74 ****
0.76 ***
0.78 ***
0.80 **
0.82 **
0.84 *


N-1 Dimensional Planes

From: andrew cooke <andrew@...>

Date: Thu, 29 Mar 2012 11:32:13 -0300

Hmmm.  So the post above came from misunderstanding what someone meant by

"Choose N numbers, X[i], that sum to S (sum(X) = S), that are uniformly
distributed" - which doesn't mean that each is indiviudally uniform, but that
the X[i], when repeatedly selected, uniformly fill an N-1 dimensional plane.

(It's not clear to me that fences do that either, but the plot above doesn't
disprove it)


Comment on this post