Where is the inverse CDF's implemented (Beta inverse cdf,normal inverse cdf etc. as menationed in [1]) . The list of functions(given at link [2]) under sub heading Statistics does not have these functions . Also i looked upon the functions there can be a lot more can be implemented . Matlab has huge number of functions .<br>
<br>One of such functions as Geometric CDF , i have planned for how to go about implenting it. Geometric CDF is the probability distribution of the number of bernoulli trials to get one success.<br><br><u>Calling Sequence</u><br>
<br>[P,Q] = cdfgeo("PQ",Xn,Pr,Ompr)<u><br><br>Arguments</u><br>P,Q,Xn,PR,Ompr<br> 5 real vecors of the same size<br><br>P,Q(Q=1-P)<br> The cumulation from 0 to Xn of the geometric distribution. <br> (Probablility
of getting one(first) success in Xn trials or less than Xn. Probability of success is Pr . ) Input range: [0,1].<br><br>Xn<br>The number of max trials .We will be looking at the geometric distribution for 1 to Xn trials and then sum them up to get geocdf.<br>
<br>Pr,Ompr<br>The probability of getting success in each trial .<br><br>In matlab it is implemented as <br><br>geocdf(3,0.5) where 3 is our Xn and Pr=Ompr=0.5<br><br>Hence now geocdf will be = (0.5^0)(0.5) + ((0.5)^1)(0.5) + ((0.5)^2)(0.5) = 0.5 + 0.25 + 0.125 <br>
Explanation of 3rd term -- to get first success, i used 3 trials(Xn) , therefore got failure in first two trials and success in the 3rd trial . <br>Explanation of 2nd Term -- to get first success , i used 2 trials (less than Xn) , therefore got failure in first trial and success in 2nd one...<br>
Explanation of 1st Term -- i got success in first trial only ... <br><br>Note -- Xn is the total number of max trials. For determining cdf we will be using fewer number tan Xn also and adding all them up .. <br><br><br>
<br>[1] <a href="http://forge.scilab.org/index.php/p/distfun/">http://forge.scilab.org/index.php/p/distfun/</a><br>[2] <a href="http://help.scilab.org/docs/5.4.0-alpha-1/en_US/index.html">http://help.scilab.org/docs/5.4.0-alpha-1/en_US/index.html</a><br>
<br><br><div class="gmail_quote">On Sat, Mar 24, 2012 at 9:50 PM, Sylvestre Ledru <span dir="ltr"><<a href="mailto:sylvestre.ledru@scilab-enterprises.com">sylvestre.ledru@scilab-enterprises.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">thanks for the forward, just don't forget me the next time ;)<br>
<br>
Le samedi 24 mars 2012 à 17:19 +0530, prateek papriwal a écrit :<br>
<div><div class="h5">><br>
><br>
> ---------- Forwarded message ----------<br>
> From: prateek papriwal <<a href="mailto:papriwalprateek@gmail.com">papriwalprateek@gmail.com</a>><br>
> Date: Sat, Mar 24, 2012 at 4:47 PM<br>
> Subject: Re: Accurate probability distribution functions (GSOC 2012)<br>
> To: <a href="mailto:michael.baudin@contrib.scilab.org">michael.baudin@contrib.scilab.org</a><br>
><br>
><br>
> ok .<br>
><br>
> yes . its true there are much more functions in matlab.<br>
> PLUS more function such as<br>
> binmial,geometric,hypergeometric,poisson,chebyshev's inequality,<br>
> markov's inequality, MGF of random variable , chi- square, bernoulli<br>
> etc.<br>
><br>
> Yes implementing these functions in scilab would be a good idea to<br>
> work upon .<br>
><br>
><br>
> On Sat, Mar 24, 2012 at 4:27 PM, <<a href="mailto:michael.baudin@contrib.scilab.org">michael.baudin@contrib.scilab.org</a>><br>
> wrote:<br>
> I think that you looked at it too quickly.<br>
> There is, by far, not all the required functions.<br>
> To see this more quickly, compare with :<br>
><br>
> <a href="http://www.mathworks.fr/help/toolbox/stats/bq_w_hm.html#bq_w_ie-8" target="_blank">http://www.mathworks.fr/help/toolbox/stats/bq_w_hm.html#bq_w_ie-8</a><br>
><br>
> There is only ~5% of the distributions in Matlab.<br>
> At least 30-40 more functions are necessary to make a complete<br>
> distribution toolbox:<br>
> just look at the "TODO" section in the readme.<br>
><br>
> As far as I am concerned, I have no interest in Queuing<br>
> theory.<br>
><br>
> Regards,<br>
><br>
> Michaël<br>
><br>
><br>
><br>
> On Mon, 19 Mar 2012 18:05:25 +0530, prateek papriwal<br>
> <<a href="mailto:papriwalprateek@gmail.com">papriwalprateek@gmail.com</a>> wrote:<br>
><br>
> @michael, thank you for ur reply .<br>
><br>
><br>
> <a href="http://forge.scilab.org/index.php/p/distfun/" target="_blank">http://forge.scilab.org/index.php/p/distfun/</a> [1] .. i<br>
> saw the<br>
><br>
> features(functions) already implemented .<br>
> i would like to have a view at the current state and<br>
> would be very<br>
> mjuch interested in implemeting concepts of Queueing<br>
> theory .<br>
><br>
><br>
> On Sun, Mar 18, 2012 at 7:24 AM, Sylvestre Ledru<br>
> wrote:<br>
><br>
> Le samedi 17 mars 2012 à 17:24 +0100,<br>
><br>
> <a href="mailto:michael.baudin@contrib.scilab.org">michael.baudin@contrib.scilab.org</a> [3]<br>
><br>
> a écrit :<br>
><br>
><br>
> > Secondly, i now work at EDF R&D, and no more in<br>
> the Scilab<br>
> > Consortium.<br>
> > How i could be involved in the GSOC 2012 has to be<br>
> clarified with<br>
> > the<br>
> > Consortium.<br>
><br>
> As long as I (or someone else in the consortium) am<br>
> in the loop, I<br>
> don't<br>
> have any issues with that.<br>
><br>
> S<br>
><br>
> --<br>
> -----------------------------<br>
> Sylvestre Ledru<br>
> Projects manager<br>
> Community manager<br>
> -----------------------------<br>
> Scilab Enterprises<br>
><br>
> <a href="http://www.scilab-enterprises.com/" target="_blank">http://www.scilab-enterprises.com/</a> [4]<br>
> <a href="http://www.scilab.org/" target="_blank">http://www.scilab.org/</a> [5]<br>
> -----------------------------<br>
><br>
><br>
><br>
> Links:<br>
> ------<br>
> [1] <a href="http://forge.scilab.org/index.php/p/distfun/" target="_blank">http://forge.scilab.org/index.php/p/distfun/</a><br>
> [2] mailto:<a href="mailto:sylvestre.ledru@scilab-enterprises.com">sylvestre.ledru@scilab-enterprises.com</a><br>
> [3] mailto:<a href="mailto:michael.baudin@contrib.scilab.org">michael.baudin@contrib.scilab.org</a><br>
> [4] <a href="http://www.scilab-enterprises.com/" target="_blank">http://www.scilab-enterprises.com/</a><br>
> [5] <a href="http://www.scilab.org/" target="_blank">http://www.scilab.org/</a><br>
><br>
><br>
><br>
<br>
</div></div>--<br>
-----------------------------<br>
Sylvestre Ledru<br>
Operation manager<br>
<div class="im">Community manager<br>
-----------------------------<br>
Scilab Enterprises<br>
<a href="http://www.scilab-enterprises.com/" target="_blank">http://www.scilab-enterprises.com/</a><br>
</div><a href="http://www.scilab.org/" target="_blank">http://www.scilab.org/</a><br>
-----------------------------<br>
<br>
</blockquote></div><br>