Submit manuscript...
eISSN: 2378-315X

Biometrics & Biostatistics International Journal

Research Article Volume 4 Issue 3

On poisson-amarendra distribution and its applications

Rama Shanker,1 Hagos Fesshaye2

1Department of Statistics, Eritrea Institute of Technology, Eritrea
2Department of Economics, College of Business and Economics, Eritrea

Correspondence: Rama Shanker, Department of Statistics, Eritrea Institute of Technology, Asmara, Eritrea

Received: July 19, 2016 | Published: August 27, 2016

Citation: Shanker R, Fesshaye H. On poisson-amarendra distribution and its applications. Biom Biostat Int J. 2016;4(3):118-126. DOI: 10.15406/bbij.2016.04.00099

Download PDF

Abstract

In this paper a simple method for obtaining moments of Poisson-Amarendra distribution (PAD) introduced by Shanker1 has been suggested and hence the first four moments about origin and the variance has been given. The applications and the goodness of fit of the PAD have been discussed with data-sets relating to ecology, genetics and thunderstorms and the fit is compared with Poisson distribution, Poisson-Lindley distribution (PLD) introduced by Sankaran2 and Poisson-Sujatha distribution introduced by Shanker3 and the fit of PAD shows satisfactory fit in most of data-sets.

Keywords: amarendra distribution poisson-amarendra distribution, poisson-lindley distribution, poisson-sujatha distribution, moments, estimation of parameter, goodness of fit

Introduction

The Poisson-Amarendra distribution (PAD) defined by its probability mass function (p.m.f.)

P(X=x)=θ4θ3+θ2+2θ+6x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)(θ+1)x+4P(X=x)=θ4θ3+θ2+2θ+6x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)(θ+1)x+4

                                 

;x=0,1,2,...,θ>0;x=0,1,2,...,θ>0 (1.1)

has been introduced by Shanker1 for modeling count data-sets. Shanker1 has shown that PAD is a Poisson mixture of Amarendra distribution introduced by Shanker4 when the  parameter λλ  of Poisson distribution follows Amarendra distribution having probability density function (p.d.f.)

f(λ;θ)=θ4θ3+θ2+2θ+6(1+λ+λ2+λ3)eθλ;λ>0,θ>0f(λ;θ)=θ4θ3+θ2+2θ+6(1+λ+λ2+λ3)eθλ;λ>0,θ>0               (1.2)

We have

P(X=x)=0eλλxx!θ4θ3+θ2+2θ+6(1+λ+λ2+λ3)eθλdλP(X=x)=0eλλxx!θ4θ3+θ2+2θ+6(1+λ+λ2+λ3)eθλdλ                     (1.3)

=θ4(θ3+θ2+2θ+6)x!0λx(1+λ+λ2+λ3)e(θ+1)λdλ=θ4(θ3+θ2+2θ+6)x!0λx(1+λ+λ2+λ3)e(θ+1)λdλ

=θ4θ3+θ2+2θ+6x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)(θ+1)x+4=θ4θ3+θ2+2θ+6x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)(θ+1)x+4 .                                   

                                                                    

;x=0,1,2,...,θ>0;x=0,1,2,...,θ>0 (1.4)

which is the Poisson-Amarendra distribution (PAD).

 It has been shown by Shanker4 that Amarendra distribution is a four component mixture of exponential (θ)(θ) distribution, a gamma (2,θ)(2,θ) distribution, a gamma (3,θ)  distribution  and a gamma (4,θ) distribution with their mixing proportions θ3θ3+θ2+2θ+6 , θ2θ3+θ2+2θ+6 , 2θθ3+θ2+2θ+6 , and 6θ3+θ2+2θ+6  respectively. Shanker4 has discussed its various mathematical and statistical properties including its shape for different values of its parameter, moment generating function, moments, skewness, kurtosis, hazard rate function, mean residual life function, stochastic orderings, mean deviations, Bonferroni and Lorenz curves, amongst others. Further, Shanker4 has also discussed the estimation of its parameter using maximum likelihood estimation and method of moments along with applications for modeling lifetime data and observed that it gives much closer fit than Akash, Shanker and Sujatha distributions introduced by Shanker,5-7 Lindley8 and exponential distributions. It would be worth mentioning that Shanker5-7 has proposed  Akash, Shanker and Sujatha, distributions along with their various mathematical and statistical properties to model lifetime data arising from engineering and biomedical sciences and showed that these distributions provide much closer fit than Lindley and exponential distributions.

