Next Article in Journal
Pairs of Positive Solutions for a Carrier p(x)-Laplacian Type Equation
Next Article in Special Issue
Mathematical Formalization and Applications to Data with Excess of Zeros and Ones of the Unit-Proportional Hazard Inflated Models
Previous Article in Journal
Explicit Parameterizations of Ortho-Symplectic Matrices in R4
Previous Article in Special Issue
On Stochastic Representations of the Zero–One-Inflated Poisson Lindley Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Review of Wrapped Distributions for Circular Data

Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(16), 2440; https://doi.org/10.3390/math12162440
Submission received: 4 July 2024 / Revised: 30 July 2024 / Accepted: 4 August 2024 / Published: 6 August 2024
(This article belongs to the Special Issue New Advances in Distribution Theory and Its Applications)

Abstract

:
The wrapped method is the most widely used method for constructing distributions for circular data. In this paper, we provide a review of all known wrapped distributions, including 45 distributions for continuous circular data and 10 distributions for discrete circular data. For each wrapped distribution, we state its nth trigonometric moment, mean direction, mean resultant length, skewness, and kurtosis. We also discuss data applications and limitations of each wrapped distribution. This review could be a useful reference and encourage the development of more wrapped distributions. We also mention an R package available for fitting all of the reviewed distributions and illustrate its applications.

1. Introduction

Circular data, also known as directional data, refer to measurements where the values are cyclical and repeat over a defined interval. This kind of data arises naturally in situations where the end of the scale reconnects with the beginning, such as angles, time, and compass directions. For example, angles measured in degrees or radians are circular because 0 is the same as 360, and similarly, 0 radians equates to 2 π radians. This cyclical nature poses unique challenges for statistical analysis because traditional linear methods are not suitable for data that wrap around.
Practical examples of circular data abound in various fields. In meteorology, wind directions recorded over time create circular datasets, as the direction can be anywhere between 0 and 360. In biology, the study of animal movement, such as the migratory patterns of birds or the rotational behavior of certain animals, often involves circular data. Time-related data are another common example; for instance, the times of day at which certain events occur (like sleep cycles or peak traffic hours) are inherently circular because they repeat every 24 h. Analyzing such data requires careful consideration of its circular nature to avoid misinterpretation and to uncover meaningful patterns and insights.
There have been many distributions proposed for circular data. Refs. [1,2,3,4,5,6] provide excellent reviews. Most of the distributions for circular data are based on the method of wrapping. Hence, we feel it is appropriate to provide a review of all known wrapped distributions to date, which is the aim of this paper. No such review is known to date. Such a review could be a useful reference for those interested in both the theory and applications of wrapped distributions. It could also enhance the development of more wrapped distributions.
The method of wrapping can be described as follows. Suppose g is a valid probability density function (PDF). If g is defined on the positive real line, then the PDF of the wrapped distribution is
f ( θ ) = k = 0 g θ + 2 k π
for 0 θ < 2 π . If g is defined on the entire real line, then the PDF of the wrapped distribution is
f ( θ ) = k = g θ + 2 k π
for 0 θ < 2 π .
On the other hand, now suppose g is a valid probability mass function (PMF). If g is defined on the positive real line, then the PMF of the wrapped distribution is
Pr Θ = 2 π r m = k = 0 g r + k m
for r = 0 , 1 , , m 1 , and m 1 . If g is defined on the entire real line, then the PMF of the wrapped distribution is
Pr Θ = 2 π r m = k = g r + k m
for r = 0 , 1 , , m 1 , and m 1 .
The most important properties of a circular random variable, say Θ , are its nth trigonometric moment denoted by m n , mean direction denoted by μ , mean resultant length denoted by ρ , skewness denoted by γ 1 , and kurtosis denoted by γ 2 , defined by
m n = E cos ( n Θ ) + i E sin ( n Θ ) ,
μ = arcsin E sin ( Θ ) E cos ( Θ ) 2 + E sin ( Θ ) 2 ,
ρ = E cos ( Θ ) 2 + E sin ( Θ ) 2 ,
γ 1 = exp ( 2 μ ) 1 ρ 3 2 E sin ( 2 Θ ) ,
and
γ 1 = exp ( 2 μ ) 1 ρ 2 E cos ( 2 Θ ) ρ 4 1 ρ 2 ,
respectively. For each reviewed distribution, we list expressions for (5)–(9) as well as expressions for the cumulative distribution function (CDF) corresponding to (1)–(4), if they are available. We also discuss data applications as well as the limitations of each reviewed distribution. Section 2 lists all known expressions corresponding to (1)–(2). Section 3 lists all known expressions corresponding to (3)–(4). The properties (5)–(9) not listed in Section 2 and Section 3 could provide readers with avenues for future work.
The expressions listed in Section 2 and Section 3 involve several special functions, including the Lerch function defined by
Φ ( z ; s , a ) = k = 0 z k ( a + k ) s ;
the gamma function defined by
Γ ( a ) = 0 t a 1 e t d t ;
the lower incomplete gamma function defined by
γ ( a , x ) = 0 x t a 1 e t d t ;
the upper incomplete gamma function defined by
Γ ( a , x ) = x t a 1 e t d t ;
the standard normal probability density function defined by
ϕ ( x ) = 1 2 π e x 2 2 ;
the standard normal cumulative distribution function defined by
Φ ( x ) = x 1 2 π e t 2 2 d t ;
and the modified Bessel function of the second kind defined by
K ν ( x ) = π csc ( π ν ) 2 I ν ( x ) I ν ( x ) , if   ν Z , lim μ ν K μ ( x ) , if   ν Z ;
the modified Bessel function of the first kind of order ν defined by
I ν ( x ) = k = 0 1 Γ ( k + ν + 1 ) k ! x 2 2 k + ν .
The properties of these special functions can be found in [7,8]. We also set i = 1 throughout.
An R software package (version 4.4.1) for fitting of all of the reviewed distributions has been created by the second author [9]. Three data applications illustrating the R package are given in Section 4.

2. A Review of Continuous Wrapped Distributions

In this section, we review wrapped Akash, wrapped Aradhana, wrapped binormal, wrapped Birnbaum–Saunders, wrapped Cauchy, wrapped chi-square, wrapped exponential, wrapped exponentiated inverted Weibull, wrapped gamma, wrapped generalized geometric stable, wrapped generalized Gompertz, wrapped generalized normal Laplace, wrapped generalized skew normal [10], wrapped generalized skew normal [11], wrapped half-logistic, wrapped half-normal, wrapped [12]’s skew Laplace, wrapped hypoexponential, wrapped Ishita, wrapped Laplace, wrapped length-biased weighted exponential, wrapped Levy, wrapped Lindley, wrapped Linnik, wrapped Lomax, wrapped modified Lindley, wrapped new Weibull–Pareto, wrapped normal, wrapped Pareto, wrapped quasi-Lindley, wrapped Rama, wrapped Richard, wrapped Shanker, wrapped skew Laplace, wrapped skew normal, wrapped stable, wrapped Student’s t, wrapped transmuted exponential, wrapped two-parameter Lindley, wrapped two-sided Lindley, wrapped variance gamma, wrapped weighted exponential, wrapped Weibull, wrapped XGamma and wrapped XLindley distributions.

2.1. Wrapped Akash Distribution

Ref. [13] took g to be the PDF of the Akash distribution to obtain the wrapped Akash distribution. Its PDF and CDF are
f ( θ ) = λ 3 e λ θ 1 + θ 2 1 e 2 π λ 2 + 4 π e 2 π λ θ + π + π e 2 π λ θ e 2 π λ λ 2 + 2 1 e 2 π λ 3
and
F ( θ ) = 4 π λ e 2 π λ 1 ( 1 + λ θ ) e λ θ 1 e 2 π λ + π λ 1 e λ θ 1 + e 2 π λ λ 2 + 2 1 e 2 π λ 3 + 1 e λ θ λ 2 + 2 λ θ 2 + λ θ e λ θ λ 2 + 2 1 e 2 π λ ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is
m n = λ 3 λ 2 + 2 λ 2 2 i n λ n 2 + 2 λ i n 3
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 3 arctan 1 λ arctan 2 λ λ 2 + 1 ,
ρ = λ 3 λ 2 + 2 λ 2 + 1 2 + 4 λ 2 λ 2 + 1 2 ,
γ 1 = λ 3 λ 2 + 2 λ 2 2 2 + 16 λ 2 λ 2 + 4 3 sin κ λ , 2 1 λ 3 λ 2 + 2 λ 2 + 1 2 + 4 λ 2 λ 2 + 1 3 3 2 ,
and
γ 2 = λ 3 λ 2 + 2 λ 2 2 2 + 16 λ 2 λ 2 + 4 3 cos ( κ λ , 2 ) λ 3 λ 2 + 2 λ 2 + 1 2 + 4 λ 2 λ 2 + 1 3 4 1 λ 3 λ 2 + 2 λ 2 + 1 2 + 4 λ 2 λ 2 + 1 3 2 ,
respectively, where κ λ , 2 = 3 arctan 2 λ arctan 4 λ λ 2 2 6 arctan 1 λ + 2 arctan 2 λ λ 2 + 1 . The wrapped Akash distribution was used on a dataset from [1] regarding the long-axis orientations of 60 feldspar laths in basalt rock. It demonstrated a better fit to the data compared to the wrapped exponential distribution and the wrapped Lindley distribution. Despite having only one parameter, the wrapped Akash distribution proved to be a flexible model. Additionally, it allows for closed form expressions for both its PDF and its CDF.

2.2. Wrapped Aradhana Distribution

Ref. [14] took g to be the PDF of the Aradhana distribution to obtain the wrapped Aradhana distribution. Its PDF and CDF are
f ( θ ) = λ 3 e λ θ λ 2 + 2 λ + 2 ( 1 + θ ) 2 1 e 2 π λ + 4 ( θ + 1 ) π e 2 π λ 1 e 2 π λ 2 + 4 π 2 e 2 π λ 1 + e 2 π λ 1 e 2 π λ 3
and
F ( θ ) = 1 e λ θ 1 e 2 π λ λ θ λ θ + 2 λ + 2 e λ θ λ 2 + 2 λ + 2 1 e 2 π λ + 4 π λ e 2 π λ ( 1 + λ ) 1 e λ θ λ θ e λ θ λ 2 + 2 λ + 2 1 e 2 π λ 2 + 4 π 2 λ 2 1 e λ θ 1 + e 2 π λ e 2 π λ λ 2 + 2 λ + 2 1 e 2 π λ 3 ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is given by
m n = λ 3 λ 2 + 2 λ n 2 + 2 2 i n ( 1 + λ ) λ 2 + 2 λ i n 3
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 3 arctan 1 λ arctan 2 ( 1 + λ ) λ 2 + 2 λ + 1 ,
ρ = λ 3 λ 2 + 2 λ + 2 λ 2 + 2 λ + 1 2 + 4 1 + λ 2 λ 2 + 1 3 ,
γ 1 = λ 3 λ 2 + 2 λ + 2 λ 2 + 2 λ 2 2 + 16 1 + λ 2 λ 2 + 4 3 sin μ 2 2 μ 1 λ 2 + 2 λ + 1 2 + 4 1 + λ 2 λ 2 + 1 3 3 2 ,
and
γ 2 = λ 3 λ 2 + 2 λ + 2 λ 2 + 2 λ 2 2 + 16 1 + λ 2 λ 2 + 4 3 cos μ 2 2 μ λ 3 λ 2 + 2 λ + 2 λ 2 + 2 λ 2 2 + 16 1 + λ 2 λ 2 + 4 3 4 1 λ 2 + 2 λ + 1 2 + 4 1 + λ 2 λ 2 + 1 3 2 ,
respectively, where μ 2 = 3 arctan 2 λ arctan 4 ( 1 + λ ) λ 2 + 2 λ 2 . The wrapped Aradhana distribution has only one parameter, which might restrict its versatility and relevance to real-world datasets. However, it does have closed form expressions for the PDF and CDF.

2.3. Wrapped Binormal Distribution

Ref. [15] took g to be the PDF of the binormal distribution to obtain the wrapped binormal distribution. Its PDF is
f ( θ ) = k = 2 π 1 σ 1 + σ 2 e θ + 2 k π μ 2 2 σ 1 2 , θ μ , k = 2 π 1 σ 1 + σ 2 e θ + 2 k π μ 2 2 σ 2 2 , θ > μ
for 0 θ < 2 π , θ > μ , and σ > 0 . The nth trigonometric moment is
m n = 2 b n cos ( p μ ) c n sin ( p μ ) + 2 i b n sin ( p μ ) + c n cos ( p μ ) ,
where
b p = 1 2 σ 1 σ 1 + σ 2 e p 2 σ 1 2 2 + σ 2 σ 1 + σ 2 e p 2 σ 2 2 2
and
c p = 1 π σ 1 σ 1 + σ 2 e p 2 σ 1 2 2 n = 1 p σ 1 2 2 n 1 ( 2 n 1 ) ( n 1 ) ! + 2 π σ 2 σ 1 + σ 2 e p 2 σ 2 2 2 n = 1 p σ 2 2 2 n 1 ( 2 n 1 ) ( n 1 ) ! .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan b 1 sin ( μ ) + c 1 cos ( μ ) b 1 cos ( μ ) c 1 sin ( μ ) ,
ρ = 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 ,
γ 1 = b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) 2 + b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) 2 sin μ 2 2 μ 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 3 2 ,
and
γ 2 = 2 b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) 2 + b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) 2 sin μ 2 2 μ 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 2 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 4 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 2 ,
respectively, where μ 2 = arctan b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) . The wrapped binormal distribution lacks practicality due to its non-closed-form PDF.

2.4. Wrapped Birnbaum–Saunders Distribution

Ref. [16] took g to be the PDF of the Birnbaum-Saunders distribution to obtain the wrapped Birnbaum–Saunders distribution. Its PDF and CDF are
f ( θ ) = e δ 2 δ + 1 4 π μ k = 0 θ + 2 k π + δ μ δ + 1 θ + 2 k π 3 2 e δ 4 θ + 2 k π ( δ + 1 ) δ μ + δ μ θ + 2 k π ( δ + 1 )
and
F ( θ ) = k = 0 Φ δ 2 θ + 2 k π ( δ + 1 ) δ μ δ μ θ + 2 k π ( δ + 1 ) k = 0 Φ δ 2 2 k π ( δ + 1 ) δ μ δ μ 2 k π ( δ + 1 ) ,
respectively, for 0 θ < 2 π , 0 μ < 2 π , and δ > 0 . The nth trigonometric moment is
m n = 1 2 1 + δ + 1 δ + 1 4 n i μ e δ δ + 1 δ + 1 4 n i μ 2 δ + 1
for n = 1 , 2 , . The wrapped Birnbaum–Saunders distribution was used on a dataset examined by [1] and originally gathered by [17]. The dataset includes 100 ant directions in reaction to a uniformly lit black target. The wrapped Birnbaum–Saunders distribution fit the data better than both symmetric and asymmetric von Mises distributions.

