beta2bs {bs}R Documentation

The coefficient of kurtosis of the BSD

Description

Estimation of the coefficient of kurtosis (CK) for the Birnbaum-Saunders distribution.

Usage

beta2bs(x)

Arguments

x Vector of observations.

Details

Methods of goodness-of-fit based on moments (β_1-β_2 chart).

Value

The function beta2bs() gives the CK value.

Author(s)

Víctor Leiva <victor.leiva@uv.cl>, Hugo Hernández <hugo.hernande@msn.com>, and Marco Riquelme <mriquelm@ucm.cl>.

References

Leiva, V., Hernández, H., and Riquelme, M. (2006). A New Package for the Birnbaum-Saunders Distribution. Rnews, 6/4, 35-40. (http://www.r-project.org)

Examples

## beta1-beta2 chart
## Example: 12 BSD samples have been simulated by using 'rbs()' for the values
## alfa = 0.2,0.4,0.5,0.6,0.8,1.0,1.2,1.4,1.8,2.0,2.2,2.4
## beta = 1.0
sample1<-rbs(30,0.2,1.0)
sample2<-rbs(30,0.4,1.0)
sample3<-rbs(30,0.6,1.0)
sample4<-rbs(30,0.8,1.0)
sample5<-rbs(30,1.0,1.0)
sample6<-rbs(30,1.2,1.0)
sample7<-rbs(30,1.4,1.0)
sample8<-rbs(30,1.6,1.0)
sample9<-rbs(30,1.8,1.0)
sample10<-rbs(30,2.0,1.0)
sample11<-rbs(30,2.2,1.0)
sample12<-rbs(30,2.4,1.0)

x1<-c(beta1bs(sample1),beta2bs(sample1))
x2<-c(beta1bs(sample2),beta2bs(sample2))
x3<-c(beta1bs(sample3),beta2bs(sample3))
x4<-c(beta1bs(sample4),beta2bs(sample4))
x5<-c(beta1bs(sample5),beta2bs(sample5))
x6<-c(beta1bs(sample6),beta2bs(sample6))
x7<-c(beta1bs(sample7),beta2bs(sample7))
x8<-c(beta1bs(sample8),beta2bs(sample8))
x9<-c(beta1bs(sample9),beta2bs(sample9))
x10<-c(beta1bs(sample10),beta2bs(sample10))
x11<-c(beta1bs(sample11),beta2bs(sample11))
x12<-c(beta1bs(sample12),beta2bs(sample12))
x<-c(x1[1],x2[1],x3[1],x4[1],x5[1],x6[1],x7[1],x8[1],x9[1],x10[1],x11[1],x12[1])
y<-c(x1[2],x2[2],x3[2],x4[2],x5[2],x6[2],x7[2],x8[2],x9[2],x10[2],x11[2],x12[2])

## Generate the graph
plot(x,y,xlab=expression(beta[1]),ylab=expression(beta[2]),col=4,lwd=2.0)
fitbetabs(0.2,2)
fitbetabs(2.4,3)

[Package bs version 1.0 Index]