Table of Contents
- -Exercice 1 : Ré-exécuter n'est pas répliquer…
-Exercice 1 : Ré-exécuter n'est pas répliquer…
+Même si la terminologie peut varier d'un auteur ou d'une communauté à l'autre, il est important de comprendre que l'on peut distinguer @@ -55,7 +45,7 @@ vous avez pu utiliser dans le module 2). Cet exercice peut donc se faire à deux niveaux :
-
-
- un niveau facile pour ceux qui repartiront du code dans le langage +
- un niveau facile pour ceux qui repartiront du code dans le langage qu'ils n'auront initialement pas utilisé et se contenteront de le ré-exécuter. Pour cela, nul besoin de maîtriser la régression logistique, il suffit de bien inspecter les sorties produites et de @@ -64,7 +54,7 @@ ceux qui ré-exécuteraient le notebook Python dans l'environnement Jupyter du MOOC, n'hésitez pas à consulter les ressources de la section 4A du module 2 qui expliquent comment y importer un notebook. -
- un niveau plus difficile pour ceux qui souhaiteront le réécrire +
- un niveau plus difficile pour ceux qui souhaiteront le réécrire
complètement (éventuellement dans un autre langage que R ou Python,
l'expérience peut être d'autant plus intéressante que nous n'avons
pas testé ces variations). Là, si les fonctions de calcul d'une
@@ -78,41 +68,41 @@ Vous pourrez alors discuter sur le forum des succès et des échecs que
vous aurez pu rencontrer. Pour cela :
-
-
- Vous publierez auparavant dans votre dépôt les différents notebooks +
- Vous publierez auparavant dans votre dépôt les différents notebooks en prenant bien soin d'enrichir votre document des informations (numéros de version, etc.) sur votre système et sur les bibliothèques installées. -
- Vous indiquerez votre résultat (que ça soit un succès ou échec à +
- Vous indiquerez votre résultat (que ça soit un succès ou échec à
obtenir les mêmes résultats) en remplissant ce tableau (vous avez
les droits d'édition donc il vous suffit d'éditer les fichiers via
l'interface GitLab). Vous vérifierez les valeurs obtenues pour :
-
-
- les coefficients de la pente et de l'intercept -
- les estimations d'erreur de ces coefficients -
- le goodness of fit -
- la figure -
- la zone de confiance +
- les coefficients de la pente et de l'intercept +
- les estimations d'erreur de ces coefficients +
- le goodness of fit +
- la figure +
- la zone de confiance
- +
Pour chacun vous indiquerez si le résultat est :
-
-
- identique -
- proche à moins de trois décimales -
- très différent -
- non fonctionnel (pas de résultat obtenu) +
- identique +
- proche à moins de trois décimales +
- très différent +
- non fonctionnel (pas de résultat obtenu)
Vous indiquerez également dans ce tableau :
-
-
- un lien vers votre espace gitlab contenant les différents notebooks -
- le nom du système d'exploitation utilisé -
- le langage utilisé et son numéro de version -
- les numéros des principales bibliothèques utilisées +
- un lien vers votre espace gitlab contenant les différents notebooks +
- le nom du système d'exploitation utilisé +
- le langage utilisé et son numéro de version +
- les numéros des principales bibliothèques utilisées
-
-
- Python : numpy, pandas, matplotlib, statsmodels… -
- R : BLAS, ggplot, dplyr si chargées +
- Python : numpy, pandas, matplotlib, statsmodels… +
- R : BLAS, ggplot, dplyr si chargées
Re-execution is not replication…
++Unfortunately terminology varies a lot between authors and +communities, but it is important to understand the distinction between +different levels of "replication". You can be satisfied with +re-running the code and get exactly the same results, but you can also +try to obtain similar results using a similar approach, changing for +example the programming language, computational method, etc. An +article we recommend on this topic is +https://arxiv.org/abs/1708.08205. +
+ ++Often the devil is in the details that one would have never thought +about, and we have had our share of surprises while preparing this +MOOC, in particular with the exercise on the Challenger catastrophe +from module 2. We therefore propose in this exercise that you re-do a +part of this analysis, following the example of Siddhartha Dallal and +co-authors almost 30 years ago in their article Risk Analysis of the +Space Shuttle: Pre-Challenger Prediction of Failure, published in the +Journal of the American Statistical Association (Vol. 84, No. 408, +Déc., 1989), but using a different language of your choosing (Python, +R, Julia, SAS…). +
+ ++Our experience shows that the estimations of slope and intercept are +generally the same, but there can be differences when looking at +variance estimators and R2 in more detail. Another source of +surprises is the final graphical presentation, depending on the +versions of the libraries that are used. +
+ ++The computations to be done are described at +https://app-learninglab.inria.fr/gitlab/moocrr-session1/moocrr-reproducibility-study/ +together with instructions for contributing. +
+ ++You will find there our replications of the computations by Dallal et +al. (in R), one in Python and one in R (very similar to what you have +used in module 2). This exercise can be done at two levels: +
+ +-
+
- an easy level at which you start from the code in the language that you did not use initially, and content yourself with re-executin it. This doesn't require mastering logistic regression, it is sufficien to inspect the outputs produced and check that they correspond to the expected values. For those who want to re-execute the Python notebook in our MOOC's Jupyter environment, check the resources for sequence 4A of module 2 that explain how to import a notebook. +
- a more difficult level at which you rewrite the analysis completely, possibly in a different language than Python or R, which makes the exercise more interesting because we have not tested such variants. If logistic regression is not already implemented for your language, you will need a good understanding of it in order to write the code yourself, which of course makes the exercise even more instructive. +
+You can discuss your successes or failures on the forum, after following these instructions: +
+-
+
- First, publish your notebooks in your repository, taking care to enrich your document with information about your system and your libraries (version numbers etc.). +
- Indicate your result by adding to this table (you have write permissions, so you can simply edit it via the GitLab interface). Check the values obtained for:
+
-
+
- the slope and intercept coefficients +
- the error estimates for these coefficients +
- the goodness of fit +
- the plot +
- the confidence region +
+ +For each of these values, specify if your result is +
+-
+
- identical +
- close, to three decimal places +
- very different +
- non functional (no result obtained) +
+Also provide in this table: +
+-
+
- a link to your GitLab workspace with your notebook(s) +
- your operating system +
- the language you used, with the version number +
- version numbers for the main libraries
+
-
+
- Python: numpy, pandas, matplotlib, statsmodels… +
- R: BLAS, ggplot, dplyr if used +
+
+
+Don't worry if these instructions seem confusing, they are reproduced above the table and you will quickly notice if something is missing when you try to add your data. +
+ ++We will compile a synthesis of the principal divergences observes and make it available at the end of the MOOC. +
+