Right-truncating a lognormal distribution

1 vue (au cours des 30 derniers jours)
george pepper
george pepper le 18 Août 2019
Modifié(e) : Matt J le 18 Août 2019
Hello,
I would like to compute the probability that P(z<a) where
z=exp(x*beta+e) and e is distributed iid N(0,sigma^2).
I would like to evaluate the CDF of z, i.e.,
P(e< log(a)-x*beta),
for various values of x*beta. However, I know that z cannot take a value greater than a certain number, b. How can I obtain the truncated distribution to evaluate the above probability?
Many thanks in advance.
  5 commentaires
Matt J
Matt J le 18 Août 2019
Possibly you want a conditional CDF?
prob(z<a | z<b)
george pepper
george pepper le 18 Août 2019
I think this is right. How can I code this?

Connectez-vous pour commenter.

Réponse acceptée

Matt J
Matt J le 18 Août 2019
Modifié(e) : Matt J le 18 Août 2019
I think this is right. How can I code this?
function out=conditionalCDF(a,b,x,beta,sigma)
pa=normcdf(log(a)-x*beta,0,sigma);
pb=normcdf(log(b)-x*beta,0,sigma);
out=min(1, pa./pb);
end
  1 commentaire
george pepper
george pepper le 18 Août 2019
Thanks so much for your time, Matt! Have a great day.

Connectez-vous pour commenter.

Plus de réponses (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by