2.5. Wrapped Cauchy Distribution

Ref. [18] took g to be the PDF of the Cauchy distribution to obtain the wrapped Cauchy distribution. Its PDF is
f ( θ ) = k = γ μ ( γ 2 + ( θ μ + 2 π n ) 2 ) = 1 2 π sinh γ cosh γ cos ( θ μ )
for 0 θ < 2 π , γ > 0 , and < μ < . The nth trigonometric moment is
m n = e i n μ n γ
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ ,
ρ = e γ ,
γ 1 = 0 ,
and
γ 2 = e 2 γ e 8 γ 1 e γ 2 ,
respectively. Notably, the wrapped Cauchy distribution has a closed-form expression for its PDF. However, it is a symmetric distribution and, as such, has limited applicability to skewed datasets.

2.6. Wrapped Chi-Square Distribution

Ref. [19] took g to be the PDF of a chi-square distribution to obtain the wrapped chi-square distribution. Its PDF is
f ( θ ) = e θ 2 π k 2 1 2 Γ k 2 Φ e π ; 1 k 2 , θ 2 π
for 0 θ < 2 π and k > 0 . The nth trigonometric moment is
m n = 1 2 i n k 2
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = k 2 arctan 2 ,
ρ = 5 k 4 ,
γ 1 = 17 k 4 sin k 2 arctan ( 4 ) k arctan ( 2 ) 1 5 k 4 3 2 ,
and
γ 2 = 17 k 4 cos k 2 arctan ( 4 ) k arctan ( 2 ) 5 k 1 5 k 4 2 ,
respectively. The wrapped chi-square distribution faces limitations due to its non-closed-form PDF, as defined in terms of the Lerch function.

2.7. Wrapped Exponential Distribution

Ref. [20] took g to be the PDF of the exponential distribution to obtain the wrapped exponential distribution. Its PDF and CDF are
f ( θ ) = λ e λ θ 1 e 2 π λ
and
F ( θ ) = 1 e λ θ 1 e 2 π λ ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is
m n = λ λ i n
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan 1 λ ,
ρ = λ 1 + λ 2 ,
γ 1 = 2 λ 1 + λ 2 1 4 4 + λ 2 1 + λ 2 λ 3 2 ,
and
γ 2 = 3 λ 2 1 + λ 2 4 + λ 2 1 + λ 2 λ 2 ,
respectively. The wrapped exponential distribution adequately fits various real-life distributions, as reported by [1]. It offers straightforward closed-form expressions for the PDF and CDF. Nevertheless, its limitation lies in being a one-parameter distribution, restricting its flexibility.

2.8. Wrapped Exponentiated Inverted Weibull Distribution

Ref. [21] took g to be the PDF of the exponentiated inverted Weibull distribution to obtain the wrapped exponential inverted Weibull distribution. Its PDF is
f ( θ ) = c k = 0 λ ( θ + 2 π k ) ( c + 1 ) e ( θ + 2 π k ) c λ
for 0 θ < 2 π , c > 0 , and λ > 0 . The nth trigonometric moment is
m n = k = 0 i t λ 1 c k k ! Γ 1 k c = b t + i c t
for n = 1 , 2 , , where
b n = k = 0 ( 1 ) k n λ 1 c 2 k ( 2 k ) ! Γ 1 2 k c
and
c n = k = 0 ( 1 ) k n λ 1 c 2 k + 1 ( 2 k + 1 ) ! Γ 1 2 k + 1 c .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan c 1 b 1 ,
ρ = b 1 2 + c 1 2 ,
γ 1 = b 2 2 + c 2 2 sin μ 2 2 μ 1 b 1 2 + c 1 2 3 2 ,
and
γ 2 = b 2 2 + c 2 2 cos μ 2 2 μ b 1 2 + c 1 2 4 1 b 1 2 + c 1 2 3 2 ,
respectively, where μ 2 = arctan c 2 b 2 . The wrapped exponential inverted Weibull distribution was used to analyze the orientation data of 76 turtles post egg-laying, as documented by [1]. This distribution exhibited a superior fit compared to the wrapped new Weibull–Pareto distribution. Nonetheless, a limitation is noted due to the absence of a closed-form PDF.

2.9. Wrapped Gamma Distribution

Ref. [22] took g to be the PDF of the gamma distribution to obtain the wrapped gamma distribution. Its PDF is
f ( θ ) = λ r Γ ( r ) e λ θ ( 2 π ) r 1 Φ e 2 λ π ; 1 r , θ 2 π
for 0 θ < 2 π , r > 0 , and λ > 0 . The nth trigonometric moment is
m n = λ r λ i n r
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = r arctan 1 λ ,
ρ = λ r λ 2 + 1 r 2 ,
γ 1 = λ r λ 2 + 4 r 2 sin μ 2 2 μ 1 λ r λ 2 + 1 r 2 3 2 ,
and
γ 2 = λ r λ 2 + 4 r 2 sin μ 2 2 μ λ r λ 2 + 1 r 2 4 1 λ r λ 2 + 1 r 2 2 ,
respectively, where μ 2 = r arctan 2 λ . The wrapped gamma distribution is constrained by its non-closed-form PDF, expressed in terms of the Lerch function.

2.10. Wrapped Generalized Geometric Stable Distribution

Ref. [23] took g to be the PDF of the generalized geometric stable distribution to obtain the wrapped generalized geometric stable distribution. Its PDF and CDF are
f ( θ ) = 1 2 π 1 + 2 k = 1 α k cos ( k θ ) + β k sin ( k θ )
and
F ( θ ) = 1 2 π 1 + 2 k = 1 α k k sin ( k θ ) + β k k β k k cos ( k θ ) ,
respectively, for 0 θ < 2 π , where
α k = 1 + σ α k α 2 + σ α k α β tan π α 2 + μ * k 2 λ 2 cos λ arctan σ α k α β tan π α 2 + μ * k 1 + σ α k α , if   α 1 , 1 + σ α k α 2 + μ * k 2 σ α k α β π log k 2 λ 2 cos λ arctan μ * k 2 σ α k α β π log k 1 + σ α k α , if   α = 1
and
β k = 1 + σ α k α 2 + σ α k α β tan π α 2 + μ * k 2 λ 2 sin λ arctan σ α k α β tan π α 2 + μ * k 1 + σ α k α , if   α 1 , 1 + σ α k α 2 + μ * k 2 σ α k α β π log k 2 λ 2 sin λ arctan μ * k 2 σ α k α β π log k 1 + σ α k α , if   α = 1
for 0 < α 2 , λ > 0 , σ > 0 , 1 β 1 , μ < , and μ * = μ mod 2 π . The nth trigonometric moment is
m n = α n + i β n
for n = 1 , 2 , . The mean direction and mean resultant length are
μ = λ arctan μ * + σ α β tan π α 2 1 + σ α   and   2 π , if   α 1 , λ arctan μ * 1 + σ α   and   2 π , if   α = 1
and
ρ = 1 + σ α 2 + μ * σ α β tan π α 2 2 λ 2 , if   α 1 , 1 + σ 2 + μ * 2 λ 2 , if   α = 1 ,
respectively. The wrapped generalized geometric stable distribution was used to analyze hourly wind direction data collected over three days at a location on Black Mountain, ACT, Australia, as reported by [24] and discussed by [1]. Although the distribution was found to fit the data effectively, its drawback lies in its high number of parameters. Simplified distributions with fewer parameters might prove more practical.

2.11. Wrapped Generalized Gompertz Distribution

Ref. [25] took g to be the PDF of the generalized Gompertz distribution to obtain the wrapped generalized Gompertz distribution. Its PDF and CDF are
f ( θ ) = 1 b Γ ( c ) k = exp c θ + 2 π k a b e θ + 2 π k a b
and
F ( θ ) = k = 1 Γ ( c ) Γ c , e 2 π k a b Γ c , e θ + 2 π k a b ,
respectively, for 0 θ < 2 π , 0 a < 2 π , b > 0 , c > 0 , and θ + 2 π k > a . The nth trigonometric moment is
m n = e i n a Γ i n b + c Γ ( c )
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = a ,
ρ = Γ ( i b + c ) Γ ( c ) ,
γ 1 = 0 ,
and
γ 2 = Γ ( 2 i b + c ) Γ ( c ) Γ ( i b + c ) Γ ( c ) 4 1 Γ ( i b + c ) Γ ( c ) 2 ,
respectively. The wrapped generalized Gompertz distribution was used to analyze the orientations of 50 noisy scrub bird nests along a creek bank, as reported by [1]. While this distribution was found to fit the data effectively, it involves several parameters. Thus, there might be alternative distributions with fewer parameters that offer a better fit to the data.

2.12. Wrapped Generalized Normal Laplace Distribution

Ref. [26] took g to be the PDF of the generalized normal Laplace distribution to obtain the wrapped generalized normal Laplace distribution. Its PDF is
f ( θ ) = 1 2 π 1 + 2 k = 1 α k cos ( k θ ) + β k sin ( k θ )
for 0 θ < 2 π , where
α k = e τ 2 k 2 1 + a 2 k 2 1 + b 2 k 2 ζ 2 cos arctan 1 + a b k 2 cos ( η k ) + ( b a ) k sin ( η k ) 1 + a b k 2 sin ( η k ) + ( b a ) k cos ( η k )
and
β k = e τ 2 k 2 1 + a 2 k 2 1 + b 2 k 2 ζ 2 sin arctan 1 + a b k 2 cos ( η k ) + ( b a ) k sin ( η k ) 1 + a b k 2 sin ( η k ) + ( b a ) k cos ( η k )
for < η < , < τ < , a > 0 , b > 0 , and ζ > 0 . The nth trigonometric moment is
m n = e i η n τ 2 n 2 2 1 i a n 1 i b n ζ
for n = 1 , 2 , . The mean direction and mean resultant length are
μ = ζ η + arctan a b 1 + a b   mod   2 π
and
ρ = e τ 2 1 + a 2 1 + b 2 ζ 2 ,
respectively. The wrapped generalized normal Laplace distribution was used to analyze a dataset containing 1827 flight headings of migrating birds, as documented by [27]. It was found that this distribution offers a superior fit compared to a five-parameter mixture of two von Mises distributions but that it is not as effective as a four-parameter mixture incorporating circular uniform and skew normal components.

2.13. Wrapped Generalized Skew Normal Distribution [10]

Ref. [10] took g to be the PDF of the generalized skew normal distribution to obtain the wrapped generalized skew normal distribution. Its PDF is
f ( θ ) = 2 σ ( α + 2 ) k = ϕ θ + 2 π k μ σ 1 + α Φ λ σ ( θ + 2 π k μ )
for 0 θ < 2 π , < μ < , σ > 0 , < λ < , and α 1 . The nth trigonometric moment is
m n = e i μ n γ α 2 β + i n 2 1 2 2 λ
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan sin μ + α α + 2 G ( δ σ ) cos μ cos μ α α + 2 G ( δ σ ) sin μ ,
ρ = η e σ 2 2 ,
γ 1 = e p 2 σ 2 2 sin p 2 μ ω + α α + 2 G ( p δ σ ) cos p 2 μ ω 1 η e σ 2 2 3 2 ,
and
γ 2 = e p 2 σ 2 2 cos p ω α α + 2 G ( δ σ p ) sin p ω η 4 e 2 σ 2 1 η e σ 2 2 2 ,
respectively, where η = 1 + α 2 + α G ( δ σ ) , G ( d ) = 2 π n = 0 d 2 n + 1 2 n n ! ( 2 n + 1 ) for d real, and ω = μ sin μ + α α + 2 G ( p δ σ ) cos μ cos μ α α + 2 G ( p δ σ ) sin μ . The wrapped generalized skew normal distribution is hindered by its non-closed-form PDF and excessive parameters, potentially outperformed by simpler models.

2.14. Wrapped Generalized Skew Normal Distribution [11]

Ref. [11] took g to be the PDF of the generalized skew normal distribution to obtain the wrapped generalized skew normal distribution. Its PDF is
f ( θ ) = 2 w k = ϕ θ + 2 π k μ w Φ α θ + 2 π k μ w + β θ + 2 π k μ w 3
for 0 θ < 2 π , < μ < , w > 0 , < α < , and < β < . No properties were derived for this distribution. The wrapped generalized skew normal distribution was utilized to analyze wind direction data from a meteorological station in Villena, Alicante, Spain. The data, collected in June 2009 using an Oregon Scientific WMR928NX automatic weather station, were split into sea breeze and mountain breeze subsets. The wrapped generalized skew normal distribution appeared to be the most suitable fit for both the entire dataset and the mountain breeze period. However, for the sea breeze period, the Jones and Pewsey sine-skewed distribution demonstrated the best fit.

2.15. Wrapped Half-Logistic Distribution

Ref. [15] took g to be the PDF of the half-logistic distribution to obtain the wrapped half-logistic distribution. Its PDF is
f ( θ ) = k = 1 2 σ sech 2 θ + 2 k π μ 2 σ
for 0 θ < 2 π , 0 μ < 2 π , and σ > 0 . The nth trigonometric moment is
m n = 2 b n cos ( p μ ) c n sin ( p μ ) + 2 i b n sin ( p μ ) + c n cos ( p μ ) ,
where
b n = 0 a cos ( p σ y ) e y 1 + e y 2 d y + n = 1 ( 1 ) n 1 n e n a n 2 + p 2 σ 2 n cos ( p σ a ) p σ sin ( p σ a )
and
c n = 0 a sin ( p σ y ) e y 1 + e y 2 d y + n = 1 ( 1 ) n 1 n e n a n 2 + p 2 σ 2 n sin ( p σ a ) + p σ cos ( p σ a ) .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan b 1 sin ( μ ) + c 1 cos ( μ ) b 1 cos ( μ ) c 1 sin ( μ ) ,
ρ = 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 ,
γ 1 = 2 b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) 2 + b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) 2 cos μ 2 2 μ 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 3 2 ,
and
γ 2 = 2 b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) 2 + b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) 2 sin μ 2 2 μ 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 2 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 4 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 2 ,
respectively, where μ 2 = arctan b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) . The utility of the wrapped half-logistic distribution is restricted by its non-closed-form PDF.

