La serie en cuestión es esta (abridla en otra pestaña para verla en la plenitud de su definición): Con tiene todas (explicar por qué no son todas sería complicado, pero se puede dar el cuantificador casi por bueno) las defunciones (diarias) en España desde la fecha indicada. Los datos brutos están en la figura superior. Las tres siguientes tienen la descomposición estacional, la tendencia y los residuos tal como los estima stl. ...

Hoy he visto aquí y he escrito m <- matrix(c(74, 15, 10, 1, 11, 50, 38, 1, 5, 4, 90, 1, 17, 4, 19, 60), 4, 4, byrow = TRUE) m <- m / 100 luego m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m%*% m%*% m%*% m%*% m%*% m%*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m %*% m%*% m%*% m%*% m%*% m%*% m%*% m # [,1] [,2] [,3] [,4] #[1,] 0.1926676 0.1133218 0.6696203 0.02439024 #[2,] 0.1926647 0.1133206 0.6696245 0.02439024 #[3,] 0.1926638 0.1133202 0.6696258 0.02439024 #[4,] 0.1926675 0.1133218 0.6696205 0.02439025 y finalmente res <- eigen(t(m)) res$vectors[,1] / sum(res$vectors[,1]) #[1] 0.19266473 0.11332059 0.66962444 0.02439024 Aquí dice por qué.

Sirve esta entrada para hacer saber lo fundamental del trabajo de fin de master (TFM en lo que sigue) de Susana Huedo (que busca trabajo y es una chica muy sabida, aplicada y espabilada) en el CIFF. Los TFM que propongo y acabo supervisando jamás tienen vocación de criogénesis anaquélica. A Susana le sugerí un tema muy punk y con recorrido: [tratar de] predecir a los predictores. Fundamentalmente, para joder. Los chefs de encuestas electorales tienen dos discursos —uno previo y otro posterior a la publicación de los resultados—, una serie de recetas y datos que solo excepcionalmente publican. Dirán que se ciñen a una metodología científica, etc. Literatura. ...

Primera entrada del nuevo curso. Que lo es para recordaros que las XI Jornadas de Usuarios de R están en marcha. Y que buscan asistentes y ponentes. Así que, ¿nos veremos por allá?

Leo (¡en prensa!): Lo ha dicho el Instituto Nacional de Estadística: desde el inicio de la presunta crisis económica —por otros llamada fraude— el porcentaje de suicidios en España ha crecido un 20%. Nosotros somos cinco hermanos: me toca uno. Algo estamos haciendo fatal.

Esto es una página del DRAE en papel de toda la vida: Esta es una de esas definiciones en mi pantalla. No obstante, a pesar de todo el espacio blanco sobrante, el DRAE se empeña en usar abreviaturas. Porque está todavía sujeto a las servidumbres del papel preveintiunero. Podéis o no estar de acuerdo conmigo. Pero en el primer caso, cuidado con manifestar demasiado abiertamente vuestra postura o

Hoy, gelmaneo así: bar <- function(n, reps = 1e4){ foo <- function(n){ x <- rnorm(n) tmp <- t.test(x) c(tmp$p.value, abs(mean(x))) } res <- replicate(reps, foo(n)) tmp <- t(res) tmp <- tmp[tmp[,1] < 0.05,] tmp[,2] } res <- lapply(c(3, 10, 20, 50, 100), bar) sapply(res, mean) #[1] 0.8662636 0.6583157 0.4934551 0.3240322 0.2337086 Resumo: Fabrico un montón de errores de tipo I. Recuérdese: error de tipo I implica artículo publicado. Hago variar el número de sujetos (3, 10, etc.), n. Mido el tamaño (promedio) del efecto, E; el estudio de su distribución, ejercicio para el lector. Y efectivamente, E es función decreciente de n. ...

Hay una cifra que me intriga. Ya hablé sobre ella hace años. Buscando el otro día (infructosísimamente) información sobre la renta (y su distribución) de los autónomos, volví a encontrármela actualizada: En resumen: el en 2015, 3.5 millones de empresarios y autónomos declararon ingresos netos (descontando gastos) de poco más de 8000 euros por barba. Hipótesis que barajo: Que ese número no signifique lo que y pienso que significa. Que en el denominador hay autónomos inactivos, autónomos a tiempo ultraparcial, autónomos que se dieron de alta un ratico para cobrar una factura pírrica, etc. Que alguien que ni es autónomo ni lo ha sido; que ni ha creado un puesto de trabajo ni es probable que lo cree nunca (ni siquiera para sí mismo) escribirá un comentario que contenga la palabra fraude rodeado del habitual blablablá.

Hace tiempo escribí: El porcentaje de mosquitos que pueden transmitir enfermedades es pequeño, muy pequeño. Sin embargo, decimos mosquitos transmiten enfermedades sin empacho. Un porcentaje mucho mayor de los libros tienen tapas blandas. Sin embargo, no decimos que los libros tienen tapas blandas. (Las frases anteriores no son estrictamente mías: en el enlace hay otros enlaces a sus fuentes y autores). Hoy leo en el periódico un artículo que lleva por título ¿Por qué los hombres matan a las mujeres? ...

Vienen a ser la versión geo de las sparklines. Por ejemplo, Notas: El gráfico anterior no es mío. El código (y datos) con el que se generó tampoco. Son de Susana Huedo, exalumna. Está basado (todo hay que decirlo) en código de terceros y debería acordarme de cuál de ellos. Pero no es el caso.