S-Expo CDF
pSe = function(x,p)
{
p = 0.10780
cos((pi/2)*exp(-p*x))
}
ks1=ad1=cvm1=NULL
ks1=ks.test(y,pSe)
ad1=ad.test(y,pSe)
cvm1=cvm.test(y,pSe)
result1=c(ks1$statistic,ks1$p.value,ad1$statistic,ad1$p.value,
cvm1$statistic,cvm1$p.value)
#IL distribution
pIL = function(x,p)
{
p = 4.9096
(1 + (p/((1+p)*x)))*exp(-p/x)
}
ks1=ad1=cvm1=NULL
ks1=ks.test(y,pIL)
ad1=ad.test(y,pIL)
cvm1=cvm.test(y,pIL)
result1=c(ks1$statistic,ks1$p.value,ad1$statistic,ad1$p.value,
cvm1$statistic,cvm1$p.value)