[R-es] MUESTRA ESTRATIFICADA AFIJ OPTIMA

J. Miguel Marin jmmarin en est-econ.uc3m.es
Lun Feb 8 00:01:02 CET 2010


Hola,

Si sigues el capítulo de Scheaffer et al. (1990) sobre Técnicas de 
Muestreo (capítulo de M. Estratificado), un programa que hace afijación 
óptima podría ser el siguiente (para entenderlo bien deberías consultar 
el libro donde viene la teoría en la que se basa el programa):

strat.sample <- function(N.vec, sigma.vec, B, proportion=F) {
# Muestreo de forma proporcional o de forma uniforme por estrato
N.tot <- sum(N.vec)
if(proportion==F){
wt.vec <- rep(1/length(N.vec),length(N.vec))}
else {wt.vec <- N.vec/N.tot}
numer <- sum(((N.vec*sigma.vec)^2)/wt.vec)
den <- sum(N.vec*(sigma.vec^2))
D <- B^2/4
n <- ceiling(numer/(D*(N.tot^2) + den))
n.vec <- round(n*wt.vec,1)
cat(" Total n =",n,"\n", " n por estrato =",n.vec,"\n")
}

# Ejempo de Aplicacion
N.vec <- c(155,62,93)
sigma.vec <- c(5,15,10)
# O
strat.sample(N.vec,sigma.vec,2)
# O bien
strat.sample(N.vec,sigma.vec,2,T)


Si quieres aplicar los métodos clásicos de muesteo necesitas saber cual 
es el tamaño poblacional de cada estrato previamente. Si no, tendrías 
que aplicar postestratificación, lo cual es más complicado.

Saludos


> Buenos días con todos.
>
> Tengo una base de datos de reembolsos (62.000 casos aprox) que hizo 
> una compañía de seguros a sus asegurados.
> Mi objetivo es obtener una muestra estratificada con afijacion óptima 
> en base a esta variable(V.PAGADO) de tal manera que el tamaño de 
> muestra n sea mínimo.
>
> Este es un mini ejemplo de los datos que tengo:
> V.PAGADO <- c(247.38, 55.59, 124.05, 84.02, 44.51, 9.34, 36.99, 
> 133.12, 17.00, 72.00, 1005.05, 58.55, 1031.25, 88.79, 58.14, 47.61, 
> 201.41, 1703.27, 69.51, 37.02, 27.35, 36.14, 70.00, 110.29, 713.30, 
> 63.18, 51.77, 44.24, 1250.00, 31.50, 1180.50, 31.17, 148.28, 22.50, 
> 21.60, 19.00, 12.00, 9.00, 9.00, 9.00, 9.00, 7.00, 6.24, 17.71, 
> 110.39, 183.71, 80.00, 31.70, 57.88, 40.00)
>
>
> Hay algún paquete que me facilite este proceso y si lo hay como lo 
> aplicaría a estos datos???
>
> Gracias de antemano.
>
> _________________________________________________________________
>
>
> 	[[alternative HTML version deleted]]
>
>




jm~

_______________________________

        J. Miguel Marin

http://www.est.uc3m.es/jmmarin

    Dep. of Statistics
University Carlos III of Madrid
        Spain (E.U.)



Más información sobre la lista de distribución R-help-es