View source for Softmax function

{{About|the smooth approximation of arg max|the smooth approximation of max|LogSumExp}}
In [[mathematics]], the '''softmax function,''' also known as '''softargmax'''{{sfn|Goodfellow|Bengio|Courville|2016|p=184}} or '''normalized exponential function''',<ref name="bishop" />{{rp|198}} is a function that takes as input a vector of ''K'' real numbers, and normalizes it into a [[Probability distribution|probability distribution]] consisting of ''K'' probabilities proportional to the exponentials of the input numbers. That is, prior to applying softmax, some vector components could be negative, or greater than one; and might not sum to 1; but after applying softmax, each component will be in the [[Interval (mathematics)|interval]] <math>(0, 1)</math>, and the components will add up to 1, so that they can be interpreted as probabilities. Furthermore, the larger input components will correspond to larger probabilities. Softmax is often used in [[Artificial neural network|neural networks]], to map the non-normalized output of a network to a probability distribution over predicted output classes.

The standard (unit) softmax function <math>\sigma : \R^K\to\R^K</math>is defined by the formula

:<math>\sigma(\mathbf{z})_i = \frac{e^{z_i}}{\sum_{j=1}^K e^{z_j}} \text{ for } i = 1, \dotsc , K \text{ and } \mathbf z=(z_1,\dotsc,z_K) \in\R^K</math>

In words: we apply the standard [[exponential function]] to each element <math>z_i</math> of the input vector <math>\mathbf z</math> and normalize these values by dividing by the sum of all these exponentials; this normalization ensures that the sum of the components of the output vector <math>\sigma(\mathbf z)</math> is 1.

Instead of {{mvar|''e''}}, a different [[Base (exponentiation)|base]] {{mvar|''b''}}&nbsp;>&nbsp;0 can be used; choosing a larger value of {{mvar|''b''}} will create a probability distribution that is more concentrated around the positions of the largest input values. Writing <math>b = e^\beta</math> or <math>b = e^{-\beta}</math>{{efn|Positive {{mvar|''β''}} corresponds to the maximum convention, and is usual in machine learning, corresponding to the highest score having highest probability. The negative {{math|−''β''}} corresponds to the minimum convention, and is conventional in thermodynamics, corresponding to the lowest energy state having the highest probability; this matches the convention in the [[Gibbs distribution]], interpreting {{mvar|''β''}} as [[coldness]].}} (for real {{mvar|''β''}}){{efn|1=The notation {{mvar|''β''}} is for the [[thermodynamic beta]], which is inverse [[temperature]]: <math>\beta = 1/T</math>, <math>T = 1/\beta.</math>}} yields the expressions:{{efn|For <math>\beta = 0</math> ([[coldness]] zero, infinite temperature), <math>b = e^\beta = e^0 = 1</math>, and this becomes the constant function {{tmath|(1/n, \dots, 1/n)}}, corresponding to the [[discrete uniform distribution]].}}

:<math>\sigma(\mathbf{z})_i = \frac{e^{\beta z_i}}{\sum_{j=1}^K e^{\beta z_j}} \text{ or } \sigma(\mathbf{z})_i = \frac{e^{-\beta z_i}}{\sum_{j=1}^K e^{-\beta z_j}} \text{ for } i = 1,\dotsc , K </math>.

In some fields, the base is fixed, corresponding to a fixed scale,{{efn|In statistical mechanics, fixing {{mvar|''β''}} is interpreted as having coldness and temperature of 1.}} while in others the parameter {{mvar|''β''}} is varied.

== Interpretations ==
=== Smooth arg max ===
{{See also|Arg max}}
The name "softmax" is misleading; the function is not a [[smooth maximum]] (a [[smooth approximation]] to the [[maximum]] function), but is rather a smooth approximation to the [[arg max]] function: the function whose value is ''which'' index has the maximum. In fact, the term "softmax" is also used for the closely related [[LogSumExp]] function, which is a smooth maximum. For this reason, some prefer the more accurate term "softargmax", but the term "softmax" is conventional in machine learning.<ref name="sako2018"/>{{sfn|Goodfellow|Bengio|Courville|2016|pp=183–184|ps=: The name "softmax" can be somewhat confusing. The function is more closely related to the arg max function than the max function. The term "soft" derives from the fact that the softmax function is continuous and differentiable. The arg max function, with its result represented as a one-hot vector, is not continuous or differentiable. The softmax function thus provides a "softened" version of the arg max. The corresponding soft version of the maximum function is <math>\operatorname{softmax}(\mathbf{z})^\top \mathbf{z}</math>. It would perhaps be better to call the softmax function "softargmax," but the current name is an entrenched convention.}} For this section, the term "softargmax" is used to emphasize this interpretation.