The Poisson-Lindley distribution (PLD) defined by its p.m.f.

P(X=x)=θ2(x+θ+2)(θ+1)x+3;

    x = 0, 1, 2,…,  θ  > 0. (1.5)

has been introduced by Sankaran2 to model count data. The PLD is a Poisson mixture of Lindley8 distribution when the parameter λ of Poisson distribution follows Lindley8 distribution with its p.d.f.

f(λ,θ)=θ2θ+1(1+λ)eθλ  ;    x>0,θ>0 (1.6)

Shanker & Hagos9 has done detailed study about applications of Poisson-Lindley distribution for modeling count data from biological sciences and showed that it gives better fit than Poisson-distribution. Shanker et al.10 discussed the comparative study of zero-truncated Poisson and Poisson-Lindley distributions and observed that in majority of data sets zero-truncated Poisson-Lindley distribution gives better fit.

Shanker11 obtained Poisson-Sujatha distribution (PSD) having p.m.f.

P(X=x)=θ3θ2+θ+2x2+(θ+4)x+(θ2+3θ+4)(θ+1)x+3;x=0,1,2,...,θ>0 (1.7)

by compounding Poisson distribution with Sujatha distribution, introduced by Shanker7 having p.d.f.

f(x;θ)=θ3θ2+θ+2(1+x+x2)eθx;x>0,θ>0                              (1.8)

Sujatha distribution introduced by Shanker7 is a better model than exponential and Lindley distributions for modeling lifetime data from biomedical science and engineering. Further, Shanker & Hagos11 has detailed study about applications of Poisson-Sujatha distribution (PSD) for modeling count data from biological science and observed that it gives better fit than Poisson-Lindley (PLD) and Poisson-distribution. Shanker & Hagos12,13 have obtained the size-biased Poisson-Sujatha distribution (SBPSD) and zero-truncated Poisson-Sujatha distribution (ZTPSD) and discussed their various mathematical and statistical properties, estimation of their parameter and applications. Further, Shanker & Hagos14 have detailed study regarding applications of zero-truncated Poisson (ZTPD), zero-truncated Poisson-Lindley distribution (ZTPLD), and zero-truncated Poisson-Sujatha distribution (ZTPSD) for modeling data-sets excluding zero counts from demography and biological sciences and concluded that in majority of data-sets ZTPSD gives better fit than ZTPD and ZTPLD.

In this paper a simple method for obtaining moments of Poisson-Amarendra distribution (PAD) introduced by Shanker1 has been suggested and hence the first four moments about origin and the variance has been presented. It seems that not much work has been done on the applications of PAD so far.  The PAD has been fitted to some data-sets relating to ecology ,genetics and thunderstorms and its goodness of fit has been compared with Poisson distribution (PD), Poisson-Lindley distribution (PLD), Poisson-Sujatha distribution (PSD).

Moments of pad

Using (1.3), the r th moment about origin of PAD (1.1) can be obtained as

μr=E[E(Xr|λ)]=θ4θ3+θ2+2θ+60[x=0xreλλxx!](1+λ+λ2+λ3)eθλdλ (2.1)

clearly that the expression under the bracket in (2.1) is the r th moment about origin of the Poisson distribution. Taking r=1  in (2.1) and using the mean of the Poisson distribution, the mean of the PAD (1.1) can be obtained as

μ1=θ4θ3+θ2+2θ+60λ(1+λ+λ2+λ3)eθλdλ=θ3+2θ2+6θ+24θ(θ3+θ2+2θ+6)

Again taking   r=2  in (2.1) and using the second moment about origin of the Poisson distribution, the second moment about origin of the PAD (1.1) can be obtained as

