Since my first representation of the rank statistic as paired was incorrect, here is the histogram produced by the simulation

perm=sample(1:20)
saple[t]=sum(abs(sort(perm[1:10])-sort(perm[11:20])))

when . It is obviously much closer to zero than previously.

An interesting change is that the regression of the log-mean on produces

> lm(log(memean)~log(enn))
Call:
lm(formula = log(memean) ~ log(enn))
Coefficients:
(Intercept) log(enn)
-1.162 1.499

meaning that the mean is in rather than in or :

> summary(lm(memean~eth-1))
Coefficients:
Estimate Std. Error t value Pr(>|t|)
eth 0.3117990 0.0002719 1147 <2e-16 ***

with a very good fit.

Filed under: R, Statistics Tagged: Le Monde, puzzle, R, Spearman rank test

*Related*

To

**leave a comment** for the author, please follow the link and comment on their blog:

** Xi'an's Og » R**.

R-bloggers.com offers

**daily e-mail updates** about

R news and

tutorials on topics such as:

Data science,

Big Data, R jobs, visualization (

ggplot2,

Boxplots,

maps,

animation), programming (

RStudio,

Sweave,

LaTeX,

SQL,

Eclipse,

git,

hadoop,

Web Scraping) statistics (

regression,

PCA,

time series,

trading) and more...

If you got this far, why not

__subscribe for updates__ from the site? Choose your flavor:

e-mail,

twitter,

RSS, or

facebook...

**Tags:** Le Monde, puzzle, R, Spearman rank test, statistics