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

C[omp]ute

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

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; Article on Didion; Cost of Living by City; British Slavery; Derrida on Metaphor; African SciFi; Traits in Julia; Alternative Japanese Lit; Pulic Key as Address (Snow); Why Information Grows; The Blindness Of The Chilean Elite; Some Victoriagate Links; This Is Why I Left StackOverflow; New TLS Implementation; Maths for Physicists; How I Am 8; 1000 Word Philosophy; Cyberpunk Reading List; Detailed Discussion of Message Dispatch in ParserCombinator Library for Julia; FizzBuzz in Julia w Dependent Types; kokko - Design Shop in Osaka; Summary of Greece, Currently; LLVM and GPUs; See Also; Schoolgirl Groyps (Maths); Japanese Lit; Another Example - Modular Arithmetic; Music from United; Python 2 and 3 compatible alternative.; Read Agatha Christie for the Plot; A Constructive Look at TempleOS; Music Thread w Many Recommendations; Fixed Version; A Useful Julia Macro To Define Equality And Hash; k3b cdrom access, OpenSuse 13.1; Week 2; From outside, the UK looks less than stellar; Huge Fonts in VirtualBox; Keen - Complex Emergencies; The Fallen of World War II; Some Spanish Fiction; Calling C From Fortran 95; Bjork DJ Set; Z3 Example With Python; Week 1; Useful Guide To Starting With IJulia; UK Election + Media; Review: Reinventing Organizations; Inline Assembly With Julia / LLVM; Against the definition of types; Dumb Crypto Paper; The Search For Quasi-Periodicity...; Is There An Alternative To Processing?; CARDIAC (CARDboard Illustrative Aid to Computation); The Bolivian Case Against Chile At The Hague; Clear, Cogent Economic Arguments For Immigration; A Program To Say If I Am Working; Decent Cards For Ill People; New Photo; Luksic And Barrick Gold; President Bachelet's Speech; Baltimore Primer; libxml2 Parsing Stream; configure.ac Recipe For Library Path; The Davalos Affair For Idiots; Not The Onion: Google Fireside Chat w Kissinger; Bicycle Wheels, Inertia, and Energy; Another Tax Fraud; Google's Borg; A Verion That Redirects To Local HTTP Server; Spanish Accents For Idiots; Aluminium Cans; Advice on Spray Painting; Female View of Online Chat From a Male; UX Reading List; S4 Subgroups - Geometric Interpretation; Fucking Email; The SQM Affair For Idiots; Using Kolmogorov Complexity; Oblique Strategies in bash; Curses Tools; Markov Chain Monte Carlo Without all the Bullshit; Email Para Matias Godoy Mercado; The Penta Affair For Idiots; Example Code To Create numpy Array in C; Good Article on Bias in Graphic Design (NYTimes); Do You Backup github?

© 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 *
0.86 
0.88 
0.90 
0.92 
0.94 
0.96 

Andrew

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
"uniform".

"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)

Andrew

Comment on this post