Next Article in Journal
Wiener Complexity versus the Eccentric Complexity
Next Article in Special Issue
Microscopically Reversible Pathways with Memory
Previous Article in Journal
Multiple Comparisons for Exponential Median Lifetimes with the Control Based on Doubly Censored Samples
Previous Article in Special Issue
Study of the Vibrational Predissociation of the NeBr2 Complex by Computational Simulation Using the Trajectory Surface Hopping Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monte Carlo Simulation of a Modified Chi Distribution with Unequal Variances in the Generating Gaussians. A Discrete Methodology to Study Collective Response Times

by
Juan Carlos Castro-Palacio
1,†,
J. M. Isidro
1,
Esperanza Navarro-Pardo
2,*,
Luisberis Velázquez-Abad
3 and
Pedro Fernández-de-Córdoba
1
1
Grupo de Modelización Interdisciplinar, Instituto Universitario de Matemática Pura y Aplicada, InterTech, Universitat Politècnica de València, E-46022 Valencia, Spain
2
Grupo de Modelización Interdisciplinar, Departamento de Psicología Evolutiva y de la Educación, InterTech, Universitat de València, E-46010 Valencia, Spain
3
Departamento de Física, Universidad Católica del Norte, Antofagasta 0610, Chile
*
Author to whom correspondence should be addressed.
Current affiliation: Department of Electrical Engineering, Electronics, Automation, and Applied Physics, Technical University of Madrid. Ronda de Valencia, 3, 28012 Madrid, Spain.
Mathematics 2021, 9(1), 77; https://doi.org/10.3390/math9010077
Submission received: 28 October 2020 / Revised: 26 December 2020 / Accepted: 28 December 2020 / Published: 31 December 2020

Abstract

:
The Chi distribution is a continuous probability distribution of a random variable obtained from the positive square root of the sum of k squared variables, each coming from a standard Normal distribution (mean = 0 and variance = 1). The variable k indicates the degrees of freedom. The usual expression for the Chi distribution can be generalised to include a parameter which is the variance (which can take any value) of the generating Gaussians. For instance, for k = 3, we have the case of the Maxwell-Boltzmann (MB) distribution of the particle velocities in the Ideal Gas model of Physics. In this work, we analyse the case of unequal variances in the generating Gaussians whose distribution we will still represent approximately in terms of a Chi distribution. We perform a Monte Carlo simulation to generate a random variable which is obtained from the positive square root of the sum of k squared variables, but this time coming from non-standard Normal distributions, where the variances can take any positive value. Then, we determine the boundaries of what to expect when we start from a set of unequal variances in the generating Gaussians. In the second part of the article, we present a discrete model to calculate the parameter of the Chi distribution in an approximate way for this case (unequal variances). We also comment on the application of this simple discrete model to calculate the parameter of the MB distribution (Chi of k = 3) when it is used to represent the reaction times to visual stimuli of a collective of individuals in the framework of a Physics inspired model we have published in a previous work.

1. Introduction

