{ "cells": [ { "cell_type": "markdown", "id": "8ea1197c", "metadata": {}, "source": [ "# Constrained Bayesian model: measurement noise\n" ] }, { "cell_type": "code", "execution_count": 1, "id": "f1a0c7ed", "metadata": { "execution": { "iopub.execute_input": "2026-09-14T14:53:03.064830Z", "iopub.status.busy": "2026-09-14T14:53:03.063787Z", "iopub.status.idle": "2026-09-14T14:53:10.135442Z", "shell.execute_reply": "2026-09-14T14:53:10.133401Z" } }, "outputs": [], "source": [ "from pathlib import Path\n", "import sys\n", "\n", "PROJECT_ROOT = Path.cwd().resolve()\n", "while not (PROJECT_ROOT / 'data.py').exists():\n", " if PROJECT_ROOT.parent == PROJECT_ROOT:\n", " raise FileNotFoundError('Open this notebook from inside organized_uncertainty_analysis.')\n", " PROJECT_ROOT = PROJECT_ROOT.parent\n", "if str(PROJECT_ROOT.parent) not in sys.path:\n", " sys.path.insert(0, str(PROJECT_ROOT.parent))\n", "\n", "import os\n", "os.chdir(PROJECT_ROOT)\n", "from organized_uncertainty_analysis.bayesian import noise_measurement as noise\n" ] }, { "cell_type": "code", "execution_count": 2, "id": "992c99d7", "metadata": { "execution": { "iopub.execute_input": "2026-09-14T14:53:10.142493Z", "iopub.status.busy": "2026-09-14T14:53:10.142493Z", "iopub.status.idle": "2026-09-14T15:02:39.895660Z", "shell.execute_reply": "2026-09-14T15:02:39.893114Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Loading SUV from: .\\suv_percentilesSLOthenUWM.mat\n", "Loading FLAGS from: .\\flags_combined.mat\n", "Loading SUV from: .\\suv_percentilesSLOthenUWM.mat\n", "Loading FLAGS from: .\\flags_combined.mat\n" ] }, { "data": { "text/html": [ "
| \n", " | Dataset | \n", "NoiseType | \n", "Sigma | \n", "x50 | \n", "Delta_refit_at_x50 | \n", "Delta_tta_at_x50 | \n", "N_refit_success | \n", "N_tta | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "FULL | \n", "multiplicative | \n", "0.129 | \n", "1.791702 | \n", "0.344217 | \n", "0.753202 | \n", "150 | \n", "300 | \n", "
| 1 | \n", "FULL | \n", "additive | \n", "0.144 | \n", "1.791702 | \n", "0.366052 | \n", "0.716211 | \n", "150 | \n", "300 | \n", "
| 2 | \n", "TRIM | \n", "multiplicative | \n", "0.129 | \n", "1.754496 | \n", "0.644678 | \n", "0.994051 | \n", "150 | \n", "300 | \n", "
| 3 | \n", "TRIM | \n", "additive | \n", "0.144 | \n", "1.754496 | \n", "1.000000 | \n", "0.941688 | \n", "150 | \n", "300 | \n", "