Visualizando densidades de población con R

Cuando se dispone de la densidad de población de una región en formato matricial, además de la visualización obvia de la misma pueden hacerse muchos juegos numéricos para representarla de mil formas creativas.

Bouncing ball con el becario motivado: ChatGPT

Apareció hace unas semanas un test al que se sometía a varios LLM's, consistente en pedir que creasen una rutina en Python que modelase una pelotilla rebotando por dentro de las paredes de un cuadrado que rota. Por lo visto solo Grok 3 dio una solución correcta, los demás modelos no detectaban bien las colisiones y la pelota se salía del cuadrado o hacía cosas raras.

La transformada de Box-Muller con R

La transformada de Box-Muller es un método de generación de pares de números aleatorios independientes con distribución normal de media nula y varianza 1, a partir de pares de números aleatorios con distribución uniforme entre 0 y 1. La transformación consiste en mapear pares de coordenadas (U1,U2) pertenecientes al cuadrado unidad definido por (0,0) y (1,1), sobre puntos (Z0,Z1) cuyo dominio se extiende a cualquier punto del plano.

Random walk browniano con R

Un random walk es un proceso estocástico que describe un recorrido formado por una sucesión de pasos aleatorios en algún espacio matemático. Cuando este movimiento es continuo, es decir cuando todas las ubicaciones del espacio son factibles, hablamos de movimiento browniano, estudiado en Física para modelar el comportamiento de partículas suspendidas en un medio líquido o gaseoso. El patrón de movimiento consiste en fluctuaciones aleatorias en la posición de una partícula dentro de un subdominio del fluido, seguido de una reubicación a otro subdominio.

Radiografía de Tenerife con R

En este ejercicio además de crear una animación que combine varios de los mapas que hemos venido creando (hillshade, curvas de nivel, proyección de sombras,...), vamos a implementarlo de forma vectorizada, para poder superponer con una línea de código sobre un mapa base información sumarizada más concreta, utilizando una flexible máscara con transparencia al más puro estilo de las máscaras de capa de Photoshop.

Cuantificando el efecto del filtro de luz azul de unas gafas graduadas

Algunas gafas graduadas (entre ellas las mías) incluyen un filtro azul que atenúa las frecuencias de luz más altas, supuestamente reduciendo la fatiga ocular y protegiendo el ojo a largo plazo frente al uso continuado de pantallas.

Teoría de colas y fiestas de pueblo

La Teoría de colas es la disciplina que estudia los procesos de demanda de Servicio por parte de una serie de Clientes sobre uno o más Recursos o Servidores. Ejemplos de colas son las que se forman en las cajas del supermercado, en la parada del bus, y también en cosas más incorpóreas como las aplicaciones que se ejecutan en un ordenador solicitando tiempo de CPU.

Calculando el promedio ponderado de una serie histórica con un único campo

En ocasiones se precisa sumarizar para cada fila de una tabla el valor que ha tenido una variable a lo largo del tiempo. La forma obvia aunque costosa de lograrlo es almacenando individualmente en columnas los últimos n valores del histórico y operar sobre ellos para obtener algún tipo de promediado (media, media ponderada, mediana,...). Vamos a ver una forma mucho más económica de calcular un promedio ponderado por antigüedad de los valores de la serie, solo necesitando el propio campo donde se almacena el dato.

Modelo geométrico de piedra de playa con R

Cuando voy a una playa de piedras me gusta fijarme en la forma ovalada que adquieren las piedras, labradas por las olas a lo largo de miles o millones de años. Vamos a programar una rutina que, a partir del contorno de una piedra de este tipo, la modele en base a 4 elipses. De cada una de ellas se utiliza 1/4 de la forma garantizando la continuidad de las derivadas, y por tanto la suavidad del contorno modelado, en los cuatro puntos de unión.

Anaglifos de un mapa de elevación con R

En 'Proyección de sombras sobre un DEM con R' partiendo de un mapa de elevaciones obteníamos una visión pseudo 3D coloreada de su hillshade, añadiéndole además una proyección de sombras. Le vamos a dar una nueva vuelta de tuerca para construir, con el mismo método de slicing que usamos en aquella ocasión, un anaglifo que permita percibir de forma real el relieve 3D del mapa.