The Chi distribution is a widely known continuous probability distribution in Statistics. It represents the distribution of a random variable which is calculated as the Euclidean norm over k variables each following a standard Normal distribution [1], where k specifies the degrees of freedom. This is the case of two familiar examples in Physics, namely, the Rayleigh distribution (Chi of k = 2 degrees of freedom) and the Maxwell-Boltzmann distribution (MB) (Chi of k = 3 degrees of freedom). The first appears when the wind velocity is analysed in two dimensions each normally distributed with equal variance, and zero mean. The second represents the distribution of the particle velocities in the Ideal Gas model of Physics, an ideal system of independent particles in thermal equilibrium with a thermostat [2].
In this work, we first generalise the expressions of the Chi distribution for the case of equal (and different from 1) variances in the generating Gaussians. We also analyse the case of unequal variances in the generating Gaussians whose distribution can be still approximated to a Chi distribution. For this purpose, we carry out Monte Carlo simulations of the generating Gaussians with unequal variances. Then, by means of non-linear fittings, we could obtain the parameter of the Chi distribution for a given k, which stands for the degrees of freedom. For this, we perform several simulations starting from different sets of unequal variances in the generating Gaussians.
In the second part of the article, we present a discrete model to calculate the parameter of the Chi distribution in an approximate way for the case of unequal variances. We also comment on the application of this simple discrete model to calculate the parameter of the MB distribution when it is used to represent the reaction times to visual stimuli of a collective of individuals. Reaction time (RT) is a key variable on many psychological experiments as, for instance, it correlates with cognitive disorders [3]. RTs are commonly represented in Experimental Psychology by means of the ex-Gaussian probability distribution [3,4], which provides three interpretable parameters, denoted as μ , σ , and τ . These parameters represent the mean and variance ( μ , σ ) of the Gaussian function, and the decay constant ( τ ) of the exponential function, which are convoluted to give rise to the ex-Gaussian function. RT distributions are asymmetric functions [4,5], which cannot be fully described by the mean and standard deviations [6]. The ex-Gaussian parameters obtained from fitting reaction time distributions to visual stimuli are interpretable in relation to cognitive disorders as certain correlations are evidenced to exist [7,8,9,10,11,12,13]. This is the case of Attention Deficit and Hyperactivity Disorder (ADHD), the most usual diagnosed cognitive disorder in school-aged children.
In our previous work [14,15], we integrated the response times of a group of individuals of similar age in a Maxwell-Boltzmann (MB) distribution in order to model the behaviour of a collective of individuals. Relevant correlation has been found between the RTs series to visual stimuli of individuals who perform experiments in a short time scale (hundreds of milliseconds) when cultural elements are not expected to have influence and without exchanging information. Similarly, relevant frequencies are found when the response times along the sequence of stimuli are analysed by Fast Fourier Transform [16,17]. These facts support the existence of a collective behaviour rather than individuals acting separately. The MB is a very simple distribution depending on one parameter only, which is proportional to the variance. One major advantage of having a probability distribution of a collective is that it allows for a rigorous and more natural way of classifying the individuals in relation to the reference group they belong to and without need for an external reference. With our discrete model, the MB distribution parameter (case of a Chi distribution of k = 3 degrees of freedom and unequal variances in the generating Gaussians) can be calculated in a simple and practical way. The results of this work are particularly relevant for the interpretation of reaction times in different contexts.
The outline of the article is as follows. In Section 2, the general aspects of the Chi distribution and its generalization for variances taking any positive value in the generating Gaussians are commented. The latter includes the case of the MB distribution (Chi of k = 3) of the Ideal Gas particles in Physics. In Section 3, we carry on developing Monte Carlo simulations of an approximate variant of the Chi distribution, namely, when the generating Gaussians have unequal variances. In Section 4, a discrete model is proposed to calculate the parameter, B, of the MB distribution when used to represent response time data of a collective of individuals.

2. The χ Distribution and Its Generalization for Variances Different from 1

The Chi distribution ( χ ) is a continuous probability distribution of a random variable obtained as:
χ =   i = 1 k X i 2 ,
where X i are random variables each following a standard Normal distribution (mean = 0 and variance = 1). The resulting probability density is expressed as:
f ( x ; k ) = 2 1 k 2   [ Γ ( k 2 ) ] 1   x k 1   exp ( x 2 2   ) .
If the variances of the generating Gaussians can take any positive value, f ( x ; k ) can be generalised as [6]:
f ( x ; B ,   k ) = 2 1 k 2   B k 2   [ Γ ( k 2 ) ] 1 x k 1   exp ( x 2 2   B ) ,
where k is the number of degrees of freedom of the distribution, and B a parameter related to its variance. The corresponding cumulative distribution function (CDF) is expressed as:
F ( x ; B ,   k ) =   1 [ Γ ( k 2 ) ] 1 Γ ( k 2 , x 2 2   B )   ,
where Γ ( k 2 , x 2 2 B ) is the upper incomplete gamma function, and Γ ( k 2 ) , the gamma function. The χ distribution as written in Equation (1) describes a multidimensional ideal gas of particles (where the k variances of the Gaussian components are equal among them).

