detecta_cluster: Detección de agrupaciones de mortalidad a revisar manualmente

In fisabio/medear: Poblaciones y Cartografia por Seccion Censal INE

Description

Esta función es útil a la hora de comprobar la geocodificación de la mortalidad, pues devuelve aquellos puntos con un exceso de fallecimientos en relación a la media de sus vecinos más próximos. Con esta dinámica es fácil identificar los centros residenciales o errores sistemáticos en la geocodificación, donde es de esperar una mayor aglomeración de defunciones.

Usage

detecta_cluster(datos, epsg = 4326, vecinos = 10, cartografia = NULL,
  secciones = NULL, limite = c(1e-10, 1e-15, 1e-20), devuelve_datos = TRUE)

Arguments

datos	Base de datos con las coordenadas que ubican cada uno de los fallecimientos. Debe contener, al menos, 3 columnas: BOD.direccion, lat y lng, las cuales deben tener exactamente esos nombres. De forma paralela la función tratará de buscar otras seis columnas fruto del protocolo de geocodificación (province, muni, tip_via, address, portalNumber y postalCode) (nuevamente deben tener exactamente esos nombres, que son los que resultan del protocolo de geocodificación), con el propósito de mostrar, junto a las direcciones del BOD, las direcciones devueltas por los servicios de geocodificado. Si la base de datos tuviera otros nombres, seria trabajo del usuario cambiárselos como paso previo al uso de la función.
epsg	Numérico con longitud igual a 1: código EPSG con la proyección de las coordenadas de la base de datos. Por defecto se usa el EPSG 4326 (longlat WGS 84).
vecinos	Numérico con longitud igual a 1: número de vecinos más próximos con los que comparar la mortalidad de cada punto. Por defecto 10.
cartografia	Objeto de clase SpatialPolygons. Con ella se representa el seccionado oficial de 2011 a modo de referencia con la que juzgar la geocodificación de un conjunto de defunciones anómalas. En principio tiene valor nulo, lo que implica usar la cartografía asociada al paquete para las ciudades MEDEA3. Si el usuario desea consultar otros municipios puede hacer uso de este argumento para tener el seccionado de fondo.
secciones	Cadena de caracteres: secciones censales anómalas (información facilitada por el nodo coordinador), para resaltarlas en el mapa. Por defecto ninguna.
limite	Numérico con longitud 1 <= x <= 9: límites de probabilidad con los que se identifican las agrupaciones sospechosas. Por defecto viene fijado a 1e-10, 1e-15 y 1e-20.
devuelve_datos	Valor lógico, por defecto TRUE. ¿La función debe devolver un data.frame con los datos identificativos de los clústeres?

Details

La función comienza calculando el número de fallecimientos en cada par único de coordenadas, identifica los n vecinos más próximos a cada punto y calcula la media de fallecimientos en los mismos. Considerando que los fallecimientos en cada par de coordenadas siguen una distribución de Poisson cuya media es la media de los fallecimientos en los puntos más cercanos, se calcula la probabilidad de superar esa media.

Value

Representa un objeto de clase leaflet en el que se marcan los puntos a revisar. Por defecto los puntos se dividen en tres colores en función de la mortalidad acontecida en las coordenadas vecinas: verde si la probabilidad de obtener ese resultado es inferior a 1e-10, azul si es inferior a 1e-15 y rojo si es inferior a 1e-20 (el caso más evidente). Haciendo clic en cada uno de los puntos se puede consultar las direcciones, tanto del BOD como las obtenidas en la geocodificación (se puede cambiar entre una y otra en el menú del extremo superior derecho), asociadas a cada punto. Al comienzo de cada dirección se indica el número de fila al que hace referencia cada dirección en los datos que se facilitaron, de forma que pueda recuperarse fácilmente dicha información para explorar los datos en profundidad.

Respecto al menú anteriormente mencionado, también permite cambiar la capa de visualización de fondo, utilizando OpenStreetMap, Google Maps o Google Satellite (de forma que pueda disponerse de información extra sin salir de la aplicación).

Si se desea efectuar el cálculo de distancias o áreas, en la esquina inferior izquierda se dispone de un menú para realizarlas (en metros y en metros cuadrados).

La función también devuelve los datos con las cooredenadas (campos lng y lat), número de defunciones (campo N), límite dentro del cual recae (campo limite, útil para filtrar resultados) y direcciones asociadas a cada punto (campo detalle, de tipo lista donde cada elemento contiene un data.frame con la información detallada).

Examples

## Not run: 
  library(medear)
  revisar <- detecta_cluster(datosmort)
  revisar

## End(Not run)

fisabio/medear documentation built on Aug. 2, 2021, 2:15 p.m.