Abstract

This work presents the neutrosophic Maxwell distribution (NMD) as a novel probability distribution. The proposed model represents a generalized design of Maxwell distribution that provides more analytical flexibility for data, including all imprecise observations or some degree of vagueness within the dataset. Important reliability characteristics and distributional properties of NMD are developed under the notion of neutrosophy. The neutrosophic forms of some commonly used functions in applied statistics such as mean, variance, moment generating function, and shape coefficients are explored. In view of uncertainties involved in the processing data and indeterminacy in the defined parameters, an estimation framework using the maximum likelihood approach is established. Additionally, the quantile function is developed to validate the distributional properties of NMD. The efficiency of the neutrosophic estimate has been studied through a Monte Carlo simulation. Finally, real data on the incubation period of COVID-19 are considered for numerical illustration, and further extensions of the NMD for future research works are discussed.

1. Introduction

Maxwell [1] derived a mathematical formulation of Maxwell distribution to illustrate the distribution of particles in thermal equilibrium. The unimodel and leptokurtic structure of the Maxwell curve is yielded due to the fact that all molecules do not travel with identical velocity, and some move faster than others [2]. A lot of fundamental characteristics of gases are explained by the Maxwell model, which impacts kinetic energy [3]. The Maxwell distribution is also known as the model of momenta, velocity, the degree of momenta, and energy of particles [4]. The Maxwell distribution is well-known for its applications in astronomy, chemistry, and engineering [5]. The failure rate of the Maxwell distribution rises steadily over time, which makes it effective in reliability and in life-testing experiments when the assumption of fixed failure rates, like that of an exponential distribution, is impractical [6]. Among many notable contributions that recommended this distribution in real-world testing investigations are [710]. For the first time, the Maxwell model was treated as a lifetime distribution and described the Bayesian estimators for the reliability function and the parameter in the study [11]. Chaturvedi and Rani [12] extended the Maxwell model by adding more parameters and provided Bayesian and classical estimators. The empirical Bayesian estimator using the Maxwell model is investigated by [13]. In these works, authors considered that known information about studied variables and distributional parameters are completely specified. The earlier studies devoted a mathematical treatment to the Maxwell model but do not throw light on neutrosophic aspects of the Maxwell probability function, especially in applied statistical research.

We encounter continuous variables in everyday life, and the observations are defined in exact numbers. There are hundreds of continuous univariate distributions found in the applied statistical applications [14]. Over the last several decades, various distributions have been extensively utilized to model data in a wide range of areas such as environmental sciences, actuarial and engineering, medical sciences, biological research, insurance, and economics. However, the common distributions do not suit well to every data in various circumstances [15]. On the other hand, measurements of any continuous variables always have some degree of imprecision [16]. It is self-evident in the scientific realm of measurements that terms like “equality” and “exact” should be forbidden since continuous processes cannot be measured precisely [17]. In addition to continuous variables, there are circumstances when accurate measurements are not feasible due to the irregular nature of studied variables [1821]. As an illustration, it is impossible to estimate a river’s depth precisely due to the wave nature of water. A similar problem arises when attempting to differentiate between a healthy and an unhealthy individual, a bad and good student, or between a normal and cool body temperature. This results in the single measurement being distorted. These considerations lead to the conclusion that real measurements do contain errors or imprecision in individual observations, which is referred to as fuzziness [22]. In a classical statistical framework, modeling variability among data is usually performed without accounting for fuzziness. Smarandache [23] proposed neutrosophic logic, which is an extension of fuzzy logic in order to tackle the uncertainties in studied variables. Smarandache [24] also provided the notions of neutrosophic statistics, which are generalizations of classical statistical procedures. When indeterminacy is zero, the neutrosophic statistical method converges to the classical approach. Valuable works that promoted the idea of neutrosophic statistics and undetermined statistics can be found in [2527]. Sherwani et al. [28] recently developed neutrosophic statistics in applied distributional theory. Based on a review of the literature and to the best of our knowledge, there is no study on the neutrosophic design of the Maxwell model for dealing with data that may be imprecisely defined.

In what follows, this study presents a statistical framework of the neutrosophic Maxwell distribution with applications, particularly in applied statistical research. Neutrosophic measures under the Maxwell model have been developed for handling imprecise data.

The work done in this study is divided into sections as follows. In Section 2, the proposed NMD, its genesis, some distributional properties, related theorems, and reliability characteristics are given. Section 3 deals with the maximum likelihood estimation under the neutrosophic environment. An extensive simulation work to validate the model authenticity is conducted in Section 4. In Section 5, the applicability of NMD has been illustrated by analyzing a real dataset on COVID-19 incubation period. Last, in Section 6, a brief discussion about the work done has been concluded.

2. Neutrosophic Maxwell Distribution