μ2=θ4θ3+θ2+2θ+60(λ2+λ)(1+λ+λ2+λ3)eθλdλ=θ4+4θ3+12θ2+48θ+120θ2(θ3+θ2+2θ+6)  

Similarly, taking r=3and4 in (2.1) and using the third and fourth moments about origin of the Poisson distribution, the third and the fourth moments about origin of the PAD (1.1) are obtained as

μ3=θ5+8θ4+30θ3+120θ2+480θ+720θ3(θ3+θ2+2θ+6)                                                                     

μ4=θ6+16θ5+84θ4+360θ3+1680θ2+5040θ+5040θ4(θ3+θ2+2θ+6)

The variance of PAD (1.1) can thus be obtained as

μ2=θ7+4θ6+14θ5+58θ4+144θ3+156θ2+240θ+144θ2(θ3+θ2+2θ+6)2

Parameter estimation

Maximum likelihood estimate (MLE) of the parameter: Let (x1,x2,...,xn) be a random sample of size n from the PAD (1.1) and let fx be the observed frequency in the sample corresponding to X=x(x=1,2,3,...,k) such that kx=1fx=n , where k is the largest observed value having non-zero frequency. The likelihood function L of the PAD (1.1) is given by

L=(θ4θ3+θ2+2θ+6)n1(θ+1)kx=1fx(x+4)kx=1[x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)]fx

The log likelihood function is thus obtained as

logL=nlog(θ4θ3+θ2+2θ+6)kx=1fx(x+4)log(θ+1)+kx=1fxlog[x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)]

The first derivative of the log likelihood function is given by

dlogLdθ=4nθn(3θ2+2θ+2)(θ3+θ2+2θ+6)n(ˉx+4)θ+1+kx=1[x2+(2θ+5)x+(3θ2+8θ+7)]fx[x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)]

where ˉx  is the sample mean.

The maximum likelihood estimate (MLE), ˆθ  of θ  of PAD (1.1) is the solution of dlogLdθ=0  and is the solution of the following  non-linear equation

4nθn(3θ2+2θ+2)(θ3+θ2+2θ+6)n(ˉx+4)θ+1+kx=1[x2+(2θ+5)x+(3θ2+8θ+7)]fx[x3+(θ+7)x2+(θ2+5θ+15)x+(θ3+4θ2+7θ+10)]=0           

This non-linear equation can be solved using any numerical iteration methods such as Newton- Raphson method, Bisection method, Regula –Falsi method etc. In this paper, Newton-Raphson method has been used to solve the above equation for estimating the parameter θ .

Method of moment estimate (MOME) of the parameter: Let (x1,x2,...,xn) be a random sample of size n  from the PAD (1.1). Equating the population mean to the corresponding sample mean, the MOME ˜θ of θ  of PAD (1.1) is the solution of the following biquadratic equation

ˉxθ4+(ˉx1)θ3+2(ˉx1)θ2+6(ˉx1)θ24=0

where ˉx is the sample mean.

Goodness of fit and applications

The Poisson distribution is a suitable statistical model for the situations where events are independent and mean equals variance, which is unrealistic in most of data sets in biological science and thunderstorms. Further, the negative binomial distribution is a possible alternative to the Poisson distribution when successive events are possibly dependent Johnson et al.15 but for fitting negative binomial distribution (NBD) to the count data, mean must be less than the variance. In biological science and thunderstorms, these conditions are not fully satisfied. Generally, the count data in biological science and thunderstorms are either over-dispersed or under-dispersed. The main reason for selecting PAD, PSD, and PLD to fit data from biological science and thunderstorms are that these distributions are always over-dispersed and PAD has some flexibility over PSD and PLD.

Applications in ecology

