Principios generales

En la propuesta se pretende tomar la estructura de la varianza con parámetros desconocidos; es decir, muestreo con reemplazo de media poblacional desconocida, varianza desconocida, tamaño de muestra fijo y población base gamma.

El resultado encontrado es que para un parámetro de forma \(\alpha\) mayor o igual que 2, la variable \(\chi^2_1=\frac{(n-1)*s^2}{\sigma^2}\) es tipo gamma con parámetro de forma mínimo \(\alpha'=\frac{n-1}{5}\) y \(\beta=5\). Aquí se define con \(b=\frac{1}{\beta}\).

mu=100
a=2
n=17
b=a/mu
sigma=sqrt(a/b^2)
library(MASS)
library(moments)
x=cbind(rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b))

chi1=(n-1)*apply(x,2,var)/sigma^2
chi2=(apply(x,2,mean)-mu)^2*n/sigma^2
chi3=chi1+chi2
cc=cbind(chi1,chi2,chi3)
px1=c(pgamma(n/2, (n-1)/5, 1/5), pgamma(n/2, (n-1)/2, 1/2), pgamma(n/2, (n-1)/(2+6/a), 1/(2+6/a)))
px2=c(pgamma(n+2, (n-1)/5, 1/5), pgamma(n+2, (n-1)/2, 1/2), pgamma(n+2, (n-1)/(2+6/a), 1/(2+6/a)))
px3=c(pgamma(2*n, (n-1)/5, 1/5), pgamma(n*2, (n-1)/2, 1/2), pgamma(n*2, (n-1)/(2+6/a), 1/(2+6/a)))
fa=ecdf(chi1)
fa1 <- fa(n/2)
fa2 <- fa(n+2)
fa3 <- fa(2*n)
rbind(px1,fa1,px2,fa2,px3,fa3)
##          [,1]      [,2]      [,3]
## px1 0.2038186 0.0674302 0.2038186
## fa1 0.1964286 0.1964286 0.1964286
## px2 0.6919364 0.7313368 0.6919364
## fa2 0.7544643 0.7544643 0.7544643
## px3 0.9569133 0.9945670 0.9569133
## fa3 0.9642857 0.9642857 0.9642857
pf=4/apply(x,2,skewness)^2
fax<-ecdf(pf)
fax1=fax(2)
fax1
## [1] 0.1875
mu=100
a=3
n=17
b=a/mu
sigma=sqrt(a/b^2)
library(MASS)
library(moments)
x=cbind(rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b))

chi1=(n-1)*apply(x,2,var)/sigma^2
chi2=(apply(x,2,mean)-mu)^2*n/sigma^2
chi3=chi1+chi2
cc=cbind(chi1,chi2,chi3)
px1=c(pgamma(n/2, (n-1)/5, 1/5), pgamma(n/2, (n-1)/2, 1/2), pgamma(n/2, (n-1)/(2+6/a), 1/(2+6/a)))
px2=c(pgamma(n+2, (n-1)/5, 1/5), pgamma(n+2, (n-1)/2, 1/2), pgamma(n+2, (n-1)/(2+6/a), 1/(2+6/a)))
px3=c(pgamma(2*n, (n-1)/5, 1/5), pgamma(n*2, (n-1)/2, 1/2), pgamma(n*2, (n-1)/(2+6/a), 1/(2+6/a)))
fa=ecdf(chi1)
fa1 <- fa(n/2)
fa2 <- fa(n+2)
fa3 <- fa(2*n)
rbind(px1,fa1,px2,fa2,px3,fa3)
##          [,1]      [,2]      [,3]
## px1 0.2038186 0.0674302 0.1661075
## fa1 0.1428571 0.1428571 0.1428571
## px2 0.6919364 0.7313368 0.6981144
## fa2 0.7053571 0.7053571 0.7053571
## px3 0.9569133 0.9945670 0.9698909
## fa3 0.9642857 0.9642857 0.9642857
pf=4/apply(x,2,skewness)^2
fax<-ecdf(pf)
fax1=fax(2)
fax1
## [1] 0.1428571
mu=100
a=4
n=17
b=a/mu
sigma=sqrt(a/b^2)
library(MASS)
library(moments)
x=cbind(rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b))

chi1=(n-1)*apply(x,2,var)/sigma^2
chi2=(apply(x,2,mean)-mu)^2*n/sigma^2
chi3=chi1+chi2
cc=cbind(chi1,chi2,chi3)
px1=c(pgamma(n/2, (n-1)/5, 1/5), pgamma(n/2, (n-1)/2, 1/2), pgamma(n/2, (n-1)/(2+6/a), 1/(2+6/a)))
px2=c(pgamma(n+2, (n-1)/5, 1/5), pgamma(n+2, (n-1)/2, 1/2), pgamma(n+2, (n-1)/(2+6/a), 1/(2+6/a)))
px3=c(pgamma(2*n, (n-1)/5, 1/5), pgamma(n*2, (n-1)/2, 1/2), pgamma(n*2, (n-1)/(2+6/a), 1/(2+6/a)))
fa=ecdf(chi1)
fa1 <- fa(n/2)
fa2 <- fa(n+2)
fa3 <- fa(2*n)
rbind(px1,fa1,px2,fa2,px3,fa3)
##          [,1]      [,2]      [,3]
## px1 0.2038186 0.0674302 0.1445620
## fa1 0.1160714 0.1160714 0.1160714
## px2 0.6919364 0.7313368 0.7028622
## fa2 0.6875000 0.6875000 0.6875000
## px3 0.9569133 0.9945670 0.9765936
## fa3 0.9687500 0.9687500 0.9687500
pf=4/apply(x,2,skewness)^2
fax<-ecdf(pf)
fax1=fax(2)
fax1
## [1] 0.09821429
mu=100
a=5
n=17
b=a/mu
sigma=sqrt(a/b^2)
library(MASS)
library(moments)
x=cbind(rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b))