The neutrosophic random variable is followed by the NMD model with probability density function () aswhere , is an indicator function, is a nonnegative variable, and subscripts and , respectively, indicate the lower value and upper value of an indeterminate parameter. The neutrosophic density function () indicates the chance of events happening between two values under the sturdy curve. In the framework of neutrosophic calculus, it is defined as the integral of the variable density over a specified range. The neutrosophic parameter denotes simply the scale factor whose different values result in a variety of neutrosophic curves of the proposed distribution. The graph of for a continuous random variable with different neutrosophic parameters is shown in Figure 1.

Figure 1 shows that the densities are asymmetric and skewed toward the right. In the neutrosophic framework, the density curve is represented by a thick layer rather than a single curve. The layer thickness (shaded region) corresponds to an indeterminacy part and total area under the sturdy curve equal to one due to completeness of . Several other supporting properties of NMD may be constructed in the following versions of the theorems.

Theorem 1. The distribution function of NMD is .

Proof. Distribution function is conventionally defined asUsing the transformation and in equation (2) yieldedFurther simplification of equation (3) resultedwhere is the incomplete gamma function. The function has the counter domain , and the domain is a real line with nondecreasing and monotone behavior. Like , it is also represented by a thick curve, as shown in Figure 2.
Figure 2 shows the overall behavior of which is right continuous and varies in the interval .

Theorem 2. The neutrosophic mean of the NMD is .

Proof. is the required result.

Theorem 3. The variance of NMD is .

Proof. Definition of variance iswhereandThus, equation (6) becomes.
Now, .
After simplifying, we get

Theorem 4. The rth moment about the origin of NMD is .

Proof. The rth moment by definition is given by

Theorem 5. The skewness of the NMD is .
The skewness coefficient is defined by , whereandSo, equation (11) becomesThus, further simplification of equation (13) yielded

Theorem 6. The kurtosis of NMD is given by .

Proof. Definition kurtosis coefficient is given bywhereandThus, using equations (9) and (17) in equation (16) yieldedorIn this section, we also considered several to illustrate the theoretical results of the proposed model.

Example 1. If the time taken by the service station employee for wheels alignment is a random variable that follows NMD with parameter , determine the probability that the attendant would take more than 15 minutes to balance four tires.

Solution: . Let be a random variable that represents the number of times (in minutes) that a service station employee takes to balance tires; then, the required probability would beUsing equation (2), we can writeBy using , equation (20) can be further simplified as

Example 2. Research published in the November 1990 edition of Chemical Engineering Progress (CEP) magazine addressed the percentage purity of neon gas from a specific supplier. Assume that the percentage purity distribution was NMD with a mean of percent. What proportion of the purity values is expected to be between 97 and 99?

Solution: . Let be a random variable that denotes the percentage purity of neon gas.
Now, .
Using the result given in Theorem 2, we can findUsing Theorem 1, we can writeHence,

Example 3. Certain pollutant concentrations in parts per billion generated by the chemical plants have typically been shown to exhibit the neutrosophic Maxwell model with distribution parameter . Find the probability that the concentrations surpass 8 parts per billion?

Solution: . Let the random variable denotes the concentrations in parts per billion; then, the required probability can be obtained as

3. Neutrosophic Estimation Procedure

In this section, we discuss the analytical procedure for estimating the parameter of NMD using the most common method of maximum likelihood (ML). A new procedure using neutrosophic statistics for estimating the parameter of NMD is introduced. Let be the observed interval values sample from NMD with density function . Assuming the parameter is unknown in the defined distribution, then be the joint probability of the observed sample.

Taking the logarithm of the product provides log-likelihood as

The ML estimate of the unknown is the value that maximizes , i.e.,

The ML estimates, namely, can be obtained by using the neutrosophic calculus aswhere and

Simplification of equation (29) provides

Equating equation (31) to yieldshence proved.

Following the work of Sultan et al. [29], the effectiveness of the neutrosophic ML estimate is also evaluated for an unknown parameter of NMD. A Monte Carlo simulation is run with varying sample sizes and parameter settings to see how the results vary. There are simulation runs performed, each one using a sample size of and with true parameter equal to . In this simulated experiment, the following two metrics are used as performance indicators for each of our simulated data samples:where is the total number of runs in the simulation experiment.

For the sake of simulation, all computations are carried out using R software. The metrics and have been calculated for the abovementioned sample sizes and are given in Table 1.

We can see from the findings in Table 1 that the biases diminish as the sample size rises. It can infer from the results that the neutrosophic estimate provides more reliable results when the sample size is increased.

4. Simulation Study

In this section, simulated data of 100000 observations from the NMD have been generated with the inverse transformation approach. Random numbers from a given distribution are created using the inverse transformation technique by solving the following equation for at predetermined values of distributional parameter:where denotes the inverse of a particular distribution function, , and is the quantile point of the distribution.

In the case of NMD, equation (35) can be expressed as

