BCE Maths approfondies HEC/ESCP ECS 2019, épreuve 2Sujet et corrigé

Epreuve de maths approfondies - ECS 2019

Pas encore noté

Explorer les annales

BanqueBCE Année2019 FilièreECS MatièreMathématiques

Annales BCE ECS Mathématiques BCE ECS 2019 ECS Mathématiques Mathématiques 2019

Téléchargements

Télécharger le sujet (PDF)

Corrigé

Lire en ligne LaTeX Word

Rapport du jury : non disponible

Source du document : annales-prepa.fr

Description

Annale de maths approfondies BCE HEC/ESCP pour la filiere ECS, session 2019.

Ces sujets peuvent vous intéresser

Lecture du sujet en ligne

L'énoncé complet, avec les formules et les figures, sans ouvrir le PDF.

Sujet en PDF

Conception : HEC Paris - ESCP Europe

OPTION SCIENTIFIQUE

MATHÉMATIQUES II

Jeudi 2 mai 2019, de 8 h. à 12 h.

La présentation, la lisibilité, l'orthographe, la qualité de la rédaction, la clarté et la précision des raisonnements entreront pour une part importante dans l'appréciation des copies.
Les candidats sont invités à encadrer dans la mesure du possible les résultats de leurs calculs.
Aucun document n'est autorisé. L'utilisation de toute calculatrice et de tout matériel électronique est interdite. Seule l'utilisation d'une règle graduée est autorisée.
Si au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, il la signalera sur sa copie et poursuivra sa composition en expliquant les raisons des initiatives qu'il sera amené à prendre.

La régression logistique permet de modéliser l'influence qu'exercent des facteurs exogènes sur une variable binaire, c'est-à-dire une variable ne pouvant prendre que deux valeurs.
Outre son domaine d'application privilégié qui est l'apprentissage automatique (machine learning), la régression logistique est couramment utilisée aussi bien en médecine qu'en actuariat et en économétrie.

Partie I. Fonction logistique et lois logistiques

On appelle fonction logistique la fonction

Λ

définie sur

R

par:

\forallx \in R, Λ(x) = 1/(1 + e^(- x))

.
1.a) Montrer que

Λ

est une bijection de

R

sur

]0, 1[

, dont la bijection réciproque est la fonction

L

définie par :

\forallx \in ]0, 1[, L(x) = ln(x/(1 - x)).

b) Calculer la dérivée de la fonction

Λ

.
c) Justifier l'existence d'un unique réel

x_0

tel que :

Λ(x_0) = x_0

.
d) Établir pour tout

x \in R

, l'inégalité :

|Λ(x) - x| ⩽ |x - x_0|

.
2. Le script Scilab suivant, dont la ligne (1) définit la fonction

Λ

, permet de calculer une valeur approchée de

x_0

par la méthode de dichotomie.
(1) deff('

y = Lambda(x)^', ^' y = 1/(1 + exp(- x))^')

;
(2)

a = 0

;
(3)

b = 1

;
(4) eps= ........;
(5) while b-a>eps ;
(6)

c = (a + b)/2

;
(7) if Lambda

(c) > c

then ........; else

b = \dots\dots\dots

; end;
(8) end;
(9)

xO = (a + b)/2

a) Compléter la ligne (7) et justifier le choix des valeurs affectées en lignes (2) et (3) aux variables a et b.
b) Quelle valeur maximale peut-on affecter en ligne (4) à la variable eps pour être assuré que l'erreur d'approximation commise ne dépasse pas

10^(- 4)

?
c) Que peut-on dire de la valeur numérique obtenue par l'instruction (10) suivante ?
(10) Lambda ( x 0 ) -x 0
3. On note

λ

la dérivée de la fonction

Λ

.
a) Vérifier que

λ

est une densité de probabilité.
b) Préciser la parité de la fonction

λ

; donner l'allure de sa courbe représentative dans le plan rapporté à un repère orthogonal et en déterminer les points d'inflexion.

On dit qu'une variable aléatoire

Z

suit la loi logistique standard si elle admet la fonction

λ

pour densité. Pour tout couple

(r, s) \in R \times R_+^*

, on dit qu'une variable aléatoire

Y

suit la loi logistique

L(r, s)

si la variable aléatoire

Z

définie par

Z = (Y - r)/s

suit la loi logistique standard.
4.a) Justifier qu'une variable aléatoire qui suit une loi logistique

L(r, s)

admet des moments de n'importe quel ordre et en indiquer l'espérance.
b) En utilisant la méthode d'inversion, écrire le script d'une fonction Scilab, function

S = grandlogis(n, p, r, s)

