From e139f2e7dc8219240198ef896836d5970ebeba02 Mon Sep 17 00:00:00 2001 From: Stefano Covino Date: Wed, 30 Apr 2025 12:45:15 +0200 Subject: [PATCH] 300425 commit --- .../Lecture-Introduction.ipynb | 9 +- .../Lecture-Time-of-Arrival.ipynb | 45 +++-- .../Science Case - FRBs/Lecture-FRBs.ipynb | 163 ++++++++++++------ Lectures/Science Case - FRBs/Manifest.toml | 2 +- README.md | 2 +- 5 files changed, 136 insertions(+), 85 deletions(-) diff --git a/Lectures/Lecture - Introduction/Lecture-Introduction.ipynb b/Lectures/Lecture - Introduction/Lecture-Introduction.ipynb index dee834b..f289f3a 100644 --- a/Lectures/Lecture - Introduction/Lecture-Introduction.ipynb +++ b/Lectures/Lecture - Introduction/Lecture-Introduction.ipynb @@ -165,10 +165,11 @@ "13. 17/4 - Wavelet analysis\n", "14. 30/4 - Time of arrival analysis\n", "15. 7/5 - Guest lecture - Exoplanets\n", - "16. 8/4 - Non-parametric methods\n", - "17. 14/5 - Gaussian processes\n", - "18. 15/5 - Science case: GRBs\n", - "19. 21/5 - Astrostatistics final considerations" + "16. 8/5 - Science case: FRBs\n", + "17. 14/5 - Non-parametric methods\n", + "18. 15/5 - Gaussian processes\n", + "19. 21/5 - Science case: GRBs\n", + "20. 22/5 - Astrostatistics final considerations" ] }, { diff --git a/Lectures/Lecture - Time of Arrival/Lecture-Time-of-Arrival.ipynb b/Lectures/Lecture - Time of Arrival/Lecture-Time-of-Arrival.ipynb index 6fe0864..909a4af 100644 --- a/Lectures/Lecture - Time of Arrival/Lecture-Time-of-Arrival.ipynb +++ b/Lectures/Lecture - Time of Arrival/Lecture-Time-of-Arrival.ipynb @@ -25,7 +25,7 @@ }, { "cell_type": "code", - "execution_count": 6, + "execution_count": 9, "id": "4ff28a4e", "metadata": { "id": "4ff28a4e" @@ -35,7 +35,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "\u001b[32m\u001b[1m Activating\u001b[22m\u001b[39m project at `/mnt/chromeos/GoogleDrive/MyDrive/Lab/Teaching/Insubria/Docs_2024_25/Lectures/Lecture - Time of Arrival`\n" + "\u001b[32m\u001b[1m Activating\u001b[22m\u001b[39m project at `~/Library/CloudStorage/GoogleDrive-stefano.covino@inaf.it/Il mio Drive/Lab/Teaching/Insubria/Docs_2024_25/Lectures/Lecture - Time of Arrival`\n" ] } ], @@ -45,7 +45,7 @@ }, { "cell_type": "code", - "execution_count": 7, + "execution_count": 10, "id": "67a89b22", "metadata": {}, "outputs": [], @@ -55,7 +55,7 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": 3, "id": "797d93b6", "metadata": {}, "outputs": [], @@ -143,7 +143,7 @@ "### Exercize: DO-climate- events\n", "***\n", "\n", - "- In the paper ([Ditlevsen et al. 2006](https://ui.adsabs.harvard.edu/abs/2006AGUFMGC24A..07D/abstract)), a possible periodicity of about 1500 years for the Dansgaard-Oeschger (DO) events,observed in the Greenland ice cores is discussed.\n", + "- In the paper ([Ditlevsen et al. 2006](https://ui.adsabs.harvard.edu/abs/2006AGUFMGC24A..07D/abstract)), a possible periodicity of about 1500 years for the Dansgaard-Oeschger (DO) events, observed in the Greenland ice cores is discussed.\n", "\n", "- According to Wikipedia, [Dansgaard–Oeschger events](https://en.wikipedia.org/wiki/Dansgaard%E2%80%93Oeschger_event) are rapid climate fluctuations that occurred during the last glacial period. Some scientists say that the events occur quasi-periodically with a recurrence time being a multiple of 1,470 years, but this is debated.\n", "\n", @@ -166,7 +166,7 @@ }, { "cell_type": "code", - "execution_count": 9, + "execution_count": 4, "id": "141dde68", "metadata": {}, "outputs": [], @@ -184,7 +184,7 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": 5, "id": "056e4daa", "metadata": {}, "outputs": [ @@ -194,7 +194,7 @@ "z2n (generic function with 1 method)" ] }, - "execution_count": 10, + "execution_count": 5, "metadata": {}, "output_type": "execute_result" } @@ -228,7 +228,7 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": 6, "id": "2e2aadae", "metadata": {}, "outputs": [], @@ -240,7 +240,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 7, "id": "748059db", "metadata": {}, "outputs": [ @@ -248,35 +248,30 @@ "data": { "text/html": [ "
\n", - "
\n", + "
\n", "
\n", " \n", " \n", - "
\n", - " \n", - "
\n", + "
\n", "
\n", "
\n", - " \n", - " \n", + " \n", + " \n", " \n", "
\n", "
\n", - "
\n", + "
\n", "
\n", " \n", "
\n", " \n", "
\n", "
\n", "
\n", - " \n", + " \n", "
\n", " \n", "
\n", @@ -285,7 +280,7 @@ "
" ] }, - "execution_count": 12, + "execution_count": 7, "metadata": {}, "output_type": "execute_result" } @@ -317,7 +312,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": 8, "id": "e39ee398", "metadata": {}, "outputs": [ @@ -387,7 +382,7 @@ "source": [ "- The Poisson Distribution probability mass function gives the probability of observing $k$ events in a time period given the length of the period and the average events per time:\n", "\n", - "$$ P(({\\rm k\\ events\\ in\\ a\\ time\\ period}) = e^{-\\frac{\\rm event}{\\rm time} \\times {\\rm \\ time\\ period}} \\frac{\\left( \\frac{\\rm event}{\\rm time} \\times {\\rm \\ time\\ period} \\right)^k}{k!}$$\n", + "$$ P({\\rm k\\ events\\ in\\ a\\ time\\ period}) = e^{-\\frac{\\rm event}{\\rm time} \\times {\\rm \\ time\\ period}} \\frac{\\left( \\frac{\\rm event}{\\rm time} \\times {\\rm \\ time\\ period} \\right)^k}{k!}$$\n", "\n", "- ${\\rm events\\ /\\ time \\times time\\ period}$ is usually simplified into a single parameter, $λ$, the rate parameter:\n", " $$ P({\\rm k\\ events\\ in\\ interval}) = e^{-\\lambda} \\times \\frac{\\lambda^k}{k!} $$\n", diff --git a/Lectures/Science Case - FRBs/Lecture-FRBs.ipynb b/Lectures/Science Case - FRBs/Lecture-FRBs.ipynb index d33b485..236ce43 100644 --- a/Lectures/Science Case - FRBs/Lecture-FRBs.ipynb +++ b/Lectures/Science Case - FRBs/Lecture-FRBs.ipynb @@ -25,7 +25,7 @@ }, { "cell_type": "code", - "execution_count": 20, + "execution_count": 1, "id": "4ff28a4e", "metadata": { "id": "4ff28a4e" @@ -35,7 +35,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "\u001b[32m\u001b[1m Activating\u001b[22m\u001b[39m project at `/mnt/chromeos/GoogleDrive/MyDrive/Lab/Teaching/Insubria/Docs_2024_25/Lectures/Science Case - FRBs`\n" + "\u001b[32m\u001b[1m Activating\u001b[22m\u001b[39m project at `~/Library/CloudStorage/GoogleDrive-stefano.covino@inaf.it/Il mio Drive/Lab/Teaching/Insubria/Docs_2024_25/Lectures/Science Case - FRBs`\n" ] } ], @@ -45,30 +45,71 @@ }, { "cell_type": "code", - "execution_count": 21, + "execution_count": 2, "id": "3a6636b9-9c79-45f7-a3bc-945c9b67b048", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "\u001b[32m\u001b[1m Installed\u001b[22m\u001b[39m Trapz ─ v2.0.3\n", + "\u001b[92m\u001b[1mPrecompiling\u001b[22m\u001b[39m project...\n", + " 30735.3 ms\u001b[32m ✓ \u001b[39mTrapz\n", + " 6781.9 ms\u001b[32m ✓ \u001b[39m\u001b[90mSixel\u001b[39m\n", + " 29613.0 ms\u001b[32m ✓ \u001b[39m\u001b[90mPNGFiles\u001b[39m\n", + " 3372.1 ms\u001b[32m ✓ \u001b[39m\u001b[90mWebP\u001b[39m\n", + " 13978.3 ms\u001b[32m ✓ \u001b[39m\u001b[90mBonito\u001b[39m\n", + " 20658.4 ms\u001b[32m ✓ \u001b[39mCSV\n", + " 138813.9 ms\u001b[32m ✓ \u001b[39m\u001b[90mMakie\u001b[39m\n", + " 108462.4 ms\u001b[32m ✓ \u001b[39mWGLMakie\n", + " 8 dependencies successfully precompiled in 307 seconds. 295 already precompiled.\n", + " \u001b[91m1\u001b[39m dependency errored.\n", + " For a report of the errors see `julia> err`. To retry use `pkg> precompile`\n" + ] + } + ], "source": [ "Pkg.instantiate()" ] }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 18, "id": "42e23f92-3b70-40a1-83ca-03919a77bee8", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "\u001b[32m\u001b[1m Resolving\u001b[22m\u001b[39m package versions...\n", + "\u001b[32m\u001b[1m Updating\u001b[22m\u001b[39m `~/Library/CloudStorage/GoogleDrive-stefano.covino@inaf.it/Il mio Drive/Lab/Teaching/Insubria/Docs_2024_25/Lectures/Science Case - FRBs/Project.toml`\n", + " \u001b[90m[27cdad02] \u001b[39m\u001b[93m~ GregoryLoredo v0.1.0 `/mnt/chromeos/GoogleDrive/Computers/Il mio iMac/Programmi/julia/GregoryLoredo.jl` ⇒ v0.1.0 `~/Programmi/julia/GregoryLoredo.jl`\u001b[39m\n", + "\u001b[32m\u001b[1m Updating\u001b[22m\u001b[39m `~/Library/CloudStorage/GoogleDrive-stefano.covino@inaf.it/Il mio Drive/Lab/Teaching/Insubria/Docs_2024_25/Lectures/Science Case - FRBs/Manifest.toml`\n", + " \u001b[90m[27cdad02] \u001b[39m\u001b[93m~ GregoryLoredo v0.1.0 `/mnt/chromeos/GoogleDrive/Computers/Il mio iMac/Programmi/julia/GregoryLoredo.jl` ⇒ v0.1.0 `~/Programmi/julia/GregoryLoredo.jl`\u001b[39m\n" + ] + } + ], "source": [ - "#Pkg.develop(path=\"/mnt/chromeos/GoogleDrive/Computers/Il mio iMac/Programmi/julia/GregoryLoredo.jl\")" + "#Pkg.develop(path=\"/mnt/chromeos/GoogleDrive/Computers/Il mio iMac/Programmi/julia/GregoryLoredo.jl\")\n", + "Pkg.develop(path=\"/Users/covino/Programmi/julia/GregoryLoredo.jl\")" ] }, { "cell_type": "code", - "execution_count": 23, + "execution_count": 19, "id": "4674aa93-46bd-467b-8738-3e5d05a38048", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "\u001b[36m\u001b[1m[ \u001b[22m\u001b[39m\u001b[36m\u001b[1mInfo: \u001b[22m\u001b[39mPrecompiling GregoryLoredo [27cdad02-83ed-4657-a9f5-b5f94b1c509a] \n" + ] + } + ], "source": [ "using CSV\n", "using DataFrames\n", @@ -205,7 +246,7 @@ "id": "0b1fe3a1", "metadata": {}, "source": [ - "### Exercize: Analyse the FRB121102 by trhe Gregory & Loredo algorithm\n", + "### Exercize: Analyse the FRB121102 by the Gregory & Loredo algorithm\n", "***\n", "\n", "- We are going to the dataset provided by [Rajwade et al. (2020)](https://ui.adsabs.harvard.edu/abs/2020MNRAS.495.3551R/abstract):" @@ -213,7 +254,7 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 20, "id": "8faaa852-2143-41b4-9caa-ba5a3119ab74", "metadata": {}, "outputs": [], @@ -240,33 +281,47 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 21, "id": "e6680267-970b-4a9e-aef0-8e5941680d41", "metadata": {}, "outputs": [ { "data": { "text/html": [ - "
\n", - "
\n", - " \n", + "
\n", + "
\n", "
\n", - " \n", + " \n", + "
\n", + "
\n", + "
\n", + " \n", + " \n", + " \n", + "
\n", + "
\n", + "
\n", + "
\n", + " \n", + "
\n", + " \n", + "
\n", "
\n", - "
\n", - "
\n", - " \n", - "
\n", - " \n", + "
\n", + " \n", + "
\n", + " \n", + "
\n", "
\n", "
\n", "
" ] }, - "execution_count": 25, + "execution_count": 21, "metadata": {}, "output_type": "execute_result" } @@ -294,7 +349,7 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 22, "id": "16fc6c4f", "metadata": {}, "outputs": [ @@ -304,7 +359,7 @@ "logrange (generic function with 1 method)" ] }, - "execution_count": 26, + "execution_count": 22, "metadata": {}, "output_type": "execute_result" } @@ -315,7 +370,7 @@ }, { "cell_type": "code", - "execution_count": 27, + "execution_count": 23, "id": "647669d9", "metadata": {}, "outputs": [ @@ -336,7 +391,7 @@ }, { "cell_type": "code", - "execution_count": 28, + "execution_count": 24, "id": "7e6fbc97", "metadata": {}, "outputs": [], @@ -366,7 +421,7 @@ }, { "cell_type": "code", - "execution_count": 29, + "execution_count": 25, "id": "a43d10da", "metadata": {}, "outputs": [], @@ -376,7 +431,7 @@ }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 26, "id": "a4496d32", "metadata": {}, "outputs": [ @@ -408,33 +463,33 @@ }, { "cell_type": "code", - "execution_count": 31, + "execution_count": 27, "id": "0d4cf70f", "metadata": {}, "outputs": [ { "data": { "text/html": [ - "
\n", + "
\n", "
\n", " \n", "
\n", " \n", "
\n", "
\n", "
\n", - " \n", - "
\n", - " \n", + " \n", + "
\n", + " \n", "
\n", "
\n", "
" ] }, - "execution_count": 31, + "execution_count": 27, "metadata": {}, "output_type": "execute_result" } @@ -467,33 +522,33 @@ }, { "cell_type": "code", - "execution_count": 32, + "execution_count": 28, "id": "6d14ad48", "metadata": {}, "outputs": [ { "data": { "text/html": [ - "
\n", + "
\n", "
\n", " \n", "
\n", " \n", "
\n", "
\n", "
\n", - " \n", - "
\n", - " \n", + " \n", + "
\n", + " \n", "
\n", "
\n", "
" ] }, - "execution_count": 32, + "execution_count": 28, "metadata": {}, "output_type": "execute_result" } @@ -529,7 +584,7 @@ }, { "cell_type": "code", - "execution_count": 33, + "execution_count": 29, "id": "c92ce08b", "metadata": {}, "outputs": [ @@ -539,7 +594,7 @@ "Z2N (generic function with 1 method)" ] }, - "execution_count": 33, + "execution_count": 29, "metadata": {}, "output_type": "execute_result" } @@ -565,7 +620,7 @@ }, { "cell_type": "code", - "execution_count": 34, + "execution_count": 30, "id": "069bdfa7", "metadata": {}, "outputs": [], @@ -575,33 +630,33 @@ }, { "cell_type": "code", - "execution_count": 35, + "execution_count": 31, "id": "7c813af5", "metadata": {}, "outputs": [ { "data": { "text/html": [ - "
\n", + "
\n", "
\n", " \n", "
\n", " \n", "
\n", "
\n", "
\n", - " \n", - "
\n", - " \n", + " \n", + "
\n", + " \n", "
\n", "
\n", "
" ] }, - "execution_count": 35, + "execution_count": 31, "metadata": {}, "output_type": "execute_result" } @@ -622,7 +677,7 @@ }, { "cell_type": "code", - "execution_count": 36, + "execution_count": 32, "id": "d21d306f", "metadata": {}, "outputs": [ diff --git a/Lectures/Science Case - FRBs/Manifest.toml b/Lectures/Science Case - FRBs/Manifest.toml index c784216..f1eda2f 100644 --- a/Lectures/Science Case - FRBs/Manifest.toml +++ b/Lectures/Science Case - FRBs/Manifest.toml @@ -473,7 +473,7 @@ version = "1.3.14+1" [[deps.GregoryLoredo]] deps = ["Trapz"] -path = "/mnt/chromeos/GoogleDrive/Computers/Il mio iMac/Programmi/julia/GregoryLoredo.jl" +path = "/Users/covino/Programmi/julia/GregoryLoredo.jl" uuid = "27cdad02-83ed-4657-a9f5-b5f94b1c509a" version = "0.1.0" diff --git a/README.md b/README.md index b5b9506..eaec2c2 100644 --- a/README.md +++ b/README.md @@ -2,4 +2,4 @@ This is a repository with material (notebooks, papers, etc.) for the **Time Domain Astrophysics** course delivered at the *Università dell'Insubria* by Stefano Covino. -*Last update: 29 April 2025.* +*Last update: 30 April 2025.* -- GitLab