numeroteca · June 19, 2018 23:21
diff --git a/analisis.R b/analisis.R
 # load library
 library(tidyverse)

 data <- read.csv("centros-escolares-euskadi-valores-aleatorios.csv", sep = "," )

 # -----------Calcula índice de Gorard por municipio -----------

 # Unmelt data to make it easier for segregation index calculation. Based sonÑ
 # https://ojs.uv.es/index.php/RASE/article/view/10129/10853
 # http://journals.sagepub.com/doi/abs/10.1177/003803850203600405
 # Indicador de Gorard
 # 1/2 * sumatorio de
 # | n extranjeros en escuela i / total extranjeros en región estudiada - total alumnos en escuela i / total alumnos en región estudiada |
 require(reshape2)
 centros_agrupados <-dcast(data,Centro + Red + Territorio + Municipio + Código ~ NIV_NAC)
 centros_agrupados[is.na(centros_agrupados)] <- 0

 # Para primaria, código 120
 # contar número de alumnos por municipio
 muncipios_alumnado <- group_by(data[data$NIV_NAC=="120EXT"|data$NIV_NAC=="120ESP",], Municipio) %>% summarise(tot_pri = sum(Alumnado))
 # contar el número de extranjeros por municipio
 muncipios_alumnado_ext <- group_by(data[data$NIV_NAC=="120EXT",], Municipio) %>% summarise(ext_pri = sum(Alumnado))
 # crea lista de municipios que tienen extranjeros en primaria
 municipios_ext_list <- levels(factor(muncipios_alumnado_ext$Municipio))

 # merge data from municipios in centros agrupados
 centros_agrupados <- merge(centros_agrupados,muncipios_alumnado,by="Municipio", all = TRUE)
 centros_agrupados <- merge(centros_agrupados,muncipios_alumnado_ext,by="Municipio", all = TRUE)
 centros_agrupados[is.na(centros_agrupados)] <- 0

 # Para secundaria, código 260
 # contar número de alumnos por municipio
 muncipios_alumnado_sec <- group_by(data[data$NIV_NAC=="260EXT"|data$NIV_NAC=="260ESP",], Municipio) %>% summarise(tot_sec = sum(Alumnado))
 # contar el número de extranjeros por municipio
 muncipios_alumnado_sec_ext <- group_by(data[data$NIV_NAC=="260EXT",], Municipio) %>% summarise(ext_sec = sum(Alumnado))
 # crea lista de municipios que tienen extranjeros en secundaria
 municipios_ext_list_sec <- levels(factor(muncipios_alumnado_sec_ext$Municipio))

 # merge data from municipios in centros agrupados
 centros_agrupados <- merge(centros_agrupados,muncipios_alumnado_sec,by="Municipio", all = TRUE)
 centros_agrupados <- merge(centros_agrupados,muncipios_alumnado_sec_ext,by="Municipio", all = TRUE)
 centros_agrupados[is.na(centros_agrupados)] <- 0

 # calcula número de extranjeros por centro
 centros_agrupados$sum_ext_centro <- centros_agrupados$`120EXT` + centros_agrupados$`260EXT`
 # calcula número de alumnos por centro
 centros_agrupados$sum_alum_centro <- centros_agrupados$`120EXT` + centros_agrupados$`120ESP` + centros_agrupados$`260EXT` + centros_agrupados$`260ESP`

 # calcula número de extranjeros para cada municipio
 centros_agrupados$ext_total_municipio <- centros_agrupados$ext_pri + centros_agrupados$ext_sec
 # calcula número de alumnos para cada municipio
 centros_agrupados$alum_total_municipio <- centros_agrupados$ext_pri + centros_agrupados$tot_pri + centros_agrupados$ext_sec + centros_agrupados$tot_sec

 # calcula Gorard para cada centro 
 centros_agrupados$gorard <- abs( 
                                (centros_agrupados$sum_ext_centro/centros_agrupados$ext_total_municipio ) - 
                                 ( centros_agrupados$sum_alum_centro /centros_agrupados$alum_total_municipio)
                                )


 # Combine the two list of municipios that have primaria and secundaria
 unique(c(municipios_ext_list, municipios_ext_list_sec))
 municipios <- unique(c(municipios_ext_list, municipios_ext_list_sec))

 # calcular el índice de gorard para cada municipio
 gorard_muni <- data.frame(matrix(ncol = 4,nrow = length(municipios)))
 names(gorard_muni) <- c("municipio","indice_gorard","ig_pub","ig_pri")

 # Crear loop para iterar portodos los municipios y sumar los valores
 i <- 0
 for (municipio in municipios_ext_list) {
  i<- i + 1
  print(municipio)
  indice_gorard <- sum(centros_agrupados[centros_agrupados$Municipio==municipio,"gorard"])
  ig_pub <- sum(centros_agrupados[centros_agrupados$Municipio==municipio & centros_agrupados$Red == "Pública","gorard"]) 
  ig_pri <- sum(centros_agrupados[centros_agrupados$Municipio==municipio & centros_agrupados$Red == "Privada","gorard"]) 
  gorard_muni[i,1] <- municipio
  gorard_muni[i,2] <- 0.5*indice_gorard
  gorard_muni[i,3] <- 0.5*ig_pub
  gorard_muni[i,4] <- 0.5*ig_pri
 }

 write.csv(gorard_muni, file = "municipios_primaria-secundaria_indice-gorard.csv")
