Skip to content

Instantly share code, notes, and snippets.

View Tillsten's full-sized avatar

Till Stensitzki Tillsten

11 followers · 10 following
  • Universität Potsdam
  • Berlin, Germany
View GitHub Profile
@Tillsten
Tillsten / exp_sum_pymc.py
Created October 29, 2012 02:56
Fit to an expontial sum in pymc
# Make test data
t=np.linspace(0,20,200)
decays = np.array([2.3,9])
x=np.linspace(-1,1, 32)
amp1 = 10*np.exp(-(x+0.5)**2/0.3)
amp2 = -8*np.exp(-(x-0.5)**2/0.4)
amps=np.column_stack((amp1,amp2))
d=amps.dot(np.exp(-t/decays[:,None]))
#Add Noise
@Tillsten
Tillsten / stacked_steamgraph_mpl.py
Created November 9, 2012 19:54
Steamgraph in mpl
import numpy as np
import random as r
import matplotlib.pyplot as plt
def layers(n, m):
'function to generate some data'
def bump(a):
x = 1 / (.1 + r.random())
y = 2 * r.random() - .5
z = 10 / (.1+ r.random())
for i in range(m):
@Tillsten
Tillsten / compartment_model.py
Last active December 12, 2015 08:39
Solving compartment models analytical and build a cython function from the solution.
# -*- coding: utf-8 -*-
"""
This modul build a kinetic matrix from given transtions between
compartments.
"""
import sympy, cython
import numpy as np
class Transition(object):
@Tillsten
Tillsten / gist:5057677
Last active December 14, 2015 08:28
{
"metadata": {
"name": "SVD on data with a peak-shift"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
@Tillsten
Tillsten / advanced version.py
Last active August 29, 2015 14:07
most basic proof of concept for using cassowary with mpl.
# -*- coding: utf-8 -*-
"""
Created on Tue Sep 30 14:20:54 2014
@author: tillsten
"""
import cassowary # can be installed via pip, pure python.
import matplotlib.pyplot as plt
import matplotlib.textpath
@Tillsten
Tillsten / emcee.ipynb
Last active August 29, 2015 14:08
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@Tillsten
Tillsten / varpro.ipynb
Last active August 29, 2015 14:10
varpro
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@Tillsten
Tillsten / gist:f109497fb605990a4b42
Created February 17, 2015 22:31
<!DOCTYPE html>
<html>
<head>
<meta charset="utf-8" />
<title>mpl_anamoaly</title>
<script src="https://cdnjs.cloudflare.com/ajax/libs/require.js/2.1.10/require.min.js"></script>
<script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/2.0.3/jquery.min.js"></script>
@Tillsten
Tillsten / mpl_inline
Created February 17, 2015 22:43
mpl inline
{
"cells": [
{
"cell_type": "code",
"execution_count": 28,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
@Tillsten
Tillsten / fixing_sigma.ipynb
Created September 4, 2015 12:55
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.