Received 5 December 2015; accepted 1 March 2016; published 4 March 2016
At present one of the major issues to know various traffic flows is in self similar nature to study and design some performance metric as that of Ethernet traffic etc. Until recently Poison approach has been used to model the road traffic irrespective of traffic intensity  . This was similar to the practice in the cases of Ethernet, LAN, WAN, and WWW traffic. But seminal studies   reveal that IP packet traffic in supposed networks tends to be bursty in nature on many time scales. This burstiness of traffic can be characterized mathematically as self- similar or long-range dependence (LRD). It is clear from the work agreed  that Poisson process could not emulate the self-similar network traffic. Markovian arrival process (MAP) emulating self-similar traffic is fitted over desired time scales by equating second-order statistics of the counts  - .
The idea of this paper is, we examine whether web users traffic data has the self similar property. This is to enhancement earlier results using a real time data  . This kind of research is useful for future studies to know the performance metrics and continuous improvement of web centers at various private and in public sector organizations. The rest of the paper has been organized as follows: Definition of self-similarity or long range dependence is given in Section II. Materials and methods are placed in Section III. In Section IV, Hurst parameter is computed using various methods. Finally, conclusions are given in Section V.
2. Long-Range Dependence and Self-Similarity
In this section we give a short description of the mathematical basis for second order self-similar processes (long-range dependence).
Exact Second-Order Self-Similar Process
The exact second-order self-similar process is defined as follows. Arrival instants are modeled as point process. Divide the time axis into disjoint intervals of unit length and let be the number of points (arrival) in the interval. Let X be a second order stationary process with variance and the autocorrelation function is given by
For each let a new time series is obtained averaging the original time series X over non-overlapping blocks of size m. That is
This new series, for each m, is also a second order stationary process with autocorrelation function.
Definition 1: The process “X” is said to be exactly second order self-similar with Hurst parameter
and variance if
Definition 2: The process “X” is said to be asymptotically second order self-similar with Hurst parameter
and variance if
In terms of variance, self-similar process is defined as follows:
Definition 3: The process “X” is said to be exactly second order self-similar with Hurst parameter
and variance if
Now we shall differentiate long range dependence (LRD) and short range dependence (SRD) processes. For from the Equation (2.3), we can see that as, and we have
The series is divergent if or otherwise they are convergent, being a posi-
tive term series. Accordingly the left hand series is divergent if or, otherwise
they are convergent. That is, for, the autocorrelation functions decays slowly, that is hyperbolically. In this case, the process x is called Long Range Dependent (LRD). The process x is Short Range Dependent (SRD) if and the autocorrelation function is summable (finite).
3. Materials and Methods
As discussed in the introduction, we are primarily interested collecting data from various sources. Real time web users data has been considered. The sample number of users logged on to an Internet server each minute over 100-minutes (see Appendix). In the study web users data can be treated as traffic and verify it is self-similar or not.
4. Methods for Estimating Hurst Parameter of Self-Similar Process
The intensity of self-similarity is given by Hurst parameter, H. The parameter H was named after the hydrologist H.E. Hurst who spent many years to investigate the problem of water storage and also to determine the level patterns of the Nile river. The parameter H has range. Estimation of H is a difficult task. Several methods are available to estimate degree of self-similarity in a time-series. We also present the three basic methods to calculate the Hurst parameter: Periodogram analysis, Correlogram method, R/S analysis, Variance-time analysis etc. Here is a method based on percentiles is applied and validated with the said methods.
4.1. Periodogram Analysis
In the frequency domain, analysis of time series is merely the analysis of a stationary process by means of its spectral representation. The periodogram  is given by
where is the Fourier frequency, n is the number of terms in the time series and is the data of the given series. To estimate h, first, one has to calculate this periodogram. Since is an good sample estimator of the spectral density, a series with long-range dependence should have a periodogram, which is proportional to close to the origin. Then a regression of the logarithm of the periodogram on the logarithm of the frequency λ should give a coefficient of. The slope of the fitted straight line is the estimate of. Using this method the H value is computed for the data given in section III. The obtained value of H in this case is 0.763.
4.2. Correlogram Method
In time series analysis  , plot of ACF (autocorrelation function) is known as correlogram where the estimated correlation can be given in terms of auto-covariance function
It has already been observed that slow decay of correlation, which is proportional to for
indicates the long-memory process. Therefore, the plot of the sample autocorrelation should exhibit this property. A much better plot for the handling of long-range dependence is the plot of ACF in logarithmic scale. If the asymptotic decay of the correlation is hyperbolic, then the points in the plot should be approximately scattered around a straight line with a negative slope of for the long memory processes but for short memory, the points should tend to diverse to minus infinity at an exponential rate. If the time series is long enough or if the series has strong long-range dependence, then this log-log correlogram is useful. Correlogram is useful as a preliminary heuristic approach to the data. Some pitfalls of sample correlation which are less known can be found in Mandelbrot  - . Even though it is neither widely used nor attractive method for estimation, still H, the self-similarity parameter, can be estimated by this method deriving an equation of the form
Using this method, the obtained value of H in this case is 0.79.
4.3. Percentile Method
In statistical methodologies, a percentile (or centile) is the value of a variable below which a certain percent of observations fall, like partition values of a process such as quartiles and deciles. There is no exact definition of percentile  , however all definitions yield similar results when the number of observations is very large. One definition of percentile, often given in texts, is that the percentile of N ordered values is obtained by first calculating the rank.
Given data set or time series. First we can find the percentiles for a given time series using
percentile, this a special type of average such as partition values in descriptive statistics like quartiles. Draw a scattered Plot percentile number against percentiles on log scales. A linear equation (say) is obtained with the slope. The Hurst parameter is then computed by
Using this method, the H value is computed for the data. The pertaining scattered data and trend line with the slope are depicted in Figure 1. The obtained value of H in this case is 0.762. One relevant paper  , which explained how the 95-percentile depends on the aggregation window size, and how this phenomenon justifies the mathematical definition of self similarity. The advantages of this method are: This method is matter of a simple empirical formula, unlike other two methods. Data however large it may be is divided into 100 parts (partition values) and the plotting involves only 100 points (percentile versus percentile number).
Figure 1. Percentile versus percentile number.
4.4. Math Lab Implementation Code is Linked with the Percentile Method
5. Some Conclusions
In this paper, real time web user’s data has been considered as traffic from various web centers and it has been proved to be self-similar. Various methods to test the self-similarity have been used. The obtained values of Hurst parameter h are reasonably close to each other. This kind of research is useful for future studies to know the performance metrics at web centers.
The number of users logged on to an Internet server each minute over 100-minutes.
 Perati, M.R., Raghavendra, K., Koppula, H.K.R., Doodipala, M.R. and Dasari, R. (2012) Self-Similar Behavior of Highway Road Traffic and Performance Analysis at Toll Plazas. Journal of Transportation Engineering, 138, 1233-1238.
 Leland, W.E., Taqqu, M.S., Willinger, W. and Wilson, D.V. (1994) On the Self-Similar Nature of Ethernet Traffic (Extended Version). IEEE/ACM Transactions on Networking, 2, 1-15.
 Anderson, A.T. and Nielsen, B.F. (1998) A Markovian Approach for Modeling Packet Traffic with Long Range Dependence. IEEE Journal on Selected Areas in Communications, 16, 719-732.
 Yoshihara, T., Kasahara, S. and Takahashi, Y. (2001) Practical Time-Scale Fitting of Self-Similar Traffic with Markov-Modulated Poisson Process. Telecommunication Systems, 17, 185-211.
 Shao, S.K., Perati, M.R., Tsai, M.G., Tsao, H.W. and Wu, J. (2005) Generalized Variance-Based Markovian Fitting for Self-Similar Traffic Modeling. IEICE Transactions on Communications, E88-B, 4659-4663.
 Doodipala, M.R. (2013) Second Order Statistics of Time Series of Various Real Time Problems in Conjunction with Periodogram Technique. International Journal of Latest Trends in Engineering and Technology (IJLTET), 3.
 Sarker, M.M.A. (2007) Estimation of the Self-Similarity Parameter in Long Memory Processes. Journal of Mechanical Engineering, ME38, 32-37. Transaction of the Mech. Eng. Div., The Institution of Engineers, Bangladesh.