Skip to content

Instantly share code, notes, and snippets.

View tupui's full-sized avatar
⚙️
Building Tansu

Pamphile Roy tupui

⚙️
Building Tansu
View GitHub Profile
@tupui
tupui / halton.py
Last active October 9, 2022 12:20
Halton Sequence in python
"""Halton low discrepancy sequence.
This snippet implements the Halton sequence following the generalization of
a sequence of *Van der Corput* in n-dimensions.
---------------------------
MIT License
Copyright (c) 2017 Pamphile Tupui ROY
@tupui
tupui / discerpancy.py
Created October 27, 2017 10:18
Discrepancy of a sample
"""Discrepancy of a sample.
Compute the centered discrepancy on a given sample.
It is a measure of the uniformity of the points in the parameter space.
The lower the value is, the better the coverage of the parameter space is.
---------------------------
MIT License
Copyright (c) 2017 Pamphile Tupui ROY
Permission is hereby granted, free of charge, to any person obtaining a copy
@tupui
tupui / crypto_wallet_status.py
Last active January 16, 2018 15:55
Gather information about your crypto wallets and the market
"""Get info about Cryptocurrency wallets.
Gather information about the market and you own wallet. Using transactions
file, it also compare your investment with the available capital.
Transactions are to be stored in JSON files formated as:
{
"2018_01_15": {
"rate": 1.505,
@tupui
tupui / quantile_dotplot.py
Last active May 12, 2021 13:49
Quantile dotplot in python
"""Quantile dotplot.
Based on R code from https://github.com/mjskay/when-ish-is-my-bus/blob/master/quantile-dotplots.md
Reference:
Matthew Kay, Tara Kola, Jessica Hullman, Sean Munson. When (ish) is My Bus?
User-centered Visualizations of Uncertainty in Everyday, Mobile Predictive Systems.
CHI 2016. DOI: 10.1145/2858036.2858558
@tupui
tupui / pod.py
Created June 20, 2018 12:25
Playing with Proper Orthogonal Decomposition in python using numpy
"""Proper Orthogonal Decomposition.
Demonstrate how to use it.
---------------------------
MIT License
Copyright (c) 2018 Pamphile Tupui ROY
@tupui
tupui / sobol_saltelli.py
Last active January 24, 2024 13:55
Sobol' variance based sensitivity indices based on Saltelli2010 in python
"""Sobol' indices.
Compute Sobol' variance based sensitivity indices.
Use formulations from Saltelli2010.
Reference:
Saltelli et al. Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index,
Computer Physics Communications, 2010. DOI: 10.1016/j.cpc.2009.09.018
@tupui
tupui / minimum_spanning_tree.py
Created July 28, 2018 00:11
Minimum Spanning Tree applied to Design of Experiments
"""Minimum Spanning Tree.
MST is used here as a discrepancy criterion.
Comparing two different designs: the higher the mean, the better the design is
in terms of space filling.
---------------------------
MIT License
@tupui
tupui / cosi.py
Created September 4, 2018 11:52
Cosine transformation sensitivity indices
"""Cosine transformation sensitivity indices.
Using a Discrete Cosine Transformation (DCT), it is possible to compute
first order sensitivity indices.
Reference:
Plischke E., How to compute variance-based sensitivity indicators with your
spreadsheet software, Environmental Modelling & Software,
2012. DOI: 10.1016/j.envsoft.2012.03.004
@tupui
tupui / sobol_visually.py
Created December 7, 2018 16:03
Visual explanation of first order Sobol' indices
r"""Visual explanation of Sobol' indices.
Sobol' indices are metrics to express sensitivity of the output from
perturbations comming from input parameters. First order indices write
.. math:: S_{x_i} = \frac{\mathbb{V}_i(Y)}{\mathbb{V}[Y]} =
\frac{\mathbb{\mathbb{V}}[\mathbb{E}(Y|x_i)]}{\mathbb{V}[Y]}
The following is using the Ishigami function
@tupui
tupui / moment_independent_visually.py
Created December 7, 2018 16:11
Visual explanation of Moment independent sensitivity analysis
r"""Visual explanation of Moment independent sensitivity analysis.
Moment-based method are based on the whole PDF to mitigate these
issues (Borgonovo2007). Based on the unconditional PDF, a conditional PDF per
parameter is computed. The more the conditional PDF deviates from the
unconditional PDF, the more the parameter has an impact on the quantity of
interest. The same procedure can be done using the Empirical Cumulative
Density Function (ECDF), respectively with the unconditional ECDF.
This visually shows this procedure. Bins of samples (red circles) are used to
compute a conditional PDF of the output. This PDF is compared to the