beta1bs {bs} | R Documentation |
Estimation of the square of skewness coefficient (SSC) for the Birnbaum-Saunders distribution.
beta1bs(x)
x |
Vector of observations. |
Methods of goodness-of-fit based on moments (β_1-β_2 chart).
The function beta1bs()
gives the SSC value.
Víctor Leiva <victor.leiva@uv.cl>, Hugo Hernández <hugo.hernande@msn.com>, and Marco Riquelme <mriquelm@ucm.cl>.
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)
## 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)