Formally, instead of considering the arg max as a function with categorical output <math>1, \dots, n</math> (corresponding to the index), consider the arg max function with [[one-hot]] representation of the output (assuming there is a unique max arg):
:<math>\operatorname{arg\,max}(z_1, \dots, z_n) = (y_1, \dots, y_n) = (0, \dots, 0, 1, 0, \dots, 0),</math>
where the output coordinate <math>y_i = 1</math> if and only if <math>i</math> is the arg max of <math>(z_1, \dots, z_n)</math>, meaning <math>z_i</math> is the unique maximum value of <math>(z_1, \dots, z_n)</math>. For example, in this encoding <math>\operatorname{arg\,max}(1, 5, 10) = (0, 0, 1),</math> since the third argument is the maximum.

This can be generalized to multiple arg max values (multiple equal <math>z_i</math> being the maximum) by dividing the 1 between all max args; formally {{math|1/''k''}} where {{mvar|k}} is the number of arguments assuming the maximum. For example, <math>\operatorname{arg\,max}(1, 5, 5) = (0, 1/2, 1/2),</math> since the second and third argument are both the maximum. In case all arguments are equal, this is simply <math>\operatorname{arg\,max}(z, \dots, z) = (1/n, \dots, 1/n).</math> Points {{mvar|'''z'''}} with multiple arg max values are [[Singular point of an algebraic variety|singular points]] (or singularities, and form the singular set) – these are the points where arg max is discontinuous (with a [[jump discontinuity]]) – while points with a single arg max are known as non-singular or regular points.

With this representation, softargmax is now a smooth approximation of arg max: as {{tmath|\beta \to \infty}}, softargmax converges to arg max. There are various notions of convergence of a function; softargmax converges to arg max [[pointwise convergence|pointwise]], meaning for each fixed input {{math|'''z'''}} as {{tmath|\beta \to \infty}}, <math>\sigma_\beta(\mathbf{z}) \to \operatorname{arg\,max}(\mathbf{z}).</math> However, softargmax does not [[uniform convergence|converge uniformly]] to arg max, meaning intuitively that different points converge at different rates, and may converge arbitrarily slowly. In fact, softargmax is continuous, but arg max is not continuous at the singular set where two coordinates are equal, while the uniform limit of continuous functions is continuous. The failure to converge uniformly is because for inputs where two coordinates are almost equal (and one is the maximum), the arg max is the index of one or the other, so a small change in input yields a large change in output. For example,
<math>\sigma_\beta(1, 1.0001) \to (0, 1),</math> but <math>\sigma_\beta(1, 0.9999) \to (1, 0),</math> and <math>\sigma_\beta(1, 1) = 1/2</math> for all inputs: the closer the points are to the singular set <math>(x, x)</math>, the slower they converge. However, softargmax does [[compact convergence|converge compactly]] on the non-singular set.

Conversely, as {{tmath|\beta \to -\infty}}, softargmax converges to arg min in the same way, where here the singular set is points with two arg ''min'' values. In the language of [[tropical analysis]], the softmax is a [[deformation theory|deformation]] or "quantization" of arg max and arg min, corresponding to using the [[log semiring]] instead of the [[max-plus semiring]] (respectively [[min-plus semiring]]), and recovering the arg max or arg min by taking the limit is called "tropicalization" or "dequantization".

It is also the case that, for any fixed {{mvar|''β''}}, if one input {{tmath|z_i}} is much larger than the others ''relative'' to the temperature, {{tmath|T = 1/β}}, the output is approximately the arg max. For example, a difference of 10 is large relative to a temperature of 1:
:<math>\sigma(0, 10) := \sigma_1(0, 10) = \left(1/(1 + e^{10}), e^{10}/(1 + e^{10})\right) \approx (0.00005, 0.99995)</math>
However, if the difference is small relative to the temperature, the value is not close to the arg max. For example, a difference of 10 is small relative to a temperature of 100:
:<math>\sigma_{1/100}(0, 10) \left(1/(1 + e^{1/10}), e^{1/10}/(1 + e^{1/10})\right) \approx (0.475, 0.525).</math>
As {{tmath|\beta \to \infty}}, temperature goes to zero, <math>T = 1/\beta \to 0</math>, so eventually all differences become large (relative to a shrinking temperature), which gives another interpretation for the limit behavior.

