| Andrew Cooke | Contents | Latest | RSS | Twitter | 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

Lepl parser for Python.

Colorless Green.

Photography around Santiago.

SVG experiment.

Professional Portfolio

Calibration of seismometers.

Data access via web services.

Cache rewrite.

Extending OpenSSH.

C-ORM: docs, API.

Last 100 entries

Small Success With Go!; Re: Quick message - This link is broken; Adding Reverb To The Echo Chamber; Sox Audio Tools; Would This Have Been OK?; Honesty only important economically before institutions develop; Stegangraphy via PS4; OpenCL Mess; More Book Recommendations; Good Explanation of Difference Between Majority + Minority; Musical Chairs - Who's The Privileged White Guy; I can see straight men watching this conversation and laffing; When it's Actually a Source of Indignation and Disgust; Meta Thread Defending POC Causes POC To Close Account; Indigenous People Of Chile; Curry Recipe; Interesting Link On Marginality; A Nuclear Launch Ordered, 1962; More Book Recs (Better Person); It's Nuanced, And I Tried, So Back Off; Marx; The Negative Of Positive; Jenny Holzer Rocks; Huge Article on Cultural Evolution and More; "Ignoring language theory"; Negative Finger Counting; Week 12; Communication Via Telecomm Bids; Finding Suspects Via Relatives' DNA From Non-Crime Databases; Statistics and Information Theory; Ice OK in USA; On The Other Hand; (Current Understanding Of) Chilean Taxes / Contributions; M John Harrison; Playing Games on a Cloud GPU; China Gamifies Real Life; Can't Help Thinking It's Thoughtcrime; Mefi Quotes; Spray Painting Bike Frame; Weeks 10 + 11; Change: No Longer Possible To Merge Metadata; Books on Old Age; Health Tree Maps; MRA - Men's Rights Activists; Writing Good C++14; Risk Assessment - Fukushima; The Future of Advertising and Surveillance; Travelling With Betaferon; I think I know what I dislike so much about Metafilter; Weeks 8 + 9; More; Pastamore - Bad Italian in Vitacura; History Books; Iraq + The (UK) Governing Elite; Answering Some Hard Questions; Pinochet: The Dictator's Shadow; An Outsider's Guide To Julia Packages; Nobody gives a shit; Lepton Decay Irregularity; An Easier Way; Julia's BinDeps (aka How To Install Cairo); Good Example Of Good Police Work (And Anonymity Being Hard); Best Santiago Burgers; Also; Michael Emmerich (Vibrator Translator) Interview (Japanese Books); Clarice Lispector (Brazillian Writer); Books On Evolution; Looks like Ara (Modular Phone) is dead; Index - Translations From Chile; More Emotion in Chilean Wines; Week 7; Aeon Magazine (Science-ish); QM, Deutsch, Constructor Theory; Interesting Talk Transcripts; Interesting Suggestion Of Election Fraud; "Hard" Books; Articles or Papers on depolarizing the US; Textbook for "QM as complex probabilities"; SFO Get Libor Trader (14 years); Why Are There Still So Many Jobs?; Navier Stokes Incomplete; More on Benford; FBI Claimed Vandalism; Architectural Tessellation; Also: Go, Blake's 7; Delusions of Gender (book); Crypto AG DID work with NSA / GCHQ; UNUMS (Universal Number Format); MOOCs (Massive Open Online Courses); Interesting Looking Game; Euler's Theorem for Polynomials; Weeks 3-6; Reddit Comment; Differential Cryptanalysis For Dummies; Japanese Graphic Design; Books To Be Re-Read; And Today I Learned Bugs Need Clear Examples; Factoring a 67 bit prime in your head; Islamic Geometric Art; Useful Julia Backtraces from Tasks; Nothing, however, is lost with less discomfort than that which, when lost, cannot be missed

© 2006-2015 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