, fournissant pour tout couple (

n, p

) d'entiers strictement positifs, une matrice

S

n

lignes et

p

colonnes dont les coefficients sont des simulations de variables aléatoires indépendantes suivant la loi logistique

L(r, s)

.
c) Décrire un procédé permettant de calculer une valeur approchée de la variance de la loi logistique standard à l'aide de la fonction grandlogis.
5. Soit

U_1

U_2

deux variables aléatoires indépendantes suivant chacune la loi exponentielle de paramètre 1 .
a) Montrer que la variable aléatoire

Z = ln((U_1)/(U_2))

suit la loi logistique standard (on pourra utiliser un changement de variable exponentiel, c'est-à-dire de la forme

t = e^x

).
b) En déduire un nouveau script Scilab permettant de simuler une variable aléatoire suivant la loi logistique standard à l'aide de la fonction grand.

Partie II. Variance de la loi logistique standard

Pour tout couple $(a, b) \in R^2$ , on note $Im(z)$ la partie imaginaire $b$ du nombre complexe $z = a + ib$ .
Pour tout polynôme $P = \sum_(k = 0)^d a_k X^k \in R[X]$ de degré $d \in N$ , les termes non nuls $a_k X^k$ sont appelés les monômes de $P$ et les $a_k$ leurs coefficients.
Dans la factorisation $P = a_d\prod_(k = 1)^d(X - z_k)$ de $P$ dans $C[X]$ (lorsque $d \neq 0$ ), la somme $\sum_(k = 1)^d z_k$ est appelée la somme des racines complexes de $P$ , que les nombres complexes $z_1, z_2, \dots, z_d$ soient distincts ou non.

Pour tout

n \in N

, on pose :

P_n = \sum_(k = 0)^n(- 1)^k((2n + 1)/(2k + 1))(X - 1)^(n - k)

.
6.a) Expliciter les polynômes

P_0

P_1

.
b) Pour tout

n \in N^*

, préciser le degré du polynôme

P_n

et donner les coefficients de ses deux monômes de plus hauts degrés.
c) Utiliser le résultat précédent pour montrer que pour tout

n \in N^*

, la somme des racines complexes de

P_n

est égale à

(2n(n + 1))/3

.
7. Soit

x \in R

n \in N

.
a) Justifier les égalités suivantes :

sin((2n + 1)x) = Im((cos(x) + isin(x))^(2n + 1)) = \sum_(k = 0)^n(- 1)^k((2n + 1)/(2k + 1))cos^(2(n - k))(x) \times sin^(2k + 1)(x)

b) En déduire, pour tout

x \in ]0, π[

, la relation :

(sin((2n + 1)x))/(sin^(2n + 1)(x)) = P_n(1/(sin^2(x)))

.
c) À l'aide du résultat de la question 6.c), montrer que pour tout

n \in N^*

, on a :

\sum_(k = 1)^n 1/(sin^2((kπ)/(2n + 1))) = (2n(n + 1))/3

Soit $x \in ]0, π/2[$ .
a) Justifier les inégalités suivantes: $sin(x) ⩽ x ⩽ tan(x)$ et $1/(sin^2(x)) - 1 ⩽ 1/(x^2) ⩽ 1/(sin^2(x))$ .
b) En utilisant le résultat de la question 7.c), en déduire, pour tout $n \in N^*$ , l'encadrement :

(n(2n - 1))/3 ⩽ ((2n + 1)^2)/(π^2)\sum_(k = 1)^n 1/(k^2) ⩽ (2n(n + 1))/3.

c) Établir l'égalité :

\sum_(k = 1)^(+ \infty)1/(k^2) = (π^2)/6

.
9. Soit

Z

une variable aléatoire suivant la loi logistique standard.
a) À l'aide d'une intégration par parties, justifier que la variance de

Z

, notée

V(Z)

, vérifie l'égalité :

V(Z) = 4\int_0^(+ \infty)(xe^(- x))/(1 + e^(- x)) dx

b) Établir pour tout

n \in N

, l'égalité :

ù \int_0^(+ \infty)(xe^(- x))/(1 + e^(- x)) dx = \sum_(k = 0)^n(- 1)^k\int_0^(+ \infty)xe^(- (k + 1)x) dx + I_n, où I_n = (- 1)^(n + 1)\int_0^(+ \infty)(xe^(- (n + 2)x))/(1 + e^(- x)) dx

c) Montrer que l'intégrale

I_n

tend vers 0 lorsque

n

tend vers

+ \infty

et en déduire l'égalité :

\int_0^(+ \infty)(xe^(- x))/(1 + e^(- x)) dx = \sum_(k = 0)^(+ \infty)((- 1)^k)/((k + 1)^2)

