Estoy intentando sacar de forma aleatoria una muestra y he establecido un bucle para ello que me seleccione de cada grupo indicado x individuos por provincia, sexo y edad. Para ello, primero calculo la cantidad representativa de cada grupo y luego, con el bucle, selecciono la muestra de un panel. El problema es que hay ciertos grupos que no tienen, del panel de donde se saca, ningún individuo representativo y cuando llega a este me da error y se para el bucle con el siguiente error: Error in slice_sample(): ! n must be positive
. Mi duda es si es posible indicar que cuando se encuentre un grupo que no tiene representación, pase al siguiente, sin que se pare.
El código que he utilizado hasta ahora es el siguiente:
Calculo lo que corresponde a cada grupo por provincia, edad y sexo
muestra <- 4000
cuotas <- datos.ine %>% mutate(total=porc*n_muestra/100)
cuotas <- cuotas %>% mutate(e18_28=porc*total[[51]]/100,
e25_34=porc*total[[52]]/100)
cuotas <- cuotas %>% mutate(h18_28=e18_28*porc[[53]]/100,
h25_34=e25_34*porc[[53]]/100,
m18_28=e18_28*porc[[54]]/100,
h25_34=e25_34*porc[[54]]/100,
Atribuyo a targets lo que he calculado en cuotas y lo coloco en cuatro filas por provincia, sexo, edad y el total que necesito de cada grupo y lo ordeno todo por las tres primeras variables
targets <- cuotas %>% pivot_longer(c("h18_28", "h25_34
"m18_28" "m25_34"),
names_to="grupo", values_to="total")
targets <- targets %>% separate(grupo, c("Sexo", "Edad"), sep=1)
targets <- targets %>% mutate(Sexo=case_when(Sexo=="h" ~ "Hombre",
Sexo=="m" ~ "Mujer"))
targets<- targets %>% arrange(provincia, Sexo, Edad)
Establezco:
1. Primero, un valor llamado n, que refleja el número de filas de nuestras sentencias.
2. Segundo, una lista vacía llamada muestra_list. En esta lista se irán guardando las bases de datos de cada grupo por el que segmentamos.
3. El bucle para que me coja la sentencias establecidas en targets del panel. Y que aquí es donde me da el error al no encontrar en el panel a mínimo un individuo de alguno de los grupos
n = nrow(targets)
muestra_list <- list()
for(i in 1:n){
muestra<- panel %>% filter(Provincia==targets$provincia[[i]],
Sexo==targets$Sexo[[i]],
Grupoedad==targets$Edad[[i]])
cantidad <- ifelse(targets$total[[i]] > nrow(muestra), nrow(muestra), targets$total[[i]])
muestra_list[[i]] <- slice_sample(muestra, n = cantidad)
}
Muchas gracias¡!
if (cantidad > 0) muestra_list[[i]] <- slice_sample(muestra, n = cantidad)
no es lo que buscas?