Ecology is the branch of biology which deals with the relations and interactions between organisms and their environment, including their organisms. The organisms and their environment in the nature are complex, dynamic, interdependent, mutually reactive and interrelated. Ecology deals with various principles which govern such relationship between organisms and their environment. Fisher et al.16 firstly discussed the applications of Logarithmic series distribution (LSD) to model count data in ecology. Later, Kempton17 fitted the generalized form of Fisher’s Logarithmic series distribution (LSD) to model insect data and concluded that it gives a superior fit as compared to ordinary Logarithmic series distribution (LSD). He also concluded that it gives better explanation for the data having exceptionally long tail. Tripathi & Gupta18 proposed another generalization of the Logarithmic series distribution (LSD) which is flexible to describe short-tailed as well as long-tailed data and fitted it to insect data and found that it gives better fit as compared to ordinary Logarithmic series distribution. Mishra & Shanker19 have discussed applications of generalized logarithmic series distributions (GLSD) to models data in ecology. Shanker & Hagos10 have tried to fit PLD for data relating to ecology and observed that PLD gives satisfactory fit. Shanker & Hagos11 has discussed applications of PSD to model count data from biological science and concluded that PSD gives superior fit than PLD in majority of data.

In this section an attempt has been made to fit Poisson distribution (PD), Poisson -Lindley distribution (PLD), Poisson-Sujatha distribution (PSD) and Poisson-Amarendra distribution (PAD) to many count data from biological sciences using maximum likelihood estimates. The data were on haemocytometer yeast cell counts per square, on European red mites on apple leaves and European corn borers per plant.

It is obvious from above tables that in Table 1, PD gives better fit than PLD, PSD and PAD; in Table 2 PAD gives better fit than PD, PLD, and PSD while in Table 3, PSD gives better fit than PD, PLD and PAD.

Number of yeast cells per square

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

213

202.1

234.0

233.7

233.7

1

128

138.0

99.4

99.6

98.4

2

37

47.1

40.5

41.0

41.0

3

18

10.71.80.20.1}

16.06.22.41.5}

16.36.72.30.9}

16.76.52.41.3}

4

3

5

1

6

0

Total

400.0

400.0

400.0

400.0

ML Estimate

ˆθ=0.6825

ˆθ=1.950236

ˆθ=2.373052

ˆθ=2.759978

χ2

10.08

11.04

10.86

12.01

d.f.

2

2

2

2

p-value

0.0065

0.004

0.0044

0.0025

Table 1 Observed and expected number of Haemocytometer yeast cell counts per square observed by ‘Student’ 1907

Number mites per leaf

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

38

25.3

35.8

35.3

35.3

1

17

29.1

20.7

20.9

20.8

2

10

16.7

11.4

11.6

11.7

3

9

6.41.80.40.20.1}

6.0
3.11.60.80.6}

6.1
3.11.50.70.8}

6.2
3.11.50.80.6}

4

3

5

2

6

1

7+

0

Total

80

80.0

80.0

80.0

80.0

ML Estimate

ˆθ=1.15

ˆθ=1.255891

ˆθ=1.264683

ˆθ=2.04047

χ2

18.27

2.47

2.52

2.41

d.f.

2

3

3

3

p-value

0.0001

0.4807

0.4719

0.4918

Table 2 Observed and expected number of red mites on Apple leaves

Applications in genetics

Genetics is the branch of biological science which deals with heredity and variation. Heredity includes those traits or characteristics which are transmitted from generation to generation, and is therefore fixed for a particular individual. Variation, on the other hand, is mainly of two types, namely hereditary and environmental. Hereditary variation refers to differences in inherited traits whereas environmental variations are those which are mainly due to environment. Much quantitative works seem to be done to model count data in genetics but so far no works has been done on fitting of PAD to data in genetics. The segregation of chromosomes has been studied using statistical tool, mainly chi-square ( χ2 ).  In the analysis of data observed on chemically induced chromosome aberrations in cultures of human leukocytes, Loeschke & Kohler20 suggested the negative binomial distribution while Janardan & Schaeffer21 suggested modified Poisson distribution. Mishra & Shanker19 have discussed applications of generalized Logarithmic series distributions (GLSD) to model data in mortality, ecology and genetics. Shanker & Hagos9 have detailed study on the applications of PLD to model data from genetics. Shanker & Hagos11 has detailed study on modeling of count data in genetics using PSD.  In this section an attempt has been made to fit to PAD, PSD, PLD and PD to data from genetics using maximum likelihood estimate. Also an attempt has been made to fit PAD, PSD, PLD, and PD to the data of Catcheside et al.22 in Table 5-7.