d) En utilisant la formule établie en 8.c), déduire de l'égalité précédente que la variance de

Z

est égale à

(π^2)/3

.
10.a) Établir la convergence des deux intégrales

\int_0^(+ \infty)ln(x)e^(- x) dx

\int_0^(+ \infty)(ln(x))^2 e^(- x) dx

.
b) On pose

I = \int_0^(+ \infty)ln(x)e^(- x) dx

J = \int_0^(+ \infty)(ln(x))^2 e^(- x) dx

En utilisant le résultat de la question 5.a), calculer

J - I^2

Partie III. Estimation à partir de données binaires

Dans cette partie,

θ

est un paramètre réel inconnu et

F

désigne la fonction de répartition d'une variable aléatoire à densité dont une densité

f

est continue et strictement positive sur

R

.
Soit

(Y_n)_(n \in N^*)

une suite de variables aléatoires indépendantes définies sur un espace probabilisé

(Ω, A, P_θ)

suivant chacune la loi de Bernoulli de paramètre

F(θ)

.
11. Justifier que

F

est une bijection de

R

sur

]0, 1[

. On note

F^(- 1)

sa bijection réciproque.
12. Pour tout

n \in N^*

, on pose :

Y¯_n = 1/n\sum_(j = 1)^n Y_j

Montrer que la suite

(\sqrtn(Y¯_n - F(θ)))_(n \in N^*)

converge en loi vers une variable aléatoire suivant une loi normale centrée dont on précisera la variance.
13. Pour tout

n \in N^*

et tout

ω \in Ω

, on pose :

T_n(ω) = {F^(− 1)(Y¯_n(ω)), si 0 < Y¯_n(ω) < 1; 0, sinon

De plus, pour tout

n \in N^*

, on note

E_n

l'événement

[0 < Y¯_n < 1]

.
a) Calculer

P_θ(E_n)

et trouver la limite de cette probabilité lorsque

n

tend vers

+ \infty

.
b) Soit

x \in R

n \in N^*

.
(i) Établir l'égalité ensembliste

{ω ∈ E_n/T_n(ω) ⩽ x} = [Y¯_n ⩽ F(x)] ∩ E_n

et montrer que

[T_n ⩽ x]

est un élément de la tribu

A

.
(ii) Justifier l'encadrement :

P_θ([Y¯_n ⩽ F(x)] \cap E_n) ⩽ P_θ([T_n ⩽ x]) ⩽ P_θ([Y¯_n ⩽ F(x)] \cap E_n) + 1 - P_θ(E_n)

c) Montrer que pour tout

x \neq θ

, on a :

lim_(n → + ∞)P_θ([T_n ⩽ x]) = {0, si x < θ; 1, si x > θ

.
d) En déduire que

(T_n)_(n \in N^*)

est une suite convergente d'estimateurs du paramètre

θ

.
14. Pour tout

n \in N^*

et tout

ω \in Ω

, on pose :

U_n(ω) = {(T_n(ω) − θ)/(Y¯_n(ω) − F(θ)), si Y¯_n(ω) ≠ F(θ); 1/(f(θ)), si Y¯_n(ω) = F(θ)

On admet sans démonstration que pour tout

n \in N^*, U_n

est une variable aléatoire sur (

Ω, A, P_θ

).
a) Soit

ε > 0

Pour tout

n \in N^*

, on note

B_n(ε)

l'événement

[|U_n - 1/(f(θ))| ⩽ ε]

.
(i) Établir l'existence d'un réel

α > 0

tel que :

\forallx \in [θ - α, θ + α], |1/(f(x)) - 1/(f(θ))| ⩽ ε

.
(ii) Pour un tel

α

, justifier l'inclusion :

[|T_n - θ| ⩽ α] \cap E_n \subset B_n(ε)

, où

E_n

a été défini dans la question 13.
b) Montrer que la suite

(U_n)_(n \in N^*)

converge en probabilité vers

1/(f(θ))

.
c) En déduire que la suite

(\sqrtn(T_n - θ))_(n \in N^*)

converge en loi vers une variable aléatoire suivant une loi normale centrée dont on précisera la variance.

Partie IV. Régression logistique

Dans toute cette partie, $p$ désigne un entier supérieur ou égal à 2 .
Pour tout couple $(n, m) \in (N^*)^2$ , on note $M_(n, m)(R)$ l'ensemble des matrices à $n$ lignes et $m$ colonnes à coefficients réels et $^t M$ la transposée de toute matrice $M \in M_(n, m)(R)$ .
Pour tout $m \in N^*$ , le produit scalaire usuel de deux vecteurs $u$ et $v$ de $R^m$ est noté $⟨u, v⟩$ . Si $U$ et $V$ sont les matrices colonnes représentant $u$ et $v$ dans la base canonique, le produit scalaire $⟨u, v⟩$ est donc l'unique coefficient de la matrice $^t UV$ .
On rappelle que les fonctions $Λ$ et $L$ ont été définies dans la partie I.