diff --git a/centros-escolares-euskadi-valores-aleatorios.csv b/centros-escolares-euskadi-valores-aleatorios.csv
	# load library
	library(tidyverse)

	data <- read.csv("centros-escolares-euskadi-valores-aleatorios.csv", sep = "," )

	# -----------Calcula índice de Gorard por municipio -----------

	# Unmelt data to make it easier for segregation index calculation. Based sonÑ
	# https://ojs.uv.es/index.php/RASE/article/view/10129/10853
	# http://journals.sagepub.com/doi/abs/10.1177/003803850203600405
	# Indicador de Gorard
	# 1/2 * sumatorio de
	# \| n extranjeros en escuela i / total extranjeros en región estudiada - total alumnos en escuela i / total alumnos en región estudiada \|
	require(reshape2)
	centros_agrupados <-dcast(data,Centro + Red + Territorio + Municipio + Código ~ NIV_NAC)
	centros_agrupados[is.na(centros_agrupados)] <- 0

	# Para primaria, código 120
	# contar número de alumnos por municipio
	muncipios_alumnado <- group_by(data[data$NIV_NAC=="120EXT"\|data$NIV_NAC=="120ESP",], Municipio) %>% summarise(tot_pri = sum(Alumnado))
	# contar el número de extranjeros por municipio
	muncipios_alumnado_ext <- group_by(data[data$NIV_NAC=="120EXT",], Municipio) %>% summarise(ext_pri = sum(Alumnado))
	# crea lista de municipios que tienen extranjeros en primaria
	municipios_ext_list <- levels(factor(muncipios_alumnado_ext$Municipio))

	# merge data from municipios in centros agrupados
	centros_agrupados <- merge(centros_agrupados,muncipios_alumnado,by="Municipio", all = TRUE)
	centros_agrupados <- merge(centros_agrupados,muncipios_alumnado_ext,by="Municipio", all = TRUE)
	centros_agrupados[is.na(centros_agrupados)] <- 0

	# Para secundaria, código 260
	# contar número de alumnos por municipio
	muncipios_alumnado_sec <- group_by(data[data$NIV_NAC=="260EXT"\|data$NIV_NAC=="260ESP",], Municipio) %>% summarise(tot_sec = sum(Alumnado))
	# contar el número de extranjeros por municipio
	muncipios_alumnado_sec_ext <- group_by(data[data$NIV_NAC=="260EXT",], Municipio) %>% summarise(ext_sec = sum(Alumnado))
	# crea lista de municipios que tienen extranjeros en secundaria
	municipios_ext_list_sec <- levels(factor(muncipios_alumnado_sec_ext$Municipio))

	# merge data from municipios in centros agrupados
	centros_agrupados <- merge(centros_agrupados,muncipios_alumnado_sec,by="Municipio", all = TRUE)
	centros_agrupados <- merge(centros_agrupados,muncipios_alumnado_sec_ext,by="Municipio", all = TRUE)
	centros_agrupados[is.na(centros_agrupados)] <- 0

	# calcula número de extranjeros por centro
	centros_agrupados$sum_ext_centro <- centros_agrupados$`120EXT` + centros_agrupados$`260EXT`
	# calcula número de alumnos por centro
	centros_agrupados$sum_alum_centro <- centros_agrupados$`120EXT` + centros_agrupados$`120ESP` + centros_agrupados$`260EXT` + centros_agrupados$`260ESP`

	# calcula número de extranjeros para cada municipio
	centros_agrupados$ext_total_municipio <- centros_agrupados$ext_pri + centros_agrupados$ext_sec
	# calcula número de alumnos para cada municipio
	centros_agrupados$alum_total_municipio <- centros_agrupados$ext_pri + centros_agrupados$tot_pri + centros_agrupados$ext_sec + centros_agrupados$tot_sec

	# calcula Gorard para cada centro
	centros_agrupados$gorard <- abs(
	(centros_agrupados$sum_ext_centro/centros_agrupados$ext_total_municipio ) -
	( centros_agrupados$sum_alum_centro /centros_agrupados$alum_total_municipio)
	)


	# Combine the two list of municipios that have primaria and secundaria
	unique(c(municipios_ext_list, municipios_ext_list_sec))
	municipios <- unique(c(municipios_ext_list, municipios_ext_list_sec))

	# calcular el índice de gorard para cada municipio
	gorard_muni <- data.frame(matrix(ncol = 4,nrow = length(municipios)))
	names(gorard_muni) <- c("municipio","indice_gorard","ig_pub","ig_pri")

	# Crear loop para iterar portodos los municipios y sumar los valores
	i <- 0
	for (municipio in municipios_ext_list) {
	i<- i + 1
	print(municipio)
	indice_gorard <- sum(centros_agrupados[centros_agrupados$Municipio==municipio,"gorard"])
	ig_pub <- sum(centros_agrupados[centros_agrupados$Municipio==municipio & centros_agrupados$Red == "Pública","gorard"])
	ig_pri <- sum(centros_agrupados[centros_agrupados$Municipio==municipio & centros_agrupados$Red == "Privada","gorard"])
	gorard_muni[i,1] <- municipio
	gorard_muni[i,2] <- 0.5*indice_gorard
	gorard_muni[i,3] <- 0.5*ig_pub
	gorard_muni[i,4] <- 0.5*ig_pri
	}

	write.csv(gorard_muni, file = "municipios_primaria-secundaria_indice-gorard.csv")