chi1=(n-1)*apply(x,2,var)/sigma^2
chi2=(apply(x,2,mean)-mu)^2*n/sigma^2
chi3=chi1+chi2
cc=cbind(chi1,chi2,chi3)
px1=c(pgamma(n/2, (n-1)/5, 1/5), pgamma(n/2, (n-1)/2, 1/2), pgamma(n/2, (n-1)/(2+6/a), 1/(2+6/a)))
px2=c(pgamma(n+2, (n-1)/5, 1/5), pgamma(n+2, (n-1)/2, 1/2), pgamma(n+2, (n-1)/(2+6/a), 1/(2+6/a)))
px3=c(pgamma(2*n, (n-1)/5, 1/5), pgamma(n*2, (n-1)/2, 1/2), pgamma(n*2, (n-1)/(2+6/a), 1/(2+6/a)))
fa=ecdf(chi1)
fa1 <- fa(n/2)
fa2 <- fa(n+2)
fa3 <- fa(2*n)
rbind(px1,fa1,px2,fa2,px3,fa3)
##           [,1]       [,2]       [,3]
## px1 0.20381860 0.06743020 0.13065068
## fa1 0.05803571 0.05803571 0.05803571
## px2 0.69193640 0.73133682 0.70649058
## fa2 0.74107143 0.74107143 0.74107143
## px3 0.95691327 0.99456698 0.98058381
## fa3 0.96875000 0.96875000 0.96875000
pf=4/apply(x,2,skewness)^2
fax<-ecdf(pf)
fax1=fax(2)
fax1
## [1] 0.08928571
mu=100
a=10
n=17
b=a/mu
sigma=sqrt(a/b^2)
library(MASS)
library(moments)
x=cbind(rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b))

chi1=(n-1)*apply(x,2,var)/sigma^2
chi2=(apply(x,2,mean)-mu)^2*n/sigma^2
chi3=chi1+chi2
cc=cbind(chi1,chi2,chi3)
px1=c(pgamma(n/2, (n-1)/5, 1/5), pgamma(n/2, (n-1)/2, 1/2), pgamma(n/2, (n-1)/(2+6/a), 1/(2+6/a)))
px2=c(pgamma(n+2, (n-1)/5, 1/5), pgamma(n+2, (n-1)/2, 1/2), pgamma(n+2, (n-1)/(2+6/a), 1/(2+6/a)))
px3=c(pgamma(2*n, (n-1)/5, 1/5), pgamma(n*2, (n-1)/2, 1/2), pgamma(n*2, (n-1)/(2+6/a), 1/(2+6/a)))
fa=ecdf(chi1)
fa1 <- fa(n/2)
fa2 <- fa(n+2)
fa3 <- fa(2*n)
rbind(px1,fa1,px2,fa2,px3,fa3)
##          [,1]      [,2]      [,3]
## px1 0.2038186 0.0674302 0.1004820
## fa1 0.0937500 0.0937500 0.0937500
## px2 0.6919364 0.7313368 0.7162644
## fa2 0.7633929 0.7633929 0.7633929
## px3 0.9569133 0.9945670 0.9881785
## fa3 0.9866071 0.9866071 0.9866071
pf=4/apply(x,2,skewness)^2
fax<-ecdf(pf)
fax1=fax(2)
fax1
## [1] 0.04017857
mu=100
a=20
n=17
b=a/mu
sigma=sqrt(a/b^2)
library(MASS)
library(moments)
x=cbind(rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b),rgamma(n,a,b))

chi1=(n-1)*apply(x,2,var)/sigma^2
chi2=(apply(x,2,mean)-mu)^2*n/sigma^2
chi3=chi1+chi2
cc=cbind(chi1,chi2,chi3)
px1=c(pgamma(n/2, (n-1)/5, 1/5), pgamma(n/2, (n-1)/2, 1/2), pgamma(n/2, (n-1)/(2+6/a), 1/(2+6/a)))
px2=c(pgamma(n+2, (n-1)/5, 1/5), pgamma(n+2, (n-1)/2, 1/2), pgamma(n+2, (n-1)/(2+6/a), 1/(2+6/a)))
px3=c(pgamma(2*n, (n-1)/5, 1/5), pgamma(n*2, (n-1)/2, 1/2), pgamma(n*2, (n-1)/(2+6/a), 1/(2+6/a)))
fa=ecdf(chi1)
fa1 <- fa(n/2)
fa2 <- fa(n+2)
fa3 <- fa(2*n)
rbind(px1,fa1,px2,fa2,px3,fa3)
##           [,1]       [,2]       [,3]
## px1 0.20381860 0.06743020 0.08426000
## fa1 0.07142857 0.07142857 0.07142857
## px2 0.69193640 0.73133682 0.72293076
## fa2 0.72321429 0.72321429 0.72321429
## px3 0.95691327 0.99456698 0.99158719
## fa3 0.99553571 0.99553571 0.99553571
pf=4/apply(x,2,skewness)^2
fax<-ecdf(pf)
fax1=fax(2)
fax1
## [1] 0.03125

Cuando pf es mayor o igual que 2 se continúa con el proceso y se encuentra que la probabilidad \(px3\) cae entre \(px1\) y \(px2\) y se forma un intervalo de probabilidad.