Dans cette partie, on note

Y

une variable aléatoire de Bernoulli, dite variable endogène, dont la loi dépend du niveau de

p

facteurs exogènes.
L'influence de ces facteurs sur la loi de

Y

est résumée par la fonction

b

qui associe à un vecteur

x \in R^p

, la probabilité

b(x)

que

Y

soit égale à 1 lorsque les niveaux des facteurs sont donnés par les composantes du vecteur

x

Dans le modèle de régression logistique envisagé dans cette partie, la fonction

b

est supposée de la forme:

b : x ⟼ Λ(⟨α, x⟩)

où

α = (α_1, α_2, \dots, α_p)

est un vecteur de

R^p

dont les composantes

α_1, α_2, \dots, α_p

sont des paramètres inconnus qui représentent les degrés d'influence des divers facteurs exogènes sur la variable endogène

Y

.
Pour estimer les paramètres du modèle, on dispose de

k

vecteurs

x^((1)), x^((2)), \dots, x^((k))

R^p(k \in N^*)

et pour tout

i \in [ [1, k] ]

, d'une suite

(Y_(i, n))_(n \in N^*)

de variables aléatoires indépendantes suivant chacune la loi de Bernoulli de paramètre

b(x^((i))) = Λ(⟨α, x^((i))⟩)

.
Pour chaque indice fixé

i

et pour tout

n \in N^*

, les variables aléatoires

Y_(i, 1), Y_(i, 2), \dots, Y_(i, n)

définissent donc un

n

-échantillon associé à la loi de la variable endogène lorsque les niveaux des facteurs exogènes sont les composantes

x_1^((i)), x_2^((i)), \dots, x_p^((i))

du vecteur

x^((i))

dans la base canonique de

R^p

.
15. On note respectivement

A

M

la matrice du vecteur

α

et la matrice de la famille (

x^((1)), x^((2)), \dots, x^((k))

) dans la base canonique de

R^p

A = (α_1; α_2; ⋮; α_p) \in M_(p, 1)(R) et M = (x_1^((1)), \dots, x_1^((k)); ⋮, ⋮; x_p^((1)), \dots, x_p^((k))) \in M_(p, k)(R)

On suppose que le rang de la matrice

M

est égal à

p

.
a) Montrer que la matrice

M^t M

est inversible.
b) Montrer que pour toute matrice

H \in M_(k, 1)(R)

, la matrice

U \in M_(p, 1)(R)

pour laquelle l'unique coefficient de la matrice

^t(^t MU - H)(^t MU - H)

est le plus petit possible, est la matrice

(M^t M)^(- 1)MH

.
c) Expliquer pourquoi les lois des variables aléatoires

Y_(i, n)

ne suffiraient pas à définir le vecteur

α

si le rang de

M

n'était pas égal à

p

.
16. Pour tout

n \in N^*

et tout

i \in [ [1, k] ]

, on pose

Y¯_(i, n) = 1/n\sum_(j = 1)^n Y_(i, j)

et pour tout

ω \in Ω

T_(i, n) = {L(Y¯_(i, n)(ω)), si 0 < Y¯_(i, n)(ω) < 1; 0, sinon

a) Soit

(c_1, c_2, \dots, c_k) \in R^k

. En utilisant les résultats de la partie III, montrer que

(\sum_(i = 1)^k c_i T_(i, n))_(n \in N^*)

est une suite convergente d'estimateurs du paramètre

\sum_(i = 1)^k c_i⟨α, x^((i))⟩

.
b) Pour tout

n \in N^*

et tout

ω \in Ω

, on pose,

H_n(ω) = (T_(1, n)(ω); T_(2, n)(ω); ⋮; T_(k, n)(ω))

(A_(1, n)(ω); A_(2, n)(ω); ⋮; A_(p, n)(ω)) = (M^t M)^(- 1)MH_n(ω)

Montrer que pour tout

j \in [ [1, p] ]

, la suite

(A_(j, n))_(n \in N^*)

est une suite convergente d'estimateurs de

α_j

FIN

Pas de description pour le moment

Commentaires• BCE Maths approfondies HEC/ESCP ECS 2019, épreuve 2

Connectez-vous pour participer aux discussions

Partagez vos avis, posez des questions et échangez avec la communauté

Chargement des commentaires...