2.16. Wrapped Half-Normal Distribution

If g is the pdf of the half-normal distribution we obtain the wrapped half-normal distribution. Its pdf is
f ( θ ) = k = 0 1 2 σ sech 2 θ + 2 k π μ 2 σ
for 0 θ < 2 π and σ > 0 . The nth trigonometric moment is
m n = 2 b n cos ( p μ ) c n sin ( p μ ) + 2 i b n sin ( p μ ) + c n cos ( p μ ) ,
where
b n = 0 a cos ( p σ y ) e y 1 + e y 2 d y + n = 1 ( 1 ) n 1 n e m n 2 + p 2 σ 2 n cos ( p σ a ) p σ sin p σ a
and
c n = 0 a cos ( p σ y ) e y 1 + e y 2 d y + n = 1 ( 1 ) n 1 n e m n 2 + p 2 σ 2 n sin ( p σ a ) p σ cos p σ a .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan b 1 sin ( μ ) + c 1 cos ( μ ) b 1 cos ( μ ) c 1 sin ( μ ) ,
ρ = 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 ,
γ 1 = 2 b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) 2 + b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) 2 sin μ 2 2 μ 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 3 2
and
γ 2 = 2 b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) 2 + b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) 2 sin μ 2 2 μ 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 2 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 4 1 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 2 ,
respectively, where μ 2 = arctan b 2 sin ( 2 μ ) + c 2 cos ( 2 μ ) b 2 cos ( 2 μ ) c 2 sin ( 2 μ ) .

2.17. Wrapped [12]’s Skew Laplace Distribution

Ref. [28] took g to be the PDF of [12]’s skew Laplace distribution to obtain the wrapped [12]’s skew Laplace distribution. Its PDF and CDF are
f ( θ ) = p c e c β + 2 π 2 β π θ + ( 1 p ) e c + ( 2 π 1 ) c β e 2 π c 1
and
F ( θ ) = p e c β + 2 π 2 β π 1 e c θ + ( 1 p ) e ( 2 π 1 ) c β e c θ 1 e 2 π c 1 ,
respectively, for 0 θ < 2 π , < β < , c > 0 , and 0 < p < 1 . The nth trigonometric moment is
m n = α n + i β n
for n = 1 , 2 , , where
α n = c 2 n 2 + c 2 n e ( 1 2 π ) β c + ( 1 n ) e ( 2 π 1 ) β c
and
β n = n c n 2 + c 2 ( n 1 ) e ( 2 π 1 ) β c n e ( 1 2 π ) β c .
The distribution was tested on the Black Mountain wind direction dataset [1]. It showed a superior fit compared to other distributions, specifically the wrapped variance gamma distribution and the generalized von Mises distribution, which have more parameters.

2.18. Wrapped Hypoexponential Distribution

Ref. [29] took g to be the PDF of the hypoexponential distribution to obtain the wrapped hypoexponential distribution. Its PDF is
f ( θ ) = λ 1 λ 2 λ 2 λ 1 e λ 1 θ k = 0 e 2 k π λ 1 e λ 2 θ k = 0 e 2 k π λ 2
for 0 θ < 2 π , λ 1 > 0 , and λ 2 > 0 . The nth trigonometric moment is
m n = λ 1 λ 2 λ 1 i n λ 2 i n
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan 1 λ 1 + 1 λ 2 ,
ρ = λ 1 λ 2 1 + λ 1 2 1 + λ 2 2 ,
γ 1 = λ 1 λ 2 4 + λ 1 2 4 + λ 2 2 sin arctan 2 λ 1 + 2 λ 2 2 arctan 1 λ 1 + 1 λ 2 1 λ 1 λ 2 1 + λ 1 2 1 + λ 2 2 3 4 ,
and
γ 2 = λ 1 λ 2 4 + λ 1 2 4 + λ 2 2 sin arctan 2 λ 1 + 2 λ 2 2 arctan 1 λ 1 + 1 λ 2 λ 1 4 λ 2 4 1 + λ 1 2 2 1 + λ 2 2 2 1 λ 1 λ 2 1 + λ 1 2 1 + λ 2 2 2 ,
respectively, The wrapped hypoexponential distribution has a PDF that can be recast into a closed form. It also offers additional flexibility compared to one-parameter models.

2.19. Wrapped Ishita Distribution

Ref. [30] took g to be the PDF of the Ishita distribution to obtain the wrapped Ishita distribution. Its PDF and CDF are
f ( θ ) = λ 3 e λ θ λ + 2 λ + θ 2 1 e 2 π λ + 4 π θ e 2 π λ 1 e 2 π λ 2 + 4 π 2 e 2 π λ 1 + e 2 π λ 1 e 2 π λ 3
and
F ( θ ) = 1 λ 3 + 2 { λ 3 + 2 λ 3 + 2 + λ θ λ θ + 2 e λ θ 1 e 2 π λ + 4 π λ 1 ( λ θ + 1 ) e λ θ e 2 λ π 1 e 2 π λ 2 + 4 π 2 λ 2 1 e λ θ 1 + e 2 π λ e 2 λ π 1 e 2 π λ 3 } ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is
m n = λ 3 λ 3 + 2 λ 3 n 2 λ + 2 2 i n λ 2 ( λ i n ) 3
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 3 arctan 1 λ arctan 2 λ 2 λ 3 λ + 2 ,
ρ = λ 3 λ 3 + 2 λ 3 λ + 2 2 + 4 λ 4 λ 2 + 1 3 ,
γ 1 = λ 3 λ 3 + 2 λ 3 4 λ + 2 2 + 16 λ 4 λ 2 + 4 3 sin κ λ , 2 1 λ 3 λ 3 + 2 λ 3 λ + 2 2 + 4 λ 4 λ 2 + 1 3 3 2 ,
and
γ 2 = λ 3 λ 3 + 2 λ 3 4 λ + 2 2 + 16 λ 4 λ 2 + 4 3 cos κ λ , 2 λ 3 λ 2 + 2 λ 2 λ + 2 2 + 4 λ 2 λ 2 + 1 3 4 1 λ 3 λ 2 + 2 λ 2 λ + 2 2 + 4 λ 2 λ 2 + 1 3 2 ,
respectively, where κ λ , 2 = 3 arctan 2 λ arctan 4 λ 2 λ 3 4 λ + 2 6 arctan 1 λ + 2 arctan 2 λ 2 λ 3 λ + 2 . The wrapped Ishita distribution, despite having closed-form expressions for PDF and CDF, might be unsuitable for some real-life datasets, as its flexibility is restricted by the fact it only has one parameter.

2.20. Wrapped Laplace Distribution

Ref. [31] took g to be the PDF of the Laplace distribution to obtain the wrapped Laplace distribution. Its PDF and CDF are
f ( θ ) = λ κ 1 + κ 2 e λ κ θ 1 e 2 π λ κ + e λ θ κ e 2 π λ κ 1
and
F ( θ ) = 1 1 + κ 2 1 e κ λ θ 1 e 2 π κ λ + κ 2 1 + κ 2 e λ θ κ 1 1 e 2 π λ κ ,
respectively, for 0 θ < 2 π , κ > 0 , and λ > 0 . The nth trigonometric moment is
m n = 1 1 i n λ κ 1 + i n λ κ
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan 1 λ κ arctan κ λ , for κ 1 , 2 π + arctan 1 λ κ arctan κ λ , for κ > 1 ,
ρ = λ 2 1 + ( λ κ ) 2 λ 2 κ 2 + 1 ,
γ 1 = 2 λ 1 + λ 2 1 4 4 + λ 2 1 + λ 2 λ 3 2 ,
and
γ 2 = 3 λ 2 1 + λ 2 4 + λ 2 1 + λ 2 | λ | 2 ,
respectively. The ant orientation data from [17] were analyzed using the wrapped Laplace distribution, which demonstrated a superior fit compared to the von Mises distribution. However, the scope of conclusions is constrained by the limited number of distributions considered. To better evaluate the goodness of fit, it is essential to compare the wrapped Laplace distribution against a broader range of distributions.

2.21. Wrapped Length-Biased Weighted Exponential Distribution

Ref. [32] took g to be the PDF of the length-biased weighted exponential distribution to obtain the wrapped length-biased weighted exponential distribution. Its PDF and CDF are
f ( θ ) = λ ( α + 1 ) 2 α ( α + 2 ) e λ β 1 1 e 2 π λ β + 2 π e 2 π λ 1 e 2 π λ e α β λ 1 e 2 π λ ( 1 + α ) β + 2 π e 2 π λ ( 1 + α ) 1 e 2 π λ ( 1 + α )
and
F ( θ ) = 1 α ( α + 2 ) { ( 1 + α ) 2 1 ( 1 + β ) e β 1 e 2 π + ( 1 + α ) 2 2 π 1 e β e 2 π 1 e 2 π 2 + e β ( 1 + α ) 1 + β ( 1 + α ) 1 1 e 2 π ( 1 + α ) + 2 π ( 1 + α ) e 2 π ( 1 + α ) e β ( 1 + α ) 1 1 e 2 π ( 1 + α ) 2 } ,
respectively, for 0 θ < 2 π , α > 0 , λ > 0 , and 0 < β 2 π . The nth trigonometric moment is
m n = ( 1 i n ) 2 1 i n 1 + α 2 1 2 i n 2 + α 1
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan ( 1 ) + arctan 1 1 + α arctan 2 2 + α ,
ρ = 1 + ( 2 + α ) 2 2 1 + ( 1 + α ) 2 ,
γ 1 = sin μ 2 2 μ 5 1 1 + ( 2 + α ) 2 1 2 2 1 + ( 1 + α ) 2 3 2 1 + 4 ( 1 + α ) 2 1 + 4 2 + α 2 ,
and
γ 2 = cos μ 2 2 μ 5 1 1 + ( 2 + α ) 2 1 2 2 1 + ( 1 + α ) 2 2 1 + 4 ( 1 + α ) 2 1 + 4 2 + α 2 1 + ( 2 + α ) 2 1 2 2 1 + ( 1 + α ) 2 4 1 1 + ( 2 + α ) 2 1 2 2 1 + ( 1 + α ) 2 2 ,
respectively, where μ 2 = arctan ( 2 ) + arctan 2 1 + α arctan 4 2 + α . The wrapped length-biased weighted exponential distribution was used to analyze the feldspar laths dataset obtained from [33] and published by [1]. Watson’s U 2 test indicated a good fit of the distribution to the data. No comparisons with other distributions were performed.

2.22. Wrapped Levy Distribution

Ref. [1] took g to be the PDF of the Levy distribution to obtain the wrapped Levy distribution. Its PDF is
f ( θ ) = k = c 2 π e c 2 ( θ + 2 π n μ ) ( θ + 2 π n μ ) 3 2
for 0 θ < 2 π , the summand is zero if θ + 2 π n μ 0 , c > 0 , and < μ < . The nth trigonometric moment is
m n = e i n μ c | n | 1 i sgn ( n ) = e i n μ c n ( 1 i )
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ + c ,
ρ = e c ,
γ 1 = e 2 c sin 2 c 2 c 1 e c 3 2 ,
and
γ 2 = e 2 c cos 2 c 2 c e 4 c 1 e c 2 ,
respectively. The practical utility of the wrapped Levy distribution is hampered by its non-closed-form PDF.

2.23. Wrapped Lindley Distribution

Ref. [34] took g to be the PDF of the Lindley distribution to obtain the wrapped Lindley distribution. Its PDF and CDF are
f ( θ ) = λ 2 1 + λ e λ θ 1 + θ 1 e 2 π λ + 2 π e 2 π λ 1 e 2 π λ 2
and
F ( θ ) = 1 1 e 2 π λ 1 e λ θ λ θ λ + 1 e λ θ 2 π λ λ + 1 e 2 π λ 1 e λ θ 1 e 2 π λ 2 ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is given by
m n = λ 2 1 + λ ( 1 + λ i n ) ( λ i n ) 2
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 2 arctan 1 λ arctan 1 λ + 1 ,
ρ = λ 2 ( λ + 1 ) 2 + 1 ( 1 + λ ) λ 2 + 1 ,
γ 1 = λ 2 4 + ( 1 + λ ) 2 sin μ 2 2 μ ( 1 + λ ) 4 + λ 2 1 λ 2 1 + ( 1 + λ ) 2 1 2 ( 1 + λ ) 1 + λ 2 3 2 ,
and
γ 2 = λ 2 4 + ( 1 + λ ) 2 cos μ 2 2 μ λ 2 1 + 1 + λ 2 ( 1 + λ ) 1 + λ 2 4 ( 1 + λ ) 4 + λ 2 1 λ 2 1 + ( 1 + λ ) 2 ( 1 + λ ) 1 + λ 2 4 ,
respectively, where μ 2 = 2 arctan 2 λ arctan 2 λ + 1 . The wrapped Lindley distribution was used to analyze a dataset concerning the orientations of 76 turtles after laying eggs, as documented by [1]. This distribution seemed to match the data effectively and performed better than the alternative distribution it was contrasted with, the wrapped exponential distribution.

2.24. Wrapped Linnik Distribution

Ref. [35] took g to be the PDF of the Linnik distribution to obtain the wrapped Linnik distribution. Its PDF and CDF are
f ( θ ) = 1 2 π 1 + 2 k = 1 cos ( k θ ) 1 + σ k α
and
F ( θ ) = 1 2 π θ + 2 k = 1 sin ( k θ ) k 1 + σ k α ,
respectively, for 0 θ < 2 π , σ > 0 , and 0 < α 2 . The nth trigonometric moment is
m n = 1 1 + σ n α
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 0 ,
ρ = 1 1 + σ ,
γ 1 = 0 ,
and
γ 2 = ( 1 + σ ) 4 1 + 2 α σ σ 2 α 1 + 2 α σ ( 1 + σ ) 2 ,
respectively. The wrapped Linnik distribution was used to analyze the frequency of traffic accidents throughout the day in Srinagar, India, in 2016, sourced from reports from the National Crime Records Bureau, India. This distribution was also applied to ant data collected by [17]. Results indicated that the wrapped Linnik distribution provided a better fit for both datasets compared to the wrapped stable distribution.

