# Noncentral beta distribution

This section contains functions for working with noncentral beta distribution. They allow to calculate density, probability, quantiles and to generate pseudo-random numbers distributed according to the corresponding law. The noncentral beta distribution is defined by the following formula:

where:

• x — value of the random variable
• a — the first parameter of beta distribution
• b — the second parameter of beta distribution
• λ — noncentrality parameter

In addition to the calculation of the individual random variables, the library also implements the ability to work with arrays of random variables.

 Function Description MathProbabilityDensityNoncentralBeta Calculates the probability density function of the noncentral beta distribution MathCumulativeDistributionNoncentralBeta Calculates the value of the noncentral beta probability distribution function MathQuantileNoncentralBeta Calculates the value of the inverse noncentral beta distribution function for the specified probability MathRandomNoncentralBeta Generates a pseudorandom variable/array of pseudorandom variables distributed according to the noncentral beta distribution law MathMomentsNoncentralBeta Calculates the theoretical numerical values of the first 4 moments of the noncentral beta distribution

Example:

 #include  #include  #include  #property script_show_inputs //--- input parameters input double a_par=2;    // the first parameter of beta distribution (shape1) input double b_par=5;    // the second parameter of beta distribution (shape2) input double l_par=1;    // noncentrality parameter (lambda) //+------------------------------------------------------------------+ //| Script program start function                                    | //+------------------------------------------------------------------+ void OnStart()   { //--- hide the price chart    ChartSetInteger(0,CHART_SHOW,false); //--- initialize the random number generator      MathSrand(GetTickCount()); //--- generate a sample of the random variable    long chart=0;    string name="GraphicNormal";    int n=1000000;       // the number of values in the sample    int ncells=53;       // the number of intervals in the histogram    double x[];          // centers of the histogram intervals    double y[];          // the number of values from the sample falling within the interval    double data[];       // sample of random values    double max,min;      // the maximum and minimum values in the sample //--- obtain a sample from the noncentral beta distribution    MathRandomNoncentralBeta(a_par,b_par,l_par,n,data); //--- calculate the data to plot the histogram    CalculateHistogramArray(data,x,y,max,min,ncells); //--- obtain the sequence boundaries and the step for plotting the theoretical curve    double step;    GetMaxMinStepValues(max,min,step);    step=MathMin(step,(max-min)/ncells); //--- obtain the theoretically calculated data at the interval of [min,max]    double x2[];    double y2[];    MathSequence(min,max,step,x2);    MathProbabilityDensityNoncentralBeta(x2,a_par,b_par,l_par,false,y2); //--- set the scale    double theor_max=y2[ArrayMaximum(y2)];    double sample_max=y[ArrayMaximum(y)];    double k=sample_max/theor_max;    for(int i=0; i=cells) ind=cells-1;       frequency[ind]++;      }    return (true);   } //+------------------------------------------------------------------+ //|  Calculates values for sequence generation                       | //+------------------------------------------------------------------+ void GetMaxMinStepValues(double &maxv,double &minv,double &stepv)   { //--- calculate the absolute range of the sequence to obtain the precision of normalization    double range=MathAbs(maxv-minv);    int degree=(int)MathRound(MathLog10(range)); //--- normalize the maximum and minimum values to the specified precision    maxv=NormalizeDouble(maxv,degree);    minv=NormalizeDouble(minv,degree); //--- sequence generation step is also set based on the specified precision    stepv=NormalizeDouble(MathPow(10,-degree),degree);    if((maxv-minv)/stepv<10)       stepv/=10.;   }