Case of the Maxwell-Boltzmann Distribution (Chi of k = 3)

In Physics, the Maxwell-Boltzmann (MB) distribution represents the probability density of the moduli of the velocities ( | v i | ) of the independent particles forming the ideal gas. This situation is given when the components of the velocities follow Gaussian distributions with the same variances and centred around zero. The MB distribution is just a Chi distribution of k = 3, where the variance in the velocity distributions is proportional to the temperature:
f M B ( | v i | ) = 2 π B 3 / 2 | v i | 2 exp ( | v i | 2 2 B ) .
In Physics, Equation (3) would represent the case of a multidimensional Ideal Gas.

3. Monte Carlo Simulations of the Case of Unequal Variances in the Generating Gaussians

As follows, we will perform Monte Carlo simulations to generate a random variable which obtained the positive square root of the sum of k squared variables each coming from non-standard Normal distributions, where the variances can take any positive value.
First, we will simulate the case of 3 generating Gaussians with zero mean and standard deviations σ 1 ,   σ 2 , and σ 3 . Then, we will evaluate the parameter B by fitting the generated data to the MB function in Equation (5) by using the non-linear fitting algorithm of Levenberg-Marquardt [18,19]. The quality of the fit is indicated by the resulting coefficient of determination R 2 .
In Table 1, we present the details of the numerical experiments, where e σ   ( % ) = 2 ×   100 %   ×   ( σ min   σ max ) / ( σ min + σ max ) is the percentage difference between the minimum and maximum variances. The parameter B fit resulting from the fit and the corresponding coefficient of determination of the fitting, R 2 are also included. The table includes the values of B calc   = [ 1 3 ( σ 1 + σ 2 + σ 3 ) ] 2 , which is proven to be a very good approximation of B fit ; e B   ( % ) is the percentage difference between B fit and B calc ; e B is the mean value over the three sets of variances taken for the same percentage difference ( e σ   ( % ) ), and SD the corresponding standard deviation. Similarly, in Table 2 we present the details of the numerical experiments and the relative errors between both evaluations of B for k = 5.
In Figure 1, the relative error between the fitted and the calculated parameters is represented as a function of the percentage difference between the minimum and maximum variances. The error bars in the figure indicate the standard deviation of this error when different possibilities are chosen for the standard deviations within the range σ min <   σ <   σ max . Results are shown for 3 and 5 Gaussian components (Table 1 and Table 2). Figure 2 depicts the same results, but for the coefficient of determination R 2 .

4. Discrete Methodology for the Calculation of the MB Parameter When Used to Represent the Reaction Times of a Collective of Individuals