2.25. Wrapped Lomax Distribution

Ref. [36] took g to be the PDF of the Lomax distribution to obtain the wrapped Lomax distribution. Its PDF and CDF are
f ( θ ) = α σ k = 0 1 + θ + 2 k π σ α 1
and
F ( θ ) = k = 0 1 + 2 k π σ α k = 0 1 + θ + 2 k π σ α ,
respectively, for 0 θ < 2 π , σ > 0 , and α > 0 . The wrapped Lomax distribution faces limitations due to its non-closed-form PDF.

2.26. Wrapped Modified Lindley Distribution

Ref. [37] took g to be the PDF of the modified Lindley distribution to obtain the wrapped modified Lindley distribution. Its PDF and CDF are
f ( θ ) = λ e λ θ 1 e 4 λ π 1 + e 2 λ π + e λ θ 1 + λ 2 λ θ 1 + 4 λ π e 4 λ π 1 e 4 λ π
and
F ( θ ) = 1 e λ θ 1 e 2 λ π + 2 λ π e 4 λ π 1 e 2 λ θ ( 1 + λ ) 1 e 4 λ θ 2 λ θ e 2 λ θ ( 1 + λ ) 1 e 4 λ θ ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is
m n = λ λ n i + λ n i ( 1 + λ ) ( 2 λ n i ) 2
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = π + arctan β 1 α 1 , if   α 1 < 0 , β 1 0 , π 2 , if   α 1 < 0 , β 1 > 0 , arctan β 1 α 1 , if   α 1 > 0 , u n d e f i n e d , if   α 1 = 0 , β 1 = 0
and
ρ = λ 2 16 λ 4 + 32 λ 3 + 24 λ 2 + 16 λ + 10 ( 1 + λ ) 2 λ 2 + 1 4 λ 2 + 1 2 ,
respectively, where
α 1 = λ 2 λ 2 + 1 + 4 λ 2 ( 1 + λ ) 4 λ 2 + 1 2
and
β 1 = λ λ 2 + 1 + 4 λ 3 ( 1 + λ ) 4 λ 2 + 1 2 λ ( 1 + λ ) 4 λ 2 + 1 2 .
The wrapped modified Lindley distribution was used on two real-life datasets from [1]. The first set involves 76 turtles’ egg-laying orientations, while the second set has 133 measurements of feldspar lath orientations in basalt. This distribution was compared with the wrapped exponential distribution, transmuted wrapped exponential distribution, and wrapped Lindley distribution for goodness of fit. The wrapped modified Lindley distribution provided a competitive fit for both datasets. It offered the best fit for the turtle dataset but did not outperform the transmuted wrapped exponential distribution for the feldspar lath dataset.

2.27. Wrapped New Weibull–Pareto Distribution

Ref. [38] took g to be the PDF of the new Weibull–Pareto distribution to obtain the wrapped new Weibull–Pareto distribution. Its PDF is
f ( θ ) = k = 1 4 σ sech 2 θ + 2 k π μ 2 σ
for 0 θ < 2 π , θ > μ , and σ > 0 . The nth trigonometric moment is
m n = π σ n sinh π σ n e i n μ
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ ,
ρ = π σ p sinh π σ p ,
γ 1 = 0 ,
and
γ 2 = 2 π σ sinh ( 2 π σ ) π σ sinh ( π σ ) 4 1 π σ sinh ( π σ ) 2 ,
respectively. The wrapped new Weibull–Pareto distribution and the wrapped exponentiated inverted Weibull distribution were tested on the turtle dataset from [1]. Both models fit the data well, but the wrapped exponentiated inverted Weibull distribution slightly outperformed the wrapped new Weibull–Pareto distribution.

2.28. Wrapped Normal Distribution

Ref. [39] took g to be the PDF of the normal distribution to obtain the wrapped normal distribution. Its PDF and CDF are
f ( θ ) = 1 2 π σ k = e ( θ μ + 2 π k ) 2 2 σ 2
and
F ( θ ) = k = Φ θ + 2 π k μ σ Φ 2 π k μ σ ,
respectively, for 0 θ < 2 π , < μ < , and σ > 0 . The nth trigonometric moment is
m n = e i n μ n 2 σ 2 2
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ ,
ρ = e σ 2 2 ,
γ 1 = 0 ,
and
γ 2 = 0 ,
respectively. The wrapped normal distribution, commonly used in literature, adequately fits real-life data but is restricted by its non-closed-form PDF.

2.29. Wrapped Pareto Distribution

Ref. [40] took g to be the PDF of a Pareto distribution to obtain the wrapped Pareto distribution. Its PDF is
f ( θ ) = α 1 α Γ 1 α 0 λ s e λ s θ 1 e 2 π λ s s 1 α 1 e s α d s
for 0 θ < 2 π , λ > 0 , and α > 0 . The nth trigonometric moment is
m n = α 1 α Γ 1 α 0 1 i n λ s 1 s 1 α 1 e s α d s
for n = 1 , 2 , . The wrapped Pareto distribution is constrained by its non-closed-form PDF.

2.30. Wrapped Quasi-Lindley Distribution

Ref. [41] took g to be the PDF of the quasi-Lindley distribution to obtain the wrapped quasi-Lindley distribution. Its PDF is
f ( θ ) = b e b θ a + 1 a + b θ 1 e 2 π b + 2 π b e 2 π b 1 e 2 π b 2
for 0 θ < 2 π , a > 0 , and b > 0 . The nth trigonometric moment is
m n = α n + i β n
for n = 1 , 2 , , where
α n = b 4 ( a + 1 ) 2 + n 2 a 2 b 2 ( a + 1 ) b 2 + n 2 cos 2 arctan n b arctan n a b ( 1 + a )
and
β n = b 4 ( a + 1 ) 2 + n 2 a 2 b 2 ( a + 1 ) b 2 + n 2 sin 2 arctan n b arctan n a b ( 1 + a ) .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = 2 arctan 1 b a arctan a b ( 1 + a ) ,
ρ = b 4 ( a + 1 ) 2 + a 2 b 2 ( a + 1 ) b 2 + 1 ,
γ 1 = β ¯ 2 ρ 3 2 ,
and
γ 2 = α ¯ 2 ρ 4 ρ 2 ,
respectively, where
α ¯ 2 = b 4 ( a + 1 ) 2 + 4 a 2 b 2 ( a + 1 ) b 2 + 4 cos 2 arctan 2 b arctan 2 a b ( 1 + a ) 4 arctan 1 b + arctan a b ( 1 + a )
and
β ¯ 2 = b 4 ( a + 1 ) 2 + 4 a 2 b 2 ( a + 1 ) b 2 + 4 sin 2 arctan 2 b arctan 2 a b ( 1 + a ) 4 arctan 1 b + arctan a b ( 1 + a ) .

2.31. Wrapped Rama Distribution

Ref. [42] took g to be the PDF of the Rama distribution to obtain the wrapped Rama distribution. Its PDF and CDF are
f ( θ ) = A e λ θ θ 3 + B θ 2 + C θ + D
and
F θ = A λ 4 γ ( 4 , λ θ ) + B λ γ ( 3 , λ θ ) + C λ 2 γ ( 2 , λ θ ) + D λ 3 γ ( 1 , λ θ )
for 0 θ < 2 π and λ > 0 , where
A = λ 4 λ 3 + 6 1 e 2 λ π ,
B = 6 π e 2 λ π 1 e 2 λ π ,
C = 2 π e 2 λ π + 1 B 1 e 2 λ π ,
and
D = 1 + 8 π 3 e 2 λ π e 4 λ π + 4 e 2 λ π + 1 1 e 2 λ π 3 .
The nth trigonometric moment is
m n = A λ 2 + n 2 4 c 0 + c 2 n 2 + c 4 n 4 + c 6 n 6 + i d 1 n + d 3 n 2 + d 5 n 5 + d 7 n 7
for n = 1 , 2 , , where
c 0 = λ 7 4 λ 6 B π e π λ 2 4 λ 7 B π 2 e 2 π λ 2 λ 7 C π e 2 π λ 2 λ 5 B e 2 π λ λ 6 C e 2 π λ 6 λ 4 e 2 π λ 12 π 2 λ 6 e 2 π λ 8 π 3 λ 7 e 2 π λ 12 π λ 5 e 2 π λ + 6 λ 4 + λ 6 C λ 7 e 2 π λ + 2 λ 5 B ,
c 2 = λ 4 C e 2 π λ + 4 λ 3 B e 2 π λ + 24 π λ 3 e 2 π λ 12 π 2 λ 4 e 2 π λ 24 π 3 λ 5 e 2 π λ + λ 4 C 3 λ 5 e 2 π λ 4 λ 3 B 36 λ 2 + 36 λ 2 e 2 π λ + 3 λ 5 4 λ 4 B π e 2 π λ 6 λ 5 C π e 2 π λ 12 λ 5 B π 2 e 2 π λ ,
c 4 = λ 2 C e 2 π λ + 6 λ B e 2 π λ + 36 π λ e 2 π λ + 12 π 2 λ 2 e 2 π λ 24 π 3 λ 3 e 2 π λ 3 λ 3 e 2 π λ λ 2 C 6 λ B 6 e 2 π λ + 6 + 3 λ 3 6 C π λ 3 e 2 π λ 12 B π 2 λ 3 e 2 π λ + 4 B π λ 2 e 2 π λ ,
c 6 = 4 B π e 2 π λ 8 π 3 λ e 2 π λ λ e 2 π λ + C e 2 π λ + 12 π 2 e 2 π λ + λ C 2 λ C π e 2 π λ 4 λ B π 2 e 2 π λ ,
d 1 = 2 λ 5 C λ 6 e 2 π λ + 6 λ 4 B 6 λ 4 B e 2 π λ 2 λ 5 C e 2 π λ 24 π 2 λ 5 e 2 π λ 8 π 3 λ 6 e 2 π λ 3 π λ 4 e 2 π λ + 24 λ 3 24 λ 3 e 2 π λ + λ 6 8 λ 5 B π e 2 π λ 4 λ 6 B π 2 e 2 π λ 2 λ 6 C π e 2 π λ ,
d 3 = 3 λ 4 e 2 π λ + 4 λ 3 C + 4 λ 2 B 4 λ 3 C e 2 π λ 4 λ 2 B e 2 π λ 24 π λ 2 e 2 π λ 48 π 2 λ 3 e 2 π λ 24 π 3 λ 4 e 2 π λ 24 λ + 24 λ e 2 π λ + 3 λ 4 6 λ 4 C π e 2 π λ 12 λ 4 B π 2 e 2 π λ 16 λ 3 B π e 2 π λ ,
d 5 = 2 λ C + 2 B e 2 π λ 3 λ 2 e 2 π λ 2 λ C e 2 π λ 24 π 2 λ e 2 π λ 24 π 3 λ 2 e 2 π λ + 12 π e 2 π λ + 3 λ 2 2 B 6 λ 2 C π e 2 π λ 12 λ 2 B π 2 e 2 π λ 8 λ B π e 2 π λ ,
and
d 7 = 1 2 C π e 2 π λ 4 B π 2 e 2 π λ 8 π 3 e 2 π λ 1 e 2 π λ .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arcsin d 1 + d 3 + d 5 + d 7 c 0 + c 2 + c 4 + c 6 2 + d 1 + d 3 + d 5 + d 7 2 ,
ρ = A λ 2 + 1 4 c 0 + c 2 + c 4 + c 6 2 + d 1 + d 3 + d 5 + d 7 2 ,
γ 1 = A e 2 μ v 3 2 λ 2 + 4 4 2 d 1 + 4 d 3 + 32 d 5 + 128 d 7 ,
and
γ 2 = A e 2 μ 1 ρ 2 λ 2 + 4 4 c 0 + 4 c 2 + 16 c 4 + 64 c 6 ρ 4 1 ρ 2 ,
respectively. The wrapped Rama distribution boasts closed-form expressions for its PDF and CDF. This feature, coupled with its good fit to various datasets, makes it a compelling choice for statistical analysis. In fact, it outperformed six other widely used distributions (with up to three parameters) when applied to two specific datasets. The first dataset, initially obtained by [33] and later published by [1], consists of long-axis orientation measurements for 60 feldspar laths in basalt, recorded in degrees. The second dataset, also reported by [1], comprises horizontal axis values for 100 outwash pebbles collected from a late Wisconsin outwash terrace near Cary, Illinois, along the Fox River.

2.32. Wrapped Richard Distribution

Ref. [43] introduced a wrapped distribution based on the Richard link function. Its PDF is
f ( θ ) = k = 0 k e k θ + 2 π n 1 m 1 1 m 1 + ( m 1 ) e k θ + 2 π n m 1 m
for 0 θ < 2 π , m 1 , and n 1 . The wrapped Richard distribution is limited by its non-closed-form PDF.

2.33. Wrapped Shanker Distribution

Ref. [44] took g to be the PDF of the Shanker distribution to obtain the wrapped Shanker distribution. Its PDF and CDF are
f ( θ ) = λ 2 λ 2 + 1 e λ θ λ + θ 1 e 2 π λ + 2 π e 2 π λ 1 e 2 π λ 2
and
F ( θ ) = 1 1 e 2 π λ 1 e λ θ λ θ e λ θ λ 2 + 1 + 2 π λ 1 e λ θ e 2 π λ λ 2 + 1 1 e 2 π λ 2 ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is
m n = λ 2 1 e 2 π λ e 2 π ( i n λ ) 1 λ ( i n λ ) 1 + 2 π λ 2 ( i n λ ) e 2 π ( i n λ ) λ 2 + 1 1 e 2 π λ 2 ( i n λ ) 2
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 2 arctan 1 λ arctan λ λ 2 + 1 ,
ρ = λ 2 λ 2 + 1 2 + λ 2 λ 2 + 1 2 ,
γ 1 = λ 2 λ 2 + 1 2 + 4 λ 2 λ 2 + 1 λ 2 + 4 sin κ 2 , λ 1 λ 2 λ 2 + 1 2 + λ 2 λ 2 + 1 2 3 2 ,
and
γ 2 = λ 2 λ 2 + 1 2 + 4 λ 2 λ 2 + 1 λ 2 + 4 cos κ 2 , λ λ 2 λ 2 + 1 2 + λ 2 λ 2 + 1 2 4 1 λ 2 λ 2 + 1 2 + λ 2 λ 2 + 1 2 2 ,
respectively, where κ λ , 2 = 2 arctan 2 λ arctan 2 λ λ 2 + 1 4 arctan 1 λ + 2 arctan λ λ 2 + 1 . The wrapped Shanker distribution has closed-form expressions for the PDF and CDF, but it might have limited utility for real-life data modeling due to the fact it only has one parameter.

