Newer
Older
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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
{
"cells": [
{
"cell_type": "markdown",
"id": "411ba7b3-7d56-45fe-b01e-205275e1988a",
"metadata": {},
"source": [
"# Des biais et des erreurs communes"
]
},
{
"cell_type": "markdown",
"id": "4e2fcf4b-d8aa-4bb2-8eab-dfe9a3210604",
"metadata": {},
"source": [
"Les exercices suivants sont destinés à vous familiariser avec les concepts appréhendés lors de l’introduction au *machine learning*. Avant toute chose, importez les librairies utiles :"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "4bbfb43b-1feb-4366-b1e7-5536f0f5aacd",
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import pandas as pd\n",
"import seaborn as sns\n",
"\n",
"sns.set_context('notebook')"
]
},
{
"cell_type": "markdown",
"id": "61c8d84f-a791-425e-ae70-306f0da93a55",
"metadata": {},
"source": [
"## Les relations à distance"
]
},
{
"cell_type": "markdown",
"id": "057d738a-a8a8-4d38-9dd2-b109d1325308",
"metadata": {},
"source": [
"Il paraît que l’univers est en expansion et que cette expansion va en s’accélérant. C’est en tout cas ce que l’étude de Wendy Freedman et al. a prouvé ([*Freedman, 2001*](../0.about-datasets.ipynb#Stellar-Objects)). Par conséquent, on s’attend à ce qu’un objet stellaire s’éloigne d’autant plus vite de nous que la distance qui nous sépare de lui est grande.\n",
"\n",
"Chargeons le jeu de données en se concentrant sur des objets proches de nous (entre 30 000 et 100 000 années-lumières) :"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "1cf3ab56-418f-46e3-bc3f-36cf0eec0dbf",
"metadata": {},
"outputs": [],
"source": [
"# load data\n",
"df = pd.read_csv(\"../files/stellar-objects.csv\", sep=\"\\t\")\n",
"\n",
"# distance: megaparsec (MPC)\n",
"# velocity: in km/s\n",
"df[\"velocity\"] = df.v_helio.fillna(df.v_flow.fillna(df.v_cmb))\n",
"\n",
"# objects close to earth, but not that close :)\n",
"data = df[(df.distance > 10) & (df.distance < 30)]"
]
},
{
"cell_type": "markdown",
"id": "f0a306e1-be3e-4431-84a3-32216340c326",
"metadata": {},
"source": [
"Affichons un nuage de points afin de vérifier la proposition de ces pontes de la NASA :"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "1fb0d73f-62bd-4777-b4e4-276554e2a599",
"metadata": {},
"outputs": [],
"source": [
"sns.scatterplot(data=data, x=\"distance\", y=\"velocity\")\n",
"\n",
"sns.despine()\n",
"\n",
"plt.title(\"Relation between distance and velocity of stellar objects\")\n",
"plt.xlabel(\"Distance (MPC)\")\n",
"plt.ylabel(\"Velocity (km/s)\")\n",
"\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"id": "fbb20849-4a22-4870-940b-8067fd06e548",
"metadata": {},
"source": [
"Rien de bien concluant à première vue, non ? Afin de déterminer visuellement s’il existe bien une relation linéaire entre la distance et la vitesse d’éloignement, affichez une droite de régression :"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "125c4241-faf9-4209-b8c6-cfc2c1b07105",
"metadata": {},
"outputs": [],
"source": [
"# your code here\n",
"\n",
"_ = sns.regplot(data=data, x=\"distance\", y=\"velocity\")"
]
},
{
"cell_type": "markdown",
"id": "aa3c4eeb-5ce4-44f2-9403-9d50a9e425e9",
"metadata": {},
"source": [
"Bon, appelez BFM TV, Wendy s’est trompée : 2/3 des points sont en dehors de l’intervalle de confiance à 95 %. Ou alors, peut-être avons-nous fait une erreur de méthodologie ?"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.6"
}
},
"nbformat": 4,
"nbformat_minor": 5
}