In previous works [14,15], we have developed a methodology to analyse the response times to visual stimuli of a group of individuals which is based on the use of the MB distribution [2] instead of the ex-Gaussian distribution [3]. The MB distribution emerges when the reaction time distributions of single individuals are combined at a higher level of hierarchy, which allows us to capture the collective behaviour in a single distribution [14,15]. The details of the computer-based experiments and for the design of the visual stimuli can be found in our previous works [4,14,15]. Although the methodology for the analysis of the reaction time data can be also found therein, we will recall in what follows the more relevant expressions.
From the first three moments (mean: M i , standard deviation: S i , and skewness: λ i ) of the reaction time distribution of each participant (denoted by the index i), the following vector v i can be defined:
v i = ( M i * ,   S i * ,   λ i * ) = ( M i M M , S i S S , λ i λ λ ) ,
where M , S , and λ are the mean values of the moments over the sample. Then, the Euclidean norm of the vector, v i , is calculated as:
| v i | =   ( M i * ) 2 + ( S i * ) 2 + (   λ i * ) 2 .
Each participant can now be represented by a single scalar, | v i | , which we have proven [1] to be distributed as a MB distribution. The following expression will be used to fit the experimental data:
Q ( | v i | ) = 2 π B 3 / 2 | v i | 2 exp ( | v i | 2 2 B ) ,
where B is the only free parameter (a measure of the dispersion of the distribution). This expression gives the probability of each individual’s time response in relation to the reference group the individual belongs to and, therefore, it allows for a classification methodology of individuals based on their reaction times to visual stimuli [14]. A classification methodology, based on a model for a collective of individuals which is characterised by a given probability distribution, allows for the classification of the behaviour of individuals in relation to the group they belong to. There is no need for an external reference.
As follows, we introduce a discrete model to evaluate the parameter B , which is consistent with the model developed in [14] and suitable for non-equal variances of the distributions of the three components of the vector v i . The methodology to estimate B develops as follows:
  • Step 1. Calculation of the moments: M i , mean; S i 2 , the variance; and λ i , the skewness for each participant of index i:
M i = 1 N T j = 1 N T t i j ,
S i 2 = 1 N T j = 1 N T ( t i j M i ) 2 ,
λ i = 1 S i 2 1 N T j = 1 N T ( t i j M i ) 3 ,
where t i j are the response times ( j = 1 ,   , N T ) of each participant ( i = 1 ,   , N P ). M i ,   S i , and λ i have time units.
  • Step 2. Evaluation of the mean values of the moments:
M = 1 N P i = 1 N P M i ;   S = 1 N P i = 1 N P S i ;   λ = 1 N P i = 1 N P λ i ,
and calculation of the mean values of the squared moments:
M 2 = 1 N S i = 1 N P M i 2 ;   S 2 = 1 N S i = 1 N P S i 2   ;   λ 2 = 1 N S i = 1 N P λ i 2 .
  • Step 3. Evaluation of the variances of the dimensionless and normalised moments.
First, we proceed to the make the moments in Equations (9)–(11) dimensionless and normalised:
M i * = M i M M   ;   S i * = S i S S   ;   λ i * = λ i λ λ .
Secondly, we evaluate the mean values of the squared dimensionless and normalised moments (D.N.M.), i.e.,
M * 2 = 1 N P i = 1 N P M i * 2 ;   S * 2 = 1 N P i = 1 N P S i * 2 ;   λ * 2 = 1 N P i = 1 N P λ i * 2 .
If we use the definitions (12) and (13), we can write:
M * 2 = 1 M 2 ( M 2 M 2 )   ;   S * 2 = 1 S 2 ( S 2 S 2 )     λ * 2 = 1 λ 2 ( λ 2 λ 2 ) .
Obviously, the mean values of the D.N.M. are zero:
M * = 1 N p i = 1 N P M i * =   0   ;   S * = 1 N p i = 1 N P S i * =   0   ;   λ * = 1 N p i = 1 N P λ i * =   0 .
Finally, we calculate the variances of the D.N.M. as:
σ M * 2 = 1 N p i = 1 N P ( M i *   M * ) 2 = 1 N p i = 1 N P M i * 2 =   M * 2     ;
σ S * 2 = 1 N p i = 1 N P ( S i *   S * ) 2 = 1 N p i = 1 N P S i * 2 = S * 2   ;  
σ λ * 2 = 1 N p i = 1 N P ( λ i *   λ * ) 2 = 1 N p i = 1 N P λ i * 2   =   λ * 2 .
Using (16), the variances and standard deviations (SD) can be expressed as:
σ M * 2 = 1 M 2 ( M 2 M 2 )   ;   σ S * 2 = 1 S 2 ( S 2 S 2 )     σ λ * 2 = 1 λ 2 ( λ 2 λ 2 ) .
σ M * = 1 M M 2 M 2   ;   σ S * = 1 S S 2 S 2   ;   σ λ * = 1 λ λ 2 λ 2   .    
  • Step 4. Estimation of B for the MB distribution.
