{{ message }}

Instantly share code, notes, and snippets.

# Mathieu Blondel mblondel

Created Aug 2, 2011
Sample variance in a single pass
View online_variance.py
 def online_mean_variance(iterable): mN = 0 mM = 0.0 mS = 0.0 for x in iterable: mN += 1 nextM = mM + (x - mM) / mN mS += (x - mM) * (x - nextM)
Last active Sep 25, 2015
Linear regression by Linear Programming
View regression_lp.py
 # (C) 2011 Mathieu Blondel # License: BSD 3 clause import numpy as np import numpy.linalg as linalg import pylab as pl from cvxopt import matrix, solvers np.random.seed(0)
Created Jul 25, 2010
Compute pi by MCMC
View mc_pi.py
 from random import random """ Find pi by the Monte-Carlo method. area of a circle = pi r^2 area of a square = (2r)^2 = 4 r^2 Perform random uniform sampling between -1 and 1. The proportion of points in the unit circle is:
Created Jul 23, 2010
Solve Number plate game by generating and interpreting Forth programs
View number_plate_solver.py
 #!/usr/bin/env python """ Find the operations needed to sum up to TARGET by using all 4 numbers in NUMBERS. """ from itertools import permutations, product NUMBERS = ["3","4","7","8"] TARGET = 10.0
Created Jul 4, 2010
Python corountine examples
View coroutines.py
 def recv_count(): try: while True: n = (yield) print "T-minus", n except GeneratorExit: print "Kaboom!" def ex1(): r = recv_count()
Created Jun 20, 2010
TCP/IP server using a main loop
View echo_server.py
 # adapted from http://roscidus.com/desktop/node/413 import socket import gobject def server(host, port): '''Initialize server and start listening.''' sock = socket.socket() sock.setsockopt(socket.SOL_SOCKET, socket.SO_REUSEADDR, 1) sock.bind((host, port))
Last active Aug 29, 2015
Local regression
View local_regression.py
 """Local regression""" # Author: Mathieu Blondel # License: BSD 3 clause import numpy as np from sklearn.base import BaseEstimator, RegressorMixin from sklearn.metrics.pairwise import pairwise_kernels from sklearn.linear_model import Ridge