diff --git a/polarimetry/figs/corner-nobound.pdf b/polarimetry/figs/corner-nobound.pdf
index ad3ac0b..5fbad5c 100644
Binary files a/polarimetry/figs/corner-nobound.pdf and b/polarimetry/figs/corner-nobound.pdf differ
diff --git a/polarimetry/figs/corner.pdf b/polarimetry/figs/corner.pdf
index 1130733..499f410 100644
Binary files a/polarimetry/figs/corner.pdf and b/polarimetry/figs/corner.pdf differ
diff --git a/polarimetry/figs/observation-fit.pdf b/polarimetry/figs/observation-fit.pdf
index 287ba2e..904519e 100644
Binary files a/polarimetry/figs/observation-fit.pdf and b/polarimetry/figs/observation-fit.pdf differ
diff --git a/polarimetry/figs/observation-mcmc-comb-fullrange.pdf b/polarimetry/figs/observation-mcmc-comb-fullrange.pdf
index 40f14af..f5a6f0f 100644
Binary files a/polarimetry/figs/observation-mcmc-comb-fullrange.pdf and b/polarimetry/figs/observation-mcmc-comb-fullrange.pdf differ
diff --git a/polarimetry/figs/observation-mcmc-comb.pdf b/polarimetry/figs/observation-mcmc-comb.pdf
index 1957446..a99ad01 100644
Binary files a/polarimetry/figs/observation-mcmc-comb.pdf and b/polarimetry/figs/observation-mcmc-comb.pdf differ
diff --git a/polarimetry/figs/trace-nobound.pdf b/polarimetry/figs/trace-nobound.pdf
index f6b314d..cd515c3 100644
Binary files a/polarimetry/figs/trace-nobound.pdf and b/polarimetry/figs/trace-nobound.pdf differ
diff --git a/polarimetry/figs/trace.pdf b/polarimetry/figs/trace.pdf
index 4ad8612..b09ff21 100644
Binary files a/polarimetry/figs/trace.pdf and b/polarimetry/figs/trace.pdf differ
diff --git a/polarimetry/pol_analyzer.ipynb b/polarimetry/pol_analyzer.ipynb
index 406d94f..cb0b243 100644
--- a/polarimetry/pol_analyzer.ipynb
+++ b/polarimetry/pol_analyzer.ipynb
@@ -101,7 +101,7 @@
"name": "stdout",
"output_type": "stream",
"text": [
- "This notebook was generated at 2021-09-03 13:54:43 (KST = GMT+0900) \n",
+ "This notebook was generated at 2021-09-08 20:38:05 (KST = GMT+0900) \n",
"0 Python 3.8.10 64bit [Clang 10.0.0 ]\n",
"1 IPython 7.25.0\n",
"2 OS macOS 10.16 x86_64 i386 64bit\n",
@@ -629,7 +629,7 @@
"name": "stderr",
"output_type": "stream",
"text": [
- "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_18676/3697028710.py:19: FutureWarning: In v4.0, pm.sample will return an `arviz.InferenceData` object instead of a `MultiTrace` by default. You can pass return_inferencedata=True or return_inferencedata=False to be safe and silence this warning.\n",
+ "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_55894/3697028710.py:19: FutureWarning: In v4.0, pm.sample will return an `arviz.InferenceData` object instead of a `MultiTrace` by default. You can pass return_inferencedata=True or return_inferencedata=False to be safe and silence this warning.\n",
" trace_b = pm.sample(20000, start=dict(h=0.1, a0=20, c1=0.1, c2=1.e-3))\n",
"Auto-assigning NUTS sampler...\n",
"Initializing NUTS using jitter+adapt_diag...\n",
@@ -655,7 +655,7 @@
" }\n",
" \n",
" \n",
- " 100.00% [84000/84000 00:48<00:00 Sampling 4 chains, 0 divergences]\n",
+ " 100.00% [84000/84000 00:42<00:00 Sampling 4 chains, 0 divergences]\n",
" \n",
" "
],
@@ -670,7 +670,7 @@
"name": "stderr",
"output_type": "stream",
"text": [
- "Sampling 4 chains for 1_000 tune and 20_000 draw iterations (4_000 + 80_000 draws total) took 59 seconds.\n"
+ "Sampling 4 chains for 1_000 tune and 20_000 draw iterations (4_000 + 80_000 draws total) took 51 seconds.\n"
]
}
],
@@ -705,7 +705,7 @@
"name": "stderr",
"output_type": "stream",
"text": [
- "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_18676/2577184737.py:1: DeprecationWarning: The function `traceplot` from PyMC3 is just an alias for `plot_trace` from ArviZ. Please switch to `pymc3.plot_trace` or `arviz.plot_trace`.\n",
+ "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_55894/2577184737.py:1: DeprecationWarning: The function `traceplot` from PyMC3 is just an alias for `plot_trace` from ArviZ. Please switch to `pymc3.plot_trace` or `arviz.plot_trace`.\n",
" pm.traceplot(trace_b)\n",
"/Users/ysbach/opt/anaconda3/envs/2005ud/lib/python3.8/site-packages/arviz/data/io_pymc3.py:96: FutureWarning: Using `from_pymc3` without the model will be deprecated in a future release. Not using the model will return less accurate and less useful results. Make sure you use the model argument or call from_pymc3 within a model context.\n",
" warnings.warn(\n"
@@ -1294,7 +1294,7 @@
"name": "stderr",
"output_type": "stream",
"text": [
- "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_18676/4211646313.py:19: FutureWarning: In v4.0, pm.sample will return an `arviz.InferenceData` object instead of a `MultiTrace` by default. You can pass return_inferencedata=True or return_inferencedata=False to be safe and silence this warning.\n",
+ "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_55894/4211646313.py:19: FutureWarning: In v4.0, pm.sample will return an `arviz.InferenceData` object instead of a `MultiTrace` by default. You can pass return_inferencedata=True or return_inferencedata=False to be safe and silence this warning.\n",
" trace_f = pm.sample(20000, start=dict(h=0.1, a0=20, c1=0.1, c2=1.e-3))\n",
"Auto-assigning NUTS sampler...\n",
"INFO:pymc3:Auto-assigning NUTS sampler...\n",
@@ -1324,7 +1324,7 @@
" }\n",
" \n",
" \n",
- " 100.00% [84000/84000 00:43<00:00 Sampling 4 chains, 0 divergences]\n",
+ " 100.00% [84000/84000 00:38<00:00 Sampling 4 chains, 0 divergences]\n",
" \n",
" "
],
@@ -1339,8 +1339,8 @@
"name": "stderr",
"output_type": "stream",
"text": [
- "Sampling 4 chains for 1_000 tune and 20_000 draw iterations (4_000 + 80_000 draws total) took 52 seconds.\n",
- "INFO:pymc3:Sampling 4 chains for 1_000 tune and 20_000 draw iterations (4_000 + 80_000 draws total) took 52 seconds.\n"
+ "Sampling 4 chains for 1_000 tune and 20_000 draw iterations (4_000 + 80_000 draws total) took 47 seconds.\n",
+ "INFO:pymc3:Sampling 4 chains for 1_000 tune and 20_000 draw iterations (4_000 + 80_000 draws total) took 47 seconds.\n"
]
}
],
@@ -1375,7 +1375,7 @@
"name": "stderr",
"output_type": "stream",
"text": [
- "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_18676/2459473891.py:1: DeprecationWarning: The function `traceplot` from PyMC3 is just an alias for `plot_trace` from ArviZ. Please switch to `pymc3.plot_trace` or `arviz.plot_trace`.\n",
+ "/var/folders/h6/tz6c2ch9771g23446ncvbczw0000gn/T/ipykernel_55894/2459473891.py:1: DeprecationWarning: The function `traceplot` from PyMC3 is just an alias for `plot_trace` from ArviZ. Please switch to `pymc3.plot_trace` or `arviz.plot_trace`.\n",
" pm.traceplot(trace_f)\n",
"/Users/ysbach/opt/anaconda3/envs/2005ud/lib/python3.8/site-packages/arviz/data/io_pymc3.py:96: FutureWarning: Using `from_pymc3` without the model will be deprecated in a future release. Not using the model will return less accurate and less useful results. Make sure you use the model argument or call from_pymc3 within a model context.\n",
" warnings.warn(\n"
@@ -1902,7 +1902,7 @@
"name": "stdout",
"output_type": "stream",
"text": [
- "This notebook was finished at 2021-09-03 13:58:50 (KST = GMT+0900) \n"
+ "This notebook was finished at 2021-09-08 20:41:48 (KST = GMT+0900) \n"
]
}
],