The participants are then represented by a single scalar | v i | distributed as a MB distribution. We evaluate its variance as:
σ 2 = 1 N p i = 1 N P   | v i | 2     μ 2 ,
where μ stands for the mean value of the distribution of | v i | . Using (7) we can write:
σ 2 = 1 N p i = 1 N P   [ ( M i * ) 2 + ( S i * ) 2 + (   λ i * ) 2 ]     μ 2 .
It is well known that the MB distribution corresponds to the distribution of the positive square root of the sum of the squares of three independent Gaussian distributed random variables, each one with zero mean and equal standard deviations ( σ M * 2 =   σ S * 2 =   σ λ * 2   A ) . The mean value of the MB distribution (3) is given by μ =   8   B π and the variance by   σ 2 = ( 3   8 π ) B . Therefore,
B =   1 3 1 N p i = 1 N P   [ ( M i * ) 2 + ( S i * ) 2 + (   λ i * ) 2 ]   =   1 3   ( σ M * 2 +   σ S * 2 +   σ λ * 2 )   =   A .
Even for the general case where the three variances of the components are not necessarily equal among them, we can still model the response time data by a MB distribution. We propose to take the parameter B from the average value of the three standard deviations (instead of taking B directly from the average value of the three variances):
B = [ 1 3 ( σ M * + σ S * + σ λ * ) ] 2 =   1 9 ( σ M * 2 + σ S * 2 + σ λ * 2 + 2   σ M * σ S * + 2   σ M * σ λ * + 2   σ S * σ λ * ) .
In Figure 3, we summarise the main steps of the proposed discrete methodology to evaluate the parameter B.
We can extend the previous procedure to the case of the distribution of the positive square roots of the sum of squares of a set of k independent random variables, each one following a Gaussian distribution with variance σ i 2   ( i = 1 , 2 , , k ). In this case, we can generalise Equation (26) as:
B = [ 1 k i = 1 k   σ i ] 2 .

5. Conclusions

Monte Carlo simulations and non-linear fittings have been carried out to study the Chi distribution when it is approximated by the square root of the sum of squared variables following non-standard Normal distributions with different values in the variance. We establish the boundaries of what to expect when we start from a set of unequal variances in the generating Gaussians. We also present a discrete model to calculate the parameter, B , of the Chi distribution in an approximate way for the case of unequal variances in the generating Gaussians. We illustrate the application of this simple discrete model to a problem within Experimental Psychology, that is, the calculation of the parameter of the MB distribution when it is used to represent the reaction times to visual stimuli of a collective of individuals in the framework of the Physics inspired model we have published in a previous work. We address two cases that involve the Chi distribution, namely, the case when the variances in the generating Gaussians are equal but different from 1 (multidimensional ideal gas in Physics), and the case when the variances are unequal (multidimensional and anisotropic ideal gas). The results from this work can find multiple applications in social sciences, where to know the behaviour of the collective, and the behaviour of an individual relative to the collective it belongs to, are of enormous importance. For instance, the probability distribution of the collective, represented in terms of a Maxwell-Boltzmann-like distribution, provides a more natural classification methodology of individuals.

Author Contributions

Conceptualization, J.C.C.-P., J.M.I., E.N.-P., L.V.-A. and P.F.-d.-C.; data curation, J.C.C.-P., J.M.I., E.N.-P., L.V.-A. and P.F.-d.-C.; formal analysis, J.C.C.-P., J.M.I., E.N.-P., L.V.-A. and P.F.-d.-C.; investigation, J.C.C.-P., J.M.I., E.N.-P., L.V.-A. and P.F.-d.-C.; methodology, J.C.C.-P., J.M.I., E.N.-P., L.V.-A. and P.F.-d.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by grant no. RTI2018-102256-B-I00 (Spain).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