It is obvious that in Table 4 & 7, PLD gives better fit than PD, PSD and PAD; in Table 5 & 6, PAD gives better fit than PD, PLD, and PSD.

Number of bores per plant

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

188

169.4

194.0

193.6

194.2

1

83

109.8

79.5

79.6

78.6

2

36

35.6

31.3

31.6

31.6

3

14

7.81.20.2}

12.04.52.7}

12.14.52.6}

12.34.72.6}

4

2

5

1

Total

324

324.0

324.0

324.0

324.0

ML Estimate

ˆθ=0.648148

ˆθ=2.043252

ˆθ=2.471701

ˆθ=2.858180

χ2

15.19

1.29

1.16

1.4

d.f.

2

2

2

2

p-value

0.0005

0.5247

0.5599

0.4966

Table 3 Observed and expected number of European corn- borer of Mc. Guire et al23

Number of aberrations

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

268

231.3

257.0

257.6

259.0

1

87

126.7

93.4

93.0

91.9

2

26

34.7

32.8

32.7

32.5

3

9

6.30.80.10.10.1}

11.2
3.81.20.40.2}

11.2
3.71.20.40.2}

11.3
3.81.30.20.0}

4

4

5

2

6

1

7+

3

Total

400

400.0

400.0

400.0

400.0

ML Estimate

ˆθ=0.5475

ˆθ=2.380442

ˆθ=2.829241

ˆθ=3.216733

χ2

38.21

6.21

6.28

6.5

d.f.

2

3

3

3

p-value

0

0.1018

0.0987

0.0897

Table 4 Distribution of number of Chromatid aberrations (0.2 g chinon 1, 24 hours)

Class/exposure
( μg|kg )

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

413

374.0

405.7

406.1

407.5

1

124

177.4

133.6

132.9

131.2

2

42

42.1

42.6

42.7

42.5

3

15

6.60.80.10.0}

13.3
4.11.20.5}

13.4
4.11.20.6}

13.6
4.31.30.6}

4

5

5

0

6

2

Total

601

601.0

601.0

601.0

601.0

ML Estimate

ˆθ=0.47421

ˆθ=2.685373

ˆθ=3.125788

ˆθ=3.492243

χ2

48.17

1.34

1.10

0.70

d.f.

2

3

3

3

p-value

0.0000

0.7196

0.7771

0.8732

Table 5 Mammalian cytogenetic dosimetry lesions in rabbit lymphoblast induced by streptonigrin (NSC-45383), Exposure -60 μg|kg

Class/exposure
( μg|kg )

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

200

172.5

191.8

192.0

192.8

1

57

95.4

70.3

70.1

69.1

2

30

26.4

24.9

24.9

24.8

3

7

4.90.70.10.0}

8.62.91.00.5}

8.62.90.90.6}

8.73.01.00.6}

4

4

5

0

6

2

Total

300

300.0

300.0

300.0

300.0

ML Estimate

ˆθ=0.55333

ˆθ=2.353339

ˆθ=2.795745

ˆθ=3.178185

χ2

29.68

3.91

3.81

3.47

d.f.

2

2

2

2

p-value

0.0000

0.1415

0.1488

0.1764

Table 6 Mammalian cytogenetic dosimetry lesions in rabbit lymphoblast induced by streptonigrin (NSC-45383), Exposure -70 μg|kg

Class/exposure
( μg|kg )

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

155

127.8

158.3

157.5

157.9

1

83

109.0

77.2

77.5

76.8

2

33

46.5

35.9

36.4

36.5

3

14

13.22.80.50.2}

16.1
7.13.12.3}

16.4
7.13.02.1}

16.6
7.23.02.0}

4

11

5

3

6

1

Total

300

300.0

300.0

300.0

300.0

ML Estimate

ˆθ=0.853333

ˆθ=2.034077

