Probly is a Python-like mini-language for probabilistic estimation. It's based on Starlark and implemented in Go.
You may use any Starlark syntax. There are only the following differences to Starlark:
math
module is imported by default, so you can directly use it e.g. math.sqrt(2)
. Probly also has a built-in sum
function not available in Starlark.This page will show you the probability distributions of all the numeric (scalar or distribution) global variables in your program (except those starting with an underscore). The values are taken at the end of the program's execution.
Name | to p10 to p90 |
pm plus/minus |
td times/divide |
Quantiles | Notes | |||||
---|---|---|---|---|---|---|---|---|---|---|
Normal |
mean |
sd |
2 | |||||||
LogNormal |
mu |
sigma |
2 | Alternatively: mean , sd |
||||||
Beta |
alpha |
beta |
||||||||
PERT |
min |
mode |
max |
[lambd] |
Like the triangular, but smoother (Wikipedia) | |||||
Uniform |
a |
b |
2 | a need not be less than b |
||||||
LogUniform |
a |
b |
2 | a need not be less than b |
||||||
Bernoulli |
p |
|||||||||
Binomial |
n |
p |
||||||||
Discrete |
x_1 |
p_1 |
x_2 |
p_2 |
... | Generic discrete distribution over any finite set of values |
These mathematical functions and constants are available in the math
module:
pow(x, y)
- Returns x raised to the power of yexp(x)
sqrt(x)
log(x, [base])
- Natural logarithm by default if base is not specifiede
pi
ceil(x)
floor(x)
fabs(x)
- Returns the absolute value of x as floatcopysign(x, y)
- Returns a value with the magnitude of x and the sign of ymod(x, y)
- Returns x
modulo y
remainder(x, y)
round(x)
- Returns the nearest integer, rounding half away from zeroacos(x)
asin(x)
atan(x)
atan2(y, x)
- Returns atan(y / x). The result is between -pi and picos(x)
sin(x)
tan(x)
degrees(x)
- Converts angle x from radians to degreesradians(x)
- Converts angle x from degrees to radiansacosh(x)
asinh(x)
atanh(x)
cosh(x)
sinh(x)
tanh(x)
hypot(x, y)
- Returns the Euclidean norm, sqrt(x^2 + y^2); the distance from the origin to (x, y)gamma(x)
- Returns the Gamma function at xThis code provides an example of the syntax of Starlark:
# Define a number
number = 18
# Define a list
numbers = [1, 2, 3, 4, 5]
# List comprehension
halves = [n / 2 for n in numbers]
# Define a function
def is_even(n):
"""Return True if n is even."""
return n % 2 == 0
# Define a dictionary
people = {
"Alice": 22,
"Bob": 40,
"Charlie": 55,
"Dave": 14,
}
names = ", ".join(people.keys()) # Alice, Bob, Charlie, Dave
# Modify a variable in a loop
sum_even_ages = 0
for age in people.values():
if is_even(age):
sum_even_ages += age
# Append to a list in a loop
over_30_names = []
for name, age in people.items():
if age > 30:
over_30_names.append(name)
If you've ever used Python, this should look very familiar. In fact, the code above is also valid Python code. Still, this short example shows most of the language. Starlark is a very small language that implements a limited subset of Python.
For our purposes, one notable difference to Python is that the exponentiation operator **
is not supported. You have to use math.pow
.
You can also look at the Starlark language specification.
Though not designed for speed, Probly is fast enough for practical purposes: around 10 milliseconds for 3,000 samples, for most examples on this page. This is due to being implemented in Go.
The time taken to return results on this page is spent overwhelmingly in web application code, not in Probly evaluation.
Interestingly, Probly is still slower than Python code that uses entirely numpy
array operations, which are very well optimised.
This should only begin to matter at very large scales, or if latency is critical.
It's not currently possible to obtain and manipulate properties of a distribution within an Probly program, like so:
x = Normal(1 to 10)
y = x.std() # Not possible
Supporting this would require some fundamental changes to the implementation of Probly, which is currently very simplistic.
The to
binary operator was inspired by Squiggle.
This example follows Section 3 of Dissolving the Fermi Paradox (Sandberg and colleagues, 2018).
From the abstract:
[The Fermi paradox] arises from the use of Drake-like equations, which implicitly assume certainty regarding highly uncertain parameters. ... When the model is recast to represent realistic distributions of uncertainty, we find a substantial ex ante probability of there being no other intelligent life in our observable universe, and thus that there should be little surprise when we fail to detect any signs of it.
The results of our simulation should approximately match Figure 2 of the paper:
The mean for N is [...] 27 million, but the median is now only 0.32 – less than one civilization per galaxy like our own. The probability of N<1 is now 52%.
Mean | 25.9 E 6 |
Std. dev. | 409 E 6 |
Variance | 167 E 15 |
Quantile | |
---|---|
0.05 | 0 |
0.25 | 4.73 E -15 |
0.50 | 1.35 |
0.75 | 8 900 |
0.95 | 15.1 E 6 |
Mean | 602 E 6 |
Std. dev. | 1.66 E 9 |
Variance | 2.74 E 18 |
Quantile | |
---|---|
0.05 | 268 |
0.25 | 10 700 |
0.50 | 1.32 E 6 |
0.75 | 128 E 6 |
0.95 | 4.40 E 9 |
Mean | 0.220 |
Std. dev. | 0.255 |
Variance | 0.065 2 |
Quantile | |
---|---|
0.05 | 0.012 6 |
0.25 | 0.031 3 |
0.50 | 0.101 |
0.75 | 0.327 |
0.95 | 0.812 |
Mean | 0.143 |
Std. dev. | 0.227 |
Variance | 0.051 7 |
Quantile | |
---|---|
0.05 | 0.001 42 |
0.25 | 0.005 73 |
0.50 | 0.030 0 |
0.75 | 0.175 |
0.95 | 0.710 |
Mean | 0.507 |
Std. dev. | 0.495 |
Variance | 0.245 |
Quantile | |
---|---|
0.05 | 0 |
0.25 | 222 E -18 |
0.50 | 0.769 |
0.75 | 1.00 |
0.95 | 1.00 |
Mean | 270 E 75 |
Std. dev. | 14.7 E 78 |
Variance | 217 E 156 |
Quantile | |
---|---|
0.05 | 9.17 E -39 |
0.25 | 237 E -18 |
0.50 | 1.47 |
0.75 | 214 E 12 |
0.95 | 5.23 E 36 |
Mean | 0.392 |
Std. dev. | 0.249 |
Variance | 0.062 1 |
Quantile | |
---|---|
0.05 | 0.113 |
0.25 | 0.180 |
0.50 | 0.314 |
0.75 | 0.566 |
0.95 | 0.893 |
Mean | 0.395 |
Std. dev. | 0.249 |
Variance | 0.062 2 |
Quantile | |
---|---|
0.05 | 0.114 |
0.25 | 0.183 |
0.50 | 0.321 |
0.75 | 0.561 |
0.95 | 0.899 |
Mean | 21.7 |
Std. dev. | 25.7 |
Variance | 659 |
Quantile | |
---|---|
0.05 | 1.26 |
0.25 | 3.11 |
0.50 | 9.58 |
0.75 | 32.6 |
0.95 | 81.9 |
star_formation | planet | habitable | r_life | ... | intelligence | detectable | longevity | n | |
---|---|---|---|---|---|---|---|---|---|
0 | 3.10 | 0.995 | 0.165 | 204 E -27 | ... | 0.008 31 | 0.261 | 104 | 0 |
1 | 1.04 | 0.138 | 0.497 | 240 E 9 | ... | 0.022 8 | 0.068 0 | 154 E 6 | 17 000 |
2 | 34.5 | 0.916 | 0.100 | 7.85 | ... | 0.386 | 0.375 | 60.5 E 6 | 27.7 E 6 |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
2997 | 1.19 | 0.212 | 0.326 | 21.9 E -21 | ... | 0.019 2 | 0.204 | 8 480 | 0 |
2998 | 15.3 | 0.719 | 0.291 | 586 E -39 | ... | 0.343 | 0.365 | 54.9 E 6 | 0 |
2999 | 35.5 | 0.661 | 0.108 | 841 E -21 | ... | 0.001 08 | 0.013 9 | 272 000 | 0 |
Get the simulation data (and more) in a machine-readable format:
/api/sim/23WJtteuWphUM6kgERUJr9/