Casos de uso de MicroDatosEs
MicroDatosEs ha sufrido algunas modificaciones recientes. En particular, un nuevo colaborador, Jorge López Pérez ha reescrito la documentación usando roxygen2.… más
r
gam/bam admiten efectos aleatorios
gam/bam admiten efectos aleatorios gam/bam admiten efectos aleatorios gam/bam admiten efectos aleatorios gam/bam admiten efectos aleatorios gam/bam admiten efectos aleatorios… más
Leaflet con capas WMS de Coreos, Catastro, etc.
Esta entrada es un subproducto del trabajo que pocería que he realizado estos días en caRtociudad. caRtociudad permite generar mapas… más
t y as.raster no conmutan; ¿por qué no conmutarán?
Creo una minimatriz, la convierto en un raster y la represento:
|
1 2 3 4 5 6 7 8 9 10 11 12 13 |
m <- matrix(c(0, 0, 0.33, 0.66, .9, .9), 2, 3) m # [,1] [,2] [,3] # [1,] 0 0.33 .9 # [2,] 0 0.66 .9 r <- as.raster(m) r # [,1] [,2] [,3] # [1,] "#000000" "#545454" "#FFFFFF" # [2,] "#000000" "#A8A8A8" "#FFFFFF" plot(r, interpolate = FALSE) |
Ahora, con la matriz traspuesta,
|
1 2 3 4 5 6 |
r_t_1 <- as.raster(t(m)) r_t_1 # [,1] [,2] # [1,] "#000000" "#000000" # [2,] "#545454" "#A8A8A8" # [3,] "#E6E6E6" "#E6E6E6" |
obtengo que… más
El malabarista de Amiga, con R (y rgl)
Un exalumno mío de unos cursos de R, Guillermo Luijk, ha creado usando R (y el paquete rgl). Toda la… más
Modelos con inflación de ceros y separación perfecta
Al estudiar problemas de conteos, la llamada inflación de ceros ocurre frecuentemente: los datos contienen más ceros de los que… más
La intrahistoria de mi libro de R
Una de las preguntas más fértiles que pueden formularse frente a algo es la del motivo de su existencia: ¿por… más
Un libro de R: mi libro de R
No quería hacerlo público aún pero alguien se ve que lo estaba leyendo por ahí. No sé si Google habrá… más
¿Un voluntario para aggiornar MicroDatosEs?
Mi paquete MicroDatosEs ya forma parte de rOpenSpain. Sin embargo, está falto de ciertas mejoras a las que aspiran los… más
¿Podría ser la solución que almas caritativas creasen viñetas espontáneamente?
Uno de los modelos más útiles potencialmente y que menos atención recibe es el de los modelos de conteos autoexcitados.… más
Documentar como el culo, no pensar en el usuario final, ser incapaz de ponerte en su situación, etc.
De vez en cuando pruebo paquetes promisorios. No es infrecuente que cosas que he intentado hace años, algún ejemplo más… más
Estructura poblacional de España: 2010-2050
Si se puede hacer para Japón, también se puede hacer para España: El código,
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 |
library(idbr) library(ggplot2) library(animation) library(ggthemes) idb_api_key("pídela en https://www.census.gov/data/developers/data-sets/international-database.html") male <- idb1('SP', 2010:2050, sex = 'male') male$SEX <- "hombres" male$POP <- -male$POP female <- idb1('SP', 2010:2050, sex = 'female') female$SEX <- "mujeres" spain <- rbind(male, female) saveGIF({ for (i in 2010:2050) { title <- as.character(i) year_data <- spain[spain$time == i, ] g1 <- ggplot(year_data, aes(x = AGE, y = POP, fill = SEX, width = 1)) + coord_fixed() + coord_flip() + annotate('text', x = 98, y = -300000, label = 'Datos: US Census Bureau IDB; idbr R package', size = 3) + geom_bar(data = subset(year_data, SEX == "mujeres"), stat = "identity") + geom_bar(data = subset(year_data, SEX == "hombres"), stat = "identity") + scale_y_continuous(breaks = seq(-300000, 300000, 150000), labels = paste0(as.character(c(seq(300, 0, -150), c(150, 300))), "k"), limits = c(min(spain$POP), max(spain$POP))) + theme_economist(base_size = 14) + scale_fill_manual(values = c('#ff9896', '#d62728')) + ggtitle(paste0('Estructura poblacional de España en, ', title)) + ylab('Población') + xlab('Edad') + theme(legend.position = "bottom", legend.title = element_blank()) + guides(fill = guide_legend(reverse = TRUE)) print(g1) } }, movie.name = 'spain_pyramid.gif', interval = 0.1, ani.width = 700, ani.height = 600) |