By solving equation (35) for at a particular value of the parameter, say and 100000 random values from the uniform distribution resulted in 100000 different neutrosophic samples from NMD. The closed-form expression for equation (36) does not exist because of the intractable nature of . The function cannot be used explicitly for generating random numbers. Alternatively, the quantile points and other useful statistics can be estimated by using R as a programming language. The basic quantile points and other useful descriptive measures for different values of estimated over 100000 runs are given in Table 2.

Table 2 provides the useful statistics of NMD estimated over simulated data in interval forms due to the fact that the parameter is not precisely determined. Notice that values of skewness and kurtosis coefficients are invariant and crisp numbers because of the fact that they are do not involve indeterminate parameter in their expressions. Additionally, results of NMD coincide with the conventional model of the Maxwell distribution if upper and lower limits become the same, as given in the last column of Table 2.

5. Real Application

In this part, we have applied our suggested model to COVID-19 mean incubation time data obtained from the source [30]. Infectious illness incubation periods or the time between infection with a microbe and disease manifestation are directly significant to prevention and control. One of the most significant indicators of disease propagation and quarantine measures is the incubation period, which is the time period between infection and the start of symptoms [31]. Explicit models of the incubation time help us better understand how diseases spread. Infectious illnesses have an incubation time that may vary from a few hours, as is usual in the case of severe food poisoning, too many decades, as is the case with HIV and TB [32]. It is helpful in clinical practice for drawing rough estimates about the origins and sources of infection of specific patients and for establishing treatment methods that allow the incubation period to be prolonged. When dealing with a point source epidemic, being aware of the incubation period model allows for statistical assessment of exposure time [33]. Furthermore, identifying the incubation time distribution permits statistical assessment of the period of exposure during a point source epidemic, as well as a testing hypothesis to assess whether the pandemic has stopped. The incubation period distribution is important in statistical techniques for reconstructing epidemic curves and forecasting future illness incidence. The incubation period is easily determined from data showing the time of disclosure, such as experimental inoculation data and case travel records, since most infection occurrences are not clearly apparent. A coronavirus that causes severe acute respiratory illness is the cause of COVID-19 [34]. It has expanded to 223 nations and has been labeled a pandemic by the World Health Organization on March 11, 2020, after spreading rapidly. The average incubation time for COVID-19 is 2.87 days; however, the variation is as broad as 17.6 days. Existing studies have demonstrated that the distribution of the incubation period varies considerably, most likely due to the study population and estimation methods used. The data on incubation periods are not precisely measured. In general, the patient self-reports the beginning date of the symptom, which is deemed accurate but very rare. As a result, determining the incubation period without clear information about the moment of exposure is often challenging. Thus, instead of exact values, the mean incubation period globally from the source is provided with uncertainties instead of crisp numbers, as given in Table 2. Uncertainties in the mean incubation period are created to the strategy designed in [35].

To see Maxwell distribution is one of the reasonable models for describing the incubation period data, an informal graphical approach has been used. The visual fit of the Maxwell distribution on actual mean incubation period data is shown in Figure 3.

The subjective examination of extensively used plots, i.e., frequency distribution, probability plot, quantile plot, and CDF clearly indicates that the Maxwell model is more likely an adequate model for incubation period data. This is inferred from Figure 3 that most real data do not much deviate from the theoretical red lines, suggesting that Maxwell distribution is a reasonable model for the meantime incubation period. More naturally, the incubation period data involving uncertainties (Table 3) cannot be analyzed using the conventional Maxwell distribution. The statistical description of the incubation period data using the density of NMD is given in Table 4.

Table 4 provides the fitted measures of mean incubation period data in the form of intervals for mean and variance due to the estimated uncertainty parameter . Note that kurtosis and skewness coefficients are crisp values because analytical results as given in equations (13) and (16) do not depend on the estimated parameter value. Thus, the proposed model provides more useful information and efficiently analyzes data involving uncertainties.

6. Conclusions

A novel extension, so-called the neutrosophic Maxwell distribution (NMD) has been presented in this study. This new generalization is designed on the notions of neutrosophic calculus. Various distributional features such as mean, variance, distribution function, and shape coefficients have been explored in the neutrosophic framework. Numerical examples considered in this study demonstrated that theoretical results of NMD are elastic and appropriate to a wide variety of data. The quantile function has been developed to simulate numerical data and investigate the statistical properties of the proposed model. The estimation approach of maximum likelihood has been established and applied to simulated data. Results from the simulation study indicate that reliable estimates from the proposed model can be obtained with a larger sample size. The COVID-19 incubation time data have been utilized to describe how the proposed model can be implemented in practice. The incubation period for COVID-19 is not precisely measured and hence can easily be analyzed using NMD. It has been shown in the application section that the suggested model can analyze classical datasets and real-world data containing imprecision, vagueness, or uncertainties.

Future studies can be designed to improve the capability of the proposed model for handling indeterminate data in estimation approaches. In addition, this study would be helpful to design improved multivariate models for estimating mean incubation periods for other infectious diseases.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work under grant number (RGP.2/4/43), received by Mohammed M. Almazah (https://www.kku.edu.sa).