2.34. Wrapped Skew Laplace Distribution

Ref. [45] took g to be the PDF of a skew Laplace distribution to obtain a wrapped skew Laplace distribution. Its PDF is
f ( θ ) = 1 2 σ e θ μ σ + e θ μ σ + e μ θ σ e 2 π σ 1 + A θ ; λ 1 , λ 2 , λ 3
for 0 θ < 2 π , < λ 1 < , λ 2 > 0 , < λ 3 < , < μ < , and σ > 0 , where
A ( θ ) = k = e θ + 2 π k μ σ 1 e g θ + 2 π k μ σ sign g θ + 2 π k μ σ
and
g ( x ) = λ 1 x + λ 3 x 3 1 + λ 2 x 2 .
The nth trigonometric moment is
m n = α n + i β n
for n = 1 , 2 , , where
α n = cos ( n μ ) + 2 e ξ sin ( n μ ) ξ n n 2 σ 2 + 1 Δ n σ sin ( n μ ) n 2 σ 2 + 1 Δ 2 C n ,
β n = sin ( n μ ) 2 e ξ cos ( n μ ) ξ n n 2 σ 2 + 1 Δ n σ cos ( n μ ) n 2 σ 2 + 1 Δ 2 S n ,
= sign λ 1 , if   λ 1 λ 3 < 0 , 0 , if   λ 1 λ 3 0 ,
Δ = 0 , if   λ 1 = λ 3 = 0 , sign λ 3 , if   λ 1 = 0 , λ 3 0 , sign λ 1 , if   λ 1 0 ,
C n = 1 2 σ k = 0 2 π cos ( n θ ) e θ + 2 π k μ σ e g θ + 2 π k μ σ sign g θ + 2 π k μ σ d θ ,
S n = 1 2 σ k = 0 2 π sin ( n θ ) e θ + 2 π k μ σ e g θ + 2 π k μ σ sign g θ + 2 π k μ σ d θ ,
ξ n = sin ( ξ n σ ) + n σ cos ( ξ n σ ) ,
and
ξ = λ 1 λ 3 , if   λ 1 λ 3 < 0 , 0 , if   λ 1 λ 3 0 .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan 1 σ Δ Δ 2 S 1 1 + σ 2 2 e ξ ξ 1 ,
ρ = σ Δ σ 2 + 1 Δ 2 S 1 2 e ξ ξ 1 σ 2 + 1 sin μ + cos μ σ 2 + 1 ,
γ 1 = 2 σ Δ Δ 2 S 2 4 σ 2 + 1 2 e ξ ξ 2 cos ( 2 μ ) sin ( 2 μ ) 4 σ 2 + 1 ( 1 ρ ) 3 2 ,
and
γ 2 = 2 σ Δ Δ 2 S 2 4 σ 2 + 1 2 e ξ ξ 2 sin ( 2 μ ) cos ( 2 μ ) 4 σ 2 + 1 ρ 4 4 σ 2 + 1 ( 1 ρ ) 2 ,
respectively. The wrapped skew Laplace distribution was used on the turtle and ant datasets discussed by [1]. It was observed that this distribution provided a better fit for the first dataset compared to the wrapped Lindley, exponential, transmuted wrapped exponential, and non-negative trigonometric sums distributions [46]. Similarly, for the second dataset, the wrapped skew Laplace distribution outperformed the symmetric wrapped Laplace and non-negative trigonometric sums distributions.

2.35. Wrapped Skew Normal Distribution

Ref. [47] took g to be the PDF of the skew normal distribution to obtain the wrapped skew normal distribution. Its PDF is
f ( θ ) = 2 η k = ϕ θ + 2 π r ξ η Φ λ θ + 2 π r ξ η
for 0 θ < 2 π , < ξ < , η > 0 , and < λ < . The nth trigonometric moment is
m n = e i n ξ 1 2 n 2 η 2 1 + i J ( δ η n )
for n = 1 , 2 , , where
J ( x ) = 0 x b e u 2 2 d u
and δ = λ 1 + λ 2 . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan sin ξ + J ( δ η ) cos ξ cos ξ J ( δ η ) sin ξ ,
ρ = e η 2 2 1 + J 2 ( δ η ) 1 2 ,
γ 1 = ω 4 J ( 2 δ η ) 1 J 2 ( δ η ) 2 J ( δ η ) 1 + J 2 ( δ η ) 1 e η 2 2 1 + J 2 ( δ η ) 1 2 3 2 ,
and
γ 2 = ω 4 1 J 2 ( δ η ) + 2 J ( δ η ) J ( 2 δ η ) 1 + J 2 ( δ η ) e η 2 2 1 + J 2 ( δ η ) 1 2 4 1 e η 2 2 1 + J 2 ( δ η ) 1 2 2 ,
respectively. The wrapped skew normal distribution was used to analyze bird heading data from the autumn migration of 1987 [27]. Although it fit the data well, its PDF does not have a simple mathematical expression, which might restrict its practical utility.

2.36. Wrapped Stable Distribution

Ref. [48] took g to be the PDF of the stable distribution to obtain the wrapped stable distribution. The nth trigonometric moment is
m n = e γ α n α 1 + i β ( γ n ) 1 α 1 tan π α 2 + i δ 0 * n , if   α 1 , e γ n 1 + 2 β i π log ( γ n ) + i δ 0 * n , if   α = 1
for n = 1 , 2 , , where 0 < α 2 , 1 β 1 , γ > 0 , < δ 0 < , and δ 0 * = δ 0 mod 2 π . The mean direction, mean resultant length, skewness, and kurtosis are
μ = δ 0 * + β γ α γ tan π α 2   mod   2 π , if   α 1 , δ 0 * 2 β γ π log γ   mod   2 π , if   α = 1 ,
ρ = e γ α ,
γ 1 = β ¯ 2 ( 1 ρ ) 3 2 ,
and
γ 2 = α ¯ 2 ρ 4 ( 1 ρ ) 2 ,
respectively, where
α ¯ 2 = ρ 2 α cos β γ α 2 α 2 tan π α 2 , if   α 1 , ρ 2 α cos 2 β γ log 2 π , if   α = 1
and
β ¯ 2 = ρ 2 α sin β γ α 2 α 2 tan π α 2 , if   α 1 , ρ 2 α sin 2 β γ log 2 π , if   α = 1 .
The wrapped stable distribution demonstrated a slightly superior fit compared to a mixture distribution incorporating circular uniform and wrapped skew normal components when analyzing the bird heading data from [27].

2.37. Wrapped Student’s t Distribution

Ref. [49] took g to be the PDF of the Student’s t distribution to obtain the wrapped Student’s t distribution. Its PDF is
f ( θ ) = c λ k = 1 + θ + 2 π k μ 0 2 λ 2 ν ν + 1 2
for 0 θ < 2 π , 0 μ < 2 π , λ > 0 , and ν > 0 , where
c = Γ ν + 1 2 Γ ν 2 π ν .
The nth trigonometric moment is
m n = K ν 2 n ν n ν ν 2 Γ ν 2 2 ν 2 1
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ ,
ρ = K ν 2 λ ν λ ν ν 2 Γ ν 2 2 ν 2 1 ,
γ 1 = e 2 c sin ( 2 c 2 c ) 1 e c 3 2 ,
and
γ 2 = e 2 c cos 2 c 2 c e 4 c 1 e c 2 ,
respectively. The wrapped Student’s t distribution was used on a dataset of 104 cross-bed measurements from Himalayan molasse in Pakistan, as discussed by [1]. The wrapped Student’s t distribution yielded a better fit to this dataset compared to the von Mises distribution.

2.38. Wrapped Transmuted Exponential Distribution

Ref. [50] took g to be the PDF of the transmuted exponential distribution to obtain the wrapped transmuted exponential distribution. Its PDF and CDF are
f ( θ ) = 2 λ Λ e λ θ e θ λ 1 c 2 λ e λ θ ( Λ + 1 ) c
and
F ( θ ) = e λ θ 1 c + Λ 1 + c e λ θ c 2 ,
respectively, for 0 θ < 2 π , λ > 0 , | Λ | 1 , and c = e 2 π λ 1 . The nth trigonometric moment is
m n = λ ( λ + i n ) ( 2 Λ + c + c Λ ) ( c + 1 ) e 2 π n i 1 c 2 λ 2 + n 2 2 λ Λ ( 2 λ + i n ) ( c + 1 ) 2 e 2 i π n 1 c 2 4 λ 2 + n 2
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan c 2 Λ + 4 λ 2 Λ + 4 c λ 2 c Λ + 2 c λ 2 Λ λ 4 c λ 2 6 Λ + c 3 c Λ ,
ρ = λ 2 2 Λ c + c Λ 2 + 4 c 2 λ 4 c 2 4 λ 2 + 1 λ 2 + 1 ,
γ 1 = λ ( λ sin 2 μ 2 cos 2 μ ) ( 2 Λ + c + c Λ ) ) c λ 2 + 4 λ Λ ( c + 2 ) ( cos 2 μ λ sin 2 μ ) c λ 2 + 1 1 λ 2 2 Λ c + c Λ 2 + 4 c 2 λ 4 c 2 4 λ 2 + 1 λ 2 + 1 3 2 ,
and
γ 2 = λ 2 2 Λ + 2 c + Λ c + 2 c ( 1 Λ ) 4 Λ c λ 4 + 5 λ 2 + 4 λ sin 2 μ + c λ 1 + λ 2 3 Λ 6 λ Λ c λ 4 + 5 λ 2 + 4 λ cos 2 μ ( 2 Λ c + Λ c ) 2 + 4 λ 4 c 2 2 c 4 4 λ 2 + 1 2 λ 2 + 1 2 1 λ 2 2 Λ c + c Λ 2 + 4 c 2 λ 4 c 2 4 λ 2 + 1 λ 2 + 1 2 ,
respectively. The transmuted wrapped exponential distribution was used on the turtle dataset from [1]. It outperformed the wrapped exponential and wrapped Lindley distributions. It offers closed-form expressions for the PDF and CDF. Additionally, it is more flexible than commonly used one-parameter distributions.

2.39. Wrapped Two-Parameter Lindley Distribution

Ref. [51] took g to be the PDF of the two-parameter Lindley distribution to obtain the wrapped two-parameter Lindley distribution. Its PDF and CDF are
f ( θ ) = ξ 2 ξ + α e ξ θ 1 + α θ 1 e 2 π ξ + 2 π α e 2 π ξ 1 e 2 π ξ 2
and
F ( θ ) = 1 1 e 2 π ξ 1 e ξ θ α ξ θ α + ξ + 2 π α ξ α + ξ 1 e ξ θ e 2 π ξ 1 e 2 π ξ 2 ,
respectively, for 0 θ < 2 π , ξ > 0 , and α > ξ . The nth trigonometric moment is
m n = ξ 2 ( ξ + α i n ) ( ξ + α ) ( ξ i n )
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 2 arctan 1 ξ arctan 1 ξ + α ,
ρ = ξ 2 ( ξ + α ) 2 + 1 1 2 ( ξ + α ) ξ 2 + 1 ,
γ 1 = ξ 2 ( α + ξ ) 2 + 4 ( α + ξ ) 4 + ξ 2 sin μ 2 2 μ 1 ξ 2 ( α + ξ ) 2 + 1 ( α + ξ ) 1 + ξ 2 3 2 ,
and
γ 2 = ξ 2 ( α + ξ ) 2 + 4 ( α + ξ ) 4 + ξ 2 cos μ 2 2 μ ξ 2 ( α + ξ ) 2 + 1 ( α + ξ ) 1 + ξ 2 4 1 ξ 2 ( α + ξ ) 2 + 1 ( α + ξ ) 1 + ξ 2 2 ,
respectively, where μ 2 = 2 arctan 2 ξ arctan 2 α + ξ . The two-parameter wrapped Lindley distribution was used for analyzing two datasets: the feldspar lath dataset from [1], Appendix B5, and a dataset on wind directions observed at Gorleston, England, during summer Sundays in 1968 [52]. It was found that this distribution provided a better fit for the first dataset compared to the wrapped exponential distribution and the wrapped Lindley distribution. Additionally, for the second dataset, it was observed that the wrapped two-parameter Lindley distribution outperformed the wrapped Lindley distribution in terms of fitting.

2.40. Wrapped Two-Sided Lindley Distribution

Ref. [53] took g to be the PDF of the two-sided Lindley distribution to obtain the wrapped two-sided Lindley distribution. Its PDF and CDF are
f ( θ ) = e α θ 2 Λ α Λ β α 2 e 2 π α Λ β ( θ + 1 ) e 2 π α θ + 2 π 1 + β 2 e θ ( α + β ) Λ α ( 2 π + 1 θ ) e 2 π β + θ 1
and
F ( θ ) = e α ( 2 π θ ) 2 Λ α ( α + 1 ) e α ( θ + 2 π ) + ( 2 π 1 ) α 1 e α θ + α ( θ 2 π + 1 ) e 2 π α ( α θ + α + 1 ) + 1 + 1 2 Λ β e β θ β ( θ 1 ) 1 + β e 2 π β 2 π β + β + 1 + e β ( θ + 2 π ) β θ + 2 π β + β + 1 + 1 ,
respectively, for 0 θ < 2 π , α > 0 , and β > 0 , where Λ α = e 2 π α 1 2 ( α + 1 ) and Λ β = e 2 π β 1 2 ( β + 1 ) . The nth trigonometric moment is
m n = α 2 n + i α 2 2 Λ α e 2 i π n e 2 π α 1 2 π α e 2 π α i n 2 i π n + 1 + α 2 2 n + i α Λ α e 2 π α e 2 i p n 2 π + e 2 π α 1 + β 2 2 Λ β β + i n 2 1 e 2 π β 1 + e 2 π β + i n 2 π β + 2 i π n 1 + β 2 2 β + i n Λ β 1 e 2 π β + i n 1 ( 1 + 2 π ) e 2 π β
for n = 1 , 2 , . The wrapped two-sided Lindley distribution has been applied to the ant dataset discussed by [1] and was shown to provide a good fit to this dataset.

