Leçon
La théorie — Laboratoire de modélisation par morceaux
Une fonction définie par morceaux n’est pas une règle avec une exception ; ce sont deux règles complètes, plus l’indication de celle qui s’applique où. Ici une règle linéaire gouverne tout ce qui est en dessous d’un point de rupture c et une seconde gouverne tout à partir de c. La question intéressante n’est jamais ce que fait chaque morceau — ce sont de simples droites — mais ce qui se passe à la couture.
Ce que signifie chaque symbole
c- le point de rupture : l’entrée à laquelle la règle change.
f(c⁻)- vers quoi tend le morceau de gauche quand
xapprochecpar en dessous. C’est une limite, pas une valeur : enc, la règle de gauche ne s’applique plus. f(c⁺)- la valeur que prend réellement le morceau de droite en
c, et donc la valeur de la fonction en ce point. jumpf(c⁺) − f(c⁻). Nul exactement quand les deux morceaux se rejoignent.
D’où vient la formule
- Écrivez la définition avec sa condition :
f(x) = m₁x + b₁pourx < c, etf(x) = m₂x + b₂pourx ≥ c. Notez à quel morceau appartient le point de rupture lui-même — ici celui de droite, ce qui explique le cercle creux sur la branche gauche et le point plein sur la droite. - En
c, la fonction a une valeur donnée par la règle de droite,m₂c + b₂. La règle de gauche n’y a pas de valeur, mais elle a une destination : lorsquexmonte versc,m₁x + b₁tend versm₁c + b₁. C’est la limite à gauche. - La continuité en un point signifie que les limites des deux côtés et la valeur coïncident. Ici la limite à droite et la valeur sont une seule et même chose : deux des trois sont donc automatiques, et toute la condition se réduit à une équation,
m₁c + b₁ = m₂c + b₂. - Une équation, et n’importe quel coefficient peut être déplacé pour la satisfaire :
b₂ = (m₁ − m₂)c + b₁. La continuité est donc toujours atteignable et jamais gratuite — il faut y consacrer un paramètre. Quand on ne le fait pas, le reliquat est exactement ce que l’affichage appelle le saut :(m₂c + b₂) − (m₁c + b₁).
Comment lire ce que vous voyez
Le segment bleu est la règle de gauche et l’orange celle de droite, chacun tracé aussi loin que sa propre condition l’autorise. Au point de rupture, la branche bleue s’achève par un cercle creux — une destination qu’elle n’occupe jamais — et l’orange commence par un point plein, la valeur réelle de f(c). Quand les deux coïncident, les cercles se superposent et la ligne du saut affiche 0. Le point turquoise est f(x₀) pour ce que vous mettez dans la case d’évaluation, et les trois dernières lignes traduisent l’image en nombres : si les morceaux se rejoignent, de combien la pente change au passage, et si la marche descendante crée une plage d’entrées qui coûtent plus cher que le point de rupture lui-même.
- Suppose
- Les deux morceaux sont des droites, il y a exactement un point de rupture, et il est situé à un endroit connu. La convention semi-ouverte compte et relève d’un choix : la règle de droite possède
x = c, doncf(c)vautm₂c + b₂. Inversez la convention et la valeur au seul pointcchange, mais pas le saut. - Ne tient plus quand
- La vraie difficulté surgit quand
cn’est pas connu mais estimé. Un point de rupture estimé à partir des mêmes données que celles servant à ajuster les droites est un paramètre choisi pour minimiser l’erreur : les erreurs types usuelles de la régression sont alors trop petites et l’ajustement paraît meilleur qu’il ne l’est — c’est précisément pour cela qu’existent les méthodes de régression segmentée. Et un modèle à deux morceaux ne dit rien d’un troisième régime : prolonger l’une des droites au-delà de la plage observée, c’est supposer une règle jamais testée.
Problème entièrement résolu
-
Un seuil tarifaire à 5 et la plage où acheter plus coûte moins cher 6 étapes
Un tarif applique 2x + 3 en dessous d'un point de rupture situé à 5 et 0,8x + 8 au-dessus. Calculez la valeur en 7 — puis trouvez l'intervalle où acheter plus coûte moins cher.
-
Deux droites et une règle déterminant celle qui s'applique. La seule précaution concerne la frontière : en x = 5 exactement, c'est la seconde branche qui s'applique.
-
Approchez le point de rupture des deux côtés. La limite à gauche et la valeur réelle diffèrent, et l'ampleur de ce désaccord constitue le saut.
-
Les pentes diffèrent également, et ce de façon indépendante. Une fonction peut être continue avec un point anguleux, ou présenter un saut sans point anguleux ; celle-ci fait les deux.
-
Passons maintenant à la question intéressante. Déterminons où la première branche a déjà dépassé la valeur initiale de la seconde branche.
-
Cette inégalité est vérifiée de 4,5 jusqu'au point de rupture, et sur tout cet intervalle, l'achat le plus grand est le moins cher.
-
La condition qui éliminerait cette fenêtre mérite d'être formulée, car c'est celle à laquelle satisfait un barème bien conçu.
Réponse
13,6, et entre 4,5 et 5, il faut toujours acheter 5. La seconde branche commence à 12 alors que la première continue de monter au-delà, de sorte que quiconque se trouve dans cette fenêtre paie plus pour avoir moins. Cette fenêtre n'est pas un artefact d'arrondi — elle représente 10% de l'ensemble de la première tranche, et c'est exactement ce qu'engendre un tarif discontinu. Les barèmes réels l'évitent en raccordant les morceaux : la continuité exige b₂ − b₁ = (m₁ − m₂)c, ce qui nécessiterait ici que 5 soit égal à 6. Les tranches d'imposition marginales sont conçues ainsi à dessein, c'est pourquoi l'idée reçue selon laquelle une augmentation peut réduire votre salaire net est fausse pour l'impôt sur le revenu et vraie pour presque toutes les remises sur quantité jamais imprimées.
-
Références (2)
- The theory section ends on what happens when the breakpoint is estimated rather than known. This is the method that handles it: V. M. R. Muggeo, "Estimating regression models with unknown break-points." Statistics in Medicine 22(19), 3055–3071, 2003.
- Insight block 3 — the piecewise-linear kink at the centre of modern networks: X. Glorot, A. Bordes and Y. Bengio, "Deep Sparse Rectifier Neural Networks." Proceedings of the 14th International Conference on Artificial Intelligence and Statistics (AISTATS), PMLR 15, 315–323, 2011.