Je suis désolé, mais je suis meilleur en informatique qu’en statistique, et j’avoue être un peu perdu.. Je vais essayer de creuser cela.
La méthode de Monte-Carlo est une méthode informatique pour estimer une valeur mathématique ; un exemple largement utilisé est celui de l’estimation du nombre Pi : source
Dans le cas présent, la valeur que nous souhaitons estimer est l’espérance mathématique de votre portefeuille ainsi que l’écart-type, afin d’obtenir une mesure du ratio rendement/risque (espérance/écart-type) de votre stratégie.
Je ne comprend toujours pas.
Ma simulation sort le PnL à la fin de la période, pour un marché aléatoire. je trace *tous* les resultats (une fois triés)
(Faux) Example avec 3 simus:
PEA +5 +6 +4 --> +4 / +5 / +6
Levier x2 +10 -2 +3 --> -3 / +3 / +10
Cela permet de voir en un coup d’oeil à quelle frequence (=probabilité) le levier fait mieux, mais aussi de combien.
Non ?
Oui et non.
Oui car votre graphique contient bien les informations importantes : probabilité/gain associé
Non car vous avez tracé une fonction de répartition qui répond à la question "combien de fois A a été mieux que B ? " et non pas "à quel point est-ce que A est mieux que B ?".
Exemple simple : le loto
Supposons qu’un ticket coute 1€. Vous perdez 99 999 fois sur 100 000 et vous gagnez 1M€ 1 fois sur 100 000.
D’après votre fonction cumulative, la stratégie "ne pas jouer" est plus souvent gagnante que "jouer" (99,999%). Pourtant, d’après les statistiques, vous auriez tout intérêt à jouer (espérance de gain de 9€). Cela n’apparaît pas au premier coup d’oeil sur la fonction de répartition.
Pour le point 1, je vais essayer quand j’aurais compris comment faire.
Pour les point 2/3/4, je pense que ma simulation est plus précise? (avec prise en compte du montant initial, des versements progressifs, des appels de marge éventuels..)
Pour une première approche de la méthode de rejet, l’explication de Wikipédia me semble pas mal.
Pour les points 2/3/4, cela ne sert pas à grand chose de simuler les VP dans la mesure où ils seront soumis à la même distribution de probabilité. C’est plutôt un genre de "martingale" inventée par les banquiers pour vous faire croire que vous augmentez vos chances de gain en achetant régulièrement ; pourtant lorsque vous avez un héritage de 100k€ ils vous conseillent une AV avec une répartition 70/30… les petits filous !
Pour s’en apercevoir, il suffit de calculer le taux de rentabilité interne de votre stratégie (je vous laisse ça en devoir à faire à la maison )
Voici quelques graphiques pour étayer mon propos :
Sur ce graphique, vous devez certainement reconnaître celui que j’ai posté dans mon dernier commentaire. Je n’ai fait que superposer l’histogramme des rendement aléatoires générés pas la méthode de rejet pour vous montrer son efficacité ; votre méthode aurait généré un histogramme proche du "gaussian fit" qui n’est pas du tout adapté.
Sur cette image, j’ai tracé l’histogramme des performances(*) de 2000 portefeuilles sur une période 9 ans pour un effet de levier (k-multiplier) allant de 1 à 4. Attention, les échelles sont logarithmiques (en y aussi). On constate (sans trop de surprises) que le pic pour k=1 se situe aux alentours de 2-3 alors que le pic pour k=4 se situe aux alentours de 8-10 (donc 4 fois plus) ; la position du pic dépend (naturellement) du levier choisi.
Mais plus important, l’histogramme "s’étale" très rapidement (et même dangereusement) lorsque l’on augmente l’effet de levier.
Avec un levier de 1, vous faites avec une bonne probabilité x2-3 en 9 ans (+100%/+200%) et à première vue, moins de 10% de chances de faire -20%.
Avec un levier de 4, vous avez des chances raisonnables de faire x100 (je n’ai pas calculé mais ça doit dépasser les 10-15% à vue d’oeil) mais également un bon 15% de chances de perdre 90% de votre investissement. Voudriez-vous jouer à la la roulette russe sur 10 000€ ?
Enfin, j’ajoute que, bien que ma distribution simulée est "proche" de la distribution historique, cela n’est pas encore optimal comparé aux modèles ARCH/GARCH qui supposent une volatilité aléatoire avec mémoire (et plus proche de la réalité). Mais ce sera (peut-être) pour une prochaine fois.
Caceray
modification : (*) performance=prix_final/prix_initial ; le 10^0 indique une performance de 1, soit valeur finale = valeur initiale.