ˆθ=2.431509

χ2

24.97

1.51

1.74

1.93

d.f.

2

3

3

3

p-value

0

0.6799

0.6281

0.5871

Table 7 Mammalian cytogenetic dosimetry lesions in rabbit lymphoblast induced by streptonigrin (NSC-45383), Exposure -90 μg|kg

Applications in thunderstorms

In thunderstorm activity, the occurrence of successive thunderstorm events (THE’s) is generally a dependent process meaning that the occurrence of a THE indicates that the atmosphere is unstable and the conditions are favorable for the formation for further thunderstorm activity. The negative binomial distribution (NBD) is a possible alternative to the Poisson distribution when successive events are possibly dependent Johnson et al.15 The theoretical and empirical justification for using the NBD to describe THE activity has been fully explained and discussed by Falls et al.24 Further, for fitting Poisson distribution to the count data equality of mean and variance must be satisfied. Similarly, for fitting NBD to the count data, mean must be less than the variance. In THE, these conditions are not fully satisfied. As a model to describe the frequencies of thunderstorms (TH’s), given an occurrence of THE, the PAD can be considered over PSD, PLD and PD because PAD, PSD and PLD are always  over-dispersed and PAD has advantage over PSD and PLD. The thunderstorms data have been considered in Tables 8-10.

No. of thunderstorms

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

187

155.6

185.3

184.8

185.4

1

77

117.0

83.5

83.6

82.7

2

40

43.9

35.9

36.3

36.3

3

17

11.02.10.30.1}

15.0
6.12.51.7}

15.2
6.12.41.6}

15.4
6.32.41.5}

4

6

5

2

6

1

Total

330

330.0

330.0

330.0

330.0

ML estimate

ˆθ=0.751515

ˆθ=1.804268

ˆθ=2.229891

ˆθ=2.625345

χ2

31.93

1.43

1.25

1.07

d.f.

2

3

3

3

p-value

0.0000

0.6985

0.741

0.7843

Table 8 Observed and expected number of days that experienced X thunderstorms events at Cape Kennedy, Florida for the 11-year period of record for the month of June, January 1957 to December 1967, Falls et al24

No. of thunderstorms

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

177

142.3

177.7

176.5

176.7

1

80

124.4

88.0

88.4

87.6

2

47

54.3

41.5

42.2

42.3

3

26

15.83.50.7}

18.9
8.46.5}

19.2
8.56.2}

19.5
8.66.3}

4

9

5

2

Total

341

341.0

341.0

341.0

341.0

ML estimate

ˆθ=0.809384

ˆθ=1.583536

ˆθ=1.995806

ˆθ=2.390474

χ2

39.74

5.15

4.67

4.35

d.f.

2

3

3

3

p-value

0.0000

0.1611

0.1976

0.2261

Table 9 Observed and expected number of days that experienced X thunderstorms events at Cape Kennedy, Florida for the 11-year period of record for the month of July, January 1957 to December 1967, Falls et al24

No. of thunderstorms

Observed frequency

Expected frequency

PD

PLD

PSD

PAD

0

185

151.8

184.8

184.1

184.7

1

89

122.9

87.2

87.5

86.6

2

30

49.7

39.3

39.8

39.8

3

24

13.42.70.5}

17.1
7.35.3}

17.3
7.35.0}

17.6
7.44.9}

4

10

5

3

Total

341

341.0

341.0

341.0

341.0

ML estimate

ˆθ=0.809384

ˆθ=0.809384

ˆθ=2.114545

ˆθ=2.511962

χ2

49.49

5.03

5.06

4.83

d.f.

2

3

3

3

p-value

0.0000

0.1696

0.1674

0.1847

Table 10 Observed and expected number of days that experienced X thunderstorms events at Cape Kennedy, Florida for the 11-year period of record for the month of August, January 1957 to December 1967, Falls et al.24

Again it is obvious from fitting of PAD to thunderstorms data that PAD gives better fit than PD, PLD, and PSD in all data .

Acknowledgments

None.

Conflicts of interest

None.