It is a pleasure to thank F. López-Ramón for technical discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Blitzstein, J.K.; Hwang, J. Introduction to Probability, 1st ed.; Chapman & Hall/CRC Texts in Statistical Science Boca: Ratón, FL, USA, 2015. [Google Scholar]
  2. Tolman, R. The Principles of Statistical Mechanics; Dover Publications Inc.: New York, NY, USA, 2003. [Google Scholar]
  3. Hernaiz-Guijarro, M.; Castro-Palacio, J.C.; Navarro-Pardo, E.; Isidro, J.M.; Fernández de Córdoba, P. A Probabilistic Classification Procedure Based on Response Time Analysis towards a Quick Pre–diagnosis of Student’s Attention Deficit. Mathematics 2019, 7, 473. [Google Scholar] [CrossRef] [Green Version]
  4. Moret-Tatay, C.; Gamermann, D.; Navarro-Pardo, E.; Fernández de Córdoba, P. ExGUtils: A Python package for statistical analysis with the ex-Gaussian probability density. Front. Psychol. 2018, 9, 1–11. [Google Scholar]
  5. Luce, R.D. Response Times: Their Role in Inferring Elementary Mental Organization; Oxford University Press: New York, NY, USA, 1986. [Google Scholar]
  6. Ratcliff, R.; Murdock, B.B. Retrieval processes in recognition memory. Psychol. Rev. 1976, 83, 190–214. [Google Scholar] [CrossRef]
  7. Gmehlin, D.; Fuermaier, A.B.M.; Walther, S.; Debelak, R.; Rentrop, M.; Westermann, C.; Sharma, A.; Tucha, L.; Koerts, J.; Tucha, O.; et al. Intraindividual variability in inhibitory function in adults with ADHD. An ex-Gaussian approach. PLoS ONE 2014, 9, e112298. [Google Scholar] [CrossRef] [PubMed]
  8. Adamo, N.; Hodsoll, J.; Asherson, P.; Buitelaar, J.K.; Kuntsi, J. Ex-Gaussian, frequency and reward analyses reveal specificity of reaction time fluctuations to ADHD and not autism traits. J. Abnorm. Child. Psychol. 2019, 47, 557–567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Moret-Tatay, C.; Leth-Steensen, C.; Quarti Irigaray, T.; Argimon, I.I.L.; Gamermann, D.; Abad-Tortosa, D.; Oliveira, C.; Saíz-Mauleo’n, B.; Vázquez-Martínez, A.; Navarro-Pardo, E.; et al. The effect of corrective feedback on performance in basic cognitive tasks: An analysis of RT components. Psychol. Belg. 2016, 56, 370–381. [Google Scholar] [CrossRef]
  10. Mira-Iglesias, A.; Navarro-Pardo, E.; Conejero, J.A. Power-law distribution of natural visibility graphs from reaction times series. Symmetry 2019, 11, 563. [Google Scholar] [CrossRef] [Green Version]
  11. Moreno-Cid, A.; Moret-Tatay, C.; Quarti Irigaray, T.; Argimon, I.I.L.; Murphy, M.; Szczerbinski, M.; Martínez-Rubio, D.; Beneyto-Arrojo, M.J.; Navarro-Pardo, E.; Fernández de Córdoba, P. The role of age and emotional valence in word recognition: An ex-Gaussian analysis. Stud. Psychol. 2015, 57, 83–94. [Google Scholar] [CrossRef]
  12. Moret-Tatay, C.; Moreno-Cid, A.; Argimon, I.I.L.; Quarti Irigaray, T.; Szczerbinski, M.; Murphy, M.; Vázquez-Martínez, A.; Vázquez-Molina, J.; Saíz-Mauleón, B.; Navarro-Pardo, E.; et al. The effects of age and emotional valence on recognition memory: An ex-Gaussian components analysis. Scand. J. Psychol. 2014, 55, 420–426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Navarro-Pardo, E.; Navarro-Prados, A.B.; Gamermann, D.; Moret-Tatay, C. Differences between young and old university students on a lexical decision task: Evidence through an ex-gaussian approach. J. Gen. Psychol. 2013, 140, 251–268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Castro-Palacio, J.C.; Fernández de Córdoba, P.; Isidro, J.M.; Navarro-Pardo, E. Brain reaction times: Linking Individual and Collective Behaviour through Physics modelling. arXiv 2019, arXiv:1910.12725. [Google Scholar]
  15. Castro-Palacio, J.C.; Fernández de Córdoba, P.; Isidro, J.M.; Navarro-Pardo, E. Percentile study of Chi distribution. Mathematics 2020, 8, 514. [Google Scholar] [CrossRef] [Green Version]
  16. Iglesias-Martínez, M.E.; Castro-Palacio, J.C.; Ferna´ndez-de-Córdoba, P.; Isidro, J.M.; Navarro-Pardo, E. Spectral Analysis of the Reaction Time Dynamics over Consecutive Visual Stimuli. Mathematics 2020, 8, 1979. [Google Scholar] [CrossRef]
  17. Xavier Castellanos, F.; Sonuga-Barke, E.J.S.; Scheres, A.; Di Martino, A.; Hyde, C.; Walters, J.R. Varieties of Attention-Deficit/Hyperactivity Disorder- Related Intra-Individual Variability. Biol. Psychiatry 2005, 57, 1416–1423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Levenberg, K. A Method for the Solution of Certain Non-Linear Problems in Least Squares. Q. Appl. Math. 1944, 2, 164–168. [Google Scholar] [CrossRef] [Green Version]
  19. Marquardt, D. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  20. Ho, T.S.; Rabitz, H. A general method for constructing multidimensional molecular potential energy surfaces from ab initio calculations. J. Chem. Phys. 1996, 104, 2584–2597. [Google Scholar] [CrossRef]
  21. Unke, O.T.; Castro-Palacio, J.C.; Bemish, R.; Meuwly, M. Collision-induced rotational excitation in N+2 (2S+g,v = 0)–Ar: Comparison of computations and experiment. J. Chem. Phys. 2016, 144, 224307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Relative error in the calculation of B as a function of the dispersion of the standard deviations of the Gaussian components (for k = 3, 5). The points in the plot have been connected by means of reproducing kernel Hilbert space (RKHS) interpolations [20,21].