=== Probability theory ===
In [[probability theory]], the output of the softmax function can be used to represent a [[categorical distribution]] – that is, a [[probability distribution]] over {{mvar|K}} different possible outcomes.

=== Statistical mechanics ===
In [[statistical mechanics]], the softmax function is known as the [[Boltzmann distribution]] (or [[Gibbs distribution]]):{{sfn|LeCun|Chopra|Hadsell|Ranzato|Huang|2006|p=7}} the index set <math>{1, \dots, k}</math> are the [[Microstate (statistical mechanics)|microstates]] of the system; the inputs <math>z_i</math> are the energies of that state; the denominator is known as the [[partition function (statistical mechanics)|partition function]], often denoted by {{mvar|''Z''}}; and the factor {{mvar|''β''}} is called the [[coldness]] (or [[thermodynamic beta]], or [[inverse temperature]]).

== Applications ==
The softmax function is used in various [[multiclass classification]] methods, such as [[multinomial logistic regression]] (also known as softmax regression)<ref name="bishop">{{cite book |first=Christopher M. |last=Bishop |year=2006 |title=Pattern Recognition and Machine Learning |publisher=Springer}}</ref>{{rp|206–209}} [http://ufldl.stanford.edu/tutorial/supervised/SoftmaxRegression/], multiclass [[linear discriminant analysis]], [[naive Bayes classifier]]s, and [[artificial neural network]]s.<ref>ai-faq [http://www.faqs.org/faqs/ai-faq/neural-nets/part2/section-12.html What is a softmax activation function?]</ref> Specifically, in multinomial logistic regression and linear discriminant analysis, the input to the function is the result of {{mvar|K}} distinct [[linear function]]s, and the predicted probability for the {{mvar|j}}'th class given a sample vector {{math|'''x'''}} and a weighting vector {{math|'''w'''}} is:

:<math>P(y=j\mid \mathbf{x}) = \frac{e^{\mathbf{x}^\mathsf{T}\mathbf{w}_j}}{\sum_{k=1}^K e^{\mathbf{x}^\mathsf{T}\mathbf{w}_k}}</math>

This can be seen as the [[function composition|composition]] of {{mvar|K}} linear functions <math>\mathbf{x} \mapsto \mathbf{x}^\mathsf{T}\mathbf{w}_1, \ldots, \mathbf{x} \mapsto \mathbf{x}^\mathsf{T}\mathbf{w}_K</math> and the softmax function (where <math>\mathbf{x}^\mathsf{T}\mathbf{w}</math> denotes the inner product of <math>\mathbf{x}</math> and <math>\mathbf{w}</math>). The operation is equivalent to applying a linear operator defined by <math>\mathbf{w}</math> to vectors <math>\mathbf{x}</math>, thus transforming the original, probably highly-dimensional, input to vectors in a {{mvar|K}}-dimensional space <math>\mathbb{R}^K</math>.

=== Neural networks ===
The softmax function is often used in the final layer of a neural network-based classifier. Such networks are commonly trained under a [[Cross entropy|log loss]] (or [[cross-entropy]]) regime, giving a non-linear variant of multinomial logistic regression.

Since the function maps a vector and a specific index ''i'' to a real value, the derivative needs to take the index into account:

:<math> \frac{\partial}{\partial q_k}\sigma(\textbf{q}, i) = \cdots =  \sigma(\textbf{q}, i)(\delta_{ik} - \sigma(\textbf{q}, k))</math>

Here, the [[Kronecker delta]] is used for simplicity (cf. the derivative of a [[sigmoid function]], being expressed via the function itself).

See [[Multinomial logit]] for a probability model which uses the softmax activation function.

=== Reinforcement learning ===
In the field of [[reinforcement learning]], a softmax function can be used to convert values into action probabilities. The function commonly used is:<ref>Sutton, R. S. and Barto A. G. ''Reinforcement Learning: An Introduction''. The MIT Press, Cambridge, MA, 1998. [http://incompleteideas.net/book/ebook/node17.html Softmax Action Selection]</ref>
:<math>
P_t(a) = \frac{\exp(q_t(a)/\tau)}{\sum_{i=1}^n\exp(q_t(i)/\tau)} \text{,}
</math>

where the action value <math>q_t(a)</math> corresponds to the expected reward of following action a and <math>\tau</math> is called a temperature parameter (in allusion to [[statistical mechanics]]). For high temperatures (<math>\tau\to \infty</math>), all actions have nearly the same probability and the lower the temperature, the more expected rewards affect the probability. For a low temperature (<math>\tau\to 0^+</math>), the probability of the action with the highest expected reward tends to 1.

== Properties ==
Geometrically the softmax function maps the [[vector space]] <math>\mathbb{R}^K</math> to the [[Interior (topology)|interior]] of the [[standard simplex|standard <math>(K-1)</math>-simplex]], cutting the dimension by one (the range is a <math>(K-1)</math>-dimensional simplex in <math>K</math>-dimensional space), due to the [[linear constraint]] that all output sum to 1 meaning it lies on a [[hyperplane]].

Along the main diagonal <math>(x, x, \dots, x),</math> softmax is just the uniform distribution on outputs, <math>(1/n, \dots, 1/n)</math>: equal scores yield equal probabilities.

More generally, softmax is invariant under translation by the same value in each coordinate: adding <math>\mathbf{c} = (c, \dots, c)</math> to the inputs <math>\mathbf{z}</math> yields <math>\sigma(\mathbf{z} + \mathbf{c}) = \sigma(\mathbf{z})</math>, because it multiplies each exponent by the same factor, <math>e^c</math> (because <math>e^{z_i + c} = e^{z_i} \cdot e^c</math>), so the ratios do not change:

:<math>\sigma(\mathbf{z} + \mathbf{c})_j = \frac{e^{z_j + c}}{\sum_{k=1}^K e^{z_k + c}} = \frac{e^{z_j} \cdot e^c}{\sum_{k=1}^K e^{z_k} \cdot e^c} = \sigma(\mathbf{z})_j.</math>

Geometrically, softmax is constant along diagonals: this is the dimension that is eliminated, and corresponds to the softmax output being independent of a translation in the input scores (a choice of 0 score). One can normalize input scores by assuming that the sum is zero (subtract the average: <math>\mathbf{c}</math> where <math display="inline">c = \frac{1}{n} \sum z_i</math>), and then the softmax takes the hyperplane of points that sum to zero, <math display="inline">\sum z_i = 0</math>, to the open simplex of positive values that sum to 1<math display="inline">\sum \sigma(\mathbf{z})_i = 1</math>, analogously to how the exponent takes 0 to 1, <math>e^0 = 1</math> and is positive.

By contrast, softmax is not invariant under scaling. For instance, <math>\sigma\bigl((0, 1)\bigr) = \bigl(1/(1 + e), e/(1+e)\bigr)</math> but <math>\sigma\bigl((0, 2)\bigr) = \bigl(1/(1 + e^2), e^2/(1+e^2)\bigr).</math>

The [[standard logistic function]] is the special case for a 1-dimensional axis in 2-dimensional space, say the ''x''-axis in the {{math|(x, y)}} plane. One variable is fixed at 0 (say <math>z_2 = 0</math>), so <math>e^0 = 1</math>, and the other variable can vary, denote it <math>z_1 = x</math>, so <math display="inline">e^{z_1}/\sum_{k=1}^2 e^{z_k} = e^x/(e^x + 1),</math> the standard logistic function, and <math display="inline">e^{z_2}/\sum_{k=1}^2 e^{z_k} = 1/(e^x + 1),</math> its complement (meaning they add up to 1). The 1-dimensional input could alternatively be expressed as the line <math>(x/2, -x/2)</math>, with outputs <math>e^{x/2}/(e^{x/2} + e^{-x/2}) = e^x/(e^x + 1)</math> and <math>e^{-x/2}/(e^{x/2} + e^{-x/2}) = 1/(e^x + 1).</math>

The softmax function is also the gradient of the [[LogSumExp]] function, a [[smooth maximum]]; defining:
:<math>\operatorname{LSE}(z_1, \dots, z_n) = \log\left( \exp(z_1)+ \cdots + \exp(z_n) \right),</math>
the partial derivatives are:

: <math> \partial_i \operatorname{LSE}(\mathbf{x}) = \frac{\exp x_i}{\sum_i \exp x_i}.</math>

Expressing the partial derivatives as a vector with the [[gradient]] yields the softmax.

== History ==
The softmax function was used in [[statistical mechanics]] as the [[Boltzmann distribution]] in the foundational paper {{harvtxt|Boltzmann|1868}}, formalized and popularized in the influential textbook {{harvtxt|Gibbs|1902}}.

The use of the softmax in [[decision theory]] is credited to {{harvtxt|Luce|1959}},{{sfn|Gao|Pavel|2017|p=1}} who used the axiom of [[independence of irrelevant alternatives]] in [[rational choice theory]] to deduce the softmax in [[Luce's choice axiom]] for relative preferences.

In machine learning, the term "softmax" is credited to John S. Bridle in two 1989 conference papers, {{harvtxt|Bridle|1990a}}:{{sfn|Gao|Pavel|2017|p=1}} and {{harvtxt|Bridle|1990b}}:<ref name="sako2018"/>
{{quote
|We are concerned with feed-forward non-linear networks (multi-layer perceptrons, or MLPs) with multiple outputs. We wish to treat the outputs of the network as probabilities of alternatives (''e.g.'' pattern classes), conditioned on the inputs. We look for appropriate output non-linearities and for appropriate criteria for adaptation of the parameters of the network (''e.g.'' weights). We explain two modifications: probability scoring, which is an alternative to squared error minimisation, and a normalised exponential ('''softmax''') multi-input generalisation of the logistic non-linearity.{{sfn|Bridle|1990a|p=227}}}}
{{quote
|For any input, the outputs must all be positive and they must sum to unity. ...

Given a set of unconstrained values, {{tmath|V_j(x)}}, we can ensure both conditions by using a Normalised Exponential transformation:
:<math>Q_j(x) = \left. e^{V_j(x)} \right/ \sum_k e^{V_k(x)} </math>
This transformation can be considered a multi-input generalisation of the logistic, operating on the whole output layer. It preserves the rank order of its input values, and is a differentiable generalisation of the ‘winner-take-all’ operation of picking the maximum value. For this reason we like to refer to it as '''softmax'''.{{sfn|Bridle|1990b|p=213}}}}

== Example ==
If we take an input of [1, 2, 3, 4, 1, 2, 3], the softmax of that is [0.024, 0.064, 0.175, 0.475, 0.024, 0.064, 0.175]. The output has most of its weight where the '4' was in the original input. This is what the function is normally used for: to highlight the largest values and suppress values which are significantly below the maximum value. But note: softmax is not scale invariant, so if the input were [0.1, 0.2, 0.3, 0.4, 0.1, 0.2, 0.3] (which sums to 1.6) the softmax would be [0.125, 0.138, 0.153, 0.169, 0.125, 0.138, 0.153]. This shows that for values between 0 and 1 softmax, in fact, de-emphasizes the maximum value (note that 0.169 is not only less than 0.475, it is also less than the initial proportion of 0.4/1.6=0.25).

Computation of this example using [[Python (programming language)|Python]] code:

<syntaxhighlight lang="python">
>>> import numpy as np
>>> a = [1.0, 2.0, 3.0, 4.0, 1.0, 2.0, 3.0]
>>> np.exp(a) / np.sum(np.exp(a)) 
array([0.02364054, 0.06426166, 0.1746813, 0.474833, 0.02364054,
       0.06426166, 0.1746813])
</syntaxhighlight>

Here is an example of [[Julia (programming language)|Julia]] code:

<syntaxhighlight lang="jlcon">
julia> A = [1.0, 2.0, 3.0, 4.0, 1.0, 2.0, 3.0];  # semicolon to suppress interactive output

julia> exp.(A) ./ sum(exp.(A))
7-element Array{Float64,1}:
 0.0236405
 0.0642617
 0.174681
 0.474833
 0.0236405
 0.0642617
 0.174681
</syntaxhighlight>

Here is an example of [[R (programming language)|R]] code:

<syntaxhighlight lang="r">
> z <- c(1.0, 2.0, 3.0, 4.0, 1.0, 2.0, 3.0)
> softmax <- exp(z)/sum(exp(z))
> softmax
[1] 0.02364054 0.06426166 0.17468130 0.47483300 0.02364054 0.06426166 0.17468130
</syntaxhighlight>

== See also ==
* [[Softplus]]
* [[Multinomial logistic regression]]
* [[Dirichlet distribution]] – an alternative way to sample categorical distributions
* [[Partition function (statistical mechanics)|Partition function]]

== Notes ==
{{notelist}}

== References ==
{{reflist|refs=
<ref name="sako2018">{{cite web|first=Yusaku |last=Sako
|title=Is the term "softmax" driving you nuts?
|url=https://medium.com/@u39kun/is-the-term-softmax-driving-you-nuts-ee232ab4f6bd
|date=2018-06-02
|publisher=[[Medium (website)|Medium]]
}}</ref>
}}
{{refbegin}}
* {{cite journal |last=Boltzmann |first=Ludwig |authorlink=Ludwig Boltzmann
|year=1868
|title=Studien über das Gleichgewicht der lebendigen Kraft zwischen bewegten materiellen Punkten
|trans-title=Studies on the balance of living force between moving material points
|journal=Wiener Berichte |volume=58 |pages=517–560
}}
* {{cite book
|first=Josiah Willard |last=Gibbs |authorlink=Josiah Willard Gibbs
|year=1902
|title=[[Elementary Principles in Statistical Mechanics]]
}}
* {{cite conference
|first=John S. |last=Bridle
|year=1990a
|title=Probabilistic Interpretation of Feedforward Classification Network Outputs, with Relationships to Statistical Pattern Recognition
|pages=227–236
|conference=Neurocomputing: Algorithms, Architectures and Applications (1989)
|editor1=Soulié F.F.
|editor2=Hérault J. 
|series=NATO ASI Series (Series F: Computer and Systems Sciences)
|volume=68
|publisher=Springer
|location=Berlin, Heidelberg
|doi=10.1007/978-3-642-76153-9_28
}}
* {{cite conference
|first=John S. |last=Bridle
|year=1990b
|title=Training Stochastic Model Recognition Algorithms as Networks can Lead to Maximum Mutual Information Estimation of Parameters
|url=https://papers.nips.cc/paper/195-training-stochastic-model-recognition-algorithms-as-networks-can-lead-to-maximum-mutual-information-estimation-of-parameters
|conference=[[Advances in Neural Information Processing Systems]] 2 (1989)
|editor = D. S. Touretzky
|publisher = Morgan-Kaufmann
}}
* {{cite book
|first1=Yann |last1=LeCun
|author-link1=Yann LeCun
|first2=Sumit |last2=Chopra
|first3=Raia |last3=Hadsell
|first4=Marc’Aurelio |last4=Ranzato
|first5=Fu Jie |last5=Huang
|chapter=A Tutorial on Energy-Based Learning
|chapter-url=http://yann.lecun.com/exdb/publis/pdf/lecun-06.pdf
|isbn=978-0-26202617-8
|publisher=MIT Press
|year=2006
|title=Predicting Structured Data
|url=https://mitpress.mit.edu/books/predicting-structured-data
|series=Neural Information Processing series
|editors=Gökhan Bakır, Thomas Hofmann, Bernhard Schölkopf, Alexander J. Smola, Ben Taskar and S.V.N Vishwanathan
}}
* {{cite book |title=Deep Learning |year=2016
|first1=Ian |last1=Goodfellow |authorlink1=Ian Goodfellow
|first2=Yoshua |last2=Bengio |authorlink2=Yoshua Bengio
|first3=Aaron |last3=Courville
|publisher=MIT Press
|url=http://www.deeplearningbook.org
|isbn=978-0-26203561-3
|chapter-url=https://www.deeplearningbook.org/contents/mlp.html
|section=6.2.2.3 Softmax Units for Multinoulli Output Distributions
|pages=180–184
}}
* {{cite arXiv|eprint=1704.00805|last1=Gao|first1=Bolin|title=On the Properties of the Softmax Function with Application in Game Theory and Reinforcement Learning|last2=Pavel|first2=Lacra|class=math.OC|year=2017}}
{{refend}}

[[Category:Computational neuroscience]]
[[Category:Logistic regression]]
[[Category:Artificial neural networks]]
[[Category:Functions and mappings]]
[[Category:Articles with example Python code]]
[[Category:Exponentials]]
[[Category:Articles with example Julia code]]