Somewhere else, part 139

Some writings worth reading

and another one, how far away from the office could I get in 10 minutes? (http://flaviogortana.com/isoscope/  ht @goulu)

(including 10min if you want to walk)

and “John Oliver may have crashed the FCC’s comments system” http://avclub.com/article/…

  • “The Overwhelming Persistence of Neighborhood Poverty” http://citylab.com/politics/2014/05/ … “For all the debate, gentrification is far from the norm”
  • Interesting exercise on conditional probabilities (via @WTFViz, meaning that there is some problem here, for those who, indeed, cannot do the maths)

et un peu de lecture en français

Did I miss something?

Approches Statistiques du Risque

Il semble que l’ouvrage approches statistiques du risque, édité par Jean Jacques Droesbeke, Myriam Maumy-Bertrand, Gilbert Saporta et Christine Thomas-Agnan soit finalement paru, selon Gilbert Saporta. Ou devraient paraître – au pire – dans les jours à venir (je le vois pas encore sur le site de Technip). Les 400 pages reprennent les mini-cours que nous avions fait aux Journées d’Etudes Statistique de 2010, au CIRM, avec Patrice Bertail, Anne-Laure Fougères, Stéphane Loisel, Philippe NaveauEric Parent et Jean-Marc Tallon. Pour ma part, mes slides et mes notes de cours traînent quelque part sur le blog. Je garde un très bon souvenir de cette semaine de mini-cours, qui présentaient l’avantage d’être très complets, avec une introduction aux mesures de risques et aux lien avec la théorie de la décision, à la théorie du risque, aux notions simples d’extrêmes, de lois multivariés, et d’extrêmes multivariés, avec des applications en hydrologie, en santé, et en assurance. Avec bien entendu l’introduction historique de Jean Jacques. Même si les documents ont été écrits voilà déjà presque cinq ans, j’ai survolé mes chapitres il y a quelques semaines, et j’ai l’impression qu’ils restent d’actualité pour tous ceux qui s’intéressent au sujet (et ils sont nombreux me semble-t-il).

Box plot, Fisher’s style

In a recent issue of Significance, I discovered an interesting – and amuzing – figure, about some box & beard plot, in Dr Fisher’s casebook: Beard the statistician in his den.

In French, the box plot (introduced by John Tukey, not George Box, as discussed in a previous post) is popular under the name boîte à moustaches (box with a mustache, for a simple translation).

> set.seed(2)
> x=rnorm(500)
> boxplot(x,horizontal=TRUE,axes=FALSE)
> axis(1)

I was wondering if it was possible to reproduce that Fisher’s style of box plot, with a beard (Ronald Fisher was famous for his beard). Technically, it is not complicated,

> boxplot(x,horizontal=TRUE,xlim=c(-1,1.3),axes=FALSE)
> axis(1)
> Q=quantile(x,c(.25,.75))
> y=cut(x[(x>=Q[1])&(x<=Q[2])],seq(Q[1],Q[2],length=11))
> tb=table(y)
> u=seq(Q[1],Q[2],length=11)
> umid=(u[1:10]+u[2:11])/2
> for(i in 1:10) segments(umid[i],1-.2,umid[i],1-.2-tb[i]/20,lwd=3)

Yes, I kept the mustaches here, but ploting a simple box, one could easily draw a goatee, or a chin strip

> rect(Q[1],1-.2,Q[2],1+.2)
> segments(median(x),1-.2,median(x),1+.2,lwd=2)

The only problem, is that between the first and the third quartiles, the distribution is much more ‘uniform‘ than it might look on the figure in Significance. Unless, perhaps, if add more bars on the left of the first quartile, and on the right of the third one,

> du=diff(umid)[1]
> y=cut(x,seq(Q[1]-du*8,Q[2]+du*8,length=11+16))
> tb=table(y)
> u=seq(Q[1]-du*8,Q[2]+du*8,length=11+16)
> umid=(u[1:26]+u[2:27])/2
> for(i in 1:8) segments(umid[i],1,umid[i],1-.2-tb[i]/20,lwd=3)
> for(i in 19:26) segments(umid[i],1,umid[i],1-.2-tb[i]/20,lwd=3)

Of course, we can try to add a smile on that face, below the mustache,

> vu=seq(-2.5,2.5,by=.02)
> vv=dnorm(vu)
> lines(vu,1-vv*4,col="red",lty=2)
> d=density(x,bw=.1)
> lines(d$x,1-d$y*4,col="red")

So far, I am not convinced. I mean, it should be possible to add something on the box plot,  to add more information. But I still wonder what… maybe using colors?

Somewhere else, part 138

Some writings worth reading, here and there,

vs “Men and prostate cancer in Europe”

et un peu de lecture en français,

Did I miss something else interesting?