RE: Triangular Probability Distributions
- To: mathgroup at smc.vnet.net
- Subject: [mg30020] RE: [mg30009] Triangular Probability Distributions
- From: "tgarza01 at prodigy.net.mx" <tgarza01 at prodigy.net.mx>
- Date: Sat, 21 Jul 2001 16:16:45 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
Hello Michael, I'm not clear as to what you mean by "methodology". Still, what you may do is define the triangular probability density function and work from it. For example, if you want your density defined in the interval [a, b], In[1]:= Clear[triPdf,triCdf] In[2]:= triPdf[x_,a_,b_]:=((1+a)/(b^2-a^2))*4*(x-a)/;a<=x<=(b+a)/2; triPdf[x_,a_,b_]:=((1+a)/(b^2-a^2))*4*(b-x)/;(b+a)/2<x<=b; triPdf[x_,a_,b_]:=0/;a>x||x<b; In[3]:= Plot[triPdf[x,0,1],{x,0,1}]; Here you obtained the graph of the triangular density in [0,1]. The distribution function is then defined as In[4]:= triCdf[x_,a_,b_]:=Integrate[triPdf[y,a,b],{y,a,x}] which you plot with (it takes a little while, due to the fact that the numerical integration is slow because of the peak at x = 0.5; you may integrate one part after the other and then it runs very quickly): In[5]:= Off[NIntegrate::"ncvb"];Off[NIntegrate::"slwcon"]; In[6]:= Plot[triCdf[x,0,1],{x,0,1}]; I turned off the messages to avoid looking at them. The k-th moment is Integrate[x^k*triPdf[x,a,b],{x,a,b}] for any k. Tomas Garza Original Message: ----------------- From: loopm at yahoo.com (Michael Loop) To: mathgroup at smc.vnet.net Subject: [mg30020] [mg30009] Triangular Probability Distributions I have been looking for a methodology for using the triangular probability distribution in Mathematica. I have not found anything that allows me to do this. Has anyone found a way to use the triangular distribution? Are there any add-on packages that would include this distribution? Thank you, Michael Loop Minneapolis MN -------------------------------------------------------------------- Mail2Web - Check your email from the web at http://www.mail2web.com/ .