Figure 1. Relative error in the calculation of B as a function of the dispersion of the standard deviations of the Gaussian components (for k = 3, 5). The points in the plot have been connected by means of reproducing kernel Hilbert space (RKHS) interpolations [20,21].
Mathematics 09 00077 g001
Figure 2. Coefficient of determination R 2 as a function of the dispersion of the SD of the Gaussian components (k = 3, 5). The points in the plot have been connected by means of RKHS interpolations [9,10].
Figure 2. Coefficient of determination R 2 as a function of the dispersion of the SD of the Gaussian components (k = 3, 5). The points in the plot have been connected by means of RKHS interpolations [9,10].
Mathematics 09 00077 g002
Figure 3. Chart summarising the main steps of the discrete methodology for the calculation of the parameter B .
Figure 3. Chart summarising the main steps of the discrete methodology for the calculation of the parameter B .
Mathematics 09 00077 g003
Table 1. Numerical experiments to estimate the accuracy in the calculation of the parameter B. We have generated 106 Gaussian-distributed random numbers for each of the three Gaussian distributions.
Table 1. Numerical experiments to estimate the accuracy in the calculation of the parameter B. We have generated 106 Gaussian-distributed random numbers for each of the three Gaussian distributions.
C σ 1 σ 2 σ 3 e σ   ( % ) R 2 B fit B calc e B   ( % ) e B SD
10.260.500.7496.00.9740.2580.2503.32.443.10
20.260.380.7496.00.9610.2130.2120.6
30.260.620.7496.00.9750.3080.2925.4
40.300.500.7080.00.9860.2540.2501.61.481.50
50.300.400.7080.00.9800.2180.2180.0
60.300.600.7080.00.9860.2930.2843.0
70.320.500.6872.00.9900.2530.2501.11.011.10
80.320.410.6872.00.9870.2210.2210.1
90.320.590.6872.00.9900.2870.2812.1
100.400.500.6040.00.9990.2500.2500.10.140.20
110.400.450.6040.00.9990.2340.2340.0
120.400.550.6040.00.9990.2680.2670.3
130.450.500.5520.01.0000.2500.2500.00.040.00
140.450.480.5520.01.0000.2420.2420.0
150.450.530.5520.01.0000.2590.2580.1
Table 2. Numerical experiments to estimate the accuracy in the evaluation of the parameter B (k = 5). We have generated 106 Gaussian-distributed random numbers for each of the five Gaussian distributions.
Table 2. Numerical experiments to estimate the accuracy in the evaluation of the parameter B (k = 5). We have generated 106 Gaussian-distributed random numbers for each of the five Gaussian distributions.
C σ 1 σ 2 σ 3 σ 4 σ 5 e σ   ( % ) R 2 B fit B calc e B   ( % ) e B SD
10.280.400.510.630.7490.20.9830.2710.2604.21.714.70
20.280.300.510.690.7490.20.9680.2710.2546.6
30.280.400.510.520.7490.20.9850.2480.2403.3
40.310.410.510.600.7077.20.9900.2630.2552.90.742.90
50.310.350.510.630.7077.20.9850.2580.2493.7
60.310.480.510.520.7077.20.9940.2590.2532.2
70.340.430.510.600.6866.70.9940.2650.2602.00.512.00
80.340.370.510.610.6866.70.9910.2590.2522.5
90.340.480.510.520.6866.70.9960.2600.2561.5
100.410.460.510.550.6037.60.9990.2570.2550.60.090.70
110.410.430.510.550.6037.60.9990.2510.2490.6
120.410.420.510.570.6037.60.9990.2530.2510.8
130.420.450.490.520.5526.81.0000.2360.2350.20.050.30
140.420.450.490.540.5526.81.0000.2400.2390.3
150.420.470.490.510.5526.81.0000.2380.2370.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Castro-Palacio, J.C.; Isidro, J.M.; Navarro-Pardo, E.; Velázquez-Abad, L.; Fernández-de-Córdoba, P. Monte Carlo Simulation of a Modified Chi Distribution with Unequal Variances in the Generating Gaussians. A Discrete Methodology to Study Collective Response Times. Mathematics 2021, 9, 77. https://doi.org/10.3390/math9010077

AMA Style

Castro-Palacio JC, Isidro JM, Navarro-Pardo E, Velázquez-Abad L, Fernández-de-Córdoba P. Monte Carlo Simulation of a Modified Chi Distribution with Unequal Variances in the Generating Gaussians. A Discrete Methodology to Study Collective Response Times. Mathematics. 2021; 9(1):77. https://doi.org/10.3390/math9010077

Chicago/Turabian Style

Castro-Palacio, Juan Carlos, J. M. Isidro, Esperanza Navarro-Pardo, Luisberis Velázquez-Abad, and Pedro Fernández-de-Córdoba. 2021. "Monte Carlo Simulation of a Modified Chi Distribution with Unequal Variances in the Generating Gaussians. A Discrete Methodology to Study Collective Response Times" Mathematics 9, no. 1: 77. https://doi.org/10.3390/math9010077

APA Style

Castro-Palacio, J. C., Isidro, J. M., Navarro-Pardo, E., Velázquez-Abad, L., & Fernández-de-Córdoba, P. (2021). Monte Carlo Simulation of a Modified Chi Distribution with Unequal Variances in the Generating Gaussians. A Discrete Methodology to Study Collective Response Times. Mathematics, 9(1), 77. https://doi.org/10.3390/math9010077

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop