From bf0eaaf6199ceadbf8d6ce22256ec26fa131922c Mon Sep 17 00:00:00 2001 From: Benoit Carry Date: Mon, 17 Jul 2023 16:41:13 +0200 Subject: [PATCH] added stats on observations --- notebooks/fig_spin_valid.ipynb | 40 ++++++++++++++++++++++++++++++++++ 1 file changed, 40 insertions(+) diff --git a/notebooks/fig_spin_valid.ipynb b/notebooks/fig_spin_valid.ipynb index 3645592..5753e05 100644 --- a/notebooks/fig_spin_valid.ipynb +++ b/notebooks/fig_spin_valid.ipynb @@ -324,6 +324,46 @@ "# len(data), len(data[mask]), len(data[mask_HG1G2]), len(data[mask_HG1G2sp]), len(data[mask_sHG1G2_g]), len(data[mask_sHG1G2_r]), len(data[mask_sHG1G2])" ] }, + { + "cell_type": "markdown", + "id": "e20c6a59", + "metadata": {}, + "source": [ + "# Stat observations" + ] + }, + { + "cell_type": "code", + "execution_count": 30, + "id": "bdc0da14", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Observation $89_{-24}^{+41}$\n", + "Days $1068_{-82}^{+64}$\n", + "Observation $101_{-29}^{+44}$\n", + "Days $1074_{-70}^{+71}$\n" + ] + } + ], + "source": [ + "q_obs = data['n_obs'].quantile([0.25,0.50,0.75])\n", + "q_day = data['n_days'].quantile([0.25,0.50,0.75])\n", + "\n", + "print('Observation ${:.0f}_{{{:-.0f}}}^{{{:+.0f}}}$'.format( q_obs[0.50], q_obs[0.25]-q_obs[0.50], q_obs[0.75]-q_obs[0.50] ))\n", + "print('Days ${:.0f}_{{{:-.0f}}}^{{{:+.0f}}}$'.format( q_day[0.50], q_day[0.25]-q_day[0.50], q_day[0.75]-q_day[0.50] ))\n", + "\n", + "q_obs = data.loc[maskFINK,'n_obs'].quantile([0.25,0.50,0.75])\n", + "q_day = data.loc[maskFINK,'n_days'].quantile([0.25,0.50,0.75])\n", + "\n", + "print('Observation ${:.0f}_{{{:-.0f}}}^{{{:+.0f}}}$'.format( q_obs[0.50], q_obs[0.25]-q_obs[0.50], q_obs[0.75]-q_obs[0.50] ))\n", + "print('Days ${:.0f}_{{{:-.0f}}}^{{{:+.0f}}}$'.format( q_day[0.50], q_day[0.25]-q_day[0.50], q_day[0.75]-q_day[0.50] ))\n", + "\n" + ] + }, { "cell_type": "markdown", "id": "ed6e70f4-d774-4c4b-bd59-449c227b9108",