2.41. Wrapped Variance Gamma Distribution

Ref. [54] took g to be the PDF of the variance gamma distribution to obtain the wrapped variance gamma distribution. Its PDF is
f ( θ ) = γ 2 λ e β ( θ μ ) π Γ ( λ ) ( 2 α ) λ 1 2 m = e 2 π m β K λ 1 2 α | θ + 2 m π μ | | θ + 2 m π μ | λ 1 2
for 0 θ < 2 π , α > 0 , β > 0 , 0 | β | < α , and λ R . The nth trigonometric moment is
m n = e i μ n γ α 2 ( β + i n ) 2 2 λ
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ ,
ρ = γ α 2 ( β + i 2 ) 1 2 2 λ ,
γ 1 = 0 ,
and
γ 2 = γ ( α 2 ( β + 2 i ) 2 2 λ γ ( α 2 ( β + i ) 2 ) 8 λ 1 γ α 2 ( β + i ) 2 2 λ 2 ,
respectively. The wrapped variance gamma distribution was used to analyze the Black Mountain wind direction dataset from [1]. It fit the data well, but its practical utility is limited because it lacks a closed-form PDF or CDF.

2.42. Wrapped Weighted Exponential Distribution

Ref. [55] took g to be the PDF of the weighted exponential distribution to obtain the wrapped weighted exponential distribution. Its PDF is
f ( θ ) = α + 1 α λ e λ θ k = 0 e 2 k π λ 1 e α λ ( θ + 2 k π )
for 0 θ < 2 π , α > 0 , and λ > 0 . The nth trigonometric moment is
m n = α n + i β n ,
where
α n = λ 2 cos arctan n λ 2 + α 1 + α λ 2 + n 2 λ 2 + n 2 ( 1 + α ) 2
and
β n = λ 2 sin arctan n λ 2 + α 1 + α λ 2 + n 2 λ 2 + n 2 ( 1 + α ) 2 .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan 2 + α λ ( 1 + α ) ,
ρ = λ 2 λ 2 + 1 λ 2 + 1 ( 1 + α ) 2 ,
γ 1 = λ 2 sin arctan 2 λ 2 + α 1 + α 2 arctan 1 λ 2 + α 1 + α λ 2 + 4 λ 2 + 4 ( 1 + α ) 2 ( 1 ρ ) 3 2 ,
and
γ 2 = λ 2 cos arctan 2 λ 2 + α 1 + α 2 arctan 1 λ 2 + α 1 + α ρ 4 λ 2 + 4 λ 2 + 4 ( 1 + α ) 2 ( 1 ρ ) 2 ,
respectively. The wrapped weighted exponential distribution is a flexible model with a PDF that can be recast into a closed form.

2.43. Wrapped Weibull Distribution

Ref. [40] took g to be the PDF of the Weibull distribution to obtain the wrapped Weibull distribution. Its PDF is
f ( θ ) = k = 0 1 2 σ sech 2 θ + 2 k π μ 2 σ
for 0 θ < 2 π and c > 0 . The nth trigonometric moment is
m n = k = 0 i k n k k ! Γ 1 + k c = b n + i c n ,
where
b n = k = 0 ( 1 ) k t 2 k ( 2 k ) ! Γ 1 + 2 k c
and
c n = k = 0 ( 1 ) k t 2 k + 1 ( 2 k + 1 ) ! Γ 1 + 2 k + 1 c .
The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan c 1 b 1 ,
ρ = 2 b 1 cos μ c 1 sin μ 2 + b 1 sin μ + c 1 cos μ 2 ,
γ 1 = b 1 2 + c 1 2 sin μ 2 2 μ 1 b 1 2 + c 1 2 3 2 ,
and
γ 2 = b 1 2 + c 1 2 cos μ 2 2 μ b 1 2 + c 1 2 4 1 b 1 2 + c 1 2 2 ,
respectively, where μ 2 = arctan c 2 b 2 . The wrapped Weibull distribution is limited by the fact that it does not admit a closed-form PDF.

2.44. Wrapped XGamma Distribution

Ref. [56] took g to be the PDF of the XGamma distribution to obtain the wrapped XGamma distribution. Its PDF and CDF are
f ( θ ) = λ 2 e θ λ ( λ + 1 ) 1 e 2 π λ 1 + λ θ 2 2 + 2 π λ ( π θ ) e 2 π λ + ( θ + π ) e 2 π λ 1 e 2 π λ 2
and
F ( θ ) = 1 1 + λ ( 1 + θ ) + θ 2 λ 2 2 λ + 1 e θ λ 1 1 e 2 π λ + 2 π λ λ + 1 1 ( 1 + θ λ ) e θ λ e 2 π λ 1 e 2 π λ 2 + 2 π 2 λ 2 λ + 1 1 e θ λ e 2 π λ 1 + e 2 π λ 1 e 2 π λ 3 ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is
m n = λ 2 λ 2 + λ ( 1 2 i n ) n 2 ( λ + 1 ) λ i n 3
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = 3 arctan 1 λ arctan 2 λ λ 2 + λ 1 ,
ρ = λ 2 λ + 1 λ 2 + λ 1 2 + 4 λ 2 λ 2 + 1 3 ,
γ 1 = λ 2 λ + 1 λ 2 + λ 4 2 + 16 λ 2 λ 2 + 4 2 sin κ λ , 2 1 λ 2 λ + 1 λ 2 + λ 1 2 + 4 λ 2 λ 2 + 1 3 3 2 ,
and
γ 2 = λ 2 λ + 1 λ 2 + λ 4 2 + 16 λ 2 λ 2 + 4 2 cos κ λ , 2 λ 2 1 + λ λ 2 + λ 1 2 + 4 λ 2 λ 2 + 1 3 4 1 λ 2 λ + 1 λ 2 + λ 1 2 + 4 λ 2 λ 2 + 1 3 2 ,
respectively, where κ λ , 2 = arctan 2 λ 6 arctan 1 λ + 2 arctan 2 λ λ 2 + λ 1 arctan 4 λ λ 2 + λ 4 . The wrapped XGamma distribution was used to analyze the feldspar lath dataset from [1]. It demonstrated superior fit compared to the wrapped exponential and the wrapped Lindley distributions.

2.45. Wrapped XLindley Distribution

Ref. [57] took g to be the PDF of the XLindley distribution to obtain the wrapped XLindley distribution. Its PDF and CDF are
f ( θ ) = λ 2 e θ λ ( λ + 1 ) 2 1 e 2 π λ 2 1 e 2 π λ ( λ + θ + 2 ) + 2 π e 2 π λ
and
F ( θ ) = 1 e θ λ 1 + λ θ ( λ + 1 ) 2 1 e 2 π λ + 2 π λ e 2 π λ 1 e θ λ ( λ + 1 ) 2 1 e 2 π λ 2 ,
respectively, for 0 θ < 2 π and λ > 0 . The nth trigonometric moment is
m n = λ 2 ( λ + 1 ) ( λ i n ) + λ 2 ( 1 + λ i n ) ( λ + 1 ) 2 ( λ i n ) 2
for n = 1 , 2 , . The mean direction, mean resultant length, skewness and kurtosis are
μ = 2 arctan 1 λ arctan λ + 2 ( λ + 1 ) 2 ,
ρ = λ 2 ( λ + 1 ) 4 + ( λ + 2 ) 2 ( λ + 1 ) 2 λ 2 + 1 ,
γ 1 = λ 2 4 λ 2 + 2 2 + ( λ + 1 ) 4 sin μ 2 2 μ ( λ + 1 ) 2 λ 2 + 4 1 λ 2 ( λ + 1 ) 4 + ( λ + 2 ) 2 ( λ + 1 ) 2 λ 2 + 1 3 2 ,
and
γ 2 = λ 2 4 λ 2 + 2 2 + ( λ + 1 ) 4 cos μ 2 2 μ ( λ + 1 ) 2 λ 2 + 4 λ 8 λ + 2 2 + ( λ + 1 ) 4 2 ( λ + 1 ) 8 λ 2 + 1 4 1 λ 2 ( λ + 1 ) 4 + ( λ + 2 ) 2 ( λ + 1 ) 2 λ 2 + 1 2 ,
respectively, where μ 2 = 2 arctan 2 λ arctan 2 ( λ + 2 ) ( λ + 1 ) 2 . The wrapped XLindley distribution was used to analyze two datasets: one consisted of sun directions recorded from 50 starhead topminnows under overcast conditions [1], and the other included 349 transactions occurring between 1 January 2020 and 29 July 2020 [58]. Although the wrapped XLindley distribution outperformed some distributions based on certain goodness-of-fit measures, it was not consistently superior in all cases.

3. A Review of Discrete Wrapped Distributions

In this section, we review wrapped binomial, wrapped discrete Cauchy, wrapped discrete exponential, wrapped discrete Mittag–Leffler, wrapped discrete skew Laplace, wrapped geometric, wrapped negative binomial, wrapped Poisson, wrapped Poisson–Lindley, and wrapped zero-inflated Poisson distributions.

3.1. Wrapped Binomial Distribution

Ref. [59] took g to be the PDF of the binomial distribution to obtain the wrapped binomial distribution. Its PMF is
P Θ = 2 π r m = k = 0 n r m n r + k m p 1 r + k m q 1 n r k m
for 0 θ < 2 π , r = 0 , 1 , , m 1 , m 1 , n m 1 , 0 < p 1 < 1 , and q 1 = 1 p 1 . The nth trigonometric moment is
m n = q 1 + p 1 cos 2 n π m + i p 1 sin 2 n π m n
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = n arctan p 1 sin 2 π m q 1 + p 1 sin 2 π m ,
ρ = p 1 2 + q 1 2 + 2 p 1 q 1 sin 2 π m n 2 ,
γ 1 = β 2 cos ( 2 μ ) α 2 sin ( 2 μ ) ( 1 ρ ) 3 2 ,
and
γ 2 = α 2 cos ( 2 μ ) + β 2 sin ( 2 μ ) ρ 4 ( 1 ρ ) 2 ,
respectively, where
α 2 = p 1 2 + q 1 2 + 2 p 1 q 1 cos 4 π m n 2 cos 2 arctan p 1 sin 4 π m q 1 + p 1 cos 4 π m
and
β 2 = p 1 2 + q 1 2 + 2 p 1 q 1 cos 4 π m n 2 sin 2 arctan p 1 sin 4 π m q 1 + p 1 cos 4 π m .
The wrapped binomial distribution faces practical limitations due to its non-closed-form PMF.

3.2. Wrapped Discrete Cauchy Distribution

Ref. [60] introduced the wrapped discrete Cauchy distribution. Its PMF is
P Θ = 2 π r m = 1 a 2 1 + a 2 m 2 a m cos ( m μ ) m 1 a 2 m 1 + a 2 2 a cos ( θ μ )
for 0 < a < 1 , 0 θ < 2 π , r = 0 , 1 , , m 1 , and m is a positive integer. The nth trigonometric moment is
m n = 1 a 2 1 + a 2 m 2 a m cos ( m μ ) m 1 a 2 m r = 0 m 1 e n θ 1 + a 2 2 a cos ( θ μ )
for n = 1 , 2 , . The wrapped discrete Cauchy distribution has the advantage of having a closed-form PMF.

3.3. Wrapped Discrete Exponential Distribution

Ref. [61] introduced the wrapped discrete exponential distribution. Its PMF and CDF are
P Θ = 2 π r m = e λ θ 1 e 2 π λ 1 e 2 π λ m
and
F ( Θ k ) = 1 e 2 π λ ( k + 1 ) m 1 e 2 π λ , k = 0 , 1 , , m 1 ,
respectively, for λ > 0 , r = 0 , 1 , , m 1 , and m is a positive integer. The nth trigonometric moment is
m n = 1 e 2 π λ m 1 e 2 π λ 1 e 2 π λ + 2 π i n 1 e 2 π λ m + 2 π i n m
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan e 2 π λ m sin 2 π m e 4 π λ m sin 2 π m 1 e 2 π λ m cos 2 π m e 2 π λ m + e 4 π λ m cos 2 π m ,
ρ = 1 e 2 π λ m cos 2 π m e 2 π λ m + e 4 π λ m cos 2 π m 2 + e 2 π λ m sin 2 π m e 4 π λ m sin 2 π m 2 1 + e 4 π λ m 2 e 2 π λ m cos 2 π m 2 ,
γ 1 = 1 e 2 π λ m cos 4 π m e 2 π λ m + e 4 π λ m cos 4 π m 2 + e 2 π λ m sin 4 π m e 4 π λ m sin 4 π m 2 1 + e 4 π λ m 2 e 2 π λ m cos 4 π m 2 sin μ 2 2 μ 1 1 e 2 π λ m cos 2 π m e 2 π λ m + e 4 π λ m cos 2 π m 2 + e 2 π λ m sin 2 π m e 4 π λ m sin 2 π m 2 1 + e 4 π λ m 2 e 2 π λ m cos 2 π m 2 3 2 ,
and
γ 2 = 1 e 2 π λ m cos 4 π m e 2 π λ m + e 4 π λ m cos 4 π m 2 + e 2 π λ m sin 4 π m e 4 π λ m sin 4 π m 2 1 + e 4 π λ m 2 e 2 π λ m cos 4 π m 2 cos μ 2 2 μ 1 1 e 2 π λ m cos 2 π m e 2 π λ m + e 4 π λ m cos 2 π m 2 + e 2 π λ m sin 2 π m e 4 π λ m sin 2 π m 2 1 + e 4 π λ m 2 e 2 π λ m cos 2 π m 2 2 1 e 2 π λ m cos 2 π m e 2 π λ m + e 4 π λ m cos 2 π m 2 + e 2 π λ m sin 2 π m e 4 π λ m sin 2 π m 2 1 + e 4 π λ m 2 e 2 π λ m cos 2 π m 2 2 1 1 e 2 π λ m cos 2 π m e 2 π λ m + e 4 π λ m cos 2 π m 2 + e 2 π λ m sin 2 π m e 4 π λ m sin 2 π m 2 1 + e 4 π λ m 2 e 2 π λ m cos 2 π m 2 2 ,
respectively. The wrapped discrete exponential distribution is a simple model with a closed-form PMF. For this reason, some may consider it more practical than competing distributions.

3.4. Wrapped Discrete Mittag–Leffler Distribution

Ref. [62] introduced the discrete Mittag–Leffler distribution. Its nth trigonometric moment is given by
m n = 1 + 1 δ δ 1 e 2 π i n m α 1
for n = 1 , 2 , , where 0 < δ < 1 , 0 < α < 1 , and m = 1 , 2 , . The wrapped discrete Mittag–Leffler distribution has the disadvantage of having a non-closed-form PMF.

3.5. Wrapped Discrete Skew Laplace Distribution

Ref. [63] obtained the wrapped discrete skew Laplace distribution. Its PMF and CDF are
P Θ = 2 π r m = ( 1 p ) ( 1 q ) 1 p q q m r 1 p m + p r 1 q m 1 p m 1 q m
and
F Θ s = ( 1 p ) ( 1 q ) ( 1 p q ) 1 p m 1 q m 1 p m q m s m q q s + 1 ( p s ) m 1 q m 1 p s , s = 0 , 1 , , m 1 ,
respectively, for m N , r = 0 , 1 , , m 1 , and λ > 0 . The nth trigonometric moment is given by
m n = ( 1 p ) ( 1 q ) 1 p e i 2 π n m 1 q e 2 i π n m
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan ( p q ) sin 2 π m 1 + p q ( p + q ) cos 2 π m ,
ρ = ( 1 p ) ( 1 q ) 1 + p q ( p + q ) cos 2 π m 2 + ( p q ) sin 2 π m 2 ,
γ 1 = ( 1 p ) ( 1 q ) 1 + p q ( p + q ) cos 4 π m 2 + ( p q ) sin 4 π m 2 sin μ 2 2 μ 1 ( 1 p ) ( 1 q ) 1 + p q ( p + q ) cos 2 π m 2 + ( p q ) sin 2 π m 2 3 2 ,
and
γ 2 = ( 1 p ) ( 1 q ) 1 + p q ( p + q ) cos 4 π m 2 + ( p q ) sin 4 π m 2 cos μ 2 2 μ ( ( 1 p ) ( 1 q ) ) 2 1 + p q ( p + q ) cos 2 π n m 2 + ( p q ) sin 2 π n m 2 1 ( 1 p ) ( 1 q ) 1 + p q ( p + q ) cos 2 π m 2 + ( p q ) sin 2 π m 2 2 ,
respectively, where μ 2 = arctan ( p q ) sin 4 π m 1 + p q ( p + q ) cos 4 π m . The wrapped discrete skew Laplace distribution may be considered more practical than certain other models due to it having a closed-form PMF.

3.6. Wrapped Geometric Distribution

Ref. [62] obtained the wrapped geometric distribution. Its PMF and CDF are
P Θ = 2 π r m = δ ( 1 δ ) r 1 ( 1 δ ) m
and
F Θ y = r = 0 y δ ( 1 δ ) r 1 ( 1 δ ) m , y = 0 , 1 , , m 1 ,
respectively, for m N , r = 0 , 1 , , m 1 , and δ > 0 . The nth trigonometric moment is
m n = δ x i y ,
where x = 1 ( 1 δ ) cos 2 π n m and y = ( 1 δ ) sin 2 π n m for n 0 (mod m), n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan ( 1 δ ) sin 2 π m 1 ( 1 δ ) cos 2 π m ,
ρ = δ δ 2 + 2 ( 1 δ ) 1 cos 2 π m ,
γ 1 = δ sin μ 2 2 μ 1 ( 1 δ ) cos 4 π m 2 + ( 1 δ ) sin 4 π m 2 1 δ δ 2 + 2 ( 1 δ ) 1 cos 2 π m 3 2 ,
and
γ 2 = δ cos μ 2 2 μ 1 ( 1 δ ) cos 4 π m 2 + ( 1 δ ) sin 4 π m 2 δ 4 δ 2 + 2 ( 1 δ ) 1 cos 2 π m 2 1 δ δ 2 + 2 ( 1 δ ) 1 cos 2 π m 2 ,
respectively, where μ 2 = arctan ( 1 δ ) sin 4 π m 1 ( 1 δ ) cos 4 π m . The wrapped geometric distribution is a simple discrete circular model with the practical advantage of having a closed-form PMF.

3.7. Wrapped Negative Binomial Distribution

Ref. [64] introduced the wrapped negative binomial distribution. Its PMF and CDF are
P Θ = 2 π r m = k = 0 r + k m + n 1 n 1 p 1 n q 1 r + k m
and
F Θ y = r = 0 y r + k m + n 1 n 1 p 1 n q 1 r + k m , y = 0 , 1 , , m 1 ,
respectively, for p 1 [ 0 , 1 ] , p 1 + q 1 = 1 , and n Z + . The nth trigonometric moment is
m n = p 1 1 q 1 cos 2 π n m + i p 1 q 1 sin 2 π n m 1 + q 1 2 2 q 1 cos 2 π n m n
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = n arctan p 1 q 1 2 π m p 1 1 q 1 cos 2 π m ,
ρ = p 1 1 + q 1 2 2 q 1 cos 2 π m n ,
γ 1 = p 1 1 + q 1 2 2 q 1 cos 4 π m n sin μ 2 2 μ 1 p 1 1 + q 1 2 2 q 1 cos 2 π m n 3 2 ,
and
γ 2 = p 1 1 + q 1 2 2 q 1 cos 4 π m n cos μ 2 2 μ p 1 1 + q 1 2 2 q 1 cos 4 π m 4 n 1 p 1 1 + q 1 2 2 q 1 cos 2 π m n 2 ,
respectively, where μ 2 = n arctan p 1 q 1 sin 4 π m p 1 1 q 1 cos 4 π m . The wrapped negative binomial distribution does not have a closed-form PMF.

3.8. Wrapped Poisson Distribution

Ref. [39] obtained the wrapped Poisson distribution. Its PMF and CDF are
P Θ = 2 π r m = k = 0 e λ λ r + k m ( r + k m ) !
and
F Θ y = r = 0 y P θ = 2 π r m , y = 0 , 1 , , m 1 ,
respectively, for m N , r = 0 , 1 , , m 1 , and λ > 0 . The nth trigonometric moment is
m n = e λ 1 cos 2 π n m + i   sin λ 2 π n m
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = λ sin 2 π m ,
ρ = e λ 1 cos 2 π m ,
γ 1 = e λ 1 cos 4 π m sin λ sin 4 π m 2 μ 1 e λ 1 cos 2 π m 3 2 ,
and
γ 2 = e λ 1 cos 4 π m cos λ sin 4 π m 2 μ e 4 λ 1 cos 2 π m 1 e λ 1 cos 2 π m 2 ,
respectively. The wrapped Poisson distribution lacks a closed-form PMF, limiting its practical usefulness.

3.9. Wrapped Poisson–Lindley Distribution

Ref. [65] obtained the wrapped Poisson–Lindley distribution. Its PMF and CDF are
P Θ = 2 π r m = θ 2 ( 1 + θ ) m ( r + θ + 2 ) ( 1 + θ ) m 1 + m ( 1 + θ ) r + 3 ( 1 + θ ) m 1 2
and
F Θ y = ( 1 + θ ) m y 3 ( 1 + θ ) m 1 2 [ 1 ( 1 + θ ) m θ m + ( 3 + y + θ ) ( 1 + θ ) m 1 + ( 1 + θ ) 1 + y 1 2 θ + m θ θ 2 + ( 1 + θ ) 2 + m ] , y = 0 , 1 , , m 1 ,
respectively, for m N and r = 0 , 1 , , m 1 . The wrapped Poisson–Lindley distribution was used to analyze the turtle dataset and the Gorleston wind direction dataset discussed by [1]. Additionally, it was applied to a dataset consisting of arrival directions of low showers of cosmic rays, with declination and right ascension as the coordinate system [66]. The goodness of fit of the wrapped Poisson–Lindley distribution was compared to that of the wrapped geometric distribution. The wrapped Poisson–Lindley distribution showed the best fit for the first and third datasets, but not for the second dataset.

3.10. Wrapped Zero-Inflated Poisson Distribution

Ref. [67] introduced the wrapped zero-inflated Poisson distribution. Its PMF and CDF are
P Θ = 2 π r m = k = p ( r + k m )
and
F Θ y = r = 0 y P θ = 2 π r m , y = 0 , 1 , , m 1 ,
respectively, for m N , r = 0 , 1 , , m 1 , λ > 0 , and w [ 0 , 1 ] , where
p ( x ) = w + ( 1 w ) e λ , if x = 0 , ( 1 w ) e λ λ x x ! , if x = 1 , 2 , .
The nth trigonometric moment is
m n = w + ( 1 w ) e λ + λ e i n
for n = 1 , 2 , . The mean direction, mean resultant length, skewness, and kurtosis are
μ = arctan ( 1 w ) e λ 1 cos 2 π m sin λ sin 2 π m w + ( 1 w ) e λ 1 cos 2 π m cos λ sin 2 π m ,
ρ = w 2 + ( 1 w ) 2 e 2 λ 1 cos 2 π p m + 2 w ( 1 w ) e λ 1 cos 2 π p m cos λ sin 2 π p m ,
γ 1 = w sin 2 μ + ( 1 w ) e λ 1 cos 4 π m sin λ sin 4 π m 2 μ 1 a 1 2 + b 1 2 + 2 a 1 b 1 cos λ sin 2 π m 3 2 ,
and
γ 2 = w cos 2 μ + ( 1 w ) e λ 1 cos 4 π m cos λ sin 4 π m 2 μ a 1 2 + b 1 2 + 2 a 1 b 1 cos λ sin 2 π m 2 1 a 1 2 + b 1 2 + 2 a 1 b 1 cos λ sin 2 π m 2 ,
respectively, where a 1 = w and b 1 = ( 1 w ) e λ 1 cos 2 π m . The wrapped zero-inflated Poisson distribution, a more flexible alternative to the one-parameter wrapped Poisson distribution, is useful for directional data with a high number of zero counts.

4. Data Applications

All of the distributions reviewed in Section 2 and Section 3 can be fit by using the package Wrapped due to [9]. The package computes the PDF, CDF, quantile function, random samples, maximum likelihood estimates, standard errors, confidence intervals, and measures of goodness of fit associated with (1)–(2) for any specified g. In this section, we illustrate three data applications by fitting all of the reviewed distributions using Wrapped. In each data application, the following distributions gave the best fits: wrapped exponential, wrapped gamma, wrapped Weibull, wrapped Pareto, wrapped normal, wrapped Cauchy, wrapped Laplace and, wrapped t distributions. Hence, the tables and figures will be limited to these distributions.

4.1. Dataset 1

This data taken from [1] contain arrival times on a 24 h clock of 254 patients at an intensive care unit over a period of about 12 months. The maximized log-likelihood values as well as the corresponding values of AIC and BIC are given in Table 1.
We see that the wrapped gamma distribution gives the best fit in terms of AIC and BIC. The wrapped exponential distribution gives the worst fit in terms of AIC and BIC. These findings are confirmed by the density plots shown in Figure 1.

4.2. Dataset 2

These data taken from [1] contain measurements of the directions taken by 76 turtles after treatment. The maximized log-likelihood values as well as the corresponding values of AIC and BIC are given in Table 2.
We see that the wrapped Cauchy distribution gives the best fit in terms of AIC and BIC. The wrapped Pareto distribution gives the worst fit in terms of AIC and BIC. These findings are confirmed by the density plots shown in Figure 2.

4.3. Dataset 3

These data taken from [1] contain measurements of long-axis orientation of 133 feldspar laths in basalt. The maximized log-likelihood values as well as the corresponding values of AIC and BIC are given in Table 3.
We can see that the wrapped normal distribution gives the best fit in terms of AIC, whereas the wrapped exponential distribution gives the best fit in terms of BIC. The wrapped Pareto distribution gives the worst fit in terms of AIC and BIC. These findings are confirmed by the density plots shown in Figure 3.

5. Conclusions

In this paper, we have reviewed 45 wrapped distributions for continuous circular data and 10 wrapped distributions for discrete circular data by listing their probability density/mass functions, cumulative distribution functions, trigonometric moments, mean directions, mean resultant lengths, skewness, and kurtosis (whenever they are available). We have also discussed data applications and limitations of the reviewed distributions. This paper could be a source of reference and may encourage further developments in the area of wrapped distributions.
Future work may provide similar reviews for wrapped bivariate distributions, wrapped multivariate distributions, wrapped matrix variate distributions, and wrapped complex variate distributions. Another future possibility is to write R packages for fitting wrapped bivariate distributions, wrapped multivariate distributions, wrapped matrix variate distributions, and wrapped complex variate distributions.

Author Contributions

Investigation, W.B. and S.N.; methodology, W.B. and S.N.; formal analysis, W.B. and S.N.; software, W.B. and S.N.; writing—review and editing, W.B. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, and further inquiries can be directed to the corresponding author.

Acknowledgments

The authors would like to thank the editor and the two referees for careful reading and comments, which greatly improved the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fisher, N.I. Statistical Analysis of Circular Data; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  2. Fisher, N.I. Statistical Analysis of Circular Data; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  3. Jones, M.C. On families of distributions with shape parameters. Int. Stat. Rev. 2015, 83, 175–192. [Google Scholar] [CrossRef]
  4. Ley, C.; Verdebout, T. Modern Directional Statistics; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  5. Pewsey, A.; Garcia-Portugues, E. Recent advances in directional statistics. Test 2021, 30, 1–58. [Google Scholar] [CrossRef]
  6. Ley, C.; Babic, S.; Craens, D. Flexible models for complex data with applications. Annu. Rev. Stat. Its Appl. 2021, 8, 369–391. [Google Scholar] [CrossRef]
  7. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1986; Volumes 1–3. [Google Scholar]
  8. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 6th ed.; Academic Press: San Diego, CA, USA, 2000. [Google Scholar]
  9. Nadarajah, S.; Zhang, Y. Wrapped: An R package for circular data. PLoS ONE 2017, 12, e0188512. [Google Scholar] [CrossRef] [PubMed]
  10. Kumar, C.S.; Alex, R. Wrapped generalized skew normal distribution. Stoch. Model. Appl. 2021, 25, 1–7. [Google Scholar]
  11. Hernandez-Sanchez, E.; Scarpa, B. A wrapped flexible generalized skew-normal model for a bimodal circular distribution of wind directions. Chil. J. Stat. 2012, 3, 129–141. [Google Scholar]
  12. Holla, M.S.; Bhattacharya, S.K. On a compound Gaussian distribution. Ann. Inst. Stat. Math. 1968, 20, 331–336. [Google Scholar] [CrossRef]
  13. Al-Khazaleh, A.; Al-Meanazel, T.R. Wrapped Akash distribution. Electron. J. Appl. Stat. Anal. 2021, 14, 305–317. [Google Scholar]
  14. Ganaie, R.; Rajagopalan, V.; Rather, A. Weighted Aradhana distribution: Properties and applications. J. Inf. Technol. Res. 2019, IX, 392–406. [Google Scholar]
  15. Ramabhadra Sarma, I.; Dattatreya Rao, A.V.; Girija, S.V.S. On characteristic functions of wrapped half logistic and binormal distributions. Int. J. Stat. Syst. 2009, 4, 33–45. [Google Scholar]
  16. Saraiva, M.A.; de Azevedo Cysneiros, F.J. Wrapped Birnbaum-Saunders distribution: Definition and estimation. Nexus Math. 2022, 5, 200010. [Google Scholar] [CrossRef]
  17. Jander, R. Die optische Richtungsorientierung derroten Waldameise (Formica rufa). Z. Fur-Vergle-Ichende Physiol. 1957, 40, 162–238. [Google Scholar] [CrossRef]
  18. Levy, P. L’addition des variables aléatoires définies sur une circonférence. Bull. SociéTé MathéMatique Fr. 1939, 67, 1–41. [Google Scholar] [CrossRef]
  19. Roy, S.; Adnan, M.A.S. Wrapped Chi-Square Distribution. Available online: https://www.researchgate.net/publication/263080292_Wrapped_Chi-square_distribution (accessed on 3 August 2024).
  20. Rao, J.S.; Kozubowski, T.J. A wrapped exponential circular model. Proc. AP Acad. Sci. 2001, 5, 43–56. [Google Scholar]
  21. Subrahmanyam, P.S.; Dattatreya Rao, A.V.; Girija, S.V.S. On wrapping of exponentiated inverted Weibull distribution. Int. J. Innov. Res. Sci. Technol. 2017, 3, 18–25. [Google Scholar]
  22. Coelho, C.A. The wrapped gamma distribution and wrapped sums and linear combinations of independent gamma and Laplace distributions. J. Stat. Theory Pract. 2007, 1, 1–29. [Google Scholar] [CrossRef]
  23. Jayakumar, K.; Sajayan, T. Wrapped generalized geometric stable distributions with an application to wind direction. Far East J. Theor. Stat. 2022, 66, 147–166. [Google Scholar]
  24. Cameron, M.A. The comparison of time series recorders. Technometrics 1983, 25, 9–22. [Google Scholar] [CrossRef]
  25. Roy, S.; Adnan, M.A.S. Wrapped generalized Gompertz distribution: An application to ornithology. J. Biom. Biostat. 2012, 3, 153–156. [Google Scholar] [CrossRef]
  26. Reed, W.J.; Pewsey, A. Two nested families of skew-symmetric circular distributions. Test 2009, 18, 516–528. [Google Scholar] [CrossRef]
  27. Bruderer, B.; Jenni, L. Migration across the alps. In Bird Migration; Springer: Berlin/Heidelberg, Germany, 1990; pp. 60–77. [Google Scholar]
  28. Varghese, J.; Jose, K.K. Wrapped hb-skewed laplace distribution and its application in meteorology. Biom. Biostat. Int. J. 2018, 7, 525–530. [Google Scholar]
  29. Adnan, M.A.S.; Roy, S. Wrapped hypo-exponential distribution. J. Stat. Manag. Syst. 2013, 16, 1–11. [Google Scholar] [CrossRef]
  30. Al-Meanazel, T.R.; Al-Khazaleh, A. Wrapped Ishita distribution. J. Stat. Appl. Probab. 2021, 10, 293–299. [Google Scholar]
  31. Jammalamadaka, S.R.; Kozubowski, T. A new family of circular models: The wrapped Laplace distributions. Adv. Appl. Stat. 2003, 3, 77–103. [Google Scholar]
  32. Bhattacharjee, S.; Borah, D. Wrapped length biased weighted exponential distribution. Thail. Stat. 2019, 17, 223–234. [Google Scholar]
  33. Smith, N.M. Reconstruction of the Tertiary Drainage Systems of the Inverell Region. Unpublished B.Sc. (Honours) Thesis, Department of Geography, University of Sydney, Sydney, Australia, 1988. [Google Scholar]
  34. Joshi, S.; Jose, K.K. Wrapped Lindley distribution. Commun. Stat.-Theory Methods 2018, 47, 1013–1021. [Google Scholar] [CrossRef]
  35. Jacob, S. Study on Circular Distributions. Ph.D. Thesis, Department of Statistics, University of Calicut, Calicut, India, 2012. [Google Scholar]
  36. Subba Rao, R.; Ravindranath, V.; Dattatreya Rao, A.V.; Prasad, G.; Ravi Kishore, P. Wrapped Lomax distribution: A new circular probability model. Int. J. Eng. Technol. 2018, 7, 150–152. [Google Scholar] [CrossRef]
  37. Chesneau, C.; Tomy, L.; Jose, M. Wrapped modified Lindley distribution. J. Stat. Manag. Syst. 2021, 24, 1025–1040. [Google Scholar] [CrossRef]
  38. Subrahmanyam, P.S.; Dattatreya Rao, A.V.; Girija, S.V.S. On wrapping of new Weibull Pareto distribution. Int. J. Adv. Res. Rev. 2017, 2, 10–20. [Google Scholar]
  39. Mardia, K.V. Statistics of Directional Data. In Probability and Mathematical Statistics; Academic Press: London, UK, 1972. [Google Scholar]
  40. Jammalamadaka, S.R.; Kozubowski, T.J. A general approach for obtaining wrapped circular distributions via mixtures. Sankhyā A 2017, 79, 133–157. [Google Scholar] [CrossRef]
  41. Bhattacharjee, S. Exploration of some distributional properties, parameter estimates and applications of the wrapped quasi Lindley distribution. J. Math. Comput. Sci. 2020, 10, 2390–2407. [Google Scholar]
  42. Bell, W.; Nadarajah, S. The wrapped Rama distribution. Preprint 2023. [Google Scholar]
  43. Demirelli, A.E.; Gürcan, M. New wrapped distribution via Richard link function. AIP Conf. Proc. 2017, 1863, 350002. [Google Scholar]
  44. Al-Meanazel, T.R.; Al-Khazaleh, A. Wrapped Shanker distribution. Ital. J. Pure Appl. Math. 2021, 46, 184–194. [Google Scholar]
  45. Yilmaz, A. Wrapped flexible skew Laplace distribution. J. Turk. Stat. Assoc. 2018, 11, 53–64. [Google Scholar]
  46. Fernandez-Duran, J.J. Circular distributions based on nonnegative trigonometric sums. Biometrics 2004, 60, 499–503. [Google Scholar] [CrossRef] [PubMed]
  47. Pewsey, A. The wrapped skew-normal distribution on the circle. Commun. Stat.-Theory Methods 2000, 29, 2459–2472. [Google Scholar] [CrossRef]
  48. Pewsey, A. The wrapped stable family of distributions as a flexible model for circular data. Comput. Stat. Data Anal. 2008, 52, 1516–1523. [Google Scholar] [CrossRef]
  49. Pewsey, A.; Lewis, T.; Jones, M. The wrapped t family of circular distributions. Aust. N. Z. J. Stat. 2007, 49, 79–91. [Google Scholar]
  50. Yilmaz, A.; Bicer, C. A new wrapped exponential distribution. Math. Sci. 2018, 12, 285–293. [Google Scholar] [CrossRef]
  51. Bhattacharjee, S.; Ahmed, I.; Das, K.K. Wrapped two-parameter Lindley distribution for modelling circular data. Thail. Stat. 2020, 19, 81–94. [Google Scholar]
  52. Mardia, K.V.; Jupp, P.E. Directional Statistics; John Wiley and Sons: New York, NY, USA, 2000. [Google Scholar]
  53. Bicer, C.; Yilmaz, A. Wrapped two sided Lindley distribution. In Research and Reviews in Science and Mathematics; Iyit, N., Dogan, H.H., Akgul, H., Eds.; Gece Academy Publishing: Ankara, Turkey, 2022; pp. 1–13. [Google Scholar]
  54. Adnan, M.A.S.; Roy, S. Wrapped variance gamma distribution with an application to wind direction. J. Environ. Stat. 2014, 6, 1–10. [Google Scholar]
  55. Roy, S.; Adnan, M.A.S. Wrapped weighted exponential distributions. Stat. Probab. Lett. 2012, 82, 77–83. [Google Scholar] [CrossRef]
  56. Al-Mofleh, H.; Sen, S. The wrapped Xgamma distribution for modeling circular data appearing in geological context. arXiv 2019, arXiv:1903.00177. [Google Scholar]
  57. Zinhom, E.; Nassar, M.M.; Radwan, S.S.; Elmasry, A. The wrapped XLindley distribution. Environ. Ecol. Stat. 2023, 30, 669–686. [Google Scholar] [CrossRef]
  58. Bahnsen, A.C.; Stojanovic, A.; Aouada, D.; Ottersten, B. Improving credit card fraud detection with calibrated probabilities. In Proceedings of the 2014 SIAM International Conference on Data Mining, Philadelphia, PA, USA, 24–26 April 2014; pp. 677–685. [Google Scholar]
  59. Girija, S.V.S.; Dattatreya Rao, A.V.; Srihari, V.L.N. On wrapped binomial model characteristics. Math. Stat. 2014, 2, 231–234. [Google Scholar] [CrossRef]
  60. Girija, S.V.S.; Srihari, G.V.L.N.; Srinivas, R. On discrete wrapped Cauchy model. Math. Theory Model. 2019, 9, 15–23. [Google Scholar]
  61. Srinivas, R.; Srihari, G.V.L.N.; Girija, S.V.S. On construction of discrete wrapped Lindley distribution. i-Manag. J. Math. 2008, 7, 10–19. [Google Scholar] [CrossRef]
  62. Jacob, S.; Jayakumar, K. Wrapped geometric distribution: A new probability model for circular data. J. Stat. Theory Appl. 2013, 12, 348–355. [Google Scholar] [CrossRef]
  63. Jayakumar, K.; Jacob, S. Wrapped skew Laplace distribution on integers: A new probability model for circular data. Open J. Stat. 2012, 2, 106–114. [Google Scholar] [CrossRef]
  64. Srihari, G.V.L.N.; Girija, S.V.S.; Rao, A.V.D. On wrapped negative binomial model. Int. J. Math. Arch. 2018, 9, 148–153. [Google Scholar]
  65. Chesneau, C.; Bakouch, H.S.; Hussain, T. Modeling and applications to circular data with a wrapped Poisson-Lindley model. J. Math. Study 2022, 55, 139–157. [Google Scholar] [CrossRef]
  66. Fisher, N.I.; Lewis, T.; Embleton, B.J.J. Statistical Analysis of Spherical Data; Cambridge University Press: New York, NY, USA, 1987. [Google Scholar]
  67. Thavaneswaran, A.; Mandal, S.; Pathmanathan, D. Estimation for wrapped zero inflated Poisson and wrapped Poisson distributions. Int. J. Stat. Probab. 2016, 5, 1–8. [Google Scholar] [CrossRef]
Figure 1. Histogram of the data and fitted densities for dataset 1.
Figure 1. Histogram of the data and fitted densities for dataset 1.
Mathematics 12 02440 g001
Figure 2. Histogram of the data and fitted densities for dataset 2.
Figure 2. Histogram of the data and fitted densities for dataset 2.
Mathematics 12 02440 g002
Figure 3. Histogram of the data and fitted densities for dataset 3.
Figure 3. Histogram of the data and fitted densities for dataset 3.
Mathematics 12 02440 g003
Table 1. Fitted distributions and values of log L , AIC, and BIC for dataset 1.
Table 1. Fitted distributions and values of log L , AIC, and BIC for dataset 1.
Distribution log L AICBIC
Wrapped exponential−468.6939.2942.8
Wrapped gamma−436.9877.8884.9
Wrapped Weibull−441.0886.1893.1
Wrapped Pareto−446.1896.2903.3
Wrapped normal−438.5881.1888.1
Wrapped Cauchy−443.6891.2898.2
Wrapped Laplace−443.4890.8897.9
Wrapped t−438.5883.1893.7
Table 2. Fitted distributions and values of log L , AIC, and BIC for dataset 2.
Table 2. Fitted distributions and values of log L , AIC, and BIC for dataset 2.
Distribution log L AICBIC
Wrapped exponential−122.3246.7249.0
Wrapped gamma−119.6243.2247.9
Wrapped Weibull−120.6245.3249.9
Wrapped Pareto−135.9275.8280.5
Wrapped normal−125.1254.1258.8
Wrapped Cauchy−113.9231.9236.6
Wrapped Laplace−116.7237.3242.0
Wrapped t−113.7233.5240.5
Table 3. Fitted distributions and values of log L , AIC, and BIC for dataset 3.
Table 3. Fitted distributions and values of log L , AIC, and BIC for dataset 3.
Distribution log L AICBIC
Wrapped exponential−243.9489.8492.6
Wrapped gamma−243.9491.8497.5
Wrapped Weibull−243.9491.8497.5
Wrapped Pareto−247.6499.2505.0
Wrapped normal−242.4488.7494.5
Wrapped Cauchy−242.7489.5495.3
Wrapped Laplace−243.9491.9497.6
Wrapped t−242.4490.7499.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bell, W.; Nadarajah, S. A Review of Wrapped Distributions for Circular Data. Mathematics 2024, 12, 2440. https://doi.org/10.3390/math12162440

AMA Style

Bell W, Nadarajah S. A Review of Wrapped Distributions for Circular Data. Mathematics. 2024; 12(16):2440. https://doi.org/10.3390/math12162440

Chicago/Turabian Style

Bell, William, and Saralees Nadarajah. 2024. "A Review of Wrapped Distributions for Circular Data" Mathematics 12, no. 16: 2440. https://doi.org/10.3390/math12162440

APA Style

Bell, W., & Nadarajah, S. (2024). A Review of Wrapped Distributions for Circular Data. Mathematics, 12(16), 2440. https://doi.org/10.3390/math12162440

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