References

  1. Shanker R. The discrete Poisson‒Amarendra. International Journal of Statistical Distributions and Applications. 2016b.
  2. Sankaran M. The discrete Poisson‒Lindley distribution. Biometrics. 1970;26(1):145‒149.
  3. Shanker R. The discrete Poisson‒Sujatha distribution. International Journal of Probability and Statistics. 2016d;5(1):1‒9.
  4. Shanker R. Amarendra distribution and Its Applications. American Journal of Mathematics and Statistics. 2016a;6(1):44‒56.
  5. Shanker R. Shanker distribution and Its Applications. International Journal of Statistics and Applications. 2015a;5(6):338‒348.
  6. Shanker R. Akash distribution and Its Applications. International Journal of Probability and Statistics. 2015b;4(3):65‒75.
  7. Shanker R. Sujatha distribution and Its Applications. Statistics in Transition new Series. 2016c;17(3).
  8. Lindley DV. Fiducial distributions and Bayes theorem. Journal of Royal Statistical Society, Series. 1958;B20(1):102‒107.
  9. Shanker R, Hagos F. On Poisson‒Lindley distribution and Its applications to Biological Sciences. Biometrics & Biostatistics International Journal. 2015;2(4):1‒5.
  10. Shanker R, Hagos F, Sujatha S, et al. On zero‒truncation of Poisson and Poisson‒Lindley distributions and Their Applications. Biometrics & Biostatistics International Journal. 2015;2(6):1‒14.
  11. Shanker R, Hagos F. On Poisson‒Sujatha distribution and Its Applications to model count data from biological sciences. Biometrics & Biostatistics International Journal. 2016d;3(4):1‒7.
  12. Shanker R, Hagos F. Size‒biased Poisson‒Sujatha distribution with Applications. American Journal of Mathematics and Statistics. 2016a;6(4):145‒154.
  13. Shanker R, Hagos F. Zero‒truncated Poisson‒Sujatha distribution with Applications. Journal of Ethiopian Statistical Association. 2016b.
  14. Shanker R, Hagos F. On zero‒truncation of Poisson, Poisson‒Lindley, and Poisson‒Sujatha distribution and their Applications. Biometrics & Biostatistics International Journal. 2016c;3(5):1‒13.
  15. Johnson NL, Kotz S, Kemp AW. Univariate Discrete Distributions. 2nd edn, John Wiley & sons Inc, USA. 1992.
  16. Fisher RA, Corpet AS, Williams CB. The relation between the number of species and the number of individuals in a random sample of an animal population. Journal of Animal Ecology. 1943;12(1):42‒58.
  17. Kempton RA. A generalized form of Fisher’s logarithmic series. Biometrika. 1975;62(1):29‒38.
  18. Tripathi RC, Gupta RC. A generalization of the log‒series distribution. Comm in Stat (Theory and Methods). 1985;14(8):1779‒1799.
  19. Mishra A, Shanker R. Generalized logarithmic series distribution‒Its nature and applications. Proceedings of the Vth International Symposium on Optimization and Statistics. 2002;155‒168.
  20. Loeschke V, Kohler W. Deterministic and Stochastic models of the negative binomial distribution and the analysis of chromosomal aberrations in human leukocytes. Biometrical Journal. 1976;18(6):427‒451.
  21. Janardan KG, Schaeffer DJ. Models for the analysis of chromosomal aberrations in human leukocytes. Biometrical Journal. 1977;19 (8):599‒612.
  22. Catcheside DG, Lea DE, Thoday JM. Types of chromosome structural change induced by the irradiation on Tradescantia microspores. J Genet. 1946;47:113‒136.
  23. Mc Guire JU, Brindley TA, Bancroft TA. The distribution of European corn‒borer larvae pyrausta in field corn. Biometrics. 1957;13(1):65‒78.
  24. Catcheside DG, Lea DE, Thoday JM. The production of chromosome structural changes in Tradescantia microspores in relation to dosage, intensity and temperature. J Genet. 1946;47:137‒149.
Creative Commons Attribution License

©2016 Shanker, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.