Take the 2-minute tour ×
Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

Binomial Distribution (binomial model inspired by Merton)

Hi I want to do the following and I hope that anyone of you can help me out.

In short, I have a portfolio consisting of 10 companies, whereas their asset values are driven by a common economic background factor u and the assets in question are correlated. Now I want to calculate the probability that e.g. 3 of 10 companies default simultaneously.

In order to do so, I choose a mixed binomial Merton model. whereas the probability that a given number of companies default is given by:

enter image description here

whereas u is the common economic background factor and a standard normal random variable.

C(i) is given by:

enter image description here

whereas I substituted r with μ

Then I entered the following definitions:

V1 = 100; D1 = 70; μ1 = 0.2; σ1 = 0.1;
V2 = 110; D2 = 75; μ2 = 0.25; σ2 = 0.2;
V3 = 95; D3 = 72; μ3 = 0.1; σ3 = 0.2;
V4 = 120; D4 = 100; μ4 = 0.4; σ4 = 0.3;
V5 = 120; D5 = 110; μ5 = 0.2; σ5 = 0.2;
V6 = 100; D6 = 65; μ6 = 0.25; σ6 = 0.12;
V7 = 99; D7 = 85; μ7 = 0.3; σ7 = 0.05;
V8 = 125; D8 = 76; μ8 = 0.2; σ8 = 0.3;
V9 = 105; D9 = 100; μ9 = 0.23; σ9 = 0.02;
V10 = 100; D10 = 75; μ10 = 0.2 ; σ10 = 0.2;

T = 5; ρ = 0.5; \[ScriptU] = 0.25;

Evaluating Subscript[C, i]

C1 = (Log[V1/D1] + (μ1 - 0.5*σ1^2)*T)/(σ1*Sqrt[T]);

C2 = (Log[V2/D2] + (μ2 - 0.5*σ2^2)*T)/(σ2*Sqrt[T]);

C3 = (Log[V3/D3] + (μ3 - 0.5*σ3^2)*T)/(σ3*Sqrt[T]);

C4 = (Log[V4/D4] + (μ4 - 0.5*σ4^2)*T)/(σ4*Sqrt[T]);

C5 = (Log[V5/D5] + (μ5 - 0.5*σ5^2)*T)/(σ5*Sqrt[T]);

C6 = (Log[V6/D6] + (μ6 - 0.5*σ6^2)*T)/(σ6*Sqrt[T]);

C7 = (Log[V7/D7] + (μ7 - 0.5*σ7^2)*T)/(σ7*Sqrt[T]);

C8 = (Log[V8/D8] + (μ8 - 0.5*σ8^2)*T)/(σ8*Sqrt[T]);

C9 = (Log[V9/D9] + (μ9 - 0.5*σ9^2)*T)/(σ9*Sqrt[T]);

C10 = (Log[V10/D10] + (μ10 - 0.5*σ10^2)*T)/(σ10*Sqrt[T]);

Now i am not sure How to calculate the individual probabilities and add them up.

Would you suggest first to integrate the terms on the right hand side of the integral or is it possible to enter the whole function in Mathematica considering both the Normaldistribution and the Binomialdistribution?

I would be very thankful for any advise.

Thanks.

share|improve this question
    
Stocks defaulting?? I know the probability of bond defaults can be relatively high... but stocks?? I mean, you don't see a company bankruptcy everyday in the market... –  Rod Jul 11 '13 at 13:57
    
@ Rod Lm Sorry, let me substitute Stocks by assets so V represents the asset value of a company and D1 etc. the debt level. I will change this in the text above. thenks for reminding me :-) –  Milan Ivica